Волновая природа активизации разломов Центральной Азии на базе сейсмического мониторинга
Аннотация
Возникновение очагов землетрясений в сейсмоактивных зонах литосферы характеризуется квазихаотичностью. Вариации пространственно-временной миграции очагов в реальном (месяцы, годы) времени не могут быть объяснены изменениями напряженно-деформированного состояния среды эндогенной природы. В статье во главу угла поставлен один из определяющих структурных факторов контролирования очагов землетрясений - разломы различных рангов, современная активизация которых и сейсмичность парагенетически взаимосвязаны. Под этим углом зрения по разработанной методике изучены пространственно-временные закономерности современной активизации разломов, ее источники и механизмы реализации. Показано, что активизация разломов и ее относительно высокая частота в масштабах реального времени вызываются медленными деформационными волнами возбуждения, генераторами которых могут быть межплитные и межблоковые подвижки хрупкой литосферы. По скорости прохождения деформационных волн возбуждения активные разломы могут классифицироваться на группы, разнящиеся по геолого-геофизическим параметрам. Они позволяют оценивать направление фронтов деформационных волн возбуждения и выделять области преимущественной активизации разломов в интервалах реального (с геологической точки зрения мгновенного) времени.
Результаты исследований современной геодинамики литосферы позволили выделить группы активных разломов, интенсивность активизаций которых изменяется в интервалах короткого реального времени и практически не зависит от функционирующих геодинамических режимов регионов. Установлено, что разломы активизируются с изменяющейся интенсивностью и активизация происходит чаще, чем фиксируются изменения в тектоническом режиме и региональном поле напряжений. Вывод получен благодаря введению новых параметров численных оценок количественной и энергетической (магнитудной) характеристик разломов. Для оценки интенсивности активизации разломов в реальном времени (месяцы, годы) предложено использовать их количественный индекс сейсмической активности (КИСА) оn (км-1), под которым понимается число сейсмических событий n определенных энергетических классов K, приходящихся на единицу длины разлома L (км) при принятой ширине области его динамического влияния M (км) за заданный промежуток времени t (годы) [1]:
оn= ?n(M, K, t)/L (1)
Ширина области динамического влияния разлома М определяется по уравнению
М=bL (2)
где L - длина разломов, км; b - коэффициент пропорциональности, зависящий от L и по эмпирическим данным изменяющийся от 0.03 до 0.09 соответственно для трансрегиональных и локальных разломов [2]. В реальном масштабе времени именно КИСА характеризует нестабильность крыльев разломов и даёт основание для анализа доли участия разнорангового разломного сообщества в сейсмическом процессе.
Для определения энергетического потенциала разломов предложено применять магнитудный (энергетический) индекс сейсмической активности разломов оk , под которым понимается значение класса максимального сейсмического события Kmax (K = lg E, Дж), приходящегося на длину разлома L (км) при принятой ширине области его динамического влияния М (км). Магнитудный индекс сейсмической активности оценивается по выражению:
оk = Kmax(t) (М, K, t) (3)
где Kmax(t) - максимальный класс землетрясения (или его максимальная магнитуда) в области динамического влияния разлома М за заданный промежуток времени t [3].
Рассмотрение вариаций количественного и магнитудного индексов сейсмической активности на примерах разломной тектоники Центральной Азии свидетельствует об отсутствии ясно выраженной пространственной закономерности в активизации территориально сближенных ансамблей разломов в чрезвычайно короткие интервалы реального времени. Создается кажущееся впечатление, что активизация разломов в границах сейсмоактивных зон происходит хаотично (рис. 1, 2). Эндогенные источники развития разломов и генетически связанной с ними сейсмичности в задаваемые, с геологической точки зрения мгновенные, интервалы времени остаются пространственно и энергетически стабильными. Следовательно, пространственно-временные закономерности в возбуждении активизации разломов следует искать в энергетически слабых, но достаточных для нарушения метастабильного состояния разломно-блоковой среды литосферы, триггерных механизмах.
разлом сейсмичность межблоковый литосфера
Рис. 1. Карта активных разломов Центральной Азии по магнитудному индексу сейсмической активности (3): 1 - разломы с оk ? 12 (весьма активные); 2 - разломы с оk =10-11 (активные); 3 - разломы с оk ? 9 (слабоактивные); 4 - положение сечений на карте и их номера; 5 - номера разломов по каталогу.
Рис. 2. Периодичность изменений магнитудного индекса сейсмической активности по событиям энергетических классов 10 и 12 на примере двух удаленных друг от друга разломов Центральной Азии.
Исследуемые вопросы и методы их решения. Для исследования интегрированной пространственной выраженности активизации разноранговых разломов и поисков ее триггерных механизмов изучена тенденция временной направленности возникновения эпицентров землетрясений в зонах динамического влияния сейсмоактивных разломов.
В последние годы на примерах тектонически активных зон Центральной Азии проведено тщательное изучение сейсмоактивных разломов по данным сейсмологического мониторинга. Критерием безусловной активизации разломов является приуроченность к ним эпицентров землетрясений. Разломы растут при каждом новом сейсмическом акте. М.В. Гзовский [4] показал, что максимальный их рост за один сейсмический акт не превышает полутора десятка метров. Опубликованные в последние годы работы по изучению последовательности вспарывания разрывов и смещений по ним при землетрясениях выявили их определенную тенденцию в одном из направлений. Рассмотрим некоторые детали локализации смещений в разрывах во времени при постоянном и импульсивном нагружении.
В экспериментальной работе В.А. Санькова и К.Ж. Семинского [5,6] было убедительно показано неравномерное смещение по простиранию в разрывах, формирующихся в сдвиговом поле напряжений. Эти, трудно воспринимаемые среди геологов в те годы знания, в последующем были подтверждены наблюдениями в эпицентральных зонах землетрясений. По этому вопросу весьма аргументированные данные получены при крупномасштабных исследованиях вариаций смещений земной поверхности вдоль сейсмоактивного разлома Emerson в его различных частях (рис. 3). Они сопоставлены с еще более детальными наблюдениями по разлому Garlock (Южная Калифорния) [7]. Закономерности между амплитудами горизонтальных и вертикальных смещений и длиной разрывов по простиранию намечаются только по огибающей кривой [7]. Таким образом, в каждом конкретном случае, связь между длинами разломов, амплитудами смещений вдоль них и магнитудами землетрясений во многом зависит от расположения точек наших фактических наблюдений в природной обстановке и замеров соответствующих параметров вдоль простирания активных разломов.
Рис. 3. Вариации горизонтальных (а) и вертикальных (б) амплитуд смещений по простиранию сейсмоактивного разлома [7]. Условные обозначения: линия - средние значения; точки - единичные замеры сдвига; под кривыми графиков - изменение состава горных пород по простиранию разлома.
По разным группам наблюдений развитие и активизация сейсмоактивных разломов происходят преимущественно в одном из двух от эпицентра (очага) землетрясения направлений [8, 9; и др.]. Богатый фактический материал дали проведенные в последние годы обобщения по расположению максимальных амплитуд смещений по простиранию разрывов и эпицентров землетрясений. Авторы [10] показали, что положение максимальных смещений в сейсмоактивных разломах по отношению к эпицентрам основных толчков не совпадает с локализацией последних, а смещения максимальны на некотором удалении по простиранию активного разрыва (рис. 4). Известно, что эпицентры последующих относительно сильных землетрясений в конкретно исследуемых разрывах происходят в местах максимальных смещений, связанных с предшествующими событиями. При этом сила нового события в определенной мере пропорциональна и длине разрыва, и амплитуде смещения [11]. Отсюда, что наиболее вероятно, следующий сейсмический акт локализуется скорее всего в точке разлома с максимальным смещением, за ним во времени следующий и т.д. (рис. 5). Будет намечаться тенденция в направлении расположения последующих очагов землетрясений и, следовательно, будут определяться и преимущественные направления прорастания разрывов при их активизации.
Рис. 4. Расположение главных очагов землетрясений и максимальных амплитуд горизонтальных смещений по простиранию активных разломов [10]: 1 - Куньлунь (2001 г., Mw = 7.8); 2 - Маный (1997 г., Mw = 7.9); 3 - Дюздже (1999 г., Mw = 7.2).
Рис. 5. Две модели разрастания активных разломов и накопление амплитуд смещений во времени [10]: одновозрастные (а) и разновозрастные (б) сегменты разлома и их синхронная активизация. Стрелкой показано направление разрастания разлома. 1, 2, 3, …, n - последовательность сейсмических актов.
В основе рассматриваемых далее авторских построений лежит представление о том, что землетрясение любого класса фиксирует нарушение равновесия в зоне разлома, сопровождающееся увеличением интенсивности трещиноватости и при сильных событиях смещением крыльев. Частота сейсмических событий в зоне разлома отражает интенсивность его активизаций, а тенденция в пространственной направленности очагов вдоль оси разлома во времени воспроизводит скорость и вектор движений триггерного источника активизаций. При этом в зоне разлома происходит реализация второго [12] механизма развития крупной трещины: её продолжающееся формирование идет по предварительно уже существующей перколяционной сети более мелких трещин и скорость дискретного развития (активизации) трещины по простиранию может быть исключительно низкой, исчисляемой годами, столетиями или более продолжительным временем.
Для выяснения тенденций векторной направленности возбуждений по наиболее активным разломам были построены индивидуальные графики, на оси абсцисс которых откладывались длины разломов с соответствующими положениями эпицентров землетрясений, на оси ординат ? время этих событий (рис. 6). По физическому смыслу каждый из графиков отражает два новых дополнительных параметра разломов: наклон его линии воспроизводит вектор пространственно-временного движения очагов землетрясений вдоль разрыва (с левого фланга разлома на правый или наоборот), а тангенс угла ее наклона к оси ординат - средние скорости пространственного распространения волны возмущения, стимулирующей возникновение очагов землетрясений.
Рис. 6. Определение вектора разрастания разлома по простиранию при активизации и оценки средней скорости движения волны возбуждения.
Обсуждение результатов. По предложенной методике, используя каталог землетрясений Байкальского филиала Геофизической службы СО РАН за 1960-2000 гг., рассмотрена специфика пространственно-временного распространения очагов землетрясений по простиранию активных разломов Центральной Азии (рис. 7).
Выделено более ста разноранговых разломов с зафиксированными в областях своего динамического влияния очагами землетрясений 12-16 классов. Построены индивидуальные графики «время события - пространство» и проанализированы временные тренды сейсмических событий по отдельным активным разломам. Эпицентры землетрясений конкретных разломов на графиках образуют системы параллельных прямых, как если бы вдоль соответствующих разломов распространялись с постоянной средней скоростью серии деформационных волновых возмущений, инициирующих сейсмические события - активизации разломов. Всего таких систем параллельных прямых с равными, но в разных направлениях углами наклона выделено 7 групп, каждая из которых соответствуют одинаковым средним скоростям деформационных возмущений. Внутри группы разрывы подразделяются на две подгруппы, соответствующие разным (противоположным) векторам возмущений (рис. 8). Группировка разломов по критерию одинаковых скоростей деформационных волн возмущений свидетельствует об идентичных параметрах их активизации (табл. 1).
Дополнительно, используя известное соотношение
л = Vt (4)
по характерному времени возникновения событий t и их средней скорости V в каждой из групп оценена вероятная длина деформационных волн возмущения л.
Принимая во внимание недостаточную обеспеченность некоторых разломов минимальным количеством исходных данных по сейсмичности, из дальнейшего анализа исключены последние три группы (табл. 1).