Тензор сейсмического момента

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

История и развитие

Концепция тензора сейсмического момента была разработана в 1960-х — 1970-х годах в рамках теории упругости и сейсмологии. Основополагающий вклад внесли такие учёные, как Кейити Аки (Япония) и Хироо Канамори (Япония, США). В 1977 году Канамори предложил метод расчёта тензора сейсмического момента по записям объёмных волн, что позволило перейти от качественных оценок механизма очага к количественным. С развитием глобальных сейсмических сетей (например, Global CMT Project, основанный в 1981 году) стало возможным автоматизированное определение тензора для землетрясений магнитудой более 5,0 по всему миру.

Математическое определение

Тензор сейсмического момента M определяется как интеграл по поверхности разрыва:

\[ M_{ij} = \int_S \mu \, (u_i n_j + u_j n_i) \, dS \]

где:

• \( \mu \) — модуль сдвига горных пород;
• \( u_i \) — вектор смещения на разрыве (подвижка);
• \( n_j \) — вектор нормали к поверхности разрыва;
• \( S \) — площадь разрыва.

В изотропной среде тензор сейсмического момента является симметричным (\( M_{ij} = M_{ji} \)) и имеет шесть независимых компонент. В декартовой системе координат он записывается как матрица 3×3:

\[ \mathbf{M} = \begin{pmatrix} M_{xx} & M_{xy} & M_{xz} \\ M_{xy} & M_{yy} & M_{yz} \\ M_{xz} & M_{yz} & M_{zz} \end{pmatrix} \]

Скалярный сейсмический момент \( M_0 \) связан с тензором через норму:

\[ M_0 = \frac{1}{\sqrt{2}} \sqrt{\sum_{i,j} M_{ij}^2} \]

Компоненты и их интерпретация

Тензор сейсмического момента можно разложить на две основные составляющие:

Изотропная часть

Изотропная часть \( M_{\text{iso}} \) соответствует объёмному расширению или сжатию в очаге. Она равна \( \frac{1}{3} \text{Tr}(\mathbf{M}) \), где \( \text{Tr} \) — след матрицы. Для большинства тектонических землетрясений изотропная часть близка к нулю, так как разрыв происходит без изменения объёма пород. Значительная изотропная компонента характерна для взрывов (например, ядерных испытаний) или вулканических событий.

Девнаторная часть

Девнаторная часть \( \mathbf{M}_{\text{dev}} \) описывает чистый сдвиг без изменения объёма. Она получается вычитанием изотропной части из полного тензора:

\[ \mathbf{M}_{\text{dev}} = \mathbf{M} - \frac{1}{3} \text{Tr}(\mathbf{M}) \mathbf{I} \]

где \( \mathbf{I} \) — единичный тензор. Девнаторная часть, в свою очередь, может быть разложена на два типа источников: двойную пару сил (DC) и компенсированный линейный векторный диполь (CLVD). Двойная пара сил соответствует классическому сдвиговому разрыву по разлому, а CLVD — более сложным процессам, таким как изгиб или трещинообразование.

Определение механизма очага

По тензору сейсмического момента можно восстановить механизм очага землетрясения — ориентацию плоскости разрыва и направление подвижки. Для этого используются собственные значения и собственные векторы тензора. Собственные значения \( \lambda_1 \geq \lambda_2 \geq \lambda_3 \) определяют тип разрыва:

• Если \( \lambda_2 \approx 0 \), то механизм близок к чистому сдвигу (двойная пара сил).
• Если \( \lambda_1 \approx -\lambda_3 \) и \( \lambda_2 \approx 0 \), то разрыв является сдвиговым.
• Если все три собственных значения одинаковы по знаку, то механизм — объёмное расширение или сжатие.

Собственные векторы задают оси напряжений: ось \( T \) (наибольшее растяжение) соответствует \( \lambda_1 \), ось \( P \) (наибольшее сжатие) — \( \lambda_3 \), а ось \( B \) (промежуточная) — \( \lambda_2 \).

Методы расчёта

Расчёт тензора сейсмического момента производится по сейсмическим записям (сейсмограммам) с использованием методов инверсии. Основные подходы включают:

Инверсия объёмных волн

Используются записи продольных (P) и поперечных (S) волн на удалённых станциях. Метод требует знания скоростной модели земной коры и мантии. Позволяет оценить тензор для землетрясений с магнитудой от 4,5.

Инверсия поверхностных волн

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

Совместная инверсия

Комбинирование данных объёмных и поверхностных волн повышает точность и устойчивость решения. Современные алгоритмы (например, в проекте Global CMT) используют автоматизированную обработку с учётом шумов и неоднородностей среды.

Применение

Тензор сейсмического момента широко используется в сейсмологии и геофизике:

• Определение механизма очага: позволяет классифицировать землетрясения по типу разрыва (сдвиг, взброс, сброс) и оценить ориентацию разлома.
• Оценка сейсмической опасности: анализ тензоров для исторических землетрясений помогает выявить активные разломы и их напряжённое состояние.
• Изучение тектоники плит: распределение осей \( P \) и \( T \) по регионам указывает на направление тектонических напряжений.
• Мониторинг ядерных испытаний: изотропная компонента тензора позволяет отличать взрывы от землетрясений (например, в рамках Договора о всеобъемлющем запрещении ядерных испытаний).
• Вулканология: аномальные тензоры сейсмического момента могут указывать на движение магмы или гидротермальную активность.

Примеры

• Землетрясение в Японии 11 марта 2011 года (магнитуда 9,0–9,1): тензор сейсмического момента показал сдвиговый механизм с субгоризонтальной осью сжатия, направленной на запад, что соответствует субдукции Тихоокеанской плиты под Охотскую.
• Землетрясение в Непале 25 апреля 2015 года (магнитуда 7,8): тензор выявил взбросовый механизм с осью сжатия, ориентированной с юга на север, что типично для столкновения Индийской и Евразийской плит.
• Ядерный взрыв в КНДР 3 сентября 2017 года: тензор сейсмического момента показал значительную изотропную компоненту (около 30% от общего момента), что подтвердило искусственное происхождение события.

Ограничения и критика

Тензор сейсмического момента основан на предположении о точечном источнике, что справедливо для землетрясений с размерами очага, малыми по сравнению с длиной волны. Для крупных событий (магнитуда > 8) требуется учёт конечного размера разрыва, что приводит к использованию более сложных моделей (например, тензора момента с конечным разрывом). Кроме того, точность расчёта зависит от качества сейсмических данных и скоростной модели среды, что может вносить систематические ошибки.

Источники

• Aki, K., & Richards, P. G. (2002). Quantitative Seismology (2nd ed.). University Science Books.
• Kanamori, H. (1977). The energy release in great earthquakes. Journal of Geophysical Research, 82(20), 2981–2987.
• Dziewonski, A. M., & Woodhouse, J. H. (1983). An experiment in systematic study of global seismicity: Centroid-moment tensor solutions for 201 moderate and large earthquakes of 1981. Journal of Geophysical Research, 88(B4), 3247–3271.
• Global CMT Project. (2024). Catalog of Centroid-Moment Tensor Solutions. Columbia University.
• Stein, S., & Wysession, M. (2003). An Introduction to Seismology, Earthquakes, and Earth Structure. Blackwell Publishing.