В 2020-21 гг были выполнены 3 блока исследований:

   В серии лабораторных акустоэмиссионных (АЭ) экспериментов исследовалась динамика фрикционного скольжения разлома и сопровождающий данное скольжение акустоэмиссионный портрет. Было показано, что в условиях фрикционного скольжения процессы, протекающие на микроуровне, контролируют динамику скольжения разлома на макроуровне; а также наблюдается высокая корреляция между скейлинговым параметров b-value (закон Гуттенберга-Рихтера) и фрикционной прочностью модельного разлома. Самоорганизация среды на микроуровне сопровождается излучением большого числа импульсов АЭ. Применение алгоритма количественной категоризации, основанного на статистическом анализе волновой формы импульсов АЭ, указывает на формирование двух популяций АЭ с условно быстро и медленно нарастающей амплитудой и которые характеризуются различными закономерностями вариаций b-value. Рассмотрение b-value для двух популяций АЭ позволяет выявить закономерности эволюции как сильных, так и слабых динамических подвижек и определить переход разлома в критическую стадию лабораторного сейсмического цикла. Снижение b-value для импульсов с «медленным нарастанием волновой формы» является надежным индикатором перехода разлома в критическую стадию и указывает на подготовку и быстрых, и медленных динамических событий, но не определяет время до события. Даже лабораторный модельный разлом является достаточно сложной системой для полного описания физики происходящих на микроуровне процессов. Применение метода машинного обучения с учителем «случайный лес» позволило решить задачи прогнозирования поведения системы, не имея полного математического описания исследуемой системы – был определен момент наступления динамического события. При обучении и прогнозировании не учитывалось прошлое или будущее состояние модельного разлома, на вход подавались лишь данные в пределах определенного короткого временного интервала. Разработанный алгоритм машинного обучения позволяют на основе обработки больших многомерных массивов данных выявлять наиболее информативные параметр системы и с высокой точностью прогнозирует момент наступления динамического событий.

   Были систематизированы данные геологических изысканий Приморского разлома Байкальской рифтовой зоны – сформирован ГИС-проект, содержащий сведения о более 200 образцов горных город. Сформированный ГИС-проект состоит из геологического модуля, представленного геологическими и тектоническими картами разных масштабов (1: 200 000 и 1:1 000 000), геодезического модуля (содержит положение всех точек отбора проб в пространстве), геофизического модуля, к которому привязана таблица со значениями химического состава для каждого образца и петрографическим описанием. В ходе проведения фактических работ в 2021 г. было выполнено петрографическое изучение 214 шлифов горных пород, слагающих зону Приморского разлома, представленных как основным комплексом пород массива, так и вмещающими породами. Произведено определение химического состава этих пород. На основе сформированного банка геолого-геофизических данных разработан метод выделения структурных связей в осевой зоне регионального разлома. В основе метода лежит комплексирование результатов кластерного анализа химического состава образцов с их петрографическим описанием. На основе данных химического анализа составлялись 10-мерные вектора, каждая компонента которого соответствует доле определенного оксида. Далее строилось кластерное дерево на основе алгоритма минимального покрывающего дерева с использованием метрики Чебышева. Независимые ветви дерева формируют три кластера, каждый из которых имеет обобщенное петрографическое описание. Один из кластеров фактически характеризует зону контакта основных и кислых пород, наиболее сильно подверженных тектоническим преобразованиям, выраженным в петрографических характеристиках пород. Породы данного кластера являются наиболее деформированными среди всех рассматриваемых кластеров. Ширина зоны, соответствующая образцам данного кластера, составляет 15 м и формирует центральную зону разлома.

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

2021-2022 гг.:

1. Разработан новый подход к анализу пространственных закономерностей локализации сейсмичности, обусловленный именно строением разломной зоны. Были использованы данные высокоточного сейсмического каталога Северной Калифорнии. На первом этапе сейсмический каталог прошел процедуру декластеризации, т.е. были удалены все фор- и афтершоковые последовательности. На втором этапе для анализа пространственной структуры сейсмичности был усовершенствован алгоритм топологической фильтрации. Введение двух критических условий, основанных на «геологическом пониманиии» структуры разломов позволяет избежать экспертной оценки при реализации алгоритма и проводить локализацию, определять размер контактных пятен в автоматическом режиме. Вторым аспектом исследования высокоточных каталогов сейсмических событий был анализ статистического самоподобия контактных пятен (asperity) без попытки выделения отдельных иерархических уровней. При сделанных выше допущениях, можно определить эффективное количество «вложенных» asperities при изменении их радиуса. Для всех кластеров, сформированных в границах очагов основных землетрясений с M> 5, эффективное количество вложенных asperiities монотонно увеличивается с уменьшением радиуса по степенному закону: n(R)~R^(-d), d – показатель самооподобия.

2. Для понимания корректности полученных результатов в рамках выполнения первой задачи было проведено численное моделирование контактного взаимодействия жестких шероховатых поверхностей. Были построены фрактальные поверхности, заданные двумерной диапазонно-ограниченной функцией Вейерштрасса. Полученные в результатах численного моделирования данные о закономерностях локализации контактных пятен сравнивались с результатами выполнения задачи 1.

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

4. Для определения закономерностей деструкции литосферы в центральной зоне тектонического разлома был проведен статистический, топологический и фрактальный анализы микротрещиноватости. Для 300 шлифов по построенным цифровым образам были оценены фрактальная размерность трещиноватости, оценена энтропия Шенона и выявлена топология трещиноватости на основе классификации точек пересечений и узлов сети трещин. Большой объем полученных данных позволил провести корреляционный анализ и установить сильную связь между топологией трещиноватости и фрактальными свойствами ансамбля трещин. Работы по данному направлению требуют более детального рассмотрения на следующем этапе выполнения проекта.

5. Полевые работы в 2021 году проводились на западном побережье озера Байкал, вдоль Приморского разлома, который представлен уступом в рельефе с перепадом высот до 600м. Приморский разлом входит в Байкальскую рифтовую систему разломов и протягивается на северо-восток от устья р. Бугульдейка до деревни Зама на расстояние более 100 км. Петрографические исследования выполнены в 251 образце по единой структуре, в которой описываются петрографический тип горной породы, тип тектонитов, процент матрикса, преобладающие типы деформаций, основные составные минералы, наличие метаморфических и метасоматических изменений, флюидный привнос.

6. Серия экспериментов была проведена на уникальном экспериментальном стенде RAMA в лаборатории Деформационных процессов в земной коре ИДГ РАН. Этот стенд представляет собой установку двухостного нагружения. Модельный разлом представлял собой контакт шероховатых поверхностей блоков размером 750х120 мм, который был заполнен мелкодисперсным гранулированным материалом. Для выявления фазовых превращений был применен метод синхронного термического анализа (СТА). Метод основан на синхронной регистрации изменения массы вещества delta(m) и количества подводимого тепла delta(Q) при нагревании. Сравнение кривых нагревания образцов заполнителя до и после экспериментов позволяет делать выводы о возможных фазовых превращениях.

Публикации:

1. Гридин Г.А., Григорьева А.В., Ружич В.В., Остапчук А.А. Выделение структурных особенностей центральной части тектонического разлома Динамические процессы в геосферах, 2020. Выпуск 12. C.45-52. DOI: 10.26006/IDG.2020.51.17.018.

2. Остапчук А.А., Кочарян Г.Г., Морозова К.Г., Павлов Д.В., Гридин Г.А. Особенности формирования динамического сдвига в тонком слое гранулированного материала Физика Земли, 2021. №5, с.91-103 DOI: 10.31857/S0002333721050136

3. Ostapchuk A, Morozova K, Markov V, Pavlov D and Popov M (2021) Acoustic Emission Reveals Multiple Slip Modes on a Frictional Fault. Front. Earth Sci. 9:657487.  DOI: 10.3389/feart.2021.657487.

4. Ostapchuk A, Polyatykin V, Popov M and Kocharyan G (2022) Seismogenic Patches in a Tectonic Fault Interface. Front. Earth Sci. 10:904814. DOI: 10.3389/feart.2022.904814

5. Ustinov, S.; Ostapchuk, A.; Svecherevskiy, A.; Usachev, A.; Gridin, G.; Grigor’eva, A.; Nafigin, I. Prospects of Geoinformatics in Analyzing Spatial Heterogeneities of Microstructural Properties of a Tectonic Fault. Appl. Sci. 2022, 12, 2864. DOI: 10.3390/app12062864

6. Кочарян Г.Г., Остапчук А.А. Мезоструктура зоны скольжения тектонического разлома Физическая мезомеханика. 2022. 

при поддержке