Волновые методы геофизики (георадиолокация и сейсморазведка)

Содержание

Слайд 2

Волновые методы геофизики Радиоволновое просвечивание Сейсморазведка Георадиолокация 08.02.2016 13:41

Волновые методы геофизики

Радиоволновое
просвечивание
Сейсморазведка
Георадиолокация

08.02.2016 13:41

Слайд 3

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

Георадилокация

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

среды на излучаемое электромагнитное поле. Метод работает в диапазоне частот от 50МГц до 3 000 МГц (3ГГц).

СЕЙСМОРАЗВЕДКА

ГЕОРАДИОЛОКАЦИЯ

упругие волны (уравнения движения)
длина волны
затухание
разрешающая способность

электромагнитные волны (уравнения Максвелла)
длина волны
затухание
разрешающая способность

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

08.02.2016 13:41

Слайд 4

2 Модели Исследователь работает в рамках двух моделей – модели реальной

2 Модели

Исследователь работает в рамках двух моделей – модели реальной среды

и модели способа распространения физического поля (теории) в модели реальной среды.
Наиболее ярким подтверждением высказанного тезиса служит книга Л.М. Бреховских «Волны в слоистых средах». Уже в названии содержится указание на две эти модели.
Вопрос о степени соответствия или адекватности той и другой модели – реальности.
Любой метод будет работать для ограниченного набора ситуаций – и при том с погрешностью!
Бесконечно разнообразная и изменяющаяся реальная среда - набор типизированных ситуаций, то есть моделей.
Ряд упрощений и допущений.
При решении задач с помощью системы уравнений стараются с помощью записи граничных условий свести эту систему к решаемой тем или иным способом, например, к системе из n уравнений с n неизвестными.
Даже самые сложные модели сред и самые сложные теории распространения физических полей в них являются существенными упрощениями реальности.

08.02.2016 13:41

Слайд 5

Надо: построить 2 модели и установить связи между ними. 1) решать

Надо: построить 2 модели и установить связи между ними.
1)

решать прямую задачу
2) решать обратную задачу
Обратные задачи геофизики НЕкорректны в смысле Адамара.
Условия Адамара: существование и единственность решения.
Вместе с тем хорошо известно о третьем условии корректности по Адамару в отношении непрерывной зависимости решения от данных задачи.
Какие условия не выполняются и как?

08.02.2016 13:41

Слайд 6

Колебания и волны ti - время прихода отраженного сигнала; А- амплитуда

Колебания и волны

ti - время прихода отраженного сигнала;
А- амплитуда отраженного сигнала;
Т

– видимый период сигнала, соответствующий 1/f , где f – центральная частота спектра сигнала.
Аi – Ai-1 - затухание георадиолокационного сигнала в слое, где Аi и Ai-1 сигналы, отраженные от кровли и подошвы i-того слоя.

08.02.2016 13:41

Слайд 7

Единство кинематических моделей среды. Какие это модели? 08.02.2016 13:41

Единство кинематических моделей среды.

Какие это модели?

08.02.2016 13:41

Слайд 8

Единство кинематических моделей среды. 08.02.2016 13:41 а-однородная среда; б-однородно-слоистая среда; в-непрерывная

Единство кинематических моделей среды.

08.02.2016 13:41

а-однородная среда;
б-однородно-слоистая среда;
в-непрерывная среда;
г-трансверсально-изотропная среда

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

Электрофизическая модель среды. Представляется в виде: а) слоистой толщи с постоянными

Электрофизическая модель среды.
Представляется в виде:
а) слоистой толщи с постоянными электрофизическими свойствами

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

08.02.2016 13:41

Слайд 10

Единство структуры данных Рис. ХХХ Результаты геофизических исследований с движущегося судна

Единство структуры данных

Рис. ХХХ Результаты геофизических исследований с движущегося судна на

акваториях:
А) методом сейсмоакустики на р. Волга; Б) методом георадиолокации на р. Чикапа (Республика Конго).

Результаты геофизических исследований с движущегося судна на акваториях:
А) методом сейсмоакустики на р. Волга;
Б) методом георадиолокации на р. Чикапа (Республика Конго).

08.02.2016 13:41

Слайд 11

Трубы (дифракция) Где с/р где г/р? 08.02.2016 13:42

Трубы (дифракция)

Где с/р где г/р?

08.02.2016 13:42

Слайд 12

Единый подход к решению прямых задач. Уравнение Гамильтона: (1) (∂t/∂x)2 +

Единый подход к решению прямых задач.
Уравнение Гамильтона: (1)
(∂t/∂x)2 + (∂t/∂y)2

+ (∂t/∂z)2 = 1/V2(x,y,z).
tотр. = 2Σ(hi/Vi); (2)
Аh = А0 Ki отр.( (1-K2j отр.))e-m2h /2h, (3)
где произведение ( (1-K2j отр.) учитывает потери интенсивности волны при прохождении границ, предшествующих находящейся на глубине h, экспоненциальный множитель e-m2h учитывает затухание в среде, а 1/2h – уменьшение амплитуды за счет сферического расхождения фронта волны.
Эти формулы действуют независимо от физической природы поля
Котр.= (ε11/2 -ε21/2)/(ε11/2 +ε21/2), (4)
где ε1 и ε2 – диэлектрические проницаемости сред выше и ниже границы соответственно, а в сейсмике
Котр. = (ρ2V2 - ρ1V1)/(ρ2V2 + ρ1V1), (5)
где произведения значений плотности на значение скорости распространения
упругих волн ρ1V1 и ρ2V2 – акустические жесткости верхней и нижней среды соответственно.


08.02.2016 13:41

Слайд 13

Свёрточная модель трассы. Моделирование. f(t) = f0(t) * g(t), (6) где

Свёрточная модель трассы. Моделирование.
f(t) = f0(t) * g(t), (6)
где трассы f(t)

– результата единичного акта “посылка-прием”,
f0(t) – зондирующий сигнал, а g(t) – импульсная характеристика среды.
Или в спектральной области:
S(jω) = S0(jω)K(jω), (7)
где S(jω)- спектр трассы,
S0(jω) – спектр зондирующего сигнала, а K(jω) –частотная характеристика среды.
При этом, g(t) или K(jω) расчитываются по заданной модели среды.
Из этого следует одинаковая обработка данных

08.02.2016 13:41

Слайд 14

RADEXPLORER Обрабатывающие процедуры: Удаление постоянной составляющей Смещение нуля времени Редакция трасс

RADEXPLORER
Обрабатывающие процедуры:
Удаление постоянной составляющей
Смещение нуля времени
Редакция трасс
Пространственная интерполяция


Подавление "звона" антенны
Двумерная фильтрация
Коррекция амплитуд
Предсказывающая деконволюция
Полосовая фильтрация
F-K Миграция Столта
Мгновенная амплитуда (преобразование Гильберта)
Учет рельефа

GeoScan32.

08.02.2016 13:41

Слайд 15

Интерпретация результатов 08.02.2016 13:41

Интерпретация результатов

08.02.2016 13:41

Слайд 16

«Радаростратиграфия» 08.02.2016 13:41

«Радаростратиграфия»

08.02.2016 13:41

Слайд 17

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

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

Слайд 18

Глубинность и разрешающая способность Различия в физической природе полей в георадиолокации

Глубинность и разрешающая способность
Различия в физической природе полей в георадиолокации и

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

По глубинности методы перекрываются в области примерно 3-15 метров, По разрешающей

По глубинности методы перекрываются в области примерно 3-15 метров,
По разрешающей

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

2. Скорость распространения волн. Волновое уравнение для диэлектрика V = c/ε1/2

2. Скорость распространения волн.

Волновое уравнение для диэлектрика

V = c/ε1/2

µабс= µ0

, µотн=1
Токов проводимости «нет», т.е. они не учитываются
Такое допущение возможно при малости потерь в среде (которые зависят от токов проводимости)
Слайд 21

В сейсморазведке скорость распространения волн определяется упругими модулями среды – плотностью-ρ,

В сейсморазведке скорость распространения волн определяется упругими модулями среды – плотностью-ρ,

модулем Юнга- Е и модулем сдвига - µ,
отражательная способности границ определяются (в основном) контрастом акустической жесткости –V*ρ.
Контраст акустической жесткости и линейные размеры поверхности локального объекта, отнесенные к длине волны, определяют его способность к образованию дифрагированных волн *

Сформулировать фразу * для георадиолокации

Слайд 22

3. Причины наличия в среде контрастов акустической жесткости и диэлектрической проницаемости.

3. Причины наличия в среде контрастов акустической жесткости и диэлектрической проницаемости.

Значения

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

4. Контрастность. 5. Различия в кинематических и динамических особенностях волновых полей.

4. Контрастность.
5. Различия в кинематических и динамических особенностях волновых полей.

В с/р

– различные типы и классы волн, разные методики, в гр – отражённые + ИП – в одной точке
Слайд 24

6. Чувствительность волновых методов к одним и тем же изменениям свойств

6. Чувствительность волновых методов к одним и тем же изменениям свойств

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

§7. Технология полевых работ волновыми методами. Принципиально возможно проведение гр работ

§7. Технология полевых работ волновыми методами.

Принципиально возможно проведение гр работ

методом ОГТ, но это сильно уменьшит производительность, не увеличив количество получаемой информации.

Измерения в режиме пешеходной съемки. Антенны в виде лыжи скользят по поверхности.
Частота посылки сигнала от единиц до тысяч раз в секунду.
Плотность наблюдений от 1 до 100 трасс на 1 погонный метр профиля Производительность и объем информации на 1-3 порядка выше, чем при производстве сейсморазведочных работ.
С/р – необходим контакт источника с поверхностью.
Г/р – можно поднять антенну – обычно при скоростной съёмке.
Исключение – с/р на акваториях.
При картировании протяженной границы на глубине порядка 8-10 метров при сходных методиках наблюдения – совмещенные источник и приемник или находящиеся на фиксированном расстоянии друг от друга (методики «t0» или «optimum offset»).
За рабочую смену в 8 часов можно отработать 10-15 км георадиолокационного профиля или 500-600 метров сейсмического профиля.

Слайд 26

Полевые работы Георадиолокация: Самый мобильный метод Простое проведение полевых работ Большая

Полевые работы

Георадиолокация:
Самый мобильный метод
Простое проведение полевых работ
Большая производительность
Сейсморазведка:
Требует определённого состояния поверхности
Минимум

механических источников шума
Самый трудоёмкий метод гф
Слайд 27

8. Комплексирование волновых методов и роль каждого из методов при комплексировании.

8. Комплексирование волновых методов и роль каждого из методов при комплексировании.
2

метода - независимые сведения о строении массива
Разные свойства – ЭФ и ЭМ
СР и ГР почти не перекрываются по глубинности, самая «глубикая» ось синфазности ГР может быть привязана к самой «мелкой» на сейсмограмме.
Слайд 28

9. Круг задач, решаемый с помощью георадиолокации . 9.1. Геологические, инженерно-геологические

9. Круг задач, решаемый с помощью георадиолокации .
9.1. Геологические, инженерно-геологические и

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

+ археологические и поисковые задачи

Слайд 29

9.2. Поиск локальных объектов, обследование инженерных сооружений, нарушения штатной ситуации а)

9.2. Поиск локальных объектов, обследование инженерных сооружений, нарушения штатной ситуации
а) трубопроводов;
б)

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

10.Примеры волновых полей и результаты решения некоторых задач. Георадар PulseEkko, 100

10.Примеры волновых полей и результаты решения некоторых задач.

Георадар PulseEkko, 100 МГц.

Профилирование дюн (а,б) в Южной Африке
Картировочные признаки УГВ.
Слайд 31

Георадарное картирование озерных отложений. Канада, PulseEkko-IV, 50 МГц

Георадарное картирование озерных отложений. Канада, PulseEkko-IV, 50 МГц

Слайд 32

Мёрзлое и талое, разные частоты Аллювиальные отложения, мелководье Маккензи

Мёрзлое и талое, разные частоты

Аллювиальные отложения, мелководье Маккензи

Слайд 33

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

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

валунов скальных пород – точек дифракции. Георадар PulseEkko-IV, антенна 200 МГц
Слайд 34

Корреляция между фотографией обнажения (а), георадиолокационным временным разрезом (b) с привязкой

Корреляция между фотографией обнажения (а), георадиолокационным временным разрезом (b) с

привязкой «время-глубина» и реконструкцией границ раздела в слоистой толще
Слайд 35

Схема проведения межскважинных георадиолокационных просвечиваний.

Схема проведения межскважинных георадиолокационных просвечиваний.

Слайд 36

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

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

8,9 м в мерзлых аллювиальных отложениях.
Слайд 37

Результаты межскважинного просвечивания с антенной 50МГц при разных расстояниях между скважинами в мерзлых аллювиальных отложениях.

Результаты межскважинного просвечивания с антенной 50МГц при разных расстояниях между скважинами

в мерзлых аллювиальных отложениях.
Слайд 38

Результаты межскважинного просвечивания с антенной 100МГц при разных расстояниях между скважинами в мерзлых аллювиальных отложениях

Результаты межскважинного просвечивания с антенной 100МГц при разных расстояниях между

скважинами в мерзлых аллювиальных отложениях
Слайд 39

Преломленная волна

Преломленная волна

Слайд 40

Обследование бетонных конструкций Старая плита Новая плита

Обследование бетонных конструкций

Старая плита

Новая плита