УСЛОВИЯ ТЕКТОНИЧЕСКОГО НАГРУЖЕНИЯ И ПРОГНОЗ ПЕРСПЕКТИВНЫХ УЧАСТКОВ МЕСТОРОЖДЕНИЯ ТАЛОВЕЙС
←
→
Транскрипция содержимого страницы
Если ваш браузер не отображает страницу правильно, пожалуйста, читайте содержимое страницы ниже
Н. Ю. Васильев*, А. О. Мострюков**, В. А. Сунцов*** УСЛОВИЯ ТЕКТОНИЧЕСКОГО НАГРУЖЕНИЯ И ПРОГНОЗ ПЕРСПЕКТИВНЫХ УЧАСТКОВ МЕСТОРОЖДЕНИЯ ТАЛОВЕЙС Введение пространственно-временного взаимодействия процес- сов тектогенеза и рудогенеза и в соответствии с этим В статье излагаются результаты исследования, достоверности развиваемых представлений о рудоло- представленные авторами на заседании Рабочей кализующих условиях тектонического нагружения гео- комиссии международного проекта «FENGOT» в логических структур (Васильев, Мострюков, 2007). п. Косалма (Карелия, 25.08.2011). В тексте рассмот- Второй вывод – об эффективности применения фазо- рены процедуры реконструкции механизмов дефор- во-циклической модели развития тектонических де- мирования структуры золоторудного объекта (Кос- формаций (Васильев, Мострюков, 2001) к анализу ус- томукшский район Карелии, см. врезку на рис. 4) и ловий развития локальных обстановок рудогенеза не проанализированы результаты реконструкций. По только на фанерозойских, но и на докембрийских эта- соотношению параметров тектонического поля на- пах (циклах) геологического развития Земли. пряжений выявлены особенности условий синруд- Результаты работы подтверждены данными гео- ного нагружения объекта и локализованы участки химической съемки и поискового бурения. их дискретного проявления в структуре месторож- дения. К обоснованию прогноза обстановок рудогенеза Практическая цель работы – выделение наиболее по реконструкции механизмов перспективных участков на одной из лицензирован- тектонических деформаций ных площадей (лицензии № ПТЗ 01278 БР от 26.10.2005 и ПТЗ 14813 БП от 11.12.2009) и опреде- В ве де н ие в пр о б ле м у ление рекомендаций для проведения поисково-оце- Одна из фундаментальных закономерностей ге- ночных работ. Работа выполнена по заказу и при фи- незиса рудных сегрегаций (лат. segregatio – отделе- нансовой поддержке ООО «Карелгеоресурс» и ООО ние): гнезд, жил, линз, рудных тел, залежей – опре- «Костомукшское золото». деляется дискретностью их проявления в иерархии Исследование основано на реконструкциях текто- структурно-вещественных неоднородностей место- нофизических условий развития сдвиговых деформа- рождений, рудных полей, узлов и т. п. Устойчивый ций в структуре указанного месторождения, включая шаг дискретности, «контролирующий» локализа- и его раму. По результатам реконструкций в истории цию сегрегаций на системно ориентированных гра- развития месторождения выделено четыре последо- ницах этих неоднородностей (на берегах разрывов, вательных этапа (фазы) изменения условий его тек- трещин, разломов, тектонизированных контактах тонического нагружения (механизмов деформирова- горных пород) и конформность сегрегаций грани- ния). Установлена аналогичность синрудных усло- цам этих неоднородностей, – хорошо известные эм- вий тектонического нагружения данного объекта, пирические закономерности. Они давно использу- сформированного в докембрии, с выявленными ра- ются с учетом принципа иерархичности в практике нее в других регионах (Васильев, Мострюков, 2007) поисков и разведки месторождений. В основе этих условиями синрудного нагружения мезозойского и закономерностей лежат бесспорные факты либо па- кайнозойского рудоносных объектов. Из данной ана- рагенетических, либо генетических соотношений, логии следуют два принципиально важных вывода. связывающих процессы тектоники и рудообразова- Первый – об эффективности анализа условий де- ния. Для развиваемого подхода к прогнозу обстано- формирования земной коры для прогноза обстановок вок рудогенеза чрезвычайно важны свидетельства именно генетической связи между этими процесса- –––––––––––––––––––––––– ми. Результаты реконструкций параметров тектони- * РГГРУ-МГРИ им. С. Орджоникидзе. ** Геофизическая обсерватория Борок ИФЗ РАН. ческих полей напряжений на эндогенных месторож- *** ООО Карелгеоресурс. дениях свидетельствуют о реальности такой связи 77
(Васильев и др., 1999; Васильев, Мострюков, 2000, С тр ук т ур но - к и не м а ти ч е с к и й м е то д р е ко н - 2001, 2007). ст р ук ц и и те к то ни че ск и х по ле й н а пр я ж е ни й Тематические исследования, посвященные ана- Метод основан на доказанной Д. П. Маккензи лизу соотношений между процессами рудогенеза и (McKenzie Dan P., 1969) закономерной причинно- тектогенеза, начаты в нашей стране в 30–40 гг. про- следственной зависимости направлений сдвиговых шлого века работами Ф. И. Вольфсона, А. В. Коро- смещений (τn) в горных породах от направлений, воз- лева, А. В. Пэка, В. М. Крейтера. Результаты иссле- действующих в этот момент на породы, главных нор- дований этих ученых и их последователей воплоще- мальных напряжений (σ1,σ2,σ3). Процедурами метода ны в подробнейших классификациях структур руд- предусмотрено использование данной зависимости ных полей и месторождений, обоснованных как для решения «обратной задачи» – нахождения харак- теоретическими разработками, так и эксперимен- теристик тектонических напряжений по совокупности тальными данными. данных о кинематике сдвиговых смещений. Результа- Вместе с тем по мере появления все более новых том реконструкции является численная информация о и детальных классификаций все острее становилась параметрах тектонических полей напряжений, необ- потребность в более информативных – количест- ходимая для «количественной» характеристики разно- венных методах анализа деформационных процес- образных кинематических (сбросовых, сдвиговых, сов тектогенеза для численных оценок изменчиво- взбросовых и «комбинированных») механизмов де- сти параметров этих процессов в истории разви- формирования земной коры. Это информация об ори- тия тектонических деформаций. Другими словами, ентации осей и плоскостей действия трех главных явно требовались такие методы анализа, чтобы вы- нормальных и двух максимальных касательных на- воды, основанные на их результатах, обеспечивали пряжений (соответственно, σ1, σ2, σ3 и τmлев, τmправ) в бы возможность рассматривать формирование ме- трехмерном пространстве, а также о величине относи- сторождения «…как результат развития процессов тельных «растяжения-сжатия» земной коры, фикси- рудогенеза на фоне закономерной последователь- руемой по численному значению коэффициента Лодэ- ной смены геодинамических обстановок в истории Надаи (–1 ≤ μσ ≤ +1), (–1 ≤ με ≤ +1). Надо отметить, что геологического развития…» (Рундквист, 1995, с. 118). практически одновременно с О. И. Гущенко француз- Представляется, что изменение структуры тектони- ским исследователем Ж. Анжелье (Angelier, 1975 и ческого поля напряжений (условий тектонического др.) был независимо разработан близкий метод рекон- нагружения земной коры) следует полагать одним струкции тектонических полей напряжений, основан- из «необходимых» событий смены геодинамических ный на той же закономерной генетической связи обстановок. сдвигов с главными нормальными напряжениями. К середине прошлого века в связи с интенсифи- Преимущество структурно-кинематического ме- кацией геологических исследований уже не только тода анализа определяется, в первую очередь, тем, перед рудной отраслью, но и в целом перед геоло- что в нем для решения задачи о кинематических ус- гией возникла проблема более глубокого и деталь- ловиях (механизмах) деформирования земной коры к ного исследования процессов тектогенеза. Соот- реконструкциям привлекается кинематическая ин- ветственно, стала реальной необходимость в со- формация о направлениях тектонических движений. вершенствовании количественных методов анали- Т. е. исходной информацией для метода являются за тектонических структур. Именно тогда и начало данные не только о пространственной ориентировке формироваться новое научное направление иссле- трещин, но и о направлениях (векторах) относитель- дований, базирующееся на численных методах ного перемещения плоскостей (берегов) этих тре- изучения тектонических деформаций земной коры щин в трехмерном пространстве. С единой кинема- с использованием достижений физики и математи- тической природой и исходной информации, и ре- ки. Выдающаяся заслуга в развитии таких исследо- конструируемых параметров деформационного про- ваний в нашей стране, без сомнения, принадлежит цесса, вероятно, следует связывать подтверждаемые М. В. Гзовскому, многогранная исследовательская практикой объективность, надежность и репрезента- деятельность которого определила рождение и раз- тивность результатов структурно-кинематического витие тектонофизики. Одной из главнейших задач метода анализа тектонических полей напряжений. тектонофизики он считал реконструкцию и моде- Исходная информация о кинематике тектониче- лирование (на основе принципов подобия) текто- ских движений, фиксируемая по «штрихам и бороз- нических полей напряжений с целью выявления дам скольжения», – это данные о знаках («правых» и информации о численных характеристиках, изме- «левых») и направлениях векторов сдвиговых (сколо- няющихся во времени и пространстве, механизмов вых) смещений по трещинам и разрывам горных по- деформирования земной коры. К важнейшим ре- род. Пример записи «полевой» информации о кине- зультатам работы научной школы М. В. Гзовского матике тектонических движений дан на рис. 1. Эту в данном направлении относится создание струк- информацию собирают in situ («на месте» лат.) – в турно-кинематического метода реконструкций естественных обнажениях горных пород и в стенках тектонических полей напряжений (Гущенко, 1973, горных выработок и «привязывают» с помощью на- 1979 и др.). вигатора GPS к координатам топографической сети. 78
Рис. 1. Пример записи полевой информации В настоящее время структурно-кинематический как будет показано ниже, процедура прогноза обста- метод реконструкции тектонических полей напряже- новок рудогенеза. ний осуществляется на основе компьютерных техно- Наиболее ярко эффективность алгоритма струк- логий. Разработан алгоритм и создан пакет про- турно-кинематического метода анализа проявляет- грамм, обеспечивающий: анализ разномасштабных ся в расчетах, проводимых на основе пространст- выборочных совокупностей полевых данных (в зави- венно-временной сепарации данных о кинематике симости от иерархического ранга изучаемого объек- сколовых смещений. В результате для последова- та), высокую скорость проведения аналитических тельного ряда тектонических полей напряжений, расчетов и надежность их результатов. В результате инверсионно сменяющих друг друга, рассчитыва- обработки первичных данных получают численные ются численные характеристики параметров этих значения одиннадцати (11) тектонофизических пара- полей. Объективным свидетельствам достоверно- метров, оптимально характеризующих механизмы сти таких временных рядов являются известные деформирования земной коры (±μσ, ±με, σ1, σ2, σ3, ε1, факты разнонаправленных борозд скольжения, ε2, ε3, τmлев, τmправ, ±∆Ζε ; пояснения к ним см. выше). фиксируемых на «одной плоскости» трещины или Определяемые по форме и пространственной ориен- разрыва (см. рис. 1). тировке эллипсоидов тектонических напряжений и Таким образом, структурно-кинематический ме- деформаций в трехмерном пространстве численные тод реконструкции тектонических полей напряжений значения этих параметров, а точнее, их комбинации стал эффективным инструментом выделения инфор- являются наиболее надежной и объективной инфор- мации о закономерностях последовательного измене- мацией, детально характеризующей особенности ния механизмов деформирования земной коры в про- развития (изменения) условий деформирования гео- цессе тектогенеза (Васильев и др., 1999; Гущенко, логических структур и во времени, и в пространстве. 1999; Васильев, Мострюков, 2000, 2001). Метод ап- Собственно, комплексом численных значений этих робирован на объектах различных иерархических параметров и характеризуются особенности тектоно- уровней (рангов): например, в платиноносном Галь- физических условий развития деформационного про- моэнанском ультраосновном массиве Корякского на- цесса в массивах горных пород и обосновывается, горья (Васильев и др., 1999; Васильев, Мострюков, 79
2000), в Байкальском рифте (Васильев, Мострюков, ных напряжений (от –1 до +1) на границах фаз. Не- 2002), в Эльджуртинском гранитном массиве Боль- прерывность развития процесса деформации выра- шого Кавказа (Vasilyev et al., 2002). жается в унаследованном сохранении ориентации од- ной из осей главных напряжений (σ1, σ3) эллипсоида Ци к л ы и фа зы те к то н и че ск и х де ф о р м а ци й и сменой индекса другой оси на индекс (σ2) проме- жуточного напряжения в «соседних» фазах (рис. 2). Одна из важнейших закономерностей тектогенеза Плавность переходов между фазами подчеркивается проявляется в последовательной смене кинематиче- выравниванием величин между одним из главных ских механизмов (взбросовых, сдвиговых, сбросовых) напряжений (σ1, σ3) и промежуточным напряжением деформирования земной коры, контролирующих эта- пы и стадии развития региональных и локальных тек- (σ2) в эллипсоидах напряжений на границах фаз. Эти тонических структур. Ведущим фактором смены ме- равенства (σ3↔σ2, σ1↔σ2) являются условием появ- ханизмов процесса является непрерывно-прерыви- ления «точек инверсии» в компонентах девиатора стое, эволюционное развитие знакопеременных текто- напряжений и, соответственно, реализации инверси- нических «импульсов», вероятно отражающих сим- онных механизмов смены фаз, когда деформируемый метричную флуктуацию направленного инверсионно- объем земной коры находится в состоянии одноосно- го изменения энергетики эндогенных процессов отно- го тектонического нагружения (или «максимального сительно разномасштабных поверхностей раздела в растяжения», или «максимального сжатия»). структурно-вещественных комплексах земной коры. Итак, подчеркнем главные признаки циклическо- Скорее всего, именно с пространственно-временными го развития условий деформирования земных недр. изменениями тектонических «импульсов» следует В течение каждой фазы цикла происходит направ- связывать квазициклические изменения структуры ленное изменение формы эллипсоида деформации тектонического поля напряжений в трехмерном про- (от με = –1 до με = +1 или наоборот), а в течение пол- странстве (X, Y, Z) – на последовательных фазах цик- ного цикла – закономерный процесс шестикратно ла деформации (рис. 2, Васильев, Мострюков, 2001). проявленной инверсионной смены ориентации осей Проявление цикла (χυχλοζ – круг, кругооборот) главных нормальных напряжений (σ1↔σ2, σ3↔σ2). деформации контролируется условием длиннопери- Природу циклического развития (цикла) тектони- одной устойчивой ориентации главных осей тензора ческих деформаций, фиксируемого по шестикратной напряжений (Tσцикл = const), а фазы (φατιζ – проявле- инверсионной смене фазовых механизмов деформа- ние, появление) деформации – условием короткопе- ционного процесса (рис. 2), вероятно, следует связы- риодной устойчивой ориентации осей главных нор- вать с проявлением, как минимум, двух «противопо- мальных напряжений (׀σ1фаз, σ2фаз, σ3фаз = const) тензо- ложно» направленных тенденций в закономерностях ра во времени и пространстве относительно коорди- развития земной коры. С одной стороны, это после- нат X, Y, Z (рис. 2). Смена фаз в течение цикла де- довательное развитие силового (энергетического) формации обусловлена непрерывным изменением импульса, направленного на нарушение изостатиче- характеристик тектонического «импульса» – соотно- ского равновесия блоковой структуры земной коры, шений между величинами главных нормальных на- а с другой стороны, – высокая инерционность пряжений (σ1, σ2, σ3) в интервалах значений от (+1, («инерция массы») направленного развития тектони- условий максимального «сжатия») до (–1, условий ческих структур, препятствующая «мгновенному» максимального «растяжения»). Предельными значе- установлению такого равновесия. ниями (+1; –1) интервалов определяются условия ин- версионной переиндексации (соответственно,σ1↔σ2; Процедуры структурно-кинематического σ3↔σ2) главных осей тензора напряжений (при со- метода реконструкции тектонических хранении ориентации их, главных осей тензора Tσ = напряжений на рудоносных объектах const, в пространстве). Явлением инверсионной пере- индексации обусловлена унаследованная смена ме- Структурно-кинематическая реконструкция тек- ханизмов деформирования земной коры тонофизических параметров деформационного про- (…взброс↔сдвиг↔сброс…) на границах фаз цикла. цесса при изучении такого иерархически разномас- Непрерывно-прерывистое и эволюционное измене- штабного и генетически многофакторного природно- ние условий деформации сопровождается унаследо- го объекта, как земная кора с ее структурно-вещест- ванным развитием трещиноватости. Объективным венными неоднородностями и рудными телами, мог- свидетельством такого развития является геометри- ла быть осуществлена только с позиций системного ческое подобие («scaling» англ.) систем трещин в подхода к процедурам анализа. Поэтому выполнение горных породах и систем микротрещин в минералах процедур – и для анализа истории развития условий этих пород (Васильев, Мострюков, 2000). тектонического нагружения формирования деформа- Полный цикл деформации реализуется в шести- ционной структуры рудоносных объектов, и для вы- кратном изменении формы и ориентации эллипсои- деления синрудных фаз деформирования земной ко- дов напряжений (деформаций) как закономерном ры и на этой основе прогноза перспективных участ- следствии непрерывной пульсации соотношений ме- ков объекта – проводится с учетом нескольких важ- жду относительными величинами главных нормаль- ных для целей исследования особенностей объекта. 80
Рис. 2. Фазово-циклическая модель инверсионного изменения условий деформации геологической среды. Выделенные сектора показывают условия деформации, опти- мальные для рудогенеза
Первая особенность определяется необходимо- но ориентировок главных нормальных напряжений стью реконструкции тектонического поля напряже- (σ1рег, σ3рег) регионального ранга на примере взбросо- ний регионального и локального иерархических ран- вой фазы регионального деформирования рудоносно- гов, каждое из которых (целое и его часть) характе- го объекта. Часть поверхности верхней полусферы ризуется своими целостными (эмерджентными, Ни- сетки Вульфа, ограниченная следом поверхности колаев, 1992, с. 42) свойствами, информацию о кото- прямого шарового конуса с осью σ1рег регионального рых представляют, соответственно, данные об «ос- ранга, запретна для ориентации осей σ3лок локального новном и дополнительных (внутрислойных)» (Гзов- ранга, а часть поверхности полусферы, ограниченная ский, 1975, с. 122, 124) механизмах деформирования следом поверхности прямого конуса с осью σ3рег ре- исследуемого объекта. Необходимость «рангового» гионального ранга, запретна для ориентации осей подхода к процедурам реконструкции продиктована σ1лок локального ранга. Эти соотношения служат кри- и отмеченной выше обязательной дискретностью териями процедуры временнòй сепарации механиз- развития рудных сегрегаций на «внутренних» грани- мов деформирования локального ранга. Процедура цах неоднородностей структуры объекта. Понятно, сепарации, или определение относительного возрас- что выявление таких границ в структуре тектониче- та этих механизмов, проводится исходя из данных ского поля напряжений возможно только по инфор- соотношений – по указанному выше необходимому мации об изменчивости «внутрислойных» – локаль- соответствию ориентации осей главных нормальных ных механизмов (условий) деформирования объекта. напряжений (σ1, σ3) механизмов локального и регио- Вторая особенность определяется необходимо- нального рангов. стью получения информации о пространственной На фрагменте 2 рис. 3 приведены эллипсоиды де- дифференциации условий «уплотнения-разуплотне- формации, форма которых определяется тремя пре- ния» структуры объекта на синрудных фазах его де- дельными значениями коэффициента Лодэ-Надаи (με формирования. В основе этого требования лежит хо- = –1; με = 0; με = +1). Каждая из форм характеризует рошо известная из петрофизики горных пород и руд процесс деформации в условиях последовательно закономерная связь рудных сегрегаций с границами сменяющихся одноосного «разуплотнения», трехос- участков, характеризующихся различиями в эффек- ного нагружения и одноосного «уплотнения» блоков тивной пористости и, соответственно, проницаемо- земной коры (слева направо). сти горных пород для процессов флюидо-массопере- Первыми условиями определяется развитие про- носа (Звягинцев, 1978, с. 139). цесса деформации при максимальном разуплотнении Третья особенность обусловлена необходимо- земной коры в направлении ориентации оси ε1. Про- стью комплексных оценок условий проницаемости цесс деформации состоит в раскрытии трещин и раз- объекта в процессе его деформирования – по количе- рывов, ориентированных и по образующей прямого ственным характеристикам условий относительных шарового конуса под углом около 45°, и по нормали «сжатия-растяжения» (+1 ≥ μσ ≥ –1), «уплотнения-ра- относительно оси ε1. Условия трехосного нагруже- зуплотнения» (+1 ≥ με ≥ –1) и вертикального «уплот- ния определяют развитие процесса деформации – нения-разуплотнения» (±ΔΖ) деформационной струк- раскрытия трещин и разрывов, ориентированных в туры объекта на синрудной фазе его деформирования. соответствии с простиранием плоскостей действия Как показывает практика, только комплексные оценки максимальных касательных напряжений (τmправ, содержат достаточно надежную и объективную ин- τmлев). Условиями одноосного сжатия определяется формацию о дифференцированном в пространстве, развитие процесса деформации при максимальном дискретном проявлении рудолокализующих механиз- уплотнении земной коры в направлении ориентации мов деформационного процесса. Численно характери- оси ε3 и раскрытии кругового пояса трещин (разры- стики этих условий определяются следующими значе- вов), ориентированных по образующей прямого ша- ниями параметров взбросовых механизмов деформи- рового конуса с осью, ортогональной (под углом 90°) рования локального ранга: отрицательными значения- к ориентации главной оси ε3. ми коэффициентов Лодэ-Надаи: (–1 < μσлок < 0) «растя- Параметр деформационного процесса (+1 ≤ με ≥ –1) жения», (–0,3 < μεлок < 0) «разуплотнения» и положи- является одним из критериев оценки условий разви- тельными значениями величины (+0,14 ≥ ΔΖлок ≥ 0,12) тия обстановок эндогенного рудогенеза по условиям «вертикального» разуплотнения. общего разуплотнения среды (0 ≥ με ≥ –0,3) на взбро- На рис. 3 приведены иллюстрации, поясняющие совой фазе деформирования земной коры. При этом содержание процедур выделения потенциально-рудо- надо отметить, что «близкий» параметр (+1 > μσ > –1), носных объемов геологической среды по данным о несущий информацию об условиях «внешнего» текто- численных соотношениях между параметрами (σ1лок, нического нагружения блоков земной коры, является σ3лок, σ1рег, σ3рег, ±μεлок, ±μσлок, ±ΔΖлок) региональных и одним из трех критериев оценки развития обстановок локальных механизмов деформирования земной коры. рудогенеза в условиях растяжения (0 > μσ > –1). На фрагменте 1 рис. 3 графически показаны гра- Однако недостаточная репрезентативность его вы- ничные условия, контролирующие возможное рас- числений, обусловленная методическими причина- пределение ориентировок главных нормальных на- ми, делает информацию о параметре наименее на- пряжений (σ1лок, σ3лок) локального ранга относитель- дежной. 82
Рис. 3. Процедуры выделения рудолокализующих условий деформации геологической среды. Дополнительные пояснения в тексте
На фрагменте 3 рис. 3 показана схема, иллюстри- щений в 103 пунктах наблюдений относительно рав- рующая процесс изменения относительной величины номерно распределенных по исследованной площади вертикальной компоненты ±ΔZ деформации. Отме- (рис. 4) и «привязанных» к координатам топографи- тим, что объем эллипсоида равновелик объему еди- ческой сети. На камеральном этапе сначала были вы- ничной сферы. На верхней части схемы показан про- полнены реконструкции региональных механизмов, цесс увеличения (+ΔZ) этой компоненты, фиксируе- а затем локальных механизмов деформирования зем- мого по отрезку вертикальной оси между «образую- ной коры. Такой порядок выполнения реконструкций щими» сферы и эллипсоида, которая проходит через объясняется тем, что параметры механизмов регио- центры этих фигур. На нижней части схемы показан нального иерархического ранга являются необходи- процесс сокращения (–ΔZ) величины вертикальной мой информационной основой (критерием) для опре- компоненты деформации. Знак и величина измене- деления «возрастной» характеристики реконструи- ния вертикальной компоненты деформации являются руемых следующей процедурой механизмов локаль- объективными свидетельствами уплотнения (–ΔZ) ного ранга. Это определение, как уже отмечалось, или разуплотнения (+ΔZ) деформационной структу- проводится по соответствию кинематики локальных ры исследуемого объекта, происходящих в процессе и региональных механизмов деформации земной ко- фазового изменения структуры тектонического поля ры (см. фрагмент 1 на рис. 3). напряжений. Верхний рисунок иллюстрирует наибо- лее благоприятные для флюидо-массопереноса и ме- Р егио н а ль ны е м е ха н и з м ы пр о цес са тасоматоза условия вертикального разуплотнения де ф о р м и р о ва н и я з ем но й ко р ы геологической среды в процессе ее деформации (ΔZ Первая процедура состояла в реконструкции по- > 0). Нижний рисунок иллюстрирует, наоборот, не- следовательного ряда фазовых условий деформиро- благоприятные для флюидо-массопереноса и метасо- вания земной коры, определявших временнýю смену матоза условия уплотнения среды в процессе ее де- параметров деформационного процесса регионально- формации (ΔZ < 0). Численная характеристика (±ΔZ) го иерархического ранга. В результате реконструк- является вторым главнейшим критерием оценки ус- ции, проведенной по данным о выборочной совокуп- ловий, благоприятных для формирования обстановок ности из 1282 направлений тектонических смещений эндогенного рудогенеза. (τn = 1282), выявлено четыре инверсионно сменяю- Фрагмент 4 рис. 3 иллюстрирует условия деформи- щих друг друга фазовых механизма деформационно- рования структуры рудного объекта на примере двух го процесса (справа налево, от наиболее древней к фаз (А и В) развития деформационного процесса. Ус- наиболее молодой фазе, это фаза D → фаза C→ фаза ловия охарактеризованы соотношениями между чис- B → фаза A на рис. 4). ленными значениями параметров (±μεлок, ±ΔZлок, Из сравнительного анализа механизмов дефор- ±μσлок) на графиках, построенных в координатах (±με, мационного процесса следует, что деформирова- ±ΔZ), – отдельно для рудной (A) и безрудной (B) фаз ние земной коры в пределах объекта происходило деформационного процесса. В левом верхнем квад- в условиях сохранения диагональных (СЗ-ЮВ и ранте графика (А) черными кружками обведены ха- СВ-ЮЗ) ориентаций у двух «пологих» (субгори- рактеристики тектонофизических параметров дефор- зонтальных) главных осей тензора напряжений на мирования «продуктивных» участков рудоносных всех четырех последовательно сменяющих друг объектов. Эти характеристики содержат принципиаль- друга фазах деформационного процесса. При этом но новую количественную информацию о рудолокали- устанавливается явное соответствие СЗ-ЮВ ори- зующих условиях тектонического нагружения геоло- ентации одной из главных осей тензора простира- гических структур в развитии обстановок рудогенеза. нию линеаментов мегарегионального иерархиче- Сопоставление характеристик тектонофизических па- ского ранга (подчеркиваемому береговыми линия- раметров деформационного процесса из рудоносных ми крупнейших озер Карелии и Финляндии), а объектов разного возраста, сформированных в разных другой оси – СВ-ЮЗ простиранию береговых ли- геологических условиях и разных регионах, позволи- ний менее крупных Каменного и Заячьего озер в ло сделать вывод о феноменальной близости этих ус- районе месторождения (смотри рис. 4 и врезку на ловий, названных рудолокализующими и, вероятно, этом рис.). Соответствие ориентации главных осей имеющих более общее значение для эндогенных руд- тензора (Т σ) генеральным направлениям «унасле- ных полей (месторождений) (Васильев, Мострюков, дованного» простирания границ новейших разно- 2007). Использование рассмотренных характеристик в масштабных морфоструктур, вероятно, можно счи- качестве критериев для прогноза обстановок рудоге- тать предварительным свидетельством достовер- неза подтвердило эффективность их применения и в ности выполненных реконструкций. Результаты настоящем исследовании. дальнейшего сопоставления собственно «кинема- тики» реконструированных механизмов деформа- Результаты исследования ции с последовательностью тектонических собы- тий, установленных (Кулешевич, Фурман, 2009) в На полевом этапе работы собраны данные об ори- исследованном районе, только подтверждают этот ентации 1282 векторов сколовых тектонических сме- вывод. 84
Рис. 4. В левом верхнем углу – географическое положение района. В правом верхнем углу – схема геологического строения района по В. Фурману, 2001 г. Внизу – си- ноптические стереограммы осей и плоскостей действия главных напряжений на последовательных фазах деформации. Дополнительные пояснения в тексте
Смена условий деформирования земной коры Аз. пр. 54° СВ в т. н. 92), ориентированных по нор- происходила при последовательном изменении кине- мали к ориентации оси «сжатия» (σ3), явно подтвер- матического типа механизмов деформационного ждают достоверность выводов о кинематике условий процесса – от сдвигового типа (D) к взбросовому (С), деформации земной коры на фазе (В). далее еще к одному взбросовому (B) и опять сдвиго- На заключительной фазе (А) деформирования вому (A). земной коры кинематический тип механизма процес- Развитие древней фазы деформационного процес- са становится опять сдвиговым, но уже с другой (об- са (D), судя по коэффициенту Лодэ-Надаи (με = ратной по сравнению с фазой D) ориентацией глав- –0,16), происходило в условиях общего разуплотне- ных напряжений (σ1, σ3, рис. 4). По величине коэф- ния структуры объекта. Наиболее вероятным собы- фициента Лодэ-Надаи (με = +0,19) можно судить о тием было раскрытие границ структурно-веществен- резком ослаблении процесса уплотнения структуры ных неоднородностей объекта, ориентированных в объекта. Поэтому важным событием заключительной соответствии с СВ-ЮЗ простиранием субвертикаль- фазы (А) стало внедрение наиболее «молодых» в ной плоскости «растяжения» (σ1, рис. 4), другим со- структуре объекта даек лампроитов рифейского воз- бытием – тектонические движения по границам раста (Аз. пр. 340°, Кулешевич, Фурман, 2009), ори- структуры объекта, ориентированным в соответст- ентированных в соответствии с элементами залега- вии с С-СВ простиранием плоскости право-сдвиго- ния реконструированной плоскости действия относи- вых максимальных касательных напряжений (τmправ, тельного «растяжения» (σ1 – Аз. пад. 57° СВ∠78°, рис. 4). С такой кинематикой согласуется процесс рис. 4). Другим событием заключительного процесса интрузии и формирования контактов Центрального и деформации, отмеченным многими исследователями Факторного гранитоидных массивов с комплексами Таловейса, явилось проявление левосдвиговых тек- основных и ультраосновных горных пород. тонических движений в полном соответствии с Ведущим фактором процесса деформирования реконструированной плоскостью действия лево- земной коры на фазе (С) следует полагать особенно- сдвиговых максимальных касательных напряжений сти переходного режима изменения кинематики про- (рис. 4) – по границе Центрального и Факторного цесса. Признаками такого режима являются: 1) смена массивов гранитоидов. кинематического типа механизма деформации от Корреляция последовательных событий в текто- «сдвигового» на «взбросовый» и 2) смена процесса ническом развитии объекта с последовательной сме- разуплотнения структуры объекта процессом ее на- ной региональных механизмов (D→C→B→A) де- растающего, но далеко не максимального уплотне- формирования земной коры свидетельствует о весь- ния. Величина коэффициента Лодэ-Надаи здесь дос- ма вероятной достоверности выполненных реконст- тигает значения только με = +0,30. Смена сдвигового рукций. Данный вывод позволяет использовать ин- механизма деформации на взбросовый определяет формацию о региональных механизмах процесса для появление важнейшего условия для увеличения ам- сепарации локальных механизмов деформирования плитуды вертикальных тектонических движений земной коры по соответствию их кинематике фаз (D, (контролируемого субвертикальной ориентацией оси C, B, A), исходя из соотношений между ориентацией «растяжения» σ1; см. рис.4) и, соответственно, для осей главных нормальных напряжений (σ1лок, σ3лок и возникновения локальных зон повышенной прони- σ1рег, σ3рег) – как это было показано на фрагменте 1 цаемости в структуре объекта. Последнее, как из- рис. 3. вестно, является необходимым условием для разви- тия процессов и флюидо-массоперноса, и метасома- Ло ка ль н ые м е х а н из м ы пр о це с са тоза, и рудогенеза. С реализацией этого же условия, де ф о р м и р о ва н и я з ем но й ко р ы по-видимому, следует связывать возникновение Реконструкция механизмов деформирования зем- бòльшего числа обстановок рудогенеза на фазе (С), ной коры, определявших изменчивость локальных которую поэтому назовем Главной рудной фазой. условий ее тектонического нагружения или структу- Деформирование земной коры на фазе (В) про- ру тектонического поля напряжений локального ие- должает контролироваться взбросовым механизмом рархического ранга основана на (тех же, что и в ре- процесса и сопровождается возрастающим уплотне- конструкциях регионального ранга) данных о 1282 нием структуры объекта, характеризуемым величи- векторах тектонических смещений из 103 пунктов ной коэффициента Лодэ-Надаи (με = +0,49). Другой наблюдений. Напомним, что при реконструкции па- важнейшей особенностью процесса деформации яви- раметров тектонических напряжений регионального лась переориентация оси главного нормального на- ранга эта информация рассматривалась в форме еди- пряжения (σ3, «сжатия») с ЮЗ румбов на фазе (С) на ной «региональной» выборочной совокупности дан- ЮВ румбы на фазе (В) (рис. 4). По-видимому, и рост ных (∑Nt = 1282tn). Исходя из новой задачи, эта ин- уплотнения структуры, и новая ориентация оси формация (имеющая привязку к географическим ко- «сжатия» (σ3 – Аз. пр. 146°∠3°, рис. 4) вместе яви- ординатам) была преобразована методом «скользя- лись главной причиной появления складок в структу- щего окна» во множество «локальных» выборочных ре объекта. Наблюдения за простиранием крутопа- совокупностей, формирующее регулярную простран- дающих осевых поверхностей складок (например, ственную сеть исходных данных. Параметры такой 86
сети: шаг между узлами сети (L) = 20 м, радиус «ок- синего, зеленого и красного цветов, вторая – круго- на» (R) = 150 м, сумма векторов (Nt) ≈ 50tn не менее вой «синоптической» диаграммой (окраска осей и 50 шт. в пределах каждого узлового «окна» – были плоскостей та же, что и у параметров локального приняты, исходя из средней плотности распределе- ранга). Из сопоставления этих материалов следует, ния полевой информации на изучаемой площади. что наиболее устойчивыми азимутальными направ- Результаты реконструкций, выполненных для ка- лениями характеризуются оси сжатия (σ3лок) и проме- ждого из узлов регулярной сети, после процедуры их жуточного напряжения (σ2лок) локального ранга, в це- временнòй сепарации явились фактической основой лом близкие к ориентации соответствующих осей ре- для построения фазовых (D, C, B, A) планов измен- гионального ранга. Азимутальные же направления чивости локальных условий развития деформацион- осей растяжения (σ1лок) локального ранга менее ус- ного процесса в геологической структуре объекта. тойчивы и ориентируются по довольно разнообраз- Причем, естественно, изменчивости не только про- ным румбам. Эти факты свидетельствуют о том, что странственной, но и временнòй. Так, при сравнении развитие локальных механизмов деформации на суммарных количеств локальных механизмов дефор- Главной рудной фазе (C) происходило в условиях мации, реконструированных для разных фаз, уста- относительно небольшого регионального сжатия новлена определенная их дифференциация от фазы к (уплотнения структуры объекта), подчеркиваемого фазе (D→C→B→A). Условия деформации земной значением величины коэффициента Лодэ-Надаи коры на наиболее «древних» фазах D и C характери- (μεрег = +0,30) регионального ранга. Согласованность зуются равными количествами данных о 138 локаль- кинематики механизмов деформации локального и ных механизмах процесса. Условия деформации на регионального ранга является, по-видимому, сущест- «средней» (по возрасту) фазе B характеризуются венным аргументом в пользу достоверности не толь- данными только о 64 механизмах процесса, а на наи- ко выполненных реконструкций, но и вывода о ре- более «молодой» фазе A – данными уже о 160 меха- альном взаимодействии региональных и локальных низмах. Таким образом, количества механизмов из- условий тектонического нагружения объекта на син- меняются в ряду фаз: 138 шт.→138 шт.→64 рудной фазе его деформирования. шт.→160 шт. Сравнивая этот ряд с рядом значений На рис. 6 дана информация об изменчивости зна- коэффициента Лодэ-Надаи (με), реконструирован- чений коэффициента Лодэ-Надаи (–1 ≤ μεлок ≤ +1) в ных для фаз деформации регионального ранга структуре месторождения, основанная на данных о ре- (–0,16→+0,30→+0,49→+0,19), следует отметить об- конструкции локальных механизмов деформирования ратное соотношение экстремумов характеристик структуры месторождения на Главной рудной фазе этих рядов на фазе деформации (В). Оно заключается (С). Сразу отметим довольно большой диапазон (по- в том, что минимальному количеству локальных ме- ложительных и отрицательных) значений коэффици- ханизмов деформации (64 шт., фаза B) соответствует ента, свидетельствующий о разнообразии «синруд- максимальное значение коэффициента Лодэ-Надаи ных» условий локального уплотнения и разуплотне- (με = +0,49, фаза B). Вероятно, данное соотношение ния структуры месторождения. При этом следует об- следует связывать с закономерным сопряжением ус- ратить внимание на несомненные факты не только ловий максимального уплотнения и минимальной дискретно-зонального (направленного) развития, но и тектонической подвижности структуры объекта ис- их дифференциацию в структуре объекта. Так, если следований в результате действия взбросового меха- условиями локального разуплотнения структуры оп- низма деформации. И наоборот, максимальное коли- ределялось ее развитие преимущественно на юго-за- чество механизмов, характеризующих условия де- падном фланге объекта, то локальное уплотнение формации объекта на фазе A (160 шт.), скорее всего, структуры объекта характерно для его северо-восточ- следует связывать с резким разуплотнением структу- ного фланга. В целом зонам локального разуплотне- ры объекта (με = +0,19) и, соответственно, увеличе- ния и уплотнения присущи два генеральных направле- нием степени ее тектонической подвижности, а так- ния их простирания (СВ-ЮЗ) и (ССЗ-ЮЮВ), согла- же лучшей сохранностью векторов тектонических сующиеся с простиранием плоскостей действия, соот- смещений, возникших в фазу действия «самого мо- ветственно, промежуточного (σ2рег) напряжения и мак- лодого» сдвигового механизма деформации. симального сжатия (σ3рег) регионального иерархиче- Анализ пространственной изменчивости локаль- ского ранга. Факт пространственно-временного сопря- ных условий тектонического нагружения земной ко- жения кинематических особенностей процесса дефор- ры проведен на примере Главной рудной фазы (С), мации, фиксируемого по параметрам механизмов и где наиболее четко проявилось действие рудолокали- локального, и регионального ранга, по всей вероятно- зующих механизмов деформации. На рис. 5 для этой сти, свидетельствует о его объективной реальности. фазы дана информация об ориентации осей главных Поэтому данный факт можно считать достаточным нормальных напряжений локального ранга (σ1лок, аргументом для вывода о генетической связи процес- σ2лок, σ3лок) и ориентации осей-плоскостей действия са синрудного разуплотнения объекта с кинематикой главных напряжений (σ1, σ2, σ3, τmлев, τmправ) региональ- механизмов деформации локального и регионального ного ранга. Первая информация представлена тремя рангов на фазе C и для вывода о прогностическом зна- облаками («множествами») стрелок, соответственно, чении параметра (±μεлок). 87
Рис. 5. Ориентация главных нормальных напряжений локального ранга для Главной рудной фазы деформации (С). На синоптической стереограмме – оси и плоскости действия главных напряжений регионального ранга для Главной рудной фазы деформации (С). Дополнительные пояснения в тексте
Рис. 6. Распределение коэффициента Лодэ-Надаи (±με) на Главной рудной фазе деформации (С). Дополнительные пояснения в тексте
На рис. 7 дана информация об изменчивости вели- надежными (0,06 < ΔΖ < 0,09; –0,3 < με < 0). Из них чины условий «вертикального» уплотнения-разуплот- оранжевым цветом обозначены более надежные, а го- нения (±ΔΖлок) в деформационной структуре объекта лубым цветом – менее надежные. Такая дробная гра- на фазе C. Изменчивость этой величины, так же как и дация степени надежности выделения участков обу- величины коэффициента Лодэ-Надаи (±μεлок), характе- словлена, во-первых, практической необходимостью ризуется, во-первых, большим диапазоном ее отрица- определения критериев очередности выполнения по- тельных и положительных значений (соответственно, исково-оценочных работ, а, во-вторых, использовани- вертикального уплотнения и разуплотнения), а, во- ем в качестве таких критериев надежности расчета па- вторых, дискретно-зональной формой распределения раметров полей напряжений. в структуре объекта. Точно так же, как и по характе- Анализ простирания границ ареалов на левом фраг- ристикам коэффициента Лодэ-Надаи (+μεлок), северо- менте рис. 8 свидетельствует об их преимущественно восточный фланг объекта характеризуется преиму- широтной и меридиональной ориентации, из чего были щественно условиями вертикального уплотнения определены рекомендации к выбору диагональной (–ΔΖлок) его деформационной структуры. Однако про- (СВ-ЮЗ; СЗ-ЮВ) ориентации ствола наклонных поис- стирания границ между зонами вертикального уплот- ковых скважин. При выборе точек заложения кустов нения (–ΔΖлок) и разуплотнения (+ΔΖлок), в отличие от скважин было рекомендовано исходить в первую оче- зон с различными значениями величины коэффициен- редь из необходимости их расположения между участ- та Лодэ-Надаи (±μεлок), явно согласуются, во-первых, с ками оптимального сочетания рудолокализующих ус- простиранием (Аз. пр. 299° СЗ) плоскости растяжения ловий деформации земной коры, обозначенными крас- (σ1рег) регионального механизма деформации, а, во- ными кружками на левом фрагменте рис. 8. вторых, с изменчивостью простирания продуктивной На правом фрагменте рис. 8 приведена схема рас- зоны контакта толщ ультраосновного и гранитоидного пределения литохимических аномалий золота (в услов- составов. Выявленную картину зависимости величи- ных единицах), составленная для Центрального участка ны (±ΔΖлок) от кинематики условий тектонического месторождения Таловейс в ИМГРЭ, Машечкин, 2011 г. нагружения структуры объекта на Главной рудной фа- Из этой схемы четко следует, что распределение метал- зе (С), по-видимому, можно считать еще одним аргу- ла в структуре объекта определяется преимущественно ментом, свидетельствующим о достоверности выпол- ортогональным – меридиональным и широтным про- ненных реконструкций и прогностической информа- стиранием границ его аномалий. Этот вывод находится тивности рассматриваемого параметра (±ΔΖ). в полном согласии с нашим выводом, сделанным по ре- В ы де л е ни е пер сп е к ти в ны х уч а ст ко в зультатам тектонофизического анализа рудолокали- место рож де н и я Т а ло ве йс зующих условий деформационного процесса. Анализом численных характеристик кинематиче- Выводы ских параметров деформационного процесса установ- лено, что во время проявления Главной рудной фазы 1. Структурно-кинематический анализ условий (C) цикла деформации в структуре месторождения деформирования земной коры является эффектив- было сформировано семь участков с оптимальным со- ным методом создания информации для прогноза об- четанием рудолокализующих условий процесса, чис- становок пространственно-временного взаимодейст- ленные характеристики которых соответствуют выяв- вия процессов тектогенеза и рудогенеза. ленным для обстановок рудогенеза (Васильев, Мос- 2. Анализ тектонических деформаций земной ко- трюков, 2007). Во время проявления фаз деформации ры докембрийского возраста на основе модели фазо- (B) и (D) были сформированы: один перспективный во-циклического изменения условий деформацион- участок на фазе (В) и два участка на фазе (D). Резуль- ного процесса показал реальную возможность при- таты анализа рудолокализующих условий деформа- менения модели к анализу не только фанерозойских ции на фазах D, C и B представлены на левом фраг- этапов развития Земли. менте рис. 8. На этом фрагменте различной штрихов- 3. Результаты проведенного анализа явились еще кой показаны ареалы проявления условий разуплотне- одним подтверждением объективной реальности ру- ния деформационной структуры объекта, благоприят- долокализующих условий деформации земной коры. ные для развития процессов флюидо-массопереноса, 4. На основе анализа рудолокализующих условий метасоматоза и рудогенеза на соответствующих фазах деформации земной коры выделены перспективные (D, C, B). В пределах этих полей кружками красного и площади на месторождении Таловейс, достоверность синего цвета показаны участки с оптимальным соче- выделения которых согласуется с результатами ми- танием условий деформации для формирования обста- нералого-геохимического изучения вещественного новок рудогенеза (ΔΖ > 0,09; –0,3 < με < 0). При этом состава земной коры. красным цветом обозначены наиболее надежные уча- В заключение разрешите выразить глубокую бла- стки, а синим цветом – менее надежные. Кружками годарность нашему старшему товарищу кандидату оранжевого и голубого цвета показаны участки с соче- геол.-минер. наук Евгению Венедиктовичу Дарагану танием условий деформации, близким для выводов о за его инициативу и поддержку авторов в проведе- формировании обстановок рудогенеза, но еще менее нии данного исследования. 90
Рис. 7. Распределение тектонофизического параметра вертикального уплотнения-разуплотнения (±ΔΖ) геологической среды. Дополнительные пояснения в тексте
Рис. 8. На левом фрагменте схема выделения перспективных участков. На правом фрагменте схема литохимических аномалий золота по данным ИМГРЭ, 2011 г. Дополнительные пояснения в тексте
ЛИТЕРАТУРА Васильев Н. Ю., Корчемагин В. А., Костенко Н. П. и др. Гущенко О. И. Метод кинематического анализа струк- Этапы и стадии тектонического нагружения в эволюцион- тур разрушения при реконструкции полей тектонических ном развитии массива габбро-пироксенит-дунитовой фор- напряжений // Поля напряжений и деформаций в литосфе- мации (Ватыно-Вывенская сутура, Корякское нагорье) ре. М., 1979. С. 7–25. // Тектоника, геодинамика и процессы магматизма и мета- Гущенко О. И. Кинематический принцип относитель- морфизма: В 2-х т. М., 1999. Т. 1. С. 121–123 (Материалы ной геохронологии палеонапряжений (основной алгоритм XXXII Тектонич. совещ.; Т. 1). тектонического стресс-мониторинга) // Теоретические и Васильев Н. Ю., Мострюков А. О. Тектонофизическая региональные проблемы геодинамики. М., 1999. С. 108– реконструкция условий размещения благородных метал- 125 (Тр. ГИН РАН; Вып. 515). лов в дунитах расслоенного массива // М. В. Гзовский и Звягинцев Л. И. Деформации горных пород и эндоген- развитие тектонофизики. М., 2000. Т. 1. С. 281–295. ное рудообразование. М., 1978. 174 с. Васильев Н. Ю., Мострюков А. О. Закономерности раз- Кулешевич Л. В., Фурман В. Н. Золоторудное месторо- вития циклов деформации в процессах тектогенеза // Тек- ждение Таловейс в Костомукшской докембрийской зелено- тоника неогея: Общие и региональные аспекты: В 2-х т. каменной структуре (Карелия) // Геология рудных место- М., 2001. Т. 1. С. 90–93 (Материалы XXXIV Тектонич. со- рождений. 2009. Т. 51, № 1. С. 58–76. вещ.; Т. 1). Николаев П. Н. Методика тектоно-динамического ана- Васильев Н. Ю., Мострюков А. О. Мезо-кайнозойский лиза. М., 1992. 295 с. цикл деформации земной коры Байкальского рифта (по Петров В. А., Мострюков А. О., Васильев Н. Ю. Струк- тектонофизическим реконструкциям) // Строение и исто- тура современного поля напряжений мезо-кайнозойского рия развития платформ Евразии: Научные чтения памяти цикла деформации Байкальской рифтовой зоны // Геофизи- проф. М. В. Муратова (материалы совещания). М., 2002. ческие исследования. 2008. Т. 9, № 3. С. 39–61. С. 17–19. Рундквист Д. В. Глобальная металлогения // Смир- Васильев Н. Ю., Мострюков А. О. Особенности рудо- новский сборник-95 (Основные проблемы рудообразо- локализующих условий деформации геологической среды вания и металлогении). М., 1995. С. 92–123. в характеристиках тектонических полей напряжений Angelier J. Sur l`analyse de measures recueillies dans des // Фундаментальные проблемы геотектоники: В 2-х т. sites failles: l`utilite d`une confrontation entre les methods М., 2007. Т. 1. С. 126–130 (Материалы XL Тектонич. со- dynamiques et cinematiquues // C. R. Acad. Sci. Paris, 1975. вещ.; Т. 1). Vol. 281. P. 1805–1808. Васильев Н. Ю., Мострюков А. О. Критерии прогноза McKenzie Dan P. The relation between fault plane обстановок рудогенеза в структуре тектонического поля solutions for earthquakes and directions of the principial напряжений (на примере кварц-золоторудного проявле- stresses // Bull. Seismol. Soc. Amer. 1969. Vol. 59, N 2. ния докембрийского возраста) // Х Междунар. конф. «Но- P. 591–601. вые идеи в науках о Земле»: Докл. в 3-х т. М., 2011. Т. 1. Vasilyev N. J., Mostrjukov A. O., Sim L. A. Die Rolle der С. 185. tektonischen Spannungen bei der Entwicklung des Eldjurtinsk- Гзовский М. В. Основы тектонофизики. М., 1975. 536 с. Granitoidkörpers und seines Rahmens (Kaukasus). Гущенко О. И. Анализ ориентировок сколовых текто- Tektonophysikalische Rekonstruktion // Tectonic & Magma нических смещений и их тектонофизическая интерпрета- 2001, Special Issue in honour of Hans Cloos. Part II. ция при реконструкции палеонапряжений // Докл. АН Zeitschrift für geologische Wissenschaften, Berlin, 2002. СССР. 1973. Т. 210, № 2. С. 331–334. Vol. 30, H. 1/2. S. 131–144.
Вы также можете почитать