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

книги из ГПНТБ / Ямщиков В.С. Геоакустика. Раздел Упругие волны в неоднородном массиве [учеб. пособие]

.pdf
Скачиваний:
9
Добавлен:
23.10.2023
Размер:
4.82 Mб
Скачать

Продольные волны первого и второго рода имеют следую­ щую физическую сущность.

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

Поперечные волны в водойаоыщенной пористой (двухфазной) среда. Закономерности распространения поперечных волн в по­ ристых средах описываются системой уравнений (1 .2 7 ). Способ получения решений уравнений (1.27) в случае рассмотрения

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

иЧ>2 вида:

 

 

yj =

С exp £(СО

 

 

4>= B e x p i (< o t - A s x )

в систему

(1.27)

и приравнивания

нулю определителя получа­

ющейся при

этом с ’стемы

можно получить дисперсионное урав­

нение для

поперечных

волн [2]

. Решая это уравнение относи­

тельно неизвестного

параметра

^

,

зависящего от волнового

числа

ft?

, и учитывая

соотношение

 

 

 

л -

/

 

»

 

 

 

Ѵ 7 =

-Ц -

1 со

можно получить выражения для

скорости поперечной волны c s

и ее

коэффициента затухания

оС5 .

 

 

Как и в случае

продольных волн, поперечные волны в зна­

чительной степени зависят от частоты распространяющейся вол­ ны (или от сил фильтрационного сопротивления, которые учиты­ ваются в па.аметре X ).

чреде льном случае при очень малых частотах или очень больших силах фильтрационного сопротивления (ма^ая проницаѳ-

20

мость среды или большая вязкооть жидкости

о ) поперечная

волна распространяется без дисперсия а без

ватухания,

причем

скорость распространения ее

такая se, как

в идеально

упругой

среде, в которой коэффициент

поперечного сдвига равви

 

а плотность

равна

(I

гп а/>г , При этой вся жидкость

колеблется

выѳстѳ

со скелетон среды.

 

 

 

 

 

 

 

(1.46)

В случае волн высоких частот или же очень налах ^ил фильтрационного сопротивления скорость поперечных волн в пре­ деле ( — 00 ) такая же, как и в идеально упругой среде о тем же коэффициентом поперечного сдвига, но плотнооть которой равна (I

(1.47)

И в атом случае сдвиговые водны распространяются без дио~ Персии и затухания, но здесь жидкость как бы вообще не участ­ вует в колебаниях скелета, т .ѳ . поперечные волны на высоких частотах распространяются в насыщенных пористых орѳдах так, как будто они насыщены газом при очень низком давлении,

Наоыщениѳ пористой среды жидкостью

 

приводит к ревко­

му возраста,

ію

величины

Cs

. Этот вывод

подтверждается на

 

ез

 

практике: скорооть

cs

не претерпевает анкетного скачка на

уровне грунтовых вод, как в

слу.ае продольных волн.

Коэффициент затухания

поперечных волн

оСа при фиксиро­

ванной частоте

 

при

очень высоких или очень малых сипах

фильтрационного сопротивления стремится в нулю, а при некото­

рых промежуточных значениях имеет максимум,

как в случае про­

дольных волн.

При переменной частоте

со

фиксированных па­

раметрах спады

ocs — о при го о

в

 

Величина ватухания попѳрѳчнь' во.л невелика по сравне­ нию с. затуханием объемных вол" второго рода.

21

§ 3» Упрощенный метод оіціедедания параметров упругих волн в пористых средах

Упрощенный метод расчета в двухфазной среде. Двухфаз­ ная среда в приведенном выше случае рассматривалась' как статистически однородная среда, состоящая из упругого порис­ того скелета и жидкого или газообразного эаполнитѳлі/. Одна­ ко при использовании полученных теоретических выводов для практического расчета скоростей и затухания упругих волн возникают значительные трудности в связи с тем, что в рас­ чет входят некоторые трудно определимые экспериментальным путем константы, В связи с зтим часто для практических рас­ четов используют результаты, полученные при помощи моделей зернистых сред, упругие свойства которых при отсутствии за­ полнителя пор определяются условиями взаимодействия на кон­ тактах зерен. В частности, наиболее удобной для конкретного анализа является модель, в которой твердая компонента грунта представлена в виде кубической упаковки идеально упругих сфер одинакового радиуса, находящихся под действием собственного веса. Реальные грунты значительно отличаются от данной идеа­ лизированной модели, однако, как показывает опыт, использова­ ние этой модели позволяет достаточно просто объяснить основ­

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

Основной принцип расчета скоростей распространения упру­ гих волн в указанных моделях заключается в следующем. Ско­ рость, например, продольных во~ч можно представить по анало­ гии о однородной инотропной средой выражением

где ж

- упругий параметр двухкомпонентной зернистой среды,

который соответствует упругому

параметру

(Я + 2р ) Для изо­

тропной

однородной среды; j °

-

средняя

плотность среды,

 

р - (/-

) / \

.

22

Основная задача заключается в установлении параметра ;г , определяемого упругими свойствами самих верен, а такжо усло­ виями взаимодействия на контактах зерен.

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

 

 

зе ■

а?,

з е л

 

 

(I.W )

 

 

~Ѵ£

 

К

* з

ѵ <

 

 

 

 

где

эег - модуль

 

у

 

у

 

 

 

объемной упругости среды, заполняющей поры;

aef - модуль объемной упругости зеоен;

- объем среды, за­

полняющей поры;

^

-

объ'ѳн зерен;

зез

- упругооть

контактов;

V - суммарный

объем.

 

 

 

 

 

 

 

Упругий коэффициент зё

зависит

от

глубины погружения 2

в грунт рассматриваемой вѳрнистой системы и пропорционален

корню кубическому

из 2

,

а

также

зависит от пористости грун­

та / п

. При малой глубине

г

упругость

контактов мала, и

роль

основного

передатчика упругих

колебаний играет

объемная

упругость орѳды. С увеличением глубины упругость контактов возрастает, и член эе6 У начинает играть основную роль. Коэф­ фициент а определяется из условий ооприкооновения упругих шаров, находящихся под действием оилы Р , направление дей­ ствия которой совпадает с линией центров шаров. Скорость по­ перечных волн Cs в рассматриваемой зернистой модели опреде­ ляется совершенно аналогично. Скорость продольных волн в пѳоке, полностью насыщенном водой, можно получить в предположе­ нии, что жидкость, заполняющая поры, и зерна движутся при распространении волны виестѳ, не испытывая относительных дви­ жений, благодаря сравнительно мало равнящимся плотностям.

Рассмотрим модель среды, содержащей грунтовые воды (рис. 4 ). Пусть среда имеет зернистое строение (двухфазная).

Скелет среды будем аппроксимировать кубической упаковкой иде­

ально упругих сфер.

Поры до глубины г=- Ну

абсолютно сухи,

а н и яѳ ^ -Му) полностью заполнены водой.

Поверхность V ,

разделяющая сухую (

с ) и водонасыщенную ( â

) части среды,

соответствует уровню (зеркалу) грунтовых вод.

23

Скорости придельных и поперечных волн в сухой л ведонаоыщенной областях среды, найденные по приведенному выше мето­

ду, согласно работе Г 5 ! выражаются следующим образом. Для части пространства о абсолютно пустыми порами

скорое ли продольных

Сс и поперечных

Ст волн равны:

 

 

 

 

г е- \ !

â >7Q

V /Z H Z T

-

 

 

 

(1.50)

 

 

 

-

 

] /( У - ^ Ѵ г

 

 

 

 

 

 

„ с= , С И -

Z L і Д « * ( * - е Л [ £ * л г '

(1*51)

 

 

 

 

s

УЛ ('-

к Л ( ' - ” ) у U- е-Ѵг ’

 

 

 

 

 

 

где

Е

-

модуль Юнга;

 

сг - коэф­

 

 

 

 

фициент Пуассона; />* -

плотность;

 

 

 

 

G

-

эффективный модуль

сдвига

среды;

 

 

 

 

т - пористость, для кубичѳокой упа­

 

 

 

 

ковки

равная 0,*76.

 

 

 

 

 

 

 

 

 

 

Соотношение

этих скороотѳй

 

 

 

 

 

 

 

Л — _ £ ^ _

Л З

 

 

(1.52)

 

 

 

 

 

 

>

сc ;f

 

Ѵ Т = ^

 

 

 

 

 

 

зависит только от коэффициента Пуас­

 

 

 

 

сона.

 

Для пространства

с

полностью

іис. *. Модель двухфазной

 

 

заполненными

порами(л -> // у )

фор­

зернистой срѳдн,содержа­

 

щей грунтовые воды

мулы для скороотѳй продольных и

 

 

 

 

 

поперечных волн имеют следующий вид

 

Л[

і

+ J 7-ß\ f E l ( Я - А ) г +Лг Н „ ] '

 

 

 

 

 

У"/А*9?.... ’

У

 

«~<Уг)г _ .

 

(1.53)

 

 

 

ІМ <~' і ) ^ т

г г

 

 

 

 

 

 

 

е

( э л г d -

СГ)

Л Е

*

-

А

*

и

'

 

 

(1.5*)

г у

Л, (і~'™)+Лг ^

V

(/-&*)*

 

 

 

 

 

где

Лг

иодуль объемного

сжатия

зерен;

Л г

"

модуль

она-

тия

жидкости;

-

плотность

жидкости;

y>f ~

плотность

зе­

рен.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При

ыводѳ формул (1.53)

и (1.54)

было учтено

действие

выталкивающей силы жидкости яа частицы среды в по. j простран­

стве 8

Н

-

 

 

 

 

 

 

 

 

 

 

 

При переходе через границу і?

-

Л

скороить Ср

скач­

кообразно увеличивается,

что

можно заметить

из формулы (1 .53).

При z = M y благодаря

изотропным свойствам

модели

скорость

продольных волн Ср

равна

граничной скорооти

с г у

продольной

преломленной волны

,

распространяющейся вдоль

границы R y ,

Скоростную дифференциацию среды вблизи границы

для

про­

дольных воля

можно характеризовать

отношением

 

 

 

 

 

 

j + —£2-

Л

 

 

 

 

Cp (z

Л іху).= .

V■» ' /--гтгЛГТ

/*,

 

 

(1.55)

 

 

 

ч + п / Л х + О - ^ / Р г

 

 

 

С£

 

 

 

£ 7 0 W сгг)*

 

 

 

Скорость

,

согласно

принятоыу выше предполохѳнню,

скачкообразно

уменьшается

на

границе

/?у в

связи

о увеличе­

нием эффективной плотности. Так, для среды с пустыми порами

эффективная

плотность

j°c =

( / — m )/> t

, для влажной орѳдн

эффективная

плотность

 

( s —m )/> f

т /> я . Скоростная

дифференциация для поперечных волн при переходе через грани­

цу раздела

Ry

определяется

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

 

 

 

 

 

 

=

c s ( Z

= H v )-

V/yuJZ L . J X

 

 

 

(1.56)

 

 

c f ( г

= К J

 

 

 

Кроме

того,

при

переходе

через

границу

водораздела

R

скачкообразно

возрастает

коэффициент

к .

определяющий--------------- -

отно» -

вѳниѳ продольных

и поперечных скоростей.

 

 

 

 

 

h .

 

 

 

 

 

 

 

 

 

 

 

(1.57)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

op<#■ :

osf

Г

 

(it-c r * )*

 

 

 

 

 

Причем

зависят от

глубины залегания

грани­

цы водораздела

р

, т .е . от

величины

.

Так,

скорооть

продольной волны при любых

И

больше скорости звука

в

воде

и плавно увеличивается с увеличением

.

Увеличение

скорос­

ти продольных волн при пѳрехс ,е через границу

f-{

тем

больше,

чем иенын«

Ң и

. Граница

Р

 

не является

преломляющей гра-

25

ницѳй для поперечных воля, а скорость сл при переходе через границу Fy существенных изменений не претерпевает, но при переходе через границу происходит реэкое возрастание отноше­

ния скоростей продольных

и поперечных волн Л ,

причем это

возрастание

тем

больше,

чем меньше

 

исходя

из

описанных выше данных о скоростном

строении

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

При наблюдении на поверхности среды в первых вступлени­

ях до некоторого расстояния от источника

сс,

прослаивается

рефрагированная

продольная волна

£р

. Кажущаяся спирость

вдоль годографа

волны

£

по мере удаления от источника воз­

растает сначала быстро, а

потом

более плавно,

іі точке

яг ,

положение которой зависит

от глубины залегания

Л ѵ гоаницы

Ру .волна ±р сменяется продольной

преломленной волной ^.расп ро ­

страняющейся вдоль границы

 

. Причем

при

2

~ £ іу

благо­

даря изотропным свойствен модели ср сгу

, где

с гу

- гра­

ничная

скорость

волны

£у

. При малых Н у

благодаря

увеличе­

нию Ср

о глубиной годограф волны

 

£ у

криволинейный: каю ­

щаяся скорость плавно увеличивается с расстоянием.

При боль­

ших

(Ну

> 3 0 „ )

возрастание

ср

о глубиной происходит

медленно и не

сказывается

на форме

годографа волны

t

 

За

волнами

t p и

 

на сейсмограмме будет регистриро­

ваться

продольная волна, отраженная

от границы

F .

Эта

волна обладает большой интенсивностью как вблио« источника, так и на больших удалениях от него.

Отсутствие существенных изменений скорости поперечных

волн на границе F y

обуславливает малую интенсивность об -

манных отракѳнных и

преломленных волн,

связанных с

этой гра­

ницей.

 

F не оказывает

 

 

Наличие

границы

существенного

влияния

на характер

поверхностных волн, распространяющихся вдоль по­

верхности ореды. Скорость рѳлеевских волн очень слабо зависит от С при. постоянной скорости поперечных волн, Б силу этого дисперсия волны Рѳлѳя будет определяться только неоднородно­ стью среды.

26

Такое же малое влияние оказывает граница Fy и на волны Лява, возникающие при неосѳсиыметричном воздействии в среде, ско­ рость поперечных волн в которой возрастает с глуОиной.

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

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

Результаты расчета величины отношения скорости продоль­

ных волн во влажном песке

Сро

на поверхности

z= * o к ско­

рости

cz

в воде

в зависимости от

пористости

т . приведена

на рис. 5.

 

 

 

 

 

 

 

Как видно

из

графика, при

 

 

 

пористости т <

30% величина ср

 

 

может отличаться от скорости с,

 

 

более чем на 20%. При дальнейшем

 

 

изменении

т

скорость во влажном

 

 

песке мало изменяется. Таким об­

 

 

разом,

при достаточно высокой

по­

Рио. 5. Расчетные зави-

ристости

т .

при условии,

что

все

-имости скорости и коэф­

поры насыщены водой, величина ско­

фициента

затухания упру­

гих волн

в зернистой

рости в двухкомпонентной

среде

 

двухфазной среде от по­

близка

к скорости

в воде

и мало

ристости

 

 

изменяется

с изменением пористости. Отсюда можно сделать

вьн

вод

о том,

что по величинам скоростей практически нельзя

су­

дить

об

изменении пористости двуі :омпонѳнтной с гедн, так

как

различия

в

величинах скоростей находятся в пределах возмож­

ных ошибок

определения величины скорооти.

 

27

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

для коэффициента поглощения в двухкоыпонентной среде

[

Ъ ]х

где

£

-

вязкость

жидкости;

J? - радиус зерен; / '

-

час­

тота;

-

средняя

плотность

двухкомпонентной среды

при

Z = 0.

 

 

 

 

 

 

 

 

 

Нак видно ив (1 .5 8 ), коэффициент поглощения пропорцио­

нален квадрату

частоты

и радиуса верен.

 

 

 

Выражение

(1.58)

можно представить в виде

 

 

t (1.59)

Соглаоно (1.59) была рассчитана зависимость г;— - = ^(т ) представленная на рис. 5, из которой видно, что коэффициент поглощения более чувствителен к изменению пористости, чем окорость.

Кроне тоге, как видно из (1.58) и (1 .5 9 ), коэффициент поглощения сильно зависит от радиуса зерен двухфазной среды.

Распространение упругих волн в трехфавных средах. В за­ ключение рассмотрим случай трехкомпонентной среды, например, когда поры между зерен одинакового радиуса заполнены жидкостью и гааом.

Для нахождения скорости продольных волн в трехкомпонент­ ной среде используем описанный выше упрощенный метод расчета. Рассмотрим окорость распространения в трехкомпонентной среде,

когда

давление на рассматриваемую среду отсутствует,

т .е .ср е ­

да располагается на поверхности.

 

 

Согласно сказанному выше скорость в трѳхкоыпонентной

среде

можно представить,как и в случае двухкомпонѳнтной сре­

ды, в

следующем виде:

 

 

 

с.

X

(1.60)

 

 

 

28

гдѳ ß и ß cp - соответственно средняя сжимаемость и средняя плотность трехкомпонентной среды, связанные с вели­ чинами сжимаемостей Д , , , д и плотностей /> , Д , и д каждой из компонент следующими зависимостями:

 

 

 

Л » “

 

 

 

 

 

,

 

 

(I,6I)

где

â[ -

объемные части

каждой

из

компонент;

 

+

=

 

Величины

А

можно найти,

зная

плотность

д -

и ско­

рость распрос

ранения

с,і

в каждой компоненте

по следующей

формуле:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А ,-

=

 

 

 

 

 

 

 

(1.62)

 

Подставляя в

(1.60)

выражения ( І .6 І )

и (1 .62),

получим

выражение

для

 

 

 

 

 

 

 

 

 

 

 

си

 

 

 

 

 

 

 

 

 

 

 

 

(1.63)

 

' У Ы 4 ' w h i - т ш т ’

где <f><£ - объем пор,

заполненных жидкостью и газом,

или

суммарная пористость среды т

.

Величина

,

так

же, как

и с

-

величина скорости в

двухкомпонентной среде,

зави­

сит

от соотношения

плотностей,

скоростей

и объемов

ь яіпо-

нент. На рис. 6 приведена теоретическая зависимость окорости сСрУа / от объема пор, заполненных водоі’.^при различных значѳниях общей пористости среды т

Срцм/се*

Рис. 6. Тѳоретичѳсказависимость скорости оаспространѳния воль от объема пор, заполнен­ ных жидкостью, в трѳхковпонѳнтноя орѳде

29

Соседние файлы в папке книги из ГПНТБ