Институт земной коры СО РАН
О динамике формирования крупных разломов литосферы (результаты физического моделирования)
А.В. Черемных
Аннотация
разлом литосфера растяжение деформационный
Обсуждаются результаты физического моделирования крупных разломных зон, формирующихся в упругопластичных моделях литосферы при сжатии, сдвиге и растяжении. Эксперименты проведены на установке «Разлом» с соблюдением теории подобия и при одинаковых исходных граничных условиях, что позволило сопоставить результаты, полученные с разных моделей. Для всех экспериментально воспроизведенных морфолого-генетических типов разломных зон выполнен анализ деформационной динамики составляющих их внутреннюю структуру систем разрывов. В развитие известных представлений о стадийности формирования разломных зон показано, что деформационная динамика внутриразломных разрывных нарушений в пределах стадий и в переходные временные интервалы между стадиями существенно отличается. Выявлена временная периодичность возникновения критических динамических состояний предразрушения систем разрывов, составляющих внутреннюю структуру разломных зон. Такие состояния зафиксированы как перед межстадийными структурными перестройками, так и в пределах стадий перед сменой этапов. При равномерной скорости деформирования в зонах крупных разломов литосферы каждые 2 или 4 млн лет должны возникать особые динамические состояния с формированием разрывных диссипативных структур (РДС). Масштабно-иерархическое строение разломных зон литосферы предполагает формирование менее крупных РДС за более короткие периоды. Установленные закономерности могут быть полезными при изучении сейсмического процесса, а также структур рудных полей и месторождений.
Ключевые слова: разломные зоны, структурные перестройки, периодичность, физическое моделирование, разрывные диссипативные структуры, подобие.
Abstract
About Dynamics of the Formation of Major Faults Lithosphere (Results of Physical Modeling)
A.V. Cheremnykh
Institute of the Earth's Crust SB RAS
The article discusses the results of physical modeling of major fault zones formed in the elasto-plastic models of the lithosphere under compression, shear and tension. The experiments were carried out at the «Fault» with respect to the theory of similarity and using the same initial boundary conditions that allow to compare results obtained with different models. For all experimentally reproduced morpho-genetic types of the fault zones the analysis of the deformation dynamics of the components of the inner structure of the systems ruptures. In the development of known ideas about the stages of formation of fault zones, it is shown that the deformation dynamics internal faults within the stages and in the transient time intervals between stages is significantly different. Revealed temporal periodicity of critical dynamical states prefracture systems break of the internal structure of fault zones. Such states are noted as before interstage structural rearrangements, and within stages, before the stage change. Fault dissipative structure (FDS) should be formed in major fault zones every 2 and 4 million years. Scale-hierarchic structure of fault zones in the lithosphere involves the formation of a smaller FDS for shorter periods. The installed regularities can be useful in the study of the seismic process and structures of ore fields and deposits.
Keywords: fault zones, structural adjustment, frequency, physical modeling, fault dissipative structure, similarity.
Введение
Известно, что формирование крупного разлома характеризуется определенной последовательностью развития его внутренней структуры: заложением многочисленных мелких разрывов, их последующим избирательным разрастанием и объединением в крупный магистральный шов [1-4; 911]. Исследование динамики структурных преобразований в зонах крупных разломов литосферы затруднено из-за большой продолжительности процесса, составляющей миллионы - десятки миллионов лет. Подобные трудности преодолеваются с помощью экспериментальных методов исследования, в частности физического моделирования.
Так, на основе физического моделирования установлено, что в процессе формирования разломной зоны выделяются четыре последовательно сменяющие друг друга стадии [9; 13]: пликативная, ранняя дизъюнктивная, поздняя дизъюнктивная и стадия полного разрушения. Пликативная стадия характеризуется деформацией модельного материала в зоне будущего разлома без образования видимых разрывов. На ранней дизъюнктивной стадии в пределах широкой зоны разлома наблюдается образование и дальнейший рост не взаимодействующих между собой мелких опережающих основной разлом разрывов. Поздняя дизъюнктивная стадия наступает после слияния мелких разрывов в более протяженные опережающие дизъюнктивные дислокации, включая сегменты будущего магистрального сместителя. На стадии полного разрушения деформации возникают только в непосредственной близости от магистрального сместителя, так как связаны с перемещениями крыльев данного дизъюнктива. На этой стадии формируются разрывы, оперяющие основной шов, что связано с волнистостью сместителя. Переход от стадии к стадии происходит быстро относительно длительности формирования всей зоны разлома и в связи со значительным изменением внутреннего строения сети разрывов в пределах зоны назван «структурной перестройкой» [9; 13].
Структурные перестройки в зоне разлома могут быть выявлены на основе анализа количественных параметров участвующих в его внутреннем строении разрывов: длина максимального разрыва, плотность разрывов, амплитуда смещения по разрывам, фрактальная размерность сети разрывов и др. [1-2; 4; 8; 9-11; 13; 15]. Как правило, в моделях выделяется несколько структурных перестроек, главная из которых соответствует переходу между ранней и поздней дизъюнктивными стадиями и выражается весьма контрастным изменением перечисленных выше параметров [13].
Целью проведенного нами экспериментального исследования являлось дальнейшее развитие представлений о деформационной динамике внутренней разрывной структуры зон крупных разломов разных морфолого - генетических типов. Все эксперименты выполнены с соблюдением условий подобия при одинаковых исходных граничных условиях: вязкости модельного материала, мощности модели и скорости ее деформирования.
Основная часть
Методика моделирования и сбора количественных показателей
Моделирование выполнено в лаборатории тектонофизики Института земной коры СО РАН (г. Иркутск). В разные годы проведено семь экспериментов, отличающихся направлением относительного смещения штампов экспериментальной установки при одинаковых остальных параметрах моделей. Три эксперимента проведены для исследования динамики формирования зоны сдвига. Зоны сжатия литосферы изучены в двух опытах. Также в двух опытах воспроизведено формирование зон растяжения литосферы. Во всех опытах относительное перемещение штампов было горизонтальным и отличалось только направлением вектора подвижного штампа относительно неподвижного (рис. 1).
Рис. 1. Схемы моделирования зон сдвига (А), сжатия (Б) и растяжения (В). 1 - модельный материал; 2 - штампы экспериментальной установки; 3 - направление перемещения подвижного штампа; 4 - разрывы
Эксперименты проведены по стандартной методике моделирования на установке «Разлом» [9-11]. Глинистая паста с вязкостью 105 Па-с размещалась на рифленой (для лучшего сцепления) поверхности штампов экспериментальной установки в виде равномерного слоя мощностью 70 мм. Затем один из штампов перемещался с постоянной скоростью, равной 5 мм в минуту. При неизменной длине модели ее ширина оставалась относительно постоянной лишь в условиях сдвига, а при сжатии и растяжении, соответственно, уменьшалась или увеличивалась. Согласно теории подобия [16] при принятых масштабных коэффициентах вязкости, плотности, линейных размеров и времени: 1 мм на модели соответствует 1 км в ее природном аналоге, а 1 минута эксперимента эквивалентна ~1 млн лет природного процесса разломообразования. При таких исходных параметрах модельная скорость смещения 5 мм/мин соответствует природной скорости смещения крыла разлома 5 мм/год.
Формирование магистрального разлома в моделях при заданных параметрах происходит за 17-26 минут в зависимости от способа приложения нагрузки. За это время со средним интервалом в 30-40 секунд производилась фоторегистрация сети разрывов на поверхности моделей (рис. 2). При компьютерной обработке фотографий все имеющиеся на них разрывы переводились в цифровую форму. Полученная таким образом база данных отражала пространственно-временное распределение разрывов в формирующихся разломных зонах. Для количественной оценки изменений, происходящих в пределах моделируемых зон, анализировались длины разрывов, а также плотность и фрактальная размерность сетей разрывов. Плотность разрывов подсчитывалась круглой палеткой с диаметром, равным половине толщины модели, а центр палетки располагался непосредственно в районе формирования магистрального шва. Определение фрактальных размерностей производилось по общей методике для объектов на плоскости [17] - методом подсчета квадратов.
Рис. 2. Фотографии сети разрывов разломных зон в моделях. А - зона сдвига, Б - зона сжатия, В-зона растяжения
Для анализа деформаций на разных стадиях развития зоны, включая пластическую, на поверхность модели были нанесены маркеры в виде окружностей диаметром 14 мм. В моделях зон растяжения и сжатия литосферы такие окружности были нанесены перпендикулярно простиранию формирующихся разломов, а в опытах простого сдвига создавалась равномерная сеть маркеров как поперек, так и вдоль разломной зоны.
Следует отметить, что, в отличие от единой сдвиговой деструктивной зоны, при растяжении и сжатии в моделях формируются две крупные раз - ломные зоны: одна со стороны подвижного штампа, другая - неподвижного (см. рис. 1). При сопоставлении динамики развития зон различного кинематического типа для растяжения и сжатия были использованы данные по сети разрывов со стороны неподвижного штампа, где процесс происходит быстрее и более контрастно выражен во внутренней структуре зоны.
Результаты экспериментов
Анализ полученных с фотографий количественных параметров разрывов в моделях разломных зон литосферы позволил установить следующую специфику динамики их формирования. Временные вариации фрактальной размерности сетей разрывов в зонах растяжения, сжатия и сдвига показали скачкообразный характер графиков (рис. 3). Увеличению фрактальной размерности сети разрывов в разломной зоне, как правило, предшествует относительное постоянство данного параметра. При растяжении завершение стадии относительного постоянства наблюдается на 720-й и 841-й секундах от начала опыта (см. рис. 3, А). По двум значениям невозможно установить периодичность процесса. Изменение фрактальной размерности сети разрывов в зоне сдвига также не позволяет надежно выявить период, так как на графике фиксируется лишь одно контрастное изменение параметра, соответствующее главной структурной перестройке (см. рис. 3, Б). Имеющиеся изменения графика до и после этого момента сложны для однозначной интерпретации. Достоверно установить периодичность увеличения фрактальной размерности сети разрывов удалось для зоны сжатия (см. рис. 3, В). На графике достаточно определенно выделяются моменты, предшествующие увеличению данного параметра - 695, 890, 1090 и 1290 секунд. Таким образом, период составил 200 секунд.
Рис. 3. Графики изменения фрактальной размерности сети разрывов в моделях разломных зон растяжения (А), сдвига (Б) и сжатия (В).
Крупными значками на графике отмечены моменты времени, предшествующие контрастному увеличению параметра
Кроме того, анализ длин разрывов во времени выявил еще одну характерную особенность. Известно, что деструктивный процесс в зоне формирующегося разлома во времени идет по пути увеличения длин разрывов с уменьшением их количества [9]. Проведенные нами эксперименты показали, что перед структурными перестройками такой ход процесса нарушается. Фиксируется увеличение количества мелких разрывов как за счет образования новых, так и за счет активизации разрывов, ранее прекративших свое развитие и находившихся в пассивном состоянии в разрывной структуре зоны разлома. Во временной последовательности гистограмм, отражающих распределение внутриразломных разрывных нарушений по длине во всех трех типах разломных зон, в различные моменты времени такие ситуации выявляются по характерному строению гистограммы (рис. 4). В частности, они возникали в зоне сдвига на 660-й, 900-й и 1141-й секундах (см. рис. 4, А), в зоне сжатия на 890, 1090 и 1290 секундах (см. рис. 4, Б) и в зоне растяжения на 720-й и 920-й секундах (см. рис. 4, В). Видно, что на гистограммах, характеризующих эти временные моменты, сеть разрывов в разломных зонах подчиняется закономерной масштабной иерархии, а именно: количество разрывов в каждом последующем интервале длины меньше, чем в предыдущем. Отмеченные особенности деформационной эволюции крупных раз - ломных зон трактуются нами с позиции синергетики и связываются с периодическим проявлением процессов самоорганизации внутриразломных разрывных систем. Совокупность разноранговых разрывов, вовлеченных в самоорганизацию, мы называем разрывными диссипативными структурами (РДС) [3]. В проведенных экспериментах РДС образуются через 200 секунд в зонах сжатия и растяжения и через 240 секунд в зонах сдвига (рис. 4). При этом в зонах сдвига РДС образуются уже на 11-й минуте эксперимента, в зонах растяжения - на 12-й минуте, а в зонах сжатия - лишь на 17-й минуте опыта.
Рис. 4. Гистограммы длин разрывов в моделях разломных зон.
А - зона сдвига, Б - зона сжатия, В-зона растяжения. Цветом на гистограммах выделены моменты образования в моделях разрывных диссипативных структур
Отмеченная периодичность возникновения РДС проявляется и во временных вариациях деформационных параметров кольцевых маркеров и плотности разрывов в центральной части сдвиговой зоны. В процессе деформации кольцо превращается в эллипс и испытывает вращательное движение в направлении смещения активного крыла разлома. С моделей снимались замеры длинной оси эллипса и угол между этой осью и простиранием зоны сдвига (рис. 5). Из приведенных графиков видно характерное изменение всех параметров через каждые 240 секунд. Эти моменты отмечены на графиках пунктирной линией. Как правило, до или после этих моментов фиксируются резкие изменения угла между длинной осью эллипса и простиранием зоны (см. рис. 5, А), длины этой оси (рис. 5, Б) и плотности разрывов (рис. 5, В). Кроме этого интервала, дополнительно намечается 120 секундный интервал. Как отмечено выше, первым временным интервалом связано возникновение РДС перед межстадийными структурными перестройками. В работах [4; 14] показано, что, так же; как весь процесс формирования зоны разлома подразделяется на стадии, каждая стадия в свою очередь подразделяется на этапы. Возможно, выделяющийся временной интервал в 120 секунд характеризует продолжительность этапов внутри стадии. Смене таких этапов предшествует появление локальных РДС, распад которых сопровождается быстрыми структурными преобразованиями разрывной сети в локальном объеме зоны разлома.