Сейсмическая томография

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

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

Ранние этапы (до 1970-х)

Первые попытки использовать сейсмические волны для изучения внутреннего строения Земли относятся к началу XX века. После землетрясения в Сан-Франциско (1906) геофизики заметили, что в некоторых областях планеты сейсмические волны распространяются быстрее или медленнее, чем в среднем. В 1920-х годах Бено Гутенберг и другие учёные построили первые одномерные модели скоростей сейсмических волн, выделив кору, мантию и ядро. Однако эти модели были глобальными и не учитывали горизонтальные вариации.

Формирование метода (1970–1990)

В 1970-х годах, с развитием вычислительной техники и появлением сети сейсмических станций (включая Глобальную сейсмографическую сеть), стало возможным обрабатывать большие объёмы данных. Американский геофизик Адам Дзевоньский (Adam Dziewonski) и его коллеги в 1980-х годах разработали первые трёхмерные модели мантии Земли, используя метод томографии. В 1984 году была опубликована модель PREM (Preliminary Reference Earth Model), ставшая эталоном для сравнения. В СССР и России исследования в этом направлении велись в Институте физики Земли РАН (Москва) и в Сибирском отделении РАН (Новосибирск).

Современный этап (с 1990-х)

С 1990-х годов сейсмическая томография стала одним из основных инструментов глобальной геофизики. Благодаря увеличению числа станций (до нескольких тысяч), улучшению алгоритмов (например, метод конечных разностей, обратное проецирование) и росту вычислительных мощностей, разрешение моделей повысилось с сотен до десятков километров. В 2000-х годах начали применяться методы полноволновой инверсии (Full-Waveform Inversion, FWI), которые учитывают не только времена прихода, но и амплитуды, фазы и форму волн.

Физические основы

Сейсмическая томография базируется на двух типах волн:

• P-волны (продольные, первичные) — распространяются в твёрдых и жидких средах, скорость зависит от модуля всестороннего сжатия и плотности.
• S-волны (поперечные, вторичные) — распространяются только в твёрдых средах, скорость зависит от модуля сдвига.

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

Методология

Сбор данных

Основные источники сейсмических волн:

• Землетрясения — естественные источники, обеспечивающие глобальное покрытие. Для локальных исследований используют события с магнитудой от 3 до 7.
• Взрывы — искусственные источники (например, при геологоразведке или ядерных испытаниях), позволяющие контролировать параметры сигнала.
• Вибрационные источники — для активной сейсморазведки на малых глубинах.

Данные регистрируются сейсмическими станциями, которые измеряют время прихода волн, их амплитуду и поляризацию. В России сеть станций включает более 300 пунктов (например, станции «Обнинск», «Москва», «Новосибирск»), а также временные сети для региональных проектов.

Обработка и инверсия

Процесс включает несколько этапов:

1. Выделение вступлений — автоматическое или ручное определение времени прихода P- и S-волн.
2. Построение лучевых траекторий — расчёт пути волны от источника к приёмнику через среду с известной или предполагаемой скоростной моделью.
3. Обратная задача — решение системы уравнений, связывающих измеренные времена прихода с распределением скоростей. Используются методы:

• Линейная томография — для малых отклонений, решается методом наименьших квадратов.
• Нелинейная томография — итеративное уточнение модели, включая методы Монте-Карло и нейронные сети.

1. Визуализация — построение трёхмерных разрезов (срезов по глубине или латеральных сечений) с цветовой кодировкой аномалий (красный — медленные, синий — быстрые).

Разрешение

Разрешение метода зависит от:

• Плотности сети станций — чем больше станций, тем выше детализация.
• Частотного диапазона — высокочастотные волны (1–10 Гц) дают лучшее разрешение на малых глубинах, низкочастотные (0.01–0.1 Гц) — проникают глубже, но с меньшей детализацией.
• Количества источников — для глобальных моделей используют данные от 10 000–50 000 землетрясений.

Виды сейсмической томографии

По масштабу

• Глобальная томография — изучает мантию и ядро Земли в целом, с разрешением от 100 до 500 км. Используется для построения моделей тектоники плит, мантийной конвекции.
• Региональная томография — охватывает территории размером от 100 до 1000 км (например, Тихоокеанское огненное кольцо, Альпийско-Гималайский пояс). Разрешение — 10–50 км.
• Локальная томография — применяется для изучения верхней коры (до 10–20 км) на участках площадью от 1 до 100 км². Разрешение — до 1–5 км. Используется в геологоразведке, сейсмическом микрорайонировании, поиске полезных ископаемых.

По типу волн

• P-волновая томография — наиболее распространена, даёт информацию о плотности и упругих свойствах.
• S-волновая томография — чувствительна к сдвиговым деформациям, позволяет выявлять зоны частичного плавления.
• Поверхностно-волновая томография — использует волны Релея и Лява, распространяющиеся вдоль поверхности. Эффективна для изучения коры и верхней мантии.
• Полноволновая томография (FWI) — учитывает полную форму волнового поля, включая отражения и преломления, что даёт наиболее точные результаты, но требует огромных вычислительных ресурсов.

По методу инверсии

• Лучевая томография — основана на приближении геометрической оптики, считает волны лучами. Проста, но не учитывает дифракцию.
• Дифракционная томография — учитывает волновые эффекты, использует интегральные преобразования (например, преобразование Радона). Требует более сложных расчётов.

Применение

Геология и геофизика

• Изучение тектоники плит — томография позволила визуализировать субдуцирующие плиты, погружающиеся в мантию до глубин 700–1000 км (например, под Японией и Южной Америкой). В России исследования проводились для зоны субдукции Курило-Камчатской дуги.
• Мантийные плюмы — выявление горячих восходящих потоков, например, под Гавайями, Исландией, а также под Сибирской платформой (Сибирский плюм, связанный с пермско-триасовым вымиранием).
• Строение коры — картирование границ Мохоровичича (Moho) и Конрада, выявление разломов, интрузий, соляных куполов.

Поиск полезных ископаемых

• Нефть и газ — томография высокого разрешения (до 10 м) используется для локализации ловушек в осадочных бассейнах, например, в Западной Сибири и на шельфе Сахалина. Метод позволяет отличать газонасыщенные пласты от водонасыщенных по аномалиям скорости.
• Рудные месторождения — выявление кимберлитовых трубок (алмазы) и колчеданных залежей (медь, цинк) по контрасту скоростей с вмещающими породами. В России томография применялась при разведке месторождений в Норильском районе и на Кольском полуострове.
• Геотермальные ресурсы — картирование зон повышенной температуры и трещиноватости для выбора мест бурения скважин.

Сейсмическая безопасность

• Сейсмическое микрорайонирование — оценка сейсмической опасности на конкретных участках (например, при строительстве АЭС, плотин, мостов) путём выявления зон разжижения грунтов и усиления колебаний.
• Прогноз землетрясений — томография помогает выявлять зоны накопления напряжений, но точный прогноз пока невозможен. В России ведутся работы по мониторингу сейсмической активности на Камчатке и в Прибайкалье.

Планетология

Метод адаптирован для изучения Луны, Марса и Венеры. Например, по данным сейсмических станций миссий «Аполлон» (США) и «Инсайт» (NASA) построены томографические модели лунной коры и мантии.

Примеры значимых результатов

• Глобальная модель S40RTS (2010) — показала, что под Африкой и Тихим океаном существуют крупные области низких скоростей (LLSVP), интерпретируемые как скопления древнего, частично расплавленного материала.
• Зона субдукции под Каскадией (Северная Америка) — томография выявила погружающуюся плиту Хуан-де-Фука, что позволило уточнить прогнозы землетрясений.
• Сибирский плюм — исследования 2010-х годов показали, что под Сибирской платформой на глубинах 100–200 км сохраняются остатки мантийного плюма, связанного с пермско-триасовым магматизмом.

Критика и ограничения

• Неоднозначность решения — обратная задача томографии имеет множество решений, и выбор модели зависит от начальных условий и регуляризации.
• Ограниченное разрешение — на глубинах более 1000 км разрешение падает до 200–500 км, что не позволяет детально изучать ядро.
• Зависимость от качества данных — ошибки в определении времени прихода волн, особенно для слабых землетрясений, могут искажать результаты.
• Вычислительная сложность — полноволновая инверсия требует суперкомпьютеров и времени расчёта до нескольких недель для региональных моделей.

Источники

1. Dziewonski, A. M., & Anderson, D. L. (1981). Preliminary reference Earth model. Physics of the Earth and Planetary Interiors, 25(4), 297–356.
2. Nolet, G. (2008). A Breviary of Seismic Tomography. Cambridge University Press.
3. Рыжов, А. В. (2015). Сейсмическая томография: методы и приложения. Геофизика, 3, 12–25.
4. Kennett, B. L. N., & Widiyantoro, S. (1999). A new global model of the Earth's mantle. Science, 286(5441), 1126–1129.
5. Материалы Института физики Земли РАН (Москва, 2020) — «Сейсмическая томография в России: современное состояние».
6. Rawlinson, N., & Sambridge, M. (2003). Seismic traveltime tomography of the crust and lithosphere. Advances in Geophysics, 46, 81–198.