№ 2 (2025)
Статьи
Применение метода сейсмической интерферометрии для исследования внутреннего ядра земли
Аннотация
Для более детального описания динамических процессов и определения свойств областей внутреннего ядра, недоступных изучению традиционными методами, в работе применен метод сейсмической интерферометрии, в основе которого лежит техника кросс-корреляционного анализа временных рядов, для различных типов данных. Кросс-корреляционный анализ окна сейсмической коды с началом через три часа после сильного события и концом через 10 часов выполнен для всех возможных пар более чем 300 станций и 6 крупных землетрясений (для каждого в отдельности), произошедших в период с 2013 по 2024 гг. Синтетические кросс-коррелограммы рассчитаны для моделей с различным затуханием и дополнительной границей во внутреннем ядре. Проведено 4 разных типа исследований внутреннего ядра методом сейсмической интерферометрии: глобальное, региональное, зависящее от широты станции и от календарного времени. Продемонстрирована устойчивость волны PKIKPPKIKP на глобальных коррелограммах, возможность ее наблюдения на территориях с высокой и низкой плотностью сейсмических станций, зависимость времени пробега волны от угла между направлением волны и осью вращения Земли, а также стационарность волны для временного периода с 2013 по 2024 гг.



Особенности зон локализации сильнейших землетрясений Курило-Камчатской дуги
Аннотация
Процессы, происходящие при взаимодействии блоков земной коры в зоне субдукции, находят отражение как в особенностях сейсмического режима, так и в распределении очагов сильнейших землетрясений исследуемой области. В данной работе проводится анализ особенностей зон локализации сильнейших землетрясений Курило-Камчатской дуги на основе данных двух каталогов: регионального каталога Камчатского филиала ФИЦ ЕГС РАН и глобального международного каталога USGS NIEC геологической службы США за 1990–2024 гг. Объединение каталогов на основе модифицированного метода “ближайшего соседа” позволило удалить дубликаты — повторяющиеся землетрясения в исходных каталогах — и получить новый единый каталог, состоящий из 52574 землетрясений. Метод “ближайшего соседа” применен для выделения двух подмножеств: независимых и связанных землетрясений в пространстве и времени, которые далее использовались для анализа плотности сейсмических событий. Полученные закономерности сопоставлены с очаговыми характеристиками сильнейших землетрясений зоны субдукции. Предложенный подход позволил выделить зоны локализации потенциально возможных очагов мегаземлетрясений вдоль Курило-Камчатской дуги.



Особенности строения зоны фазовых переходов мантии в районе о. Сахалин по данным обменных волн
Аннотация
В работе представлены оценки рельефа границ зоны фазовых переходов мантии на глубинах около 410 и 660 км по данным сети сейсмических станций о. Сахалин методом функций приемника. Проанализирован представительный набор данных, состоящий из 2500 индивидуальных PRF. Выявлено, что граница 660 испытывает опускание в центральной и северной частях острова. Граница 410 демонстрирует сильное поднятие в южной части о. Сахалин, с последующим опусканием ниже стандартных значений во всем регионе исследований. Высказано предположение о связи погружения границы 410 с присутствием в зоне фазовых переходов горячего вещества нижней мантии в северной части острова.



Инстантонное представление форшок—афтершоковых последовательностей
Аннотация
В сейсмологии из практических и теоретических соображений весьма важно описание процессов форшоковой и афтершоковой активизации. Эмпирически установлена аналогия математических соотношений, описывающих характер прямого и обратного законов Омори–Утсу. Исследования обобщенной окрестности сильного землетрясения (ООСЗ) выявили еще более тесную аналогию свойств форшоков и афтершоков. Такая аналогия распространяется также на характеристики процесса активизации, в частности, на аномальные изменения наклона графика повторяемости. В качестве унифицирующей модели всего форшок–афтершокового процесса в данной работе предложено использовать кинетические уравнения, имеющие решения в виде локализованных во времени сильно выраженных максимумов, называемых инстантонами (по аналогии с солитонами — локализованными волнами). Наглядный образ инстантонного решения — график производной по времени от логистической зависимости, описывающей переходный процесс. Скорость такого процесса сначала значительно возрастает, достигает максимума, а затем асимптотически убывает до нуля.
Цель работы — продемонстрировать эффективность использования инстантонной модели, которая обобщает модель саморазвивающихся процессов (СРП), но не предусматривает развитие физически нереализуемой сингулярности, обычно моделирующей взрывной рост числа фор- и афтершоков в окрестности главного события. Сравнение новой модели с эмпирическими данными проводится на примере землетрясений в наиболее обеспеченной средствами регистрации южной части о-ва Сахалин за 2003–2023 гг. Показано удовлетворительное сходство теоретических и эмпирических временных зависимостей как для ООСЗ, построенной для территории в пределах (44.5°–50.5° с. ш., 141.5°– 43.5° в. д.), так и для отдельных сильных землетрясений на Сахалине.



Кольцевые структуры сейсмичности в районе Южной Камчатки: возможная подготовка сильнейшего землетрясения
Аннотация
Рассматриваются некоторые характеристики сейсмичности в районе Южной Камчатки. Исследованы характеристики облака афтершоков сильного землетрясения 17.08.2024 г. (Mw = 7.0). Показано, что в районе Южной Камчатки сформировались кольцевые структуры сейсмичности в трех диапазонах глубин: 0–33, 34–70 и 71–110 км. Как и в других зонах субдукции, структуры характеризуются пороговыми значениями магнитуд (соответственно Мп1, Мп2 и Мп3), а также длинами больших осей (L1, L2 и L3). Эпицентры землетрясения 17.08.2024 г. и его наиболее сильных афтершоков попадают на выделенную здесь неглубокую кольцевую структуру (Мп1 = 5.3), что подтверждает предположение о подготовке в районе Южной Камчатки сильного землетрясения.
Ранее были построены корреляционные зависимости параметров Мп1 и Мп2 от магнитуд Mw сильных землетрясений для запада Тихого океана (в диапазоне Mw = 7.0–9.0). Используя эти зависимости, мы оценили магнитуду возможного здесь сильнейшего события: Mw = 8.6 ± 0.2. Обсуждаются причины формирования кольцевых структур сейсмичности на разных глубинах в зонах субдукции.



Развитие зоны нарушенной породы в окрестности динамической подвижки по тектоническому разлому
Аннотация
В статье приведены результаты 2D-расчетов процесса формирования зоны нарушенной породы при развитии динамического сдвига по горизонтальному тектоническому разлому. Исследованы различные режимы скольжения — субрелеевский (скорость распространения разрыва Vr не превышает скорости релеевской волны в среде) и сверхсдвиговый (величина Vr выше скорости поперечных волн). Рассмотрен вклад механизмов отрыва и сдвига в развитие нарушенной зоны вблизи разлома на разной глубине. Проведена оценка степени изменения физико-механических свойств массива на разных расстояниях от разлома. Показано, что на больших глубинах литостатические напряжения полностью подавляют отрыв, и разрушение породы происходит исключительно за счет деформации сдвига. На малых же глубинах механизм разрушения отрывом становится преобладающим. Сброс напряжений, связанный с появлением трещин отрыва, приводит к резкому уменьшению зоны разрушения сдвигом, которая локализована лишь в непосредственной близости от плоскости разрыва. Увеличение прочности на отрыв приводит к увеличению размеров зоны сдвигового разрушения. У сверхсдвиговых разрывов зона разрушения может иметь сложный, неодносвязанный характер. Изменение скорости распространения продольных волн Cp более чем на 15–20% имеет место лишь в непосредственной близости от плоскости скольжения на расстоянии 10–20 м. На больших расстояниях величина не превышает 10%. На малых глубинах могут иметь место трещины отрыва, которые распространяются на значительные расстояния от плоскости скольжения.



Природа анизотропного отклика флюидонасыщенной среды на прохождение поверхностных сейсмических волн
Аннотация
Мониторинг изменения порового давления или уровня воды в наблюдательных скважинах показывает значительные колебания как при прохождении Р-волн и волн Релея, так и при прохождении S-волн и волн Лява. Недавние измерения в скважине показали азимутальную зависимость колебаний порового давления относительно ориентации напряжений и направления простирания разломной зоны. В зоне активной разломной зоны анизотропия, индуцированная трещиноватостью, соответствует преимущественной ориентации микротрещин и других дефектов сплошности среды. Настоящая работа посвящена разработке модифицированного уравнения Скемптона для количественного описания изменений порового давления в коллекторе, вызванных прохождением поверхностных волн и связанных с ориентацией и главными значениями тензора напряжений и поврежденности (трещиноватости) пород. Разработанные соотношения позволяют описать азимутальную зависимость отклика порового давления с помощью безразмерного параметра, определяемого как отношение амплитуд колебаний давлений, вызванных сдвиговой компонентой и объемной деформацией. Согласно предложенной теоретической модели, максимальный пороупругий отклик коллектора на прохождение сейсмической волны проявляется в случае субпараллельности направлений преимущественной трещиноватости пород и максимального горизонтального напряжения.
Для верификации предложенной теоретической модели использованы данные мониторинга порового давления в резервуаре, используемом под захоронение сточных вод Арбакл (Оклахома, США). Показано, что наблюдаемое разнообразие отклика порового давления в скважинах, расположенных в окрестностях разломной зоны, пересекающей резервуар, на прохождение сейсмических волн от разноудаленных сейсмических событий с высокой точностью описывается разработанной моделью.



Реконструкция параметров сдвиговых напряжений при формировании разноранговых разломов западного Прибайкалья на основе тектонофизической интерпретации линеаментов
Аннотация
С помощью авторского программного обеспечения удалось осуществить детальное выделение линеаментов на различных масштабных уровнях для района Обручевской системы разломов западного Прибайкалья. Выделенные линеаменты значительно дополняют закартированный каркас разрывных нарушений и согласуются с простиранием соответствующих по рангу разрывных структур. На основе анализа относительной удельной плотности линеаментов локального масштабного уровня, отражающих оперяющие мегатрещины крупных разломов, установлены неоднородные зоны динамического влияния региональных структур, которые разделены на относительно однородные сегменты. Для каждого выделенного сегмента и каждой структуры в целом с помощью разработанного программного обеспечения “Lineament Stress Calculator” проведена реконструкция параметров сдвиговых напряжений с использованием модели П. Л. Хэнкока. Ранее доказано, что основные черты раннепалеозойского этапа развития региона в процессе аккреции Ольхонского террейна к южной окраине Сибирского кратона сопровождались активизацией правосдвиговых смещений по ЮЗ–СВ аккреционным швам и процессами активного метаморфизма. Полученными результатами подтверждено, что магистральные разломы ЮЗ–СВ простирания, субпараллельные краевому шву Сибирской платформы, на раннем этапе развития формировались как правые сдвиги при ориентировке оси сжатия ≈ 90°. Разломы второго порядка СЗ–ЮВ ориентировки определены как левые сдвиги и, вероятно, закладывались в это время как антитетические сколы по отношению к магистральным структурам, получив свое развитие при дальнейших структурных перестройках региона.



Восстановление пространственного распределения фильтрационных свойств неоднородных геологических сред по вариациям микросейсмичности, возникающей в результате закачки жидкости
Аннотация
Определение свойств неоднородных коллекторов по данным об эволюции микросейсмичности является важной задачей при разработке месторождений. Анализ распространения микросейсмических событий, возникающих при закачке/отборе флюида, позволяет получить ценную информацию о проницаемости и напряженном состоянии пласта. В настоящей статье рассматривается обратная задача определения фильтрационных свойств коллектора по данным о распространении микросейсмических событий. Для этого исследуется влияние различных геологических факторов на распределение очагов микросейсмических событий. Для выявления корреляционных связей между параметрами геологической модели и эволюцией микросейсмичности были использованы методы машинного обучения. Из-за недостаточной вариативности натурных данных для обучения модели была создана искусственная база каталогов микросейсмических событий, содержащих координаты очагов и времена их возникновения. Для этого проведено численное моделирование закачки жидкости и генерации микросейсмических событий в синтетических моделях проницаемых сред с разным геологическим строением. Таким образом, предложен комплексный подход к восстановлению фильтрационных свойств неоднородных коллекторов по данным об эволюции микросейсмичности с использованием методов машинного обучения. Предложенная методика может быть применена для оптимизации разработки месторождений, повышения эффективности извлечения флюидов и снижения рисков, связанных с возникновением нежелательной техногенной сейсмической активности.



Акустическая эмиссия, сопровождающая подготовку динамического проскальзывания по модельному гетерогенному разлому метрового масштаба
Аннотация
Закономерности распределения участков тектонических разломов с разными фрикционными свойствами в значительной степени контролируют динамику их скольжения. Невозможность прямого изучения структуры разломных зон на сейсмогенных глубинах делает особенно актуальным разработку методов диагностики, позволяющих получать информацию о структурных особенностях областей формирования очагов землетрясений и тем самым прогнозировать динамику скольжения.
В настоящей работе представлены результаты лабораторных экспериментов, в которых исследовались закономерности излучения упругих колебаний в процессе эволюции напряженно-деформированного состояния модельного разлома, имеющего пространственно-неоднородную структуру плоскости скольжения. Модельный разлом представлял собой нагруженный контакт блоков диабаза и имел размер 750 × 120 мм2. На интерфейсе разлома были сформированы две зоны круглой формы диаметром по 100 мм, которые обладают повышенной прочностью со свойством скоростного разупрочнения, так называемые асперити. В ходе экспериментов изменялось относительное расположение таких зон.
Процесс формирования динамического проскальзывания, обусловленного разрушением асперити, сопровождался излучением большого количества акустических импульсов, регистрируемых в полосе частот 20–80 кГц. В ходе экспериментов данные о пространственном распределении импульсов позволяют выявить две отдельные контактные области только при расстояниях между этими участками более 20 мм. При этом наблюдаются различия в статистике импульсов, излученных на различных асперити.



Эволюция микротрещин в процессе деформирования горных пород: рентгеновская микротомография и метод дискретных элементов
Аннотация
В работе выполнено прямое наблюдение микротрещин, образующихся в образце горной породы под действием одноосной сжимающей нагрузки. Детектирование дефектов в объеме осуществлялось с помощью рентгеновской компьютерной микротомографии. Особенность проводимых экспериментов состоит в том, что выполнялась томографическая съемка образца, находящегося под действием нагрузки. На основе анализа томографических сечений вычислены фрактальная размерность и относительный объем микротрещин на трех этапах нагружения. Построены трехмерные модели дефектной структуры, которые иллюстрируют изменение морфологии магистральной трещины. С помощью модели дискретных элементов проведены численные эксперименты по разрушению образцов гетерогенных материалов. Исследовано изменение фрактальной размерности магистральных трещин в процессе их роста. Установлено хорошее согласие результатов компьютерного моделирования и лабораторных экспериментов, что указывает на адекватность предложенной модели и позволяет в дальнейших исследованиях использовать ее для изучения поведения локальных параметров, которые не могут быть измерены экспериментально.



О режиме вызванной сейсмичности
Аннотация
Проблема вызванной сейсмичности имеет важный как практический, так и теоретический аспект. Практический аспект связан с опасностью вызванной сейсмичности. В ряде случаев опасность сильной вызванной сейсмичности приводила к закрытию важных промышленных проектов. Теоретический аспект связан с известным парадоксом сейсмичности, о невозможности реализации обычных землетрясений на глубинах более нескольких десятков километров. Отсюда следует, что физика вызванных, как правило, приповерхностных землетрясений может отличаться от физики большинства более глубоких событий. Рассмотрены примеры ряда областей вызванной сейсмичности, представляющих собой как окрестности крупных водохранилищ, так и области интенсивной добычи углеводородного и рудного сырья. Во всех рассмотренных районах, более или менее определенно, выявляется ряд общих тенденций. После роста вызванной сейсмичности, даже и при продолжающемся сильном техногенном воздействии, наблюдается тенденция спада сейсмической активности. Также, методом анализа обобщенной окрестности сильного землетрясения (ООСЗ), для зон вызванной сейсмичности выявляется близость интенсивности фор- и афтершокового процесса, что контрастирует со случаем обычной сейсмичности, для которой активность афтершокового процесса обычно много сильнее. Можно предположить, что спад вызванной сейсмичности связан с разгрузкой исходных тектонических напряжений, а близость интенсивности фор- и афтершокового процесса указывает на отличие физического механизма вызванных приповерхностных землетрясений от обычных, более глубоких землетрясений.



К 20-летию создания малоапертурной группы “Михнево”. Мониторинг наведенной сейсмичности
Аннотация
За 20 лет непрерывной работы малоапертурной сейсмической группы (МСГ) “Михнево” был накоплен огромный опыт регистрации сверхслабых сигналов, генерируемых региональной и глобальной сейсмичностью, разработаны и применены методы обработки данных с высокой разрешающей способностью, включая направленное суммирование и кросс-корреляцию волновых форм. В рамках данного обзора результатов инструментальных наблюдений и обработки рассмотрены два подхода к снижению порога регистрации сейсмических событий при мониторинге наведенной сейсмичности — использование станций группирования и метода кросс-корреляции волновых форм (ККВФ). Эффективность подходов по отношению к обнаружению, локации и идентификации слабых сейсмических источников проиллюстрирована афтершоковой последовательностью землетрясения вблизи г. Мариуполь, произошедшего 07.08.2016 г., а также афтершоками пятого и шестого объявленных взрывов в КНДР, обнаруженных за период с 09.09.2016 по 11.09.2021 гг. Координаты землетрясения были оценены с помощью данных МСГ “Михнево” и временной МСГ ИДГ РАН “Ростов Дон”. Точность локации сопоставима с точностью, обеспечиваемой 49 трехкомпонентными (3-С) станциями ФИЦ ЕГС РАН и Международной системы мониторинга (МСМ). За пять дней после землетрясения методом ККВФ было обнаружено и лоцировано относительно главного толчка 12 афтершоков. В обнаружении и оценке параметров также участвовали групповые станции МСМ AKASG и BRTR и 3-С станция МСМ KBZ. Сеть станций ФИЦ ЕГС РАН обнаружила 5 афтершоков, а МСМ ни одного. Локация взрывов в КНДР с помощью ККВФ позволила определить их относительное местоположение с точностью до 100–200 м. Шестой взрыв не удалось точно лоцировать относительно других из-за конечного размера его источника, вносящего значимые изменения дифференциального времени пробега, зависящие от направления на станцию. Метод ККВФ использовался и для обнаружения и идентификации слабых сейсмических событий в пределах полигона КНДР Пунгери с помощью шаблонных волновых форм от взрывов и афтершоков пятого и шестого испытаний, зарегистрированных на групповых станциях МСМ KSRS и USRK. За пятилетний период наблюдений было обнаружено 89 событий. По оценкам кросс-корреляционных характеристик сигналов на обеих станциях оказалось возможным разделить общую афтершоковую последовательность на две отдельные, связанные с процессами в зонах влияния пятого и шестого взрывов.



Выделение энергии в атмосфере при падении космических тел размером 20–250 метров
Аннотация
Приведены результаты расчетов разрушения, испарения и торможения каменных метеороидов размером от 20 до 250 метров в атмосфере Земли. Подробно исследуется перераспределение тепловой и кинетической энергии между конденсированным веществом метеороида, его парами и воздухом. Показано, что при размерах ударника в несколько десятков метров испаренное вещество не тормозится сразу, а долго летит вдоль траектории, постепенно передавая энергию воздуху. В результате основное энерговыделение в атмосфере происходит на стадии торможения струи паров, уже после того, как метеороид и его фрагменты полностью испарились.


