Scientific journal
Advances in current natural sciences
ISSN 1681-7494
"Перечень" ВАК
ИФ РИНЦ = 0,775

THE COMBINED CALCULATION METHODOLOGY OF SEISMIC-WAVE FRONTS FOR THE ANISOTROPIC MEDIA BY THE EXAMPLE OF HEXAGONAL CRYSTALS

Melnikova N.A. 1 Nemirovich-Danchenko М.М. 2 Islyamova А.А. 2 Kpata Мc E.R.E. 2
1 STI National Research Nuclear University «MEPhI»
2 National research Tomsk Polytechnic University
The calculation methodology of seismic-wave fronts for the anisotropic media having elastic transversal isotropy is presented in the article. The obtaining of Christoffel’s tensor components in 2D hexagonal syngony is described. The values of density, four independent tensor components of elastic constants and two dimensionless anisotropy coefficients are given for some hexagonal monocrystals. Radial velocity indicatrixes are plotted for typical cases of hexagonal syngony when the loading axis matches the main anisotropy axis. When the loading axis doesn’t match the main anisotropy axis, waves fronts are obtained by finite and incremental modelling. The preceding methodology may be used for the analysis of laboratory data and full-scale experiments in geophysics, especially in seismic i.e. in case of plotting time seismic cuts and studying features of seismic waves distribution in the horizontally stratified medium.
seismic
wave fronts
hexagonal symmetry
transversal isotropy
Christoffel’s tensor
elastic constants
symmetry axis rotation
numerical modelling

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

Фундаментальные основы для решения задач кристаллоакустики заложены во второй половине нашего столетия Ф.И. Федоровым [1]. Удается получать и исследовать общее решение уравнений эластодинамики для широкого класса практических случаев. Эти результаты находят свое применение в физической акустике, дефектоскопии, пьезотехнике, сейсморазведке. При этом для многих задач возможно лишь численное решение, но для правильной интерпретации результатов расчетов необходимо опираться на некие общие аналитические зависимости. Авторами ниже рассмотрены несколько характерных случаев гексагональной сингонии и построены индикатрисы лучевых скоростей. При несовпадении оси нагружения с главной осью анизотропии результат получен конечно-разностным моделированием.

Тензор Кристоффеля. В трансверсально изотропных средах наибольшая анизотропия упругих свойств наблюдается в плоскости (010). Распространение упругих волн будем рассматривать в этой плоскости. Совместим орты кристаллографической системы координат с осями X, Z. Будем считать, что среда находится в условиях плоской деформации. Запишем уравнения движения в перемещениях

mel01.wmf (1)

Уравнения (1) являются линейными однородными дифференциальными уравнениями второго порядка относительно вектора смещения [2].

Решение уравнений (1) представляет собой семейство плоских волн

mel02.wmf (2)

Здесь mel03.wmf – радиус-вектор частицы среды, А – амплитуда волн, mel04.wmf – вектор поляризации, mel05.wmf – волновой вектор (вектор пространственных частот), mel07.wmf – круговая частота, где f – частота (число полных колебаний вектора смещения за одну секунду). Единичный вектор волновой нормали mel08.wmf, такой, что mel09.wmf, где λ – длина волны, в оптике называется вектором рефракции (медленность или обратная скорость). Возмущение частоты f в среде пробегает путь λ за время mel10.wmf, а величину mel11.wmf называют фазовой скоростью, где k – длина вектора mel12.wmf.

Выражение (2) для вектора смещения плоской упругой волны является комплексным, тогда как реальные смещения должны выражаться вещественными величинами. Разрешение этого противоречия заключается в том, что уравнения (1) являются линейными однородными уравнениями, причем все операции, производимые над неизвестным вектором, не содержат мнимых выражений [2].

Учитывая вышесказанное, равенство (2) запишем в виде

mel13.wmf (3)

подставляя (3) в уравнения (1), получим

mel14.wmf (4)

Здесь nj, nk – направляющие косинусы вектора волновой нормали, а свертка cijkl nj nk есть тензор второго ранга – тензор Кристоффеля Гil. Введем вместо тензора cijkl тензор mel15.wmf, который называется приведенным тензором модулей упругости. Тогда свертка ijkl nj nk есть приведенный тензор Кристоффеля mel16.wmf. Перепишем уравнение (4) в виде

mel17.wmf (5)

Векторное уравнение (5) есть уравнение Кристоффеля. Оно является основным для всей теории упругих волн в кристаллах. Собственными значениями приведенного тензора Кристоффеля mel18.wmf являются квадраты фазовых скоростей, а его собственным вектором – вектор поляризации mel19.wmf. При вычислении компонент тензора Кристоффеля нужно учитывать, что, например, c2323 = c2332 = c3232 = c3223 = c44.

В плоскости (010) матрица приведенного тензора Кристоффеля будет иметь вид

mel20.wmf

где индекс 4 соответствует сочетаниям индексов (13) и (31).

Характеристическое уравнение матрицы тензора Кристоффеля

mel21.wmf

имеет корни

mel22.wmf (6)

Трансверсально изотропная среда (в частности, гексагональные кристаллы) является единственным примером анизотропной среды, для которой уравнение Кристоффеля при любом направлении волновой нормали распадается на линейное и квадратное уравнения и поэтому допускает сравнительно простое решение [2]. Поскольку такие среды в отношении упругих свойств обладают симметрией вращения вокруг оси 6-го порядка, то все волновые поверхности будут поверхностями вращения. Чтобы их получить, достаточно найти вид соответствующей кривой в меридиональной плоскости (плоскости симметрии, проходящей через ось 6-го порядка), после чего вращение кривой вокруг оси 6-го порядка даст искомую поверхность.

Результаты исследования и их обсуждение

Индикатрисы фазовых скоростей υ1,2 и медленностей 1/υ1,2 легко строятся непосредственно по формулам (6). Для построения волновых (лучевых) поверхностей нами использован геометрический метод Postma [8]. Для иллюстрации характера этих кривых выберем гексагональные кристаллы с различной степенью анизотропии.

Упругие модули, плотности и параметры анизотропии гексагональных кристаллов

Вещество

С11

С12

С13

С33

С44

ρ

Kp

Ksv

Лёд Ih [6]

13,84

7,06

5,81

14,99

3,19

919

0,96

1,16

Be [3]

292,3

26,7

14

336,4

162,5

1848

0,932

0,96

α–C [5]

1109

139

0

38,7

4,95

2267

5,35

7,61

Апатит [4]

167

13,1

66

140

66,3

3200

1,09

0,81

Ti [3]

162,4

92

69

180,7

46,7

4506

0,95

1,05

Zn [3]

165

34,2

50,1

62

39,6

7100

1,63

0,89

Co [3]

307

165

103

358

75,5

8900

0,93

1,23

WC [7]

720

254

267

972

328

15770

0,86

0,94

В таблице приведены константы для ряда гексагональных кристаллов (компоненты cij даны в ГПа, плотность – в кг/м3). Данные расположены по возрастанию плотности монокристаллов от разновидности Ih пресного льда до карбида вольфрама WC. Даны ссылки на работы, из которых были взяты константы. Далее, кроме обычно используемых упругих модулей и плотности в таблице приводятся так называемые безразмерные эффективные параметры анизотропии [9]:

mel23.wmf – характеризует анизотропию по скорости квази-P–волн;

mel24.wmf – характеризует анизотропию по скорости квази-SV-волн.

Эти параметры удобны при моделировании в плоскости (010), так как для изотропной среды mel25.wmf.

Из таблицы выбраны три характерных среды: 1 – стабильная разновидность пресного льда Ih, 2 – цинк (Zn), 3 – стабильный графит (α–C).

На рис. 1 построены кривые фазовых скоростей (а), обратных скоростей (б) и лучевых скоростей (в) для льда, на рис. 2, a, б, в – для цинка, и на рис. 3, a, б, в – для графита, одного из предельных случаев анизотропии. В изотропной среде кривые (индикатрисы) для всех скоростей были бы концентрическими окружностями, а поверхности фазовых скоростей, медленностей и лучевых скоростей в пространственном случае – сферами.

melnik1.tif

Рис. 1. Индикатрисы медленностей (а), фазовых скоростей (б) и лучевых скоростей (в) монокристалла пресного льда Ih. Черная линия – квазипродольная волна, серая – квазипоперечная волна

melnik3.tif

Рис. 3. Индикатрисы медленностей (а), фазовых скоростей (б) и лучевых скоростей (в) монокристалла графита. Черная линия – квазипродольная волна, серая – квазипоперечная волна

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

По индикатрисам льда заметна слабая анизотропия, и областей неоднозначности не видно. Зато для монокристалла цинка хорошо видна выраженная анизотропия по продольной скорости (латеральная вытянутость, квази-эллиптичность внешней кривой рис. 2, в). Внутренняя кривая, соответствующая индикатрисе квазипоперечной скорости, обладает характерными признаками сильной анизотропии – зонами рефракции, так называемыми «петлями» на лучах ± π/2.

melnik2.tif

Рис. 2. Индикатрисы медленностей (а), фазовых скоростей (б) и лучевых скоростей (в) монокристалла цинка. Черная линия – квазипродольная волна, серая – квазипоперечная волна

Монокристалл графита является граничным случаем анизотропии в том смысле, что для него, во-первых, характерны максимальные отклонения луча от волновой нормали (ромбовидный фронт), и, во-вторых, для широкого диапазона лучей скорость одной из ветвей квазипоперечной волны практически равна скорости квазипродольной волны – индикатрисы касаются друг друга.

Для случая несовпадения главной оси анизотропии и оси нагружения было решено использовать численный метод расчета полного волнового поля, основанный на конечно-разностном моделировании [1]. Приведем для сравнения волновую картину в задаче Лэмба для монокристалла цинка конечных размеров в случае, когда ось нагружения коллинеарна оси С (рис. 4) и для случая, когда угол между ними составляет 30 ° (рис. 5).

melnik4.tif

Рис. 4. Волновые фронты в монокристалле цинка в случае коллинеарности осей анизотропии и нагружения для двух моментов времени. Буквами обозначены: P – квазипродольная волна, SV – квазипоперечная волна, LR – волна Рэлея, C – коническая волна

melnik5.tif

Рис. 5. Волновые фронты в монокристалле цинка для случая угла 30 градусов между осью нагружения (черная стрелка) и осью анизотропии (линия). Буквой R обозначена отраженная от свободной поверхности продольная волна

Каждая из волновых картин на рис. 4–5 представляет собой мгновенные снимки (скриншоты) волнового поля для двух моментов времени. Вид фронта продольной и поперечной волн на рис. 4 соответствует индикатрисам рис. 2, в. Хорошо различимы петли для квази-SV волн, на скриншоте они воспринимаются как отдельные ветви фронта. Для более раннего момента фронт продольной волны еще не достиг нижней грани, для более позднего начинается отражение. Это показано на рис. 5, б: фронт продольной волны, достигнув нижней грани (на ней задано условие свободной поверхности), начинает отражаться от неё, фронт отраженной волны отмечен буквой R.

Выводы

Изложенная методология может быть использована в сейсморазведке при построении временных сейсмических разрезов и при оценке особенностей распространения сейсмических волн в горизонтально-слоистой среде. Кроме того, методы построения волновых фронтов могут применяться на стадии планирования и обработки результатов лабораторного моделирования в материаловедении, кристаллофизике и дефектоскопии.