Моделирование и определение основных свойств волны Лява

Оглавление

Введение

Теоретическое обоснование волн Лява

Расчётная часть

Список литературы

Введение

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

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

Упругие волны, возбуждаемые упругим сейсмическим источником, находящимся на поверхности однородного упругого изотропного полупространства с упругими модулями λ, μ, ρ, впервые изучены Лэмбом. Из полученных им результатов вытекает, что в такой среде появляются поверхностные волны. Поверхностная волна Лява одна из таких волн.

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

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

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

Теоретическое обоснование волн Лява

Рассмотрим среду, состоящую из однородного полупространства, перекрытого однородным слоем мощности h. Кровля этого слоя является свободной границей. Сделаем естественное предположение, что скорость распространения поперечных волн в слое VS1 меньше скорости распространения поперечных волн VS2 в полупространстве. Предположим, что в таких условиях проводятся сейсмические работы с использованием поперечных волн типа SН. В этом случае наряду с объемными волнами SН могут существовать родственные им поверхностные волны типа Лява. Начало декартовой системы координат поместим на дневной поверхности. Ось 0Z направим вертикально вниз, ось ОХ - по направлению распространения сейсмических волн. В соответствии с предположением о структуре сейсмических волн вектор смещения в искомой волне должен иметь вид в слое:

(1)

в полупространстве:

, (2)

где ω - круговая частота; VL - скорость распространения искомой волны Лява. U и W компоненты вектора смещения, которые вводятся следующим образом:

и

Как было показано ранее, вектор смещения и все его компоненты должны удовлетворять волновому уравнению вида:

. (3)

Подставляя в это уравнение соответственно компоненты V1 (х, z, t) и V2 (х, у,t), получим, что искомые функции А (z) и В (z) должны являться решением следующих дифференциальных уравнений:

; ; (4)

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

. (5)

При этих предположениях в качестве решений, обеспечивающих убывание амплитуды колебаний в волне Лява в бесконечности (по оси OZ), следует взять:

; (6)

, (7)

где А01, А02, В01 - константы интегрирования.

Граничные условия на свободной поверхности требуют обращения в нуль всех компонент нормального вектора напряжения Рz. В силу сделанного предположения о структуре компонент вектора смещения, это условие эквивалентно следующему:

. (8)

На границе слоя и полупространства должна выполняться непрерывность компонент вектора напряжения и вектора смещения. Эти условия, как легко видеть, трансформируются в следующие два уравнения:

; (9)

. (10)

С учетом (5) компоненты смещения могут быть представлены в виде:

; (11)

, (12)

где ; .

Подставляя выражение для смещения в первом слое в граничное условие (8), получим: А01 = 0. С учетом этого результата из условия (9) легко получить уравнение

. (13)

Выполнив третье краевое условие (10), получим уравнение

. (14)

Рассматривая (13) и (14) как систему двух линейных уравнений относительно двух неизвестных А02 и B01, можно сделать заключение, что нетривиальное решение этой системы будет возможно только в случае равенства нулю ее главного определителя

. (15)

Это уравнение целесообразно представить в виде:

, (16)

где ; .

Если это уравнение записать в другом виде и учесть многозначность функции arctg, то основное характеристическое уравнение (16) решаемой задачи можно представить в виде:

, (17)

где k - целое число, определяющее порядок (номер) гармоники волны Лява. Значение k = 0 соответствует первой, основной гармонике волн Лява. Полученное характеристическое уравнение является трансцендентным. Оно определяет бесчисленное множество зависимостей (VL (ω)) k искомой скорости волн Лява от параметров разреза и частоты.

 Скачать реферат