Файл: Ямщиков В.С. Геоакустика. Раздел Упругие волны в неоднородном массиве [учеб. пособие].pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 19.07.2024

Просмотров: 137

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

среду, параметры которой отмочены индексом ( I ) , а другая (неоднородности) - однородные шарики с радиусом ѵР , распо­ ложенные в узлах пространственной решетки с периодом S , параметры которых отмечены индексом (2) (рис. 8 ). На неодно­ родную область из однородной, состоящей из вмещающей орѳды, падает плоская монохроматическая продольная волна перпенди­

кулярно к границе раздела.

Поле смещений в неоднородной срѳ-

де будет

значительно отличаться

 

от поля смещений в двухслойной

 

однородной среде ввиду много­

 

кратного

рассеяния каждой

из

 

сфер и интерференции рассеянных

 

и прошедшей волн. Выделим в по­

 

ле смещения неоднородной среды

 

только ту

часть, которая пред-

Рио. 8. ІІодѳль микроиѳод-

отавляет

ообой плоскую продолъ-

нородной среды

нул волну, Эта волна описывается параметрами, характеризую­ щими среду, в которой она распространяется. Если сравнить ее о плоской продольной отраженной волной в двухслойной од­ нородной среде, можно определить эффективные "однородные* параметры неоднородной двухкомпонентной среды.

Определение рассеянного поля для впиничной неодг.ородно_стих Продольную волну, падающую на границу раздеда, опреде­

лим потенциалом

 

 

.с о

 

 

 

У>= у>

е ‘ С/"

,

где cp t - окорость

продольных волн в среде ( I ) ,

Поле вблизи частицы, расположенной в любой точке неод­

нородной среды (

х 0 ,

У' ,

z o ), складывается из суммы сме­

щений падающего

.юля и суммы полей,

рассеянных всеми осталь­

ными частицаги. Причем в си,^ симметрии задачи полное вовбуждающѳѳ поле будет зависеть лишь от координаты s B .

Будем считать, что радиус ча..тиц неоднородностей мал в сравнении с длиной волны во вмещающей среде и мг .юм решетки.

33


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

кодѳблющимоя в

направлении

оси z .

Продольный и поперечный

потенциалы поля смещений, рассеянного сейсмическим дрполѳм,

имеют вид:

 

 

 

 

 

 

 

 

 

а

>-i c/*r г

 

 

 

 

У=А 9 е

 

c o s & —Л a z

г

 

>

(П.І)

 

 

 

 

 

а г

 

Х-ЛЛ*

 

 

 

а

 

 

г

 

 

 

Ч> = В Л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s ‘n

&= B hs>

 

 

 

 

где

J°,

 

&

-

полярные

координаты;

г , 8 , С2

-

сферические

координаты;

со/

и

cs

-

соответственно

окорооть

продольных

и поперечных волн во вмещающей среде ( I ) .

 

 

 

 

 

Поле

смещений

П

известным

образом

зависит

от

потенци­

алов

смещений

У

и

ф

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и >=g z a a l

 

w

i

<р,

 

 

 

(П.2)

причем потенциал

имеет

только

одну составляющую

вдоль

оси

Q

,

а

У

 

і.е зависит

от угла

Q .

 

 

 

 

 

Поле, рассеянное оѳиоыическин диполем, должно удовлетво­

ряя,

волновым уравнениям:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и)

у = О

 

 

 

 

 

 

 

 

 

 

 

/

Л)

ср<

 

 

)

у> — О .

 

 

(П.з)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для полного описания раосѳ: иного неоднородностью поля не­

обходимо

найти

коэффициенты

А

и

5

і (П .І),

которые опреде­

ляются из граничных условий,

заключающихся:

 

 

 

 

I) в равенстве смещений (касательных

L/s

и нормальных й х )

на границе

неоднородности;

 

 

6 г

 

 

 

 

гг

напряже­

 

2)

в

равенстве

нормальных

и касательных

ний на площадке, пѳрпѳндикулпрной

радиус-вектору

*<?

 

^ .

 

 

Поле внутри шарика с учетом конечности потенциалов смеще­

ний

в его

центре

будет

иметь

вил:

 

 

 

 

 

 


 

— л '

д

' - oV

 

 

 

 

У=А Ѳ 2

S in

-рг~r Z

 

 

 

 

 

 

 

 

 

 

 

/

â>

 

 

 

(П.4)

 

 

S i n

~C~J

 

 

Ч > ~В д/>

 

 

 

 

 

Касательные

UB

 

и нормальные

Ux

смешения,а также

нормальные <оХі

и касательные

Ttg

напряже"ия выражаются

через потѳнциалы_оледующим образом;

 

 

 

 

ие '

/ З У ,_ L p l~ £ 2 > .

 

 

*

дѲ

 

э г

(П.5)

 

 

9 !f

/

9 у

COS Ѳ

 

 

V ,

 

 

 

 

~Т~ дѲ + s i n Ѳ

Т' - А&

где

со 1

■ ГЗ*У-L

3 /

 

0¥> ) .

c a s e

( Л )7.

9*У>

Z 9 У

“>

м

 

 

1 (П.6)

<?г

г & z

 

9Т& Ѳ ~ ~ г * д ѳ ~ с /

г

 

Ли у / - коэффициенты Ияме.

Граничные условия будут иметь вид:

V и„

 

<51

(П.7)

и , = ^

х = к

 

 

Выражая ив1 ,

а г{ , е гхС

, Zx e i

черев потенциалы с

помощью выражений (П .І),(1і.4)

и (П.5)

и подставляя их значе­

ния т> граничные условия (П.7) .можно получить систему из четы­

рех уравнений относительно неизвестных коэффициентов

А В ,

А ' , В \

 

 

 

 

неоднородностей /?

 

Если предположить,

что радиуо

иного

меньше,

чем длины продольных и поперечных воля во вмещающей

среде,

или

 

 

;

- р - к «

/

,

то,

реяая подучен­

ную сиотѳму уравнений относительно А и

В

,

можно получить

следующие выражения:

 

 

^

 

 

 

 

 

 

 

ы г / _

.

у э у

 

 

tf);

 

 

 

 

 

 

 

 

 

 

(П.8)

 

со

 

 

 

О*+JV, f ) ,

 

 

 

 

 

 

 

 

 

-

сложные

где ^ = —

a*.

Г“ т~

R »

P <Qc> Q*

 

4

 

‘ /V

 

 

-

сложные

 

 

 

 

 

 

 

 

 

 

 


функции,

зависящие от

параметров

сред (I) и (2) и координат

выражения для которых

иокно найти

в работе / 97.

§ 2.

Эффективные

параметры микронеоднородной среды

Для нахождения аффективных параметров неоднородной среды определим постоянную распространения и коэффициент от­ ражения для рассматриваемой неоднородной среды.

Рассмотрим смещение в точке ( ~яа , і/0 , Z0 ) . Это смещение окладываѳтся из смещения падающей волны о потенци­ алом

У= У <?"" * £ ^ (П.9)

и смещений, получаемых от рассеяния всеми сферами за исклю­

чением находящейся в точке (я?0

, ув , z a

) . Суммарное смеще­

ние ut

в

точке

( х в , уо , z o ) вследствие

симметрии зависит

только

от

z o

и равно

 

 

 

 

 

иж~ a z B

'

 

Смещение в поле, рассеянном кавдой неоднородностью, при учете симметрии относительно азимутального угла опреде­

ляется следующим обраэом:

 

^

 

/нгс

9

 

 

10)

иг ~ А а ( г - z0)*

 

 

 

где

X= |/(JE -Хо) г+- (у-У.) +(Z

го) >

 

 

 

 

/ >~ ]j(jc - Я . ) + ( у-ІГоУ.

 

 

 

Полное смещение ut —

равно:

 

 

 

 

 

(П .Н )

где знаі JET означает, что

суммирование

производится

по

всем

частицам эа исключением

находящейся в точке

, у , лоХ

%


 

Подотавляя выражения

(П .8),

(П .9),

(П.ІО) в уравнение

(П .ІІ),

 

 

после некоторых упрощений С 9 J

можно получить интегральное

 

 

уравнение для

потенциала Ѵ’

 

 

я

 

 

 

 

 

 

9Z.

 

 

 

 

 

 

 

 

 

 

 

 

 

d r ­

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

T

L cp/" / (®'дг

 

Щ Ra r * ' 9Z* e

dS+-3d(Pée y

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3Z_ *-M / f 9SX'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(IJI)

 

 

где

_/

 

4 Я R 5

 

-

отношение

объема

частицы к

объему

 

 

a

= —j y

s

 

 

 

окоужающей среды;

V

= n s

%

п

-

число

периодов

S

от плос­

 

 

кости частицы

( х 0 , уо ,

z 0 )

до

каждой другой частицы.-

 

 

 

 

Уравнение (П .І2) определяет

 

полное

поле в точке (гс„,

 

 

уо ,

z g ) .

Первый член уравнения (П .І2) соответствует

падаю­

 

 

щей неискаженной волне,

остальные

 

члены учитывают влияние

 

 

частицы,

расположенной в точке z

= г с , и определяют отражен­

 

ную волну.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегральное уравнение (П .І2) после замены дифференци­

 

рования

по

z o

дифференцированием по V"

и после интегриро­

 

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

 

оператора:

 

 

 

. _

S’*

сОг

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

* г о

с *,

 

'

 

 

 

 

 

 

В

результате

этих

операций интегральное уравнение

 

(П.І2) можно привести к виду дифференциального уравнения от­

 

 

 

У * d [ 0

 

Ң - л/

 

 

.( М +N, ?...) ] - О..•*/. (П.ІЗ)

 

носительноѵ і я

¥Jr

 

 

.

 

g

 

Г

 

 

 

 

л

j ^ I - d (Зр ѵ-Qf

 

 

 

 

Щ * { э М Щ ^ 2 Л ? У } *

 

 

Уравнение (П .ІЗ) имеет решение следующего вида:

У = В ¥п

(П .І4)

37