
книги из ГПНТБ / Ямщиков В.С. Геоакустика. Раздел Упругие волны в неоднородном массиве [учеб. пособие]
.pdfПродольные волны первого и второго рода имеют следую щую физическую сущность.
Волны первого рода определяются сжимаемостью материала фае (волны сжимаемости), а волны второго рода - переупаков кой, т .ѳ . взаимными смещениями твердых частиц (волны пере упаковки). В экспериментах на образцах двухфазных сред наблю дается наличие продольных волн двух типов. При экспериментах в натурных условиях фиксируется только один тип продол ,ных волн, что объясняется большим затуханием волн переупаковки в силу значительных фильтрационных потерь (движение жидкости относительно скелета), вследствие чего волны переупаковки быстро становятся ненаблюдаемыми.
Поперечные волны в водойаоыщенной пористой (двухфазной) среда. Закономерности распространения поперечных волн в по ристых средах описываются системой уравнений (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 5» |
Н |
- |
|
|
|
|
|
|
|
|
|
|
|
При переходе через границу і? |
- |
Л |
скороить Ср |
скач |
|||||||
кообразно увеличивается, |
что |
можно заметить |
из формулы (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