Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Геологические основы

.pdf
Скачиваний:
10
Добавлен:
17.01.2021
Размер:
2.62 Mб
Скачать

СПБГУАП группа 4736 https://new.guap.ru

Полученные выражения сравним с функциями, описывающими смещения

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

R , а в поперечной, что описывается слагаемыми Uϕ и Uθ , ортогональны R и

направлены по касательным к вертикальным большим кругам θ.

Амплитуды колебаний на всех компонентах зависят от направлений θ и ϕ. Так амплитуда продольной волны максимальна при θ=π/2 сила, действующая вдоль осей Х и Y и при θ = 0° для силы направленной вдоль оси Z. Максимумы поперечных составляющих Uθ наблюдаются при углах θ = 0° для сил направленных по осям Х и Y и θ =π/2 при силе направленной по оси Z.

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

Всамом деле, амплитуды продольной компоненты и поперечной Uθ

максимальны при угле θ = 0°, если сила принадлежит лучевой плоскости и направлена по оси Х. В случае силы направленной по оси Y, амплитуды UR и Uθ

максимальны при ϕ =π/2, т.е. в направлениях перпендикулярных профилю наблюдений. Сила, направленная по оси Z продуцирует в среде симметричное волновое поле, которое не зависит от азимутального угла ϕ. С другой стороны, сравнивая эти выражения, а мы замечаем, что зависимость амплитуд от расстояния R полностью соответствует случаю сферических волн от точечных центров расширения и вращения.

3.2. Двойная сила без момента.

Пусть теперь одна из сосредоточенных сил приложена по оси Х, а другая, равная ей по величине, но противоположного направления в точке оси Х на расстоянии х от первой. При условии дальней зоны, когда х<<λ<<R и

77

СПБГУАП группа 4736 https://new.guap.ru

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

Пара сил, действующих по оси Х:

URxx

=

 

 

 

 

 

 

 

F0

 

 

 

 

 

 

x sin2 θ cos2 ϕT (t τ p );

 

4πρCp2 R

 

 

 

 

 

 

λp

Uθxx

=

 

 

 

 

 

 

F0

 

 

 

 

 

x sinθ cosθ cos2 ϕT (t τs );

4πρCs2 R

 

 

 

 

 

λs

Uϕxx

= −

 

 

 

 

F0

 

 

 

 

 

 

x sinθ sinϕ cosϕT (t τs )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πρCs2 R λs

 

Пара сил, приложенных по оси Y:

URyy

=

 

 

 

 

 

 

 

F0

 

 

 

 

 

y sin2 θ sin2 ϕT (t τ p );

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πρCp2 R λp

Uϕy y =

 

 

 

 

 

 

F0

 

 

 

 

 

y sinθ sinϕ cosϕT (t τs );

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πρCs2 R λs

Uθxy

= −

 

 

 

F0

 

 

 

 

 

 

y sinθ cosθ sin2 ϕT (t τs ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πCs2 ρR λs

 

Пара сил, приложенная по оси Z:

URzz

 

 

=

 

 

F0

 

 

 

 

 

z cos2 θ T (t τ p );

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πρCp2 R λp

Uϕz z

 

= 0;

 

 

 

 

 

 

 

 

Uθz z

 

= −

 

F0

 

 

z sinθ cosθ T (t τs ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πρCs2 R λs

Рассмотренные выше

 

 

комбинации из двух равных и взаимно

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

В отличие от полей амплитуд простых направленных сил в этих случаях возникает дополнительная зависимость от отношения плеча диполя ( x, y, z) к

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

78

СПБГУАП группа 4736 https://new.guap.ru

Z, амплитуда пропорциональна cos2θ и положительна при всех θ от 0 до π . Это означает, что смещение всегда происходит от начала координат в сторону возрастающих значений R.

Поле амплитуд поперечных волн содержит ту же зависимость отношения плеча диполя к длине излучаемой поперечной волны, что и продольные волны. Но в этом случае накладываются более существенные ограничения по угловым координатам. Так амплитуды поперечных колебаний пропорциональны множителю sin2θ, это означает, что как в вертикальном (ось Z), так и в горизонтальном (ось Х) направлениях величина Uθ = 0. Она достигает максимума при θ = π/4 (45°).

Зависимость амплитуды поперечных волн от азимутального угла ϕ также определяется множителем sin2ϕ . Причем она равна 0 вдоль осей X и Y и

достигает максимальных значений при ϕ = π/4 ± kπ/2, k = 0,1,2,3. Заметим, что в случае парных сил одновременное изменение направления действия их на взаимообратное, то есть приращение углов θ и ϕ на π, не приводит к перемене полярности поля амплитуд поперечных и продольных волн, поскольку функции направленности L(θ,ϕ) не меняют свои знаки.

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

 

Две силы без момента, приложенные на

Простая сила: Lx(θ,ϕ)

расстоянии Х друг от друга: Lxx(θ,

79

СПБГУАП группа 4736 https://new.guap.ru

URx =

 

 

F0

sinθ cosϕT (t τ p )URxx

=

 

Fo

x sinθ sinθ cosϕ cosϕT (t τ p );

 

4πρCp2 R

 

 

 

 

 

 

 

 

 

 

4πρCp2 R λp

Uϕ x = −

F0

 

sinϕT (t τs )Uϕx x = −

 

 

F0

 

 

x sinθ sinϕ cosϕT (t τs );

4πρCs2 R

 

 

 

 

 

 

 

 

4πρCs2 R λs

Uθx =

 

 

F0

 

cosθ cosϕT (t τs )Uθxx

=

 

F0

 

x sinθ cosθ cosϕ cosϕT (t τs ).

4πρCs2 R

 

 

 

 

 

 

 

4πρCs2 R λs

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

URxx

=URx

 

x LxRx (θ,ϕ);

LxRx (θ,ϕ)= sinθ cosϕ.

 

 

 

 

λ

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

U x

=U

ϕ x λ

x Lx

 

(θ,ϕ);

Lx

(θ,ϕ)= sinθ cosϕ.

ϕ x

 

 

ϕ x

 

 

ϕ x

 

 

 

 

 

 

 

s

 

 

 

 

 

U x

=U

θx

λ

x Lx

(θ,ϕ);

Lx

(θ,ϕ)= sinθ cosϕ.

θx

 

s

θx

 

 

θx

 

 

 

 

 

 

 

 

 

 

 

 

 

То же самое для двойных сил по оси Y и Z:

URYY

=URY

 

 

 

y LYRY (θ,ϕ);

LYRY (θ,ϕ)= sinθ sinϕ;

 

 

 

 

λ

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U Y

=U

ϕY

 

λ

y LY

(θ,ϕ);

LY

(θ,ϕ)= sinθ sinϕ;

ϕY

 

 

s

ϕY

 

ϕY

 

 

 

 

 

 

 

 

 

 

 

 

 

U Y

=U

θY

λ

y LY

(θ,ϕ);

LY

(θ,ϕ)= sinθ sinϕ;

θY

 

s

 

 

θY

 

θY

 

 

 

 

 

 

 

 

 

 

 

 

 

 

URZZ

=URZ

 

 

 

Z

LZRZ (θ,ϕ);

LZRZ (θ,ϕ)= cosθ;

 

λ

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U Z

= 0;

 

 

 

 

 

 

 

 

LZ

 

(θ,ϕ)= 0;

ϕZ

 

 

 

 

 

 

 

 

 

 

ϕZ

 

U Z

=U

θZ

 

Z

LZ

(θ,ϕ);

LZ

(θ,ϕ)= cosθ;

λ

 

 

 

 

θZ

 

s

θZ

 

θZ

 

 

 

 

 

 

 

 

 

 

 

 

3.3. Двойные силы с моментом.

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

80

СПБГУАП группа 4736 https://new.guap.ru

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

r r

M x = [Z0 Mr y = [Xr0 Mr z = [Xr0

yr]= [Y0 zr]= [Zr0 yr]= [Yr0

zr]= Z0

y sin(Z0 , y)=Y0 z sin(Y0 , z);

xr]= X 0

z sin(X 0 , z)= Z0

x sin(Z0 x);

xr]= X 0

y sin(X 0 y)=Y0

x sin(Y0 , x).

Формулы для функций смещения во фронте упругой волны генерируемой источниками указанного типа имеют вид:

Момент относительно оси Х:

M x

= Z0

 

y sin(Z0 , y);

 

 

 

 

 

 

UMRx

=

 

 

 

Z0

 

 

2π

y

sinθ cosθ sinϕTp (t τ p );

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πρCp2 R λp

 

 

 

 

 

 

U x

= −

 

Z0

 

 

 

 

 

2π y

cosθ sinϕT

(t τ

s

);

 

 

 

 

 

 

 

 

Mϕ

 

 

 

 

4πρCs2 R

λs

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

UMx θ

= 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Момент относительно оси Y.

 

 

 

M y

= Z0

 

x sin(Z0 x);

 

 

 

 

 

 

UMRy

=

 

 

 

Z0

 

 

2π x

sinθ cosθ cosϕTp (t τ p );

 

4πρCp2 R

 

 

 

 

 

 

 

 

 

 

λp

 

 

 

 

 

 

UMy ϕ

= 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U y

= −

 

Z0

 

 

 

 

 

2π

x sin2 θ cosϕT (t τ

s

).

 

 

 

 

 

 

 

Mθ

 

 

 

4πρCs2 R λs

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Момент относительно оси Z.

 

 

 

 

M Z

 

= x0

 

y sin(x0 y);

 

 

 

 

 

UMRz

= 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

UMz ϕ

= −

 

 

 

 

x0

 

 

 

2π y sinθ(t τs );

 

 

4πρCs2 R

 

 

 

 

 

 

 

 

λs

 

 

 

 

 

UMz θ

= 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

81

СПБГУАП группа 4736 https://new.guap.ru

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

Первая из них (А,ρ,R,Сp,s) состоит из комбинации характеристик механических сил, развиваемых источником, и упругой среды, в которой производится возмущение ее первоначального состояния. При этом эта функция, в сущности, не что иное, как амплитудный множитель наиболее простого источника типа центра расширения для продольных и центра вращения для поперечных волн. Главное отличие заключено в дополнительном множителе вида r/λp,s, учитывающем меру пространственной распределенности той или иной комбинации пары сил по сравнению с длиной генерируемой продольной или поперечной волн.

Вторая функция –LM(θ,ϕ) определяет ориентацию момента пары сил относительно осей выбранной системы координат. В отличие от других комбинированных источников, функция направленности II рода определяет направления, в которых излучение продольных, либо поперечных волн не происходит совсем. Так для момента относительно оси Х, помимо значений углов

θ = 0 и π/2, где излучение продольных и поперечных волн отсутствует, этот источник не генерирует поперечные волны первого типа, т.е. те, направления колебаний частиц которых перпендикулярны волновому вектору krs1 и лежат в лучевой плоскости.

Источник типа пары сил с моментом относительно оси Y напротив излучает только поперечную волну S1 и не генерирует полностью волну типа S2. При этом направление вертикально вниз, т.е. θ = 0°, характеризуется нулевой интенсивностью излучения и волны S1.

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

82

СПБГУАП группа 4736 https://new.guap.ru

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

Момент пары сил относительно оси Х:

U x

=U

RZ

2π y

 

Lx

(θ,ϕ)T

p

Lx

(θ,ϕ)= sinθ sinϕ;

 

λ

 

 

MR

 

 

p

 

MR

 

MR

 

 

 

 

 

 

 

 

 

 

 

 

U x

=U

θ Z

 

2π y

 

Lx

(θ,ϕ)T Lx

(θ,ϕ)= cosθ sinϕ;

 

λ

 

 

Mϕ

 

 

s

 

Mϕ

 

s

Mϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

Момент пары относительно оси Y:

UMRy =URZ 2πλp x LyMR (θ,ϕ)Tp LyMR (θ,ϕ)= sinθ cosϕ;

U y

=U

θ Z

2π x

Ly

(θ,ϕ)T Ly

(θ,ϕ)= sinθ cosϕ.

 

 

 

Mθ

 

 

λ

s

 

 

Mθ

 

s

Mθ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Момент пары сил относительно Z:

 

 

 

U z

=U

 

 

2π

y

Lz

 

(θ,ϕ)T Lz

 

(θ,ϕ)= sin

θ

1

.

 

λs

 

 

 

sinϕ

Mϕ

 

 

ϕ x

 

 

 

Mϕ

s

Mϕ

 

 

 

Третья функция Tp,s(t-τp,s), как и раньше, представляет собой временную зависимость колебательного процесса в упругой волне и полностью определяется механическим воздействием, размером источника и упругими характеристиками среды Сs, Сp. Кроме того, сюда входят все фазовые особенности излучения волн,

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

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

83

СПБГУАП группа 4736 https://new.guap.ru

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

Рис. 13 иллюстрирует поведение функции L(θ,ϕ) для простейшего типа воздействий.

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

84

СПБГУАП группа 4736 https://new.guap.ru

а — Z-сила, сечение в плоскости XOZ; б — X-сила, трехмерная функция для S и Р

– волн. Для построения этой же функции для Y силы достаточно развернуть диаграмму X на 90°.

3.4. Точечные источники на свободной поверхности полупространства.

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

Горизонтальная сила вдоль оси X

LRx (θ,ϕ)= L0Rx

 

(θ,ϕ)

 

 

 

cosθ 1γ 2 sin2 θ

 

 

 

 

 

;

 

(12γ 2 sin2

θ)2 + 2γ 3 sinθ sin 2θ

1γ

 

 

 

 

 

 

 

 

 

 

 

 

2 sin2 θ

L

 

(θ

,ϕ)= L0

 

 

1;

 

 

 

 

 

 

 

 

 

 

 

 

ϕ x

 

 

ϕ x

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

(θ

,ϕ)= L0

 

 

(θ,ϕ)

 

 

 

cos2θ

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θ x

 

 

θ x

 

 

 

 

 

cos2 2θ + 2sinθ sin 2θ γ 2 sin2

θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Горизонтальная сила вдоль оси Y:

 

 

 

 

 

 

 

LRy (θ,ϕ)= L0Ry

 

 

 

cosθ 1γ 2 sin2 θ

 

 

 

 

;

 

(12γ 2 sin2 θ)2 +2γ 3 sinθ sin 2θ 1γ 2 sin2 θ

 

L

 

(θ,ϕ)=

L0

(θ,ϕ) 1;

 

 

 

 

 

 

 

 

 

ϕ y

 

 

 

ϕ y

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

(θ,ϕ)=

L0

(θ,ϕ)

 

cos2θ

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

θ y

 

 

 

θ y

 

 

 

 

cos2 2θ +2sinθ sin 2θ γ 2 sin2 θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вертикальная сила вдоль оси Z:

85

СПБГУАП группа 4736 https://new.guap.ru

 

0

 

 

 

1

2γ 2 sin2 θ

 

 

 

L

(θ,ϕ)= L

(θ,ϕ)

 

 

 

 

 

 

;

 

(12γ 2 sin2 θ)2 +2γ 3 sinθ sin 2θ

 

 

Rz

Rz

 

 

 

1γ 2 sin2 θ

 

 

 

 

 

L

(θ,ϕ)= L0

 

(θ,ϕ)

= 0;

 

 

 

 

ϕ z

ϕ z

 

 

 

 

 

 

 

 

L

(θ,ϕ)= L0

 

(θ,ϕ)

 

cosθ γ 2

sin2 θ

 

.

 

 

 

 

 

 

 

θ z

θ z

 

 

 

 

cos2 2θ + 2sinθ sin 2θ γ 2 sin2

θ

 

 

 

 

 

 

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

силы Z0 определяется параметром

γ =

Cs =

μ

, т.е. основными упругими

λ + 2μ

 

 

Cp

 

постоянными среды. Если в случае безграничного пространства углы θ с нулевой интенсивностью излучения упругих продольных волн для этих же самых источников равны 0 для Х0 и Y0 сил и π/2 для Z0 силы, то, располагая источники на границе полупространства, сюда добавляются углы θ равные π/2 и arcsin1/2+nπ

(n = 0,1,2…) для сил X0 и Y0 arcsin

2

+nπ (n = 0,1,2...) для Z0.

2γ

 

 

Для компоненты поперечных колебаний в волне S1, характеризующейся функцией направленности L0j(θ,ϕ), добавляется направление с углом θ = π/4 для сил Х0 и Y0, а для Z0 силы угол θ = π/2 и arcsin 1/γ + nπ (n = 0,1,2 …), при которых интенсивность поперечных колебаний равна 0.

Только для компоненты волнового поля соответствующей волне S2 с

функцией направленности Lϕj(θ,ϕ) интенсивность поперечных колебаний не зависит от того, где располагается источник возбуждения.

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

86