CC BY
164
GEODYNAMICS & TECTONOPHYSICS
PUBLISHED BY THE INSTITUTE OF THE EARTH’S CRUST SIBERIAN BRANCH OF RUSSIAN ACADEMY OF SCIENCES
2011 VOLUME 2 ISSUE 3 PAGES 219-243 DOI:10.5800/GT-2011-2-3-0044
ISSN 2078-502X
Stress-and-strain reconstruction for the eastern SEGMENT OF THE BALTIC SHIELD L. A. Sim1, D. V. Zhirov2, A. V. Marinin1
Schmidt Institute of Physics of the Earth RAS, 123995, Moscow, Bol’shaya Gruzinskaya str., 10-1, Russia Geological Institute of the Kola Science Centre RAS, 184209, Apatity, Fersman str., 14, Russia
Abstract: The publication presents stress determinations from geological and kinematical indicators of tectonic stress fields, varying in ranks, for the Kola Peninsula. The objective is to determine possible mechanisms of formation of recent structures in the eastern segment of the Baltic shield and to forecast seismogenic and technogenic hazard of fracturing.
The study is focused on the Kovdor and Khibin Paleozoic alkaline-ultrabasic blocks. Tectonic stresses are reconstructed by M.V. Gzovsky’s method [1954; 1975] based on identification of conjugated shear systems. Neotectonic stresses are studied by the kinematic method [Gushchenko, 1979] on the basis of measurements of tectonic displacment vectors from slicken-sides (Figure 2). Local stress data processed by the method for determination of general stress fields provide for reconstruction of main normal stresses which are arbitrarily considered as regional stresses [Sim, 1980; 2000]. This study uses the method of band-pattern distribution of fracturing in fault zones [Danilovich, 1961] which determines the main displacement line on the fault plane.
The study of the Zhelezny mining site (Kovdor block, Figures 3 and 4) revealed that elements of fractures of two different ages (centroclinal fractures of the prototectonic genesis and fractures of later tectonic activation) are spatially overlapping due to rock collapse and lacking stability of benches. Numerous inversions and changes of kinematics of relative displacements were reviewed. It was observed that the south-eastern wall of the quarry collapsed due to local extension (Fig. 6B and Photos 1 and 2), and a large fault, i.e. a prototectonic fracture, was dipping towards the quarry. Based on the analyses of local stresses at various points of the quarry (Table 1), two ‘regional’ stress fields can be revealed (Figures 7 and 8). The first paleostess field was associated with reverse faults of the WNW sub-horizontal axis of compression and the steeply dipping axis of extension. The second field was related to shear faults; its axis of compression maintained the WNW orientation, while the extension axis was reoriented, and the axis of extension attained sub-meridional position and a less steep dip. The second field is younger as suggested by traces of two displacements identified on several planes, the youngest of which are shears.
From the analysis of measurements taken at 273 planes with striations, it is evident that striations are distributed in a band-shape pattern. The band of fractures is coincident with the plane of the transition axis of the young ‘regional field (Figure 9); main maximums of density of the planes with striations are symmetrically scattered in a fan-like pattern from the axis of compression and extension of this field. Generally, the striations reflect traces of younger displacements, and their consistency with the axes of the young field supports our conclusions on age relationships between the two ‘regional’ fields. Four major stages of deformation of the Kovdor block under study are distinguished by analyses of the tectonic stresses (Figure
10).
Within the limits of the Khibin block, 14 local states of stresses are reconstructed for three mining sites (Table 2, Figure 11). At the Central mining site, re-indexation of local axis of compression and extension in the fault wings give evidence of the fault activity during the neotectonic stage (Figure 13). The ‘regional’ stress field of the Khibin block is associated with a reserve fault with low-angle NNW orientation of the axis of compression (Figure 14). The tectonophysical studies conducted at the mining sites of the Kola Peninsula give grounds to conclude that activity of faults, which positions are different at the mining sites, is variable, depending on orientations of the faults against the youngest ‘regional’ main stress. From sets of indicators, a relative age of the revealed ‘regional’ fields of stresses is accepted as neotectonic and recent.
For the segments of the Kovdor block under study, four phases of deformation are distinguished, including two early phases revealed by structural indicators, and two last phases revealed from orientations of the axes of the main stresses in the reconstructed ‘regional’ fields. The reverse field of stresses of Deformation Phase 3 (which is a more ancient ‘regional’ field according to reconstructed tectonic stresses) at the Kovdor block and the reverse-fault field at the mining sites of the Khibin block may reflect a phase of brittle deformation of the rocks after the blocks were exposed to the day surface. Since then the deformation mechanisms might have been determined by two factors which controlled sub-horizontal compression: residual gravity stresses in considerably eroded magmatic bodies as ‘recollections’ of being subject to constraint environment at depth [Rebetsky, 2008], and the impact of rifting in the Northern Atlantics. The fact that the neotectonic ‘regional’ stress field of the Kovdor block is fully similar to that of the Central Karelia (Figure 1) give grounds to conclude that the mechanism of deformation of the block under study might have been controlled by both factors. The Khibin block has a lopolithic shape which
С
Ф
gradually converts into a central-type conic structure with depth. It should thus be squeezed upward by the impact of horizontal compression of any genesis, as evidenced by the sub-vertical extension axis of the general field of the Khibin area and the recent topography as the highest mountains of the Kola Peninsula are located in the block under study.
Key words: Kola Peninsula, tectonophysics, fault tectonics.
Recommended by V.A. San’kov 29 August 2011.
Citation: Sim L.A., Zhirov D.V., Marinin A.V. Stress-and-strain reconstruction for the eastern segment of the Baltic shield // Geodynamics & Tectonophysics. 2011. V. 2. № 3. P. 219-243. doi:10.5800/GT-2011-2-3-0044.
Реконструкция напряженно-деформированного состояния восточной части Балтийского щита Л. А. Сим1, Д. В. Жиров2, А. В. Маринин1
гИнститут физики Земли им. О.Ю. Шмидта РАН, 123995, Москва, ул. Большая Грузинская, д. 10, стр. 1, Россия
2Геологический институт Кольского научного центра РАН, 184209, Апатиты, ул. Ферсмана, д. 14, Россия
Аннотация: Целью исследований являлась реконструкция по геолого-кинематическим индикаторам разноранговых полей тектонических напряжений на Кольском полуострове для определения возможный механизмов формирования новейших структур восточной части Балтийского щита, а также прогноза сейсмогенной и техногенной опасности разрышнык нарушений.
Основные исследования проведены в пределах Ковдорского и Хибинского палеозойских щелочно-ультра-основных массивов. Для реконструкции тектонических напряжений в ограниченном объеме использовался метод выделения сопряженнык сколовык систем трещин, разработанный М.В. Гзовским [1954, 1975]. Основным методом изучения неотектонических напряжений быи кинематический метод О.И. Гущенко [1979], фактическим материалом для которого служили замеры векторов тектонических перемещений на зеркалах скольжения (рис. 2). По данным о локальных стресс-состояниях методом нахождения общих полей напряжений восстанавливаются главные нормальные напряжения условно регионального уровня [Сим, 1980, 2000]. В работе использовался метод поясного распределения трещиноватости в зонах разломов, предложенный В.Н. Даниловичем [1961], который позволяет определять линию главного перемещения по плоскости разлома.
На руднике «Железный» (Ковдорский массив, рис. 3, 4) установлено, что ряд обрушений и нарушений устойчивости уступов вызывает пространственное наложение и совмещение разрышнык элементов двух различнык этапов образования (центриклинальнык трещин прототектонического генезиса и разрывных нарушений более поздней тектонической активизации), а также многочисленные инверсии и изменения характера кинематики относительным смещений. Так, обрушение юго-восточной стенки карьера обусловлено наличием обстановки локального растяжения (рис. 6), а также крупного разрывного нарушения - прототектонической трещины с падением в сторону карьера. В результате анализа локальнык стресс-состояний, определенных в разных точках карьера (табл. 1), мы предполагаем наличие двух «региональным» полей напряжений (рис. 7, 8). Первое действовавшее здесь поле палеонапряжений -взбросовое с ЗСЗ ориентацией субгоризонтальной оси сжатия и крутой осью растяжения. Второе поле сдвигового геологического типа характеризуется осью сжатия, сохранившей свою ЗСЗ ориентировку, а также переиндексацией осей растяжения и промежуточной, при которой ось растяжения стала субмеридиональной и пологой. Возраст 2-го поля напряжений принимается более молодым по следам двух перемещений на отдельный плоскостях, из которык более молодые - сдвиговые.
Анализ плотности всех измеренным 273 плоскостей с бороздами скольжения показал их отчетливое поясное распределение. Пояс трещин совпадает с плоскостью действия промежуточной оси молодого «регионального» поля (рис. 9), а главные максимумы плотностей плоскостей с бороздами скольжения расходятся от оси сжатия и растяжения этого поля симметрично. Борозды скольжения обычно отражают следы наиболее молодык перемещений, такая согласованность их с осями молодого поля напряжений подтверждает справедливость определения возрастнык взаимоотношений двух «региональных» полей. На основе изучения тектонических напряжений вы1делены1 четыре основные этапа деформирования изученной части Ковдорского массива (рис. 10).
В пределах Хибинского массива на трех рудниках восстановлено 14 локальнык стресс-состояний (табл. 2, рис.
11). На руднике «Центральный» переиндексация локальнык осей сжатия и растяжения в разнык крыльях разлома свидетельствует о его активности на неотектоническом этапе (рис. 13). «Региональное» поле напряжений Хибинского массива характеризуется взбросовым типом с пологой ССЗ ориентировкой оси сжатия (рис. 14).
Проведенные на рудниках Кольского п-ова тектонофизические исследования позволили сделать несколько вы-
водов. Степень активности разломов с разными элементами залегания на рудниках различна и зависит от их ориентации по отношению к новейшим «региональным» главным напряжениям. По комплексу признаков установлен относительный возраст выделенных «региональных» полей напряжений, который принят за неотектонический и современный.
На исследованной части Ковдорского массива выделены: два ранних этапа деформирования - по структурным признакам, два последних - по ориентации осей главных напряжений восстановленных «региональных» полей. Взбросовое поле напряжений 3-го этапа деформирования (более древнее «региональное» по реконструкции тектонических напряжений) на Ковдорском массиве и взбросовое поле на рудниках Хибинского массива могут отражать, предположительно, этап хрупкого деформирования массива пород непосредственно после выведения массивов на дневную поверхность. С этого времени механизм деформирования мог определяться двумя факторами, обусловливающими субгоризонтальное сжатие: остаточными гравитационными напряжениями, сохранившимися в значительно эродированных магматических телах как «память» о пребывании в стесненных условиях на глубине [Rebetsky, 2008], и влиянием рифтогенеза в Северной Атлантике. Полная аналогия неотектонического «регионального» поля напряжений Ковдора с таковым в Центральной Карелии (рис. 1) позволяет считать, что механизм деформирования указанного массива могут определять оба источника. Хибинский массив, имеющий лополитообразную форму, переходящую ниже в коническую структуру центрального типа, под воздействием горизонтального сжатия любого генезиса должен выдавливаться вверх, что фиксируется субвертикальной осью растяжения общего поля Хибин и подтверждается современным рельефом - наиболее высокие вершины Кольского п-ова находятся именно на этом массиве.
Ключевые слова: Кольский полуостров, тектонофизика, разломная тектоника.
1. Введение
Изучение напряженного состояния в массивах горных пород по геологическим индикаторам имеет свои специфические особенности. Они связаны с неоднородной информативностью геологических данных, используемых для реконструкции тектонических напряжений, с определением их возраста, с проблемами ранжирования и установления взаимосвязи полей напряжений и разломов разных рангов, с сопоставимостью полученных результатов с данными инструментальных наблюдений и т.д. В то же время тектонофизические исследования указанного типа имеют важное теоретическое и практическое значение. В предлагаемой работе обсуждаются некоторые аспекты тектонофизических исследований, проведенных по геологическим данным о напряжениях. Работы1 были выполнены в 2009 и 2010 гг. на Ковдорском и Хибинском щелочно-ультрабазитовых массивах и в отдельных точках Мурманского блока Кольского п-ова. Они возобновляют и продолжают тектонофизические исследования, выполнявшиеся на Балтийском щите в 70-80 годах прошлого века.
Целью исследований являлась реконструкция тектонических напряжений на рудных объектах Кольского п-ова для прогноза степени опасности разрывных нарушений и исследование возможных механизмов формирования новейших структур восточной части Балтийского щита.
2. Краткая история геологического развития
РАЙОНА РАБОТ
Район исследований расположен в пределах северовосточной части Балтийского (Фенноскандинавского)
докембрийского щита. Главные структурно-вещественные комплексы и блоковая (доменная - террейно-вая) структура Кольского региона сформировались преимущественно в период неоархея - раннего протерозоя. Первые «микроконтиненты» - сиалические ядра
- начали формироваться, вероятно, ранее 3 млрд лет. Рост коры и ее эволюция осуществлялись за счет аккреции и орогении. Все составные части Кольского кратона имеют отличительные индивидуальные черты архейской истории развития, что свидетельствует о значительной пространственной разобщенности этих блоков в рассматриваемый период. В это же время активно развивались вулканогенно-осадочные пояса, глубокометаморфизованные зеленокаменные и пара-гнейсовые комплексы которых сегодня представлены в различной степени фрагментированными и микшированными реликтами архейских зеленокаменных поясов. Породы нестратифицированного основания, как первично осадочного, так и первично плутонического генезиса, неоднократно испытали в неоархее (2.8-2.9 млрд лет) метаморфизм высоких ступеней с различными по характеру и интенсивности деформациями и дислокациями [Строение..., 2006]. В масштабах всего Кольского кратона в архее преобладали тектонические движения блоков с преимущественной горизонтальной составляющей и в меньшей мере - с вертикальной.
В период последующей коллизии (2.7-2.5 млрд. лет назад) позднеархейские комплексы были дислоцированы, метаморфизованы в условиях от амфиболитовой до гранулитовой фации, а регион был преобразован в архейский интраконтинентальный Кольско-Беломор-ский коллизион [Геология., 2002]. Тектоника архея рассматривается как тектоника малых и тонких плит, подстилаемых мигматитовым слоем средней коры и анортозит-эндербит-гранулитовой нижней корой. По
Ф.П. Митрофанову, масштабным горизонтальным и субдукционным перемещениям на этой стадии препятствовали гетерофазные условия их взаимодействия [Митрофанов и др., 1997].
В раннепротерозойский период (2.5-1.7 млрд лет назад) протоплатформа имела значительные размеры и явилась «ловушкой» для поднимающегося мантийного плюма, над которым образовалось сводовое поднятие коры с возникновением интракратонных рифтовых зон (Печенга, Имандра-Варзуга и др.) и интенсивным проявлением базит-гипербазитового магматизма. Тектонические процессы сумийского рифтинга привели к образованию авлакогенов, а ятулийского и свекофен-ского - палеорифтов. На поздних этапах преобладали значительные тангенциальные движения, с которыми связано формирование многочисленных складчатых и надвиговых структур. Оно происходило зонально с различным распределением интрузий базитов и гипер-базитов в разных зонах [Геология..., 2002].
В позднем протерозое (1.65-0.65 млрд лет назад) в режиме платформы на юго-востоке и востоке региона сформировались континентальные отложения, а на севере - образования шельфа и континентального склона. С этого времени в регионе преобладали блоковые радиальные тектонические движения. При общей тенденции к воздыманию происходила неравномерная эрозия и пенепленизация региона [Геология., 2002].
Палеозойский этап характеризуется большим разнообразием продуктов магматизма. Складчато-над-виговые процессы на северо-западе щита вызвали тектоническую и эндогенную активизацию в консолидированной коре региона. В это время сформировались новые и обновились старые разломы с неоднородными движениями по ним, а также был проявлен активный щелочно-ультраосновной магматизм: произошло внедрение Ковдорского массива (378.5-372.2 млн лет), Хибинского и Ловозерского плутонов (371-362 млн лет) и других массивов. Наиболее древними палеозойскими магматическими образованиями Кольского п-ова являются кимберлиты и кимберлитоподобные породы трубок взрыва Терского побережья (465-376 млн лет).
В палеозое сформировались вулканогенные, вулканогенно-осадочные и осадочные породы, сохранившиеся в виде останцов девонских пород на Хибинском и Ловозерском массивах и каменноугольных - в Кон-тозерской кальдере.
К мезозойским образованиям относится дайка не-фелинитов Турьего мыса - 175-150 млн лет [Рухлов и др., 1999; Геология., 2002].
Спецификой неотектонического этапа Балтийского щита и Кольского региона как его составной части является практическое отсутствие новейших отложений и совмещение тектонических движений с гляциоизо-статическими. Начало новейшего этапа относится к концу олигоцена, когда море ушло с северо-запада Русской плиты [Бабак, Николаев, 1983; Карта новей-
шей тектоники.1985; Са^Бку е1 а1., 1999]. Амплитуды новейших движений на плите рассчитывались по деформациям исходной поверхности - кровли нижнего олигоцена (рюпельские слои, харьковская свита), а в пределах щита - по деформациям поверхностей выравнивания.
В связи с отсутствием палеоген-неогеновых отложений на щите кайнозойская тектоника региона лучше изучена на позднеплейстоценовом-голоценовом этапе. Некоторая определенность в вопросе о начале неотектонического этапа была получена при оценке возможной величины эрозионного сноса с Фенноскандии в мезозое - кайнозое. Для такой оценки по литологопалеогеографическим картам мира [Ронов и др., 1989] с определенными допущениями были проанализированы мощности осадков этого времени в бассейнах Баренцева, Норвежского и Северного морей, в осред-ненном виде соответствующие величине эрозионного сноса. Мощности осадков позволили оценить скорость воздымания щита: в эоцене она равна 0.022 мм/год, в олигоцене - 0.053 мм/год, в миоцене - 0.009 мм/год и в плиоцене - 0.146 мм/год [Сим, 2010]. Согласно полученным данным подтверждено общепринятое начало неотектонического этапа в олигоцене, при этом в плиоцене отмечается ярко выраженный импульс активности внутри новейшего этапа, даже если в расчете скоростей подъема щита не учтены некоторые специфические особенности осадконакопления этого времени.
В позднем плейстоцене - голоцене для разных участков Кольского региона характерны региональные и локальные движения с неоднородной вертикальной амплитудой [Кошечкин, 1979].
На последних стадиях развития региона значительная роль принадлежит оледенениям. Амплитуда опускания из-за ледникового покрова на Кольском п-ове оценивается в первые десятки метров. Суммарная амплитуда перемещений за олигоцен - плиоцен достигала 1000 м, а в голоцене - 150 м и более. Постледнико-вое расчленение рельефа обусловлено «клавишным» характером дифференцированных вертикальных блоковых движений как докембрийских, так и палеозойских структур и массивов. Процессы поднятия продолжаются и в настоящее время. Линейные зоны -границы крупных структур региона - характеризуются повышенной сейсмичностью [Геология., 2002].
К расположенной в восточной части Балтийского щита флексуре Полканова (имеющей северо-восточное простирание) ортогонально расположены шовные зоны северо-западного простирания, к которым относятся и Ладожская, и Кандалакшская горсто-грабеновые системы.
Шовные зоны сопровождаются ярко выраженными осевыми линиями зон сосдвиговых движений второго порядка, указывающих на левосдвиговые раннелопий-ские перемещения.
В новейший этап по большинству из этих шовных
Рис. 1. Схема изученности тектонических напряжений восточной части Балтийского щита по геологическим индикаторам (по Л.А. Сим [2000], с изменениями). 1-3 - стереограммы неотектониче-ских напряжений, реконструированных кинематическим методом О.И. Гущенко [1973, 1979]: 1 - I, 2 - II, 3 - III ранга; показаны плоскости действия максимальных касательных напряжений и оси: красный цвет - максимальные, синий - минимальные, черный -промежуточные главные нормальные напряжения. Стереограммы темно-желтого цвета - реконструкции 2009-2010 гг. (сетка Вульфа, верхняя полусфера); 4 - сдвиги; 5-7 - разломы: 5 - I, 6 - II, 7 -III ранга; 8-9 - оси максимальных сжимающих напряжений в горизонтальной плоскости, восстановленные структурно-геоморфологическим методом [Сим, 1991]: 8 - I, 9 - II ранга; 10-12 - геологические структуры: 10 - Фенноскандинавский (Балтийский) щит, 11 - Русская плита, 12 - байкалиды Тимана.
Fig. 1. Schematic map of the Eastern segment of the Baltic shield covered by studies of geological indicators (modified from [Сим, 2000]). 1-3 - stereograms of neotectonic stresses determined by O.I. Gushchenko’s method [Гущенко, 1973, 1979]: 1 - Rank I, 2 -Rank II, 3 - Rank III; planes of maximum tangential stresses and axes: red - maximum; blue - minimum; black - transitional main normal stresses. Dark-yellow stereograms show stresses reconstructed for 2009-2010 (The Wulff net, upper hemisphere); 4 - strike-slip faults; 5-7 - faults: 5 - Rank I, 6 - Rank II, 7 - Rank III; 8-9 - axes of maximum compression stresses in the horizontal plane that are restored by the structural-geomorphological method [Сим, 1991]: 8 - Rank I, 9 -Rank II; 10-12 - geological structures: 10 - Fennoscandian (Baltic) Shield, 11 - Russian plate, 12 - baikalids of Timan.
зон, дешифрирующихся как крупные разломы, с помощью тектонофизических исследований также восстановлены левовосдвиговые компоненты перемещений, которые либо «реанимировали» направление ло-пийских смещений, либо развиваются унаследованно [Сим и др., 2009].
3. Изученность тектонических напряжений восточной части Балтийского щита по
ГЕОЛОГИЧЕСКИМ ИНДИКАТОРАМ
С конца 70-х годов прошлого века преимущественно коллективами специалистов геологического факультета МГУ [Бабак и др., 1979, 1981; Кольская сверхглубокая, 1984; Сим и др., 1987, 1990; Сим, 1991, 1996; Николаев, 1992] в пределах восточной части Балтийского щита были инициированы по существу первые исследования тектонических напряжений по геологическим индикаторам. Наиболее детально охарактеризовано напряженное состояние Центральной Карелии в связи с изысканиями, проводившимися в 1987-1988 гг. при выборе площадки под строительство Карельской АЭС. Кинематическим методом, разработанным О.И. Гущенко [1973, 1979] по бороздам скольжения, измеренным в древних архей-протерозойских отложениях, была установлена значительная изменчивость локальных стресс-состояний на площади 100 х 100 км.
По осям напряжений локального уровня было восстановлено общее поле напряжений сдвигового типа с запад-северо-западной ориентировкой оси сжатия и меридиональной - растяжения. Острая проблема возраста восстановленных тектонических напряжений была решена путем комплексного геохимического профилирования на детальных участках, сложенных голоценовыми ледниковыми отложениями. По данным об ориентации региональных главных напряжений Центральной Карелии для опорных детальных участков, выделенных по дешифрированию, были спрогнозированы кинематические типы разломов. Геохимические исследования показали аномально высокие содержания радона и углекислого газа вблизи широтного разлома, нормального к оси растяжения, а также в секторах локального растяжения, образованного при пересечении спрогнозированных на основании восстановленных в древних породах тектонических напряжений разнонаправленных сдвигов. Таким образом, тектонические напряжения, реконструированные по бороздам скольжения, характеризуют и современное напряженное состояние, соответственно, возраст этих напряжений неотектонический и современный. Кроме результатов геохимических исследований неотектонический возраст восстановленных кинематических напряжений подтверждается новейшей структурой Кандалакшского грабена. По геологическим данным он образовался в конце мезозоя - кайнозое; ось растяже-
ния II ранга, ориентированная нормально к оси грабена, соответствует механизму его формирования [Сим и др., 1990; Сим, 1996, 2000].
В конце ХХ в. проведено обобщение практически всех геологических данных по новейшему напряженному состоянию восточной части Балтийского щита (рис. 1). Основной вывод, последовавший из этого обобщения, заключался в том, что региональные тектонические напряжения указанного региона характеризуются сдвиговым типом с запад-северо-западной ориентацией оси сжатия и субмеридиональным растяжением, что хорошо коррелируется с данными инструментальных измерений современных напряжений и с данными по механизмам очагов землетрясений вдоль северо-западного побережья Фенноскандии [^Ьаек еЬ а!., 1989].
4. Методика работ
Основные принципы анализа полей напряжений, действующих в определенных объемах земной коры, были сформулированы М.В. Гзовским, впервые начавшим систематическое изучение современных и древних тектонических напряжений. Им предложен метод выделения сопряженных пар трещин [Гзов-ский, 1954, 1975], использовавшийся нами в единичных случаях для решения частных задач. Метод основан на положениях теории прочности Кулона, согласно которой разрушение материала и образование разрыва связываются с плоскостями действия максимальных касательных напряжений. Положение сколовых трещин отличается от указанных плоскостей на угол скалывания а. Таким образом, сколовые трещины расположены от оси сжатия а3 под углом 45°-а (т.е. биссектриса острого угла между сопряженной парой сколовых трещин есть ось сжатия а3, а тупого угла, соответственно, ось растяжения а1). Промежуточная ось а2 совпадает с линией пересечения сопряженных сколов. Проверкой справедливости восстановленных по сопряженным сколам ориентировок осей главных нормальных напряжений служат разнонаправленные векторы перемещений на сопряженных трещинах (разрывах).
Наиболее широко в работах 2009-2010 гг. использовался кинематический метод О.И. Гущенко. В этом методе реконструкции направлений действия главных напряжений и типа напряженного состояния используется взаимосвязь всей совокупности касательных напряжений, действующих на различно ориентированных площадках в однородно-напряженных участках деформируемой среды, с направлениями действия главных нормальных напряжений [Гущенко, 1973, 1979; Гущенко, Сим, 1974]. Обоснование метода опирается на теорию дислокационного скольжения Бат-дорфа-Будянского и содержит три основные предпосылки: упруго-пластическое деформирование в задан-
ном однородно-осном поле напряжений в геологической среде происходит как за счет смещений по вновь образующимся поверхностям трещин (разрывов), так и за счет перемещений по уже существующим ослабленным сечениям любого возраста и генезиса; каждое смещение вызывает возмущение поля напряжений только более мелкого масштабного уровня, а исходное поле напряжений остается неизменным; направление смещения вдоль всей поверхности сместителя совпадает с направлением касательного напряжения на этой поверхности и отвечает заданному полю напряжений. Важным условием применимости кинематического метода является положение о том, что «борозды скольжения являются следами последних... наиболее молодых сдвиговых тектонических подвижек, а следы1 более ранних смещений или стираются, или имеют худшую сохранность» [Гущенко, 1979, с. 9], что согласуется с данными, приводимыми почти во всех учебниках по структурной геологии.
Для характеристики взаимосвязи всей совокупности касательных напряжений на различно ориентированных площадках в однородно-напряженных участках деформируемой среды с направлениями действия главных нормальных напряжений были рассчитаны и построены стереограммы-номограммы распределения ориентировок касательных составляющих напряжений для различных типов тензоров напряжений (для коэффициента Лоде-Надаи ца=+1 и ца=0) [Гущенко, 1979]. Согласно этим номограммам, векторы перемещений расходятся от оси сжатия и сходятся к оси растяжения, границей разнонаправленных векторов служат плоскости действия осей а1 и а3; при одноосном напряженном состоянии векторы совпадают с дугами больших кругов, расходящихся от оси сжатия а3 (ца =1) или оси растяжения а1 (ца =-1), а при трехосном напряженном состоянии они должны находиться в створе острого угла, образованного этими дугами (рис. 2).
Метод нахождения общего поля напряжений по данным о локальных стресс-состояниях и выделение рангов полей напряжений. Данные математического моделирования характера локальных полей напряжений в окрестностях разрыва при смещении по нему [Осоки-на, 1987, 1997] позволили обосновать принцип выделения общего поля напряжений, вызвавшего перемещение по указанному разрыву. Он сводится к следующему: если все оси а1 или а3 локальных напряжений являются одноранговыми и вызваны перемещением по разрыву более высокого (крупного) ранга, то они должны описываться конусами сжатия и растяжения с углом при вершине 90°. При этом в конусе сжатия не должно быть осей растяжения и, наоборот, в конусе растяжения не должно быть осей сжатия. Оси конусов взаимно перпендикулярны и являются осями сжатия и растяжения общего поля напряжений, а точки их касания - полюсами площадок действия максимальных касательных напряжений. Взаимно перпендикулярные локальные оси сжатия и растяжения в разных крыльях
Рис. 2. Стереограммы локальных стресс-состояний Приполярного Урала с различными значениями коэффициента Лоде-Надаи. Сетка Вульфа, верхняя полусфера: а, е - трехосное напряженное состояние; б-д - одноосное сжатие; ж, з - одноосное растяжение. 1-3 - оси главных нормальных напряжений и плоскости их действия: 1 - минимальных (aj, 2 - промежуточных (о2) и 3 - максимальных сжимающих напряжений (a3); 4-5 - полюса плоскостей действия касательных напряжений: 4 - максимальных, 5 - прочих; 6 - границы конусов сжатия и растяжения; 7 - векторы перемещений на зеркалах скольжения: м, н - согласующиеся с найденными осями главных нормальных напряжений: м - с установленным, н - с неопределенным знаком смещения; о - противоречащие найденным осям напряжений; 8 - дуги больших кругов, расходящиеся от осей сжатия и растяжения.
Fig. 2. Stereograms of local states of stresses of the Polar Urals with various values of the Lode-Nadai factor. The Wulff net, upper hemisphere: а, е
- three-axes state of stresses; б-д - single-axis compression; ж, з - single-axis extension. 1-3 - axes of main normal stresses and planes: 1 - minimum (a^, 2 - transitional (a2) and 3 - maximum compression stresses (a3); 4-5 - poles of tangential stress planes: 4 - maximum, 5 - other; 6 -boundaries of compression and extension cones; 7 - displacement vectors on slickensides: м, н - consistent with the determined axes of main normal stresses: м - with established displacement sign, н - with non-determined displacement sign; о - contradicting to the determined stress axes; 8 - arcs of large circles diverging from axes of compression and extension.
разрывов на их концах попадают в область между конусами сжатия и растяжения [Сим, 1980, 2000].
Определение общего поля напряжений по данным о локальных стресс-состояниях по технологии совпадает с нахождением региональных тектонических напряжений по сейсмологическим данным, предложенным в кинематическом методе [Гущенко, 1979].
Метод поясов трещин В.Н. Даниловича. Метод, предложенный В.Н. Даниловичем [1961], основан на анализе трещиноватости вблизи крупных разломов. Большое число натурных данных показало, что вблизи разломов сколового типа более мелкие трещины образуют пояс трещиноватости. Кроме этого еще один локальный максимум плотности трещиноватости часто наблюдается вблизи полюса пояса трещиноватости.
Этот метод, прежде всего, был ориентирован на определение направления движения по крупным разломам с известным залеганием плоскости сместителя. В дальнейшем С.И. Шерман [1966] развил метод В.Н. Даниловича для определения элементов залегания плоскости разлома. Сопоставление экспериментальных и натурных данных позволило использовать этот метод для реконструкции главных осей тензора напряжений. Было установлено, что ось пояса совпадает с осью промежуточного главного напряжения, а анализ ориентации трещин отрыва и скола позволяет определить ориентацию двух других главных осей.
5. Результаты реконструкции тектонических
НАПРЯЖЕНИЙ В ВОСТОЧНОЙ ЧАСТИ БАЛТИЙСКОГО
ЩИТА
В 2009 и 2010 годах на Ковдорском и Хибинском щелочно-ультрабазитовых массивах Кольского полуострова и в отдельных точках севернее них были возобновлены тектонофизические исследования, проводившиеся на Балтийском щите в 70-80 годах прошлого века. Основные исследования проводились на карьерах рудников «Железный» (Ковдорский массив), «Центральный», «Олений Ручей» и «Коашва» (Хибинский массив)
Ковдорский массив. Ковдорский массив прорывает биотитовые гнейсы и гранитогнейсы кольско-бело-морского комплекса архея, которые образуют крупную моноклинальную структуру с общим падением на северо-восток. С вмещающими породами массив имеет эруптивные контакты. Структура массива представлена сложным сочетанием цилиндрических, конических и неполнокольцевых форм, обусловленных последовательным внедрением различных по составу интрузий. Массив отличает очень разнообразный состав пород, отраженный на геологической карте (рис. 3). Помимо разного состава магматических пород значительная неоднородность минералогического состава пород обусловлена неоднократными метасоматическими процессами и длительным периодом выветривания,
которому подвергся массив на различных этапах своего развития. По данным КЬ-Бг-метода возраст оливи-нитов и пироксенитов массива (центр и северо-запад массива) колеблется в интервале 373-376 млн лет №е№еу, Ве!!, 1995]. Близко к нему находится определение возраста бадделеитов из фоскоритов и карбона-титов, полученный И-КЬ методом [Геология рудных районов... , 2002].
Тектонофизические исследования были сосредоточены в пределах карьера рудника «Железный» (рис. 4).
На верхних и средних горизонтах карьера были произведены массовые замеры трещиноватости (около 2000 записей в карьере и более 1300 - по ориентированному керну), а в 23 точках наблюдений в разных бортах были измерены векторы тектонических перемещений. Кроме того, в наших исследованиях была использована база выполненных ранее замеров ФГУП ВИОГЕМ, насчитывающая более 17000 данных, относительно равномерно покрывающих карьерное поле.
Классификация трещиноватости Ковдорского бадделеит-апатит-магнетитового месторождения
(КБАММ) проведена по трем основным признакам: парагенезизу, структурно-формационной иерархии, а также характеру и стадийности кинематики. Первый признак позволил разделить всю трещиноватость на две большие группы: прототектонический парагенезис, включающий трещиноватость этапов «собственно магматического (позднемагматического) - контракци-онного генезиса», и парагенезис наложенной тектоники посткарбонатитовых этапов, включающий трещины выветривания и разгрузки приповерхностного слоя. Более дробное рассчленение этой структуры с выделением в группах различных подтипов и стадий на данном этапе исследований нецелесообразно. По отношению к месторождению элементы трещинной тектоники классифицируются на «внутриформационные» (т.е. развитые только в пределах месторождения) и «трансформационные» соответственно. По характеру развития выделяются монокинематические - одноактного образования и мультикинематические, испытавшие минимум два этапа деформаций. Внутриформационная парагенетическая ассоциация в масштабах всего КБАММ имеет представительность 90-95 % (авторская оценка) и включает следующие основные системы пород: радиальную с углами падения 65-90°, кольцевую субвертикальную (тангенциальную) с углами падения 70-90° и две диагонально-конические: цен-триклинальную, падающую радиально к центру массива под углом 25-55°, и периклинальную, падающую радиально от центра массива под углом 0-45°. Система субгоризонтальных трещин проявлена незначительно (преимущественно в приповерхностной части).
Формирование и развитие трещиноватости месторождения разнообразно по генезису, иерархической соподчиненности парагенетических ассоциаций и отдельных элементов трещинной тектоники и происходило полистадийно. Следствием этого является пре-
Рис. 3. Геологическая карта Ковдорского массива (составлена Б.В. Афанасьевыш и др. [Мультимедийный справочник2001] с добавлениями). 1 - гнейсы, амфиболиты и мигматиты архейского беломорского комплекса; 2 - фениты; 3: а - метаоливиниты, б - ослюденелые метаоливиниты; 4: а - оливиниты, б - рудные оливиниты, в - ослюденелые рудные пироксениты; 5: а - пироксениты и рудные пироксени-ты, б - пегматоидные рудные пироксениты, в - ослюденелые рудные пироксениты, г -нефелинизированные рудные пироксениты, д -нефелин-пироксеновые породы по гипербазитам; 6: а - мельтейгиты, б -ийолиты» 7: а - турьяиты1, б - ослюденелые турьяиты; 8: а - пи-роксен-мелилитовые породы, б - скарнированные пироксен-мелилитовые породы; 9: а - монтичеллитовые породы, б - скарнированные монтичеллитовые породы, в - гранат-амфибол-диопсидовые скарноподобные породы с монтичеллитом и везувианом; 10: а - амфибол-флогопит-форстеритовые породы, б - флогопитовые руды; 11: а - апатит-флогопит-форстеритовые породы (маложелезистые апатитовые руды), б - бадделеит-апатит-магнетитовые руды; 12 - нефелиновые сиениты; 13 - кальцит-апатит-флогопитовые породы; 14 -эгирин-биотит-карбонатитовые карбонатиты; 15 - форстерит-кальцитовые карбонатиты; 16: апатит-штаффелитовые руды; 17 - четвертичные отложения. Цифры в кружках - месторождения: 1 - апатит-магнетитовое, 2 - вермикулит-флогопитовое, 3 - апатит-штаффелитовое, 4 - руд-нык оливинитов.
Fig. 3. Geological map of the Kovdor block (modified by B.V. Afanasjev from [Multi-Media Guidebook 2001]). 1 - gneisses, amphiboles and migmatites of the Archean Belomorsky series; 2 - fenites; 3: a - metaolivinites, 6 - mica metaolivinites; 4: a - olivinites, 6 - ore olivinites, b - mica ore pyroxenites, r - nephelinic ore pyroxenites, g - nepheline-pyroxene rocks as per hyperbasites; 5: a - pyroxenites and ore pyroxenites, 6 - pegma-toid ore pyroxenites, b - mica ore pyroxenites, r - nephelinized ore pyroxenites, g - nepheline-pyroxene rocks by hyperbasites; 6: a - melteigites, 6
- iyolites; 7: a - hydrohematite (turiyaites), 6 - mica hydrohematites; 8: a - pyroxene-melilitics rocks, 6 - scarn pyroxene-melilitics rocks; 9: a -montichellite rocks, 6 - scarn montichellite rocks, b - garnet-amphibole-diopside scarn-like rocks with montichellite and vesuvian; 10: a - amphi-bole-phlogopite-fosterite rocks, 6 - phlogopite ores; 11: a - apatite-phlogopite-forsterite rocks (low-ferric apatite ores), 6 - baddeleyite-apatite-magnetite ores; 12 - nephelinic syenite; 13 - calcium-apatite-phlogopite rocks; 14 - aegirine-biotite-carbonatite carbonatites; 15 - fosterite-calcite carbonatites; 16 - apatite-shtaffelite ores; 17 - Quaternary sediments. Numbers in circles show fields: 1 - apatite-magnetite, 2 - vermiculite-phlogipite, 3 - apatite-shtaffelite, 4 - ore olivinites.
Рис. 4. Схема инженерно-геологического районирования рудника «Железный». На врезке показано положение Ковдорского месторождения в пределах Кольского полуострова.
Fig. 4. Engineering geological zonation of the Zhelezny mining site. Inset shows the position of Kovdor ore deposit in the frame of Kola Peninsula.
Рис. 5. Схема локальных осей сжатия и растяжения (рудник «Железный»), Ковдорский массив. 1 - разлом - максимально активный в молодом поле напряжений; 2 - прочие разломы; 3 - однородно-нагруженные объемы и их номера; 4-5 - оси локальных главных нормальных напряжений: 4 - сжатия, 5 - растяжения; 6 -косинус угла наклона осей напряжений 0 .
Fig. 5. Scheme of local axis of compression and extension for the Zhelezny mining site of the Kovdor block. 1 - fault Ri which is most active in the young field of stresses; 2 - other faults; 3 - volumes subject to regular loading and their numbers; 4-5 - axes of local main normal stresses: 4 - compression, 5 - extension; 6 - cosine of stress axis inclination angle of 0 .
имущественно комбинаторно-наложенный тип распределения индивидуальных трещин и зон трещиноватости различных систем и парагенетических ассоциаций в массиве пород.
На Ковдорском бадделеит-апатит-магнетитовом месторождении трещинная тектоника и трещиноватость изучались с применением системы многопараметрической индексированной документации (азимутальные, морфометрические, геометрические и статистические характеристики), которая предоставляет широкие возможности дифференцированной статистической обработки с возможностью фильтрации выборки по заданному набору параметров. Применение этой методики позволило выявить закономерную изменчивость азимутальных характеристик парагенетической ассоциации систем прототектонических трещин в зависимости от положения относительно геометрического центра рудного штока. Выявленная закономерность положена в основу инженерно-структурного районирования месторождения по секториальному принципу (рис. 4). Каждый из восьми выделенных в пределах карьерного поля (см. стереограммы на рис. 4, сетка
Ламберта, нижняя полусфера) секторов имеет закономерно ориентированную по месту своего расположения прототектоническую парагенетическую ассоциацию, а также в большей или меньшей степени проявленную парагенетическую(ие) ассоциацию(ии) наложенной тектоники посткарбонатитовых этапов. Стереограммы секторов построены на выборках массовых замеров, включающих сотни структурных элементов, задокументированных авторами в 2008-2010 гг. и ранее - силами специалистов ФГУП ВИОГЕМ.
Банк данных по бороздам скольжения состоит из 273 векторов на плоскостях со следами смещений, в дальнейшем именуемых «реализованными плоскостями». Ориентация, генезис и возраст плоскостей различный: прототектонические и тектонические трещины, контакты даек и жил, плоскости рассланцевания и др. На 32 плоскостях зафиксированы по два векторы перемещений, из них в 14 случаях удалось установить последовательность их формирования (молодые / древние). Как правило, к более молодым бороздам скольжения относятся здесь тонкие следы перемещений с преобладанием сдвиговой составляющей, в то время как к древним - более мощные, грубые борозды, часто параллельные линии падения плоскости смести-теля, иногда сопровождающиеся переориентировкой минералов в планпараллельные структуры (гнейсовидные, рассланцованные, отчетливо трахитоидные и т.п.) согласно вектору перемещения. У этих грубых борозд знак смещения (взброс или сброс) часто определяется неуверенно.
Анализ векторов перемещений, согласно кинематическому методу (графический вариант), позволил выделить на исследованной площади 17 однородно-нагруженных объемов. Они ограничены по структурным признакам (разные крылья разломов, крупных трещин, бока даек карбонатитов и т.д.). При этом к одному объему относились группы близлежащих точек в случае, если замеры в этих точках не противоречили нахождению единого поля напряжений. В каждом из них восстановлена ориентация осей главных нормальных напряжений, плоскости действия максимальных касательных напряжений (табл. 1). Необходимо отметить, что в массиве пород, вскрытых карьером, широкое развитие получили волнистые трещины прототектонического генезиса, на которых значительно меняются элементы залегания. Векторы перемещений, измеренные на одной из таких плоскостей (рис. 6, А, замеры 13а, б, в), находятся в согласии с найденными осями а1 и а3, несмотря на разброс их ориентировок.
В каждом объеме восстановлена ориентация осей и плоскости действия максимальных касательных напряжений (табл. 1). Ранее было показано, что резкое изменение ориентировок осей главных нормальных напряжений происходит в разных крыльях разломов и служит признаком активности разлома в восстановленном поле напряжений [Сим, 1980, 2000]. Два объе-
Таблица 1. Ориентировка тектонических напряжений рудника «Железный» (Ковдорский массив)
Table 1. Orientation of tectonic stresses on the Zhelezny mining site (Kovdor massive)
№ объема Точки наблюдений Месторасположение точек наблюдений 01, азимут и угол падения, ° 02, азимут и угол падения, ° 03, азимут и угол падения, ° Tmax, азимут и угол падения, ° тип Tmax, азимут и угол падения, ° тип
1 1, 2 восточный борт 246Z30 91Z60 342Z10 317Z60 193Z43
lw, 34 pw, 29
2 3 восточный борт 286Z15 114Z75 194Z0 150Z80 55Z80
lw, 10 pw, 10
3 3а, 4 восточный борт 232Z45 58Z46 325Z4 106Z56 359Z64
lw, 27 psw, 32
4 5 восточный борт 16Z30 220Z60 110Z10 248Z62 150Z80
lw, 18 ps, 25
5 7а северо-восточный и 258Z25 356Z20 118Z58 272Z73 38Z27
северный борт ls, 65 ps, 16
6 8 северо-восточный и 278Z80 62Z10 153Z5 144Z52 345Z42
северный борт lw, 50 pw, 40
7 9 северо-восточный и 30Z35 148Z36 270Z35 240Z89 148Z36
северный борт lw, 85 pw, 2
8 10, 11 северо-восточный и 231Z20 345Z44 125Z40 265Z78 7Z45
северный борт ls, 40 ps, 12
9 12 северо-восточный и 15Z30 214Z60 110Z10 244Z62 145Z78
северный борт lw, 14 pw, 30
10 13 юго-восточный борт 351Z70 102Z10 197Z18 188Z63 29Z28
lw, 60 pw, 25
11 15 юго-восточный борт 161Z70 308Z18 40Z10 25Z58 237Z38
lw, 52 pw, 35
12 16 юго-восточный борт 316Z25 123Z66 223Z15 180Z76 86Z70
lw, 20 pw, 15
13 17 юго-восточный борт 302Z40 176Z40 56Z30 186Z40 88Z85
lw, 12 pw, 50
14 18 западный борт 204Z30 296Z2 27Z60 205Z76 25Z15
15 19 западный борт 217Z40 46Z51 311Z5 s 91Z59 s 346Z68
lw, 25 pw, 30
16 20 западный борт 307Z26 154Z62 50Z30 176Z63 62Z80
lw, 10 pw, 25
17 22, 23 западный борт 102Z70 312Z20 218Z20 15Z42 234Z58
lw, 32 pw, 51
* Региональное поле 220Z65 4Z20 98Z15 82Z62 302Z 38
напряжений (1) lw, 53 pw, 32
** Региональное поле 4Z20 220Z65 98Z15 140Z84 233Z66
напряжений (2) pw, 35 lw, 7
Примечание. Азимут и угол погружения осей главных нормальных напряжений: а: - ось растяжения, а2 - промежуточная ось, а3 -ось сжатия; тшах - азимут и угол падения плоскости действия максимальных касательных напряжений; кинематический тип: w - взброс, 8 - сброс, 1 - левый сдвиг, р - правый сдвиг.
Note. Azimuth and dip angle of major normal stresses: 0^ - extension axis; o2 - transition axis; o3 - compression axis; Tmax - azimuth and dip angle of the plane of maximum tangential stresses; kinematic type: w - reserve fault; s - normal fault, l - sinistral strike-slip fault; p - dextral strike-slip fault.
ма (№ 5, т.н. 7а; № 14, т.н. 18), для которых определены крутые оси сжатия, расположены максимально близко к разлому Ш и могут являться подтверждением его активности (рис. 5).
Данные о локальных стресс-состояниях позволяют сделать предварительные выводы о возможной дополнительной роли тектонических напряжений в обрушении на юго-восточном борту карьера «Железный». Выполненные в 2009 г. структурно-тектонические исследования показали наличие плоскостей крупных трещин вдоль стенки обрушения, которые падают в сторону карьера. На этих плоскостях практически отсутствуют следы тектонических перемещений. Ближайшая от обрушения точка, в которой восстановлены
оси главных нормальных напряжений, показала, что вид напряженного состояния близок к одноосному растяжению (рис. 6, В, объем № 10, точка 13), что могло способствовать обрушению в сочетании с благоприятно ориентированными плоскостями трещиноватости. В 20-30 м на ВСВ от зоны обрушения в горных породах задокументированы трещины отрыва с характерной ребристой морфологией плоскости, на которой видно перистое расхождение «ребер», кальцитовая минерализация в породе и «дезинтеграция» карбона-тита неясного генезиса до состояния дресвы. Вышеперечисленные и не свойственные для остальной части карьера особенности находятся в согласии с локальной обстановкой растяжения. К западу от стенки обруше-
Рис. 6. Локальные стресс-состояния в однородно-нагруженных объемах (сетка Вульфа, верхняя полусфера): А - № 4 (табл. 1), точка 5; B -№ 10 (табл. 1), точка 13, юго-восточный борт карьера. 1 - оси главных нормальных напряжений и плоскости их действия: а - ^ -минимальных, б - а2 - промежуточных, в - а3 - максимальных сжимающих; 2 - полюса плоскостей действия максимальных касательных напряжений; 3 - векторы перемещений в точке полюса плоскости с бороздами скольжения: а и б - соответствующие найденным осям а1 и а3, а - с определенным, б - с неопределенным знаком смещения; в - противоречащие найденным осям а1 и а3. Цифрами показаны номера замеров векторов перемещений.
Fig. 6. Local states of stresses in the volumes subject to regular loading. The Wulff net, upper hemisphere: А - № 4 (Table 1), point 5; B - № 10 (Table 1), Point 13, the SE side of the quarry. 1 - axes and planes of main normal stresses: а - a1 - minimum, б - a2 - transitional, в - a3 - maximum compression; 2 - poles of planes of maximum tangential stresses; 3 - vectors of displacements at the pole point of the plane with striations: а and б - determined axes of a1 and a3, correspondingly (а - with determined sign, б - with undetermined sign); в - in contradiction to determined axes of a1 and a3. Numbers correspond to numbering of displacement vector measurements.
ния восстановлено локальное стресс-состояние с трехосным напряженным состоянием.
Данные о локальных стресс-состояниях позволили найти два варианта общего поля напряжений. Первый вариант (№ 1) для всего карьера имеет следующие характеристики: ось растяжения (а1) крутая, азимут погружения - 220^65°; ось сжатия (а3): 98^15°, проме жуточная ось (а2): 4°Z20° (рис. 7, табл. 1). Плоскости действия максимальных касательных напряжений т^
- азимут падения 82°Z62°, левый взбросо-сдвиг и азимут падения 302°Z38°, правый взбросо-сдвиг. Вторая плоскость тщах достаточно хорошо соответствует разлому И3, имеющему азимут падения 303-305°Z35°. Этому общему полю противоречит единственное определение оси растяжения в объеме № 2 (восточная стенка карьера), но, учитывая некоторую ненадежность этого определения из-за малого числа борозд скольжения и их слабой выраженности, этим противоречием можно пренебречь.
Второй вариант общего поля напряжений (№ 2) характеризуется переиндексацией промежуточной оси и оси сжатия, соответственно, оно характеризуется следующими параметрами: ось растяжения - 4°Z20°, промежуточная ось крутая - 220°Z65°; ось сжатия (а3): 98°Z15° (рис. 8, таблица 1). Плоскости действия максимальных касательных напряжений тщах - 140°Z84°,
правый взбросо-сдвиг с примерно равными сдвиговой и взбросовой составляющими и азимут падения 233°Z66°, левый взбросо- сдвиг с преобладанием сдвиговой составляющей смещения. Вторая плоскость тщах близка по положению в пространстве упомянутому выше разлому Я1. Оси растяжения общего поля № 2 противоречат оси локального сжатия в объемах № 1 и № 2; оси сжатия в объемах № 3, 6 и 10 в принципе находятся в пределах ошибки определения осей (и измерений векторов перемещения). Обращает на себя внимание то, что все противоречащие сдвиговому общему полю определения локальных стресс-состояний находятся в восточном борту и, возможно, отвечают более древнему этапу развития.
Весьма важные уточнения относительно возрастных взаимоотношений взбросового и сдвигового полей напряжений внес анализ поясного распределения реализованных плоскостей. Произведенный расчет плотности плоскостей с векторами перемещений (рис. 9) показал, что эти плоскости образуют отчетливый пояс РР1 с полюсом Р1 (аз. пад. 4^20°) и что пояс образован двумя максимумами плотностей векторов, ограниченных изолинией 2 %, вытянутыми вдоль пояса. Внутри максимумов, в свою очередь, проявлены по два локальных максимума. Так, в 1-м максимуме с плоскостями ССЗ простирания выделяются локальные
Рис. 7. Стереограмма общего поля напряжений № 1 (р-к «Железный»). 1а-3а - оси главных нормальных напряжений в однородно-нагруженных объемах и плоскости действия осей общего поля: 1 -растяжения, 2 - промежуточная, в - сжатия; 1б и 3б - оси сжатия (1б) и растяжения (3б) однородно-нагруженных объемов, 2б - полюса плоскостей действия максимальных касательных напряжений; 4 - конусы: а - сжатия, б - промежуточных осей, в - растяжения; 5 - простирание плоскостей действия максимальных касательных напряжений и кинематические типы: а - взбросо-сдвиги, б - сдвиго-взбросы.
Fig. 7. Stereogram of the general stress field № 1 (Zhelezny mining site). 1а-3а - axes of main normal stresses in the volumes subject to regular loading, and planes of axes of the general field: 1 - extension, 2 - transitional, в - compression; 16 and 3б - axes of compression (1б) and extension (3б) of the volumes subject to regular loading; 26 - poles of planes of maximum tangential stresses; 4a - compression cone; 4б -cone of transitional stresses; 4в - extension cone; 5 - strike of planes of maximum tangential stresses and kinematic types: 5a - strike-slip faults with reverse component, 56 - reverse faults with shear component.
максимумы плоскостей с аз. пад. 220^75° и аз. пад. 264^80°, которые оконтуриваются изолиниями в 3 и 6 %; во 2-м максимуме с плоскостями СВ простирания два локальных максимума, которые оконтурены изолинией в 2.5 %, имеют: аз. пад. 144°Z80° (восемь плоскостей) и 115°Z75° (семь плоскостей). Два локальных максимума (с плотностью в 2.5 % (координаты 144°Z80°) и с плотностью 3 % (координаты 220°Z75°) близки к плоскостям действия касательных напряжений ттах на плоскости действия оси растяжения найденного взбросового поля.
Ось пояса реализованных плоскостей с бороздами скольжения (условной трещиноватости, так как векторы фиксируются и на плоскостях рассланцевания, и на контактах жил и т.д.) близка к положению промежуточной оси сдвигового поля напряжений № 2. На рис. 8 промежуточная ось сдвигового поля напряжений находится вблизи оси пояса трещин, что может указывать, согласно С.И. Шерману [Шерман, 1966], на образование пояса трещин при перемещении по раз-
Рис. 8. Стереограмма общего поля напряжений № 2. Условные обозначения на рис. 7.
Fig. 8. Stereogram of the general stress field № 2. Legend same as in Figure 7.
лому И1 именно в этом сдвиговом поле напряжений.
По методике В.Н. Даниловича соединяются точки простирания разлома И1 на внешнем круге стереограммы с точкой Р1 (полюс пояса трещин) и восстанавливается генеральное положение плоскости данного разлома (аз. пад. 152°Z75°), что согласуется с данными картирования этого разлома в бортах карьера. Расхождение вычисленных параметров разлома на 10° от данных измерения плоскости сместителя в карьере допустимо из-за его сложного строения с разветвлением плоскости в приповерхностной части и сложной морфологии. По кинематическому типу он является правым взбросо-сдвигом с явным преобладанием сдвиговой компоненты и с поднятым юго-восточным крылом. Точка К1, лежащая на пересечении главной кинематической плоскости (или плоскости действия оси промежуточных напряжений) с плоскостью разлома, показывает выход векторы К1О - направление перемещения по плоскости разлома в поле напряжений № 2 - на верхнюю полусферу; отрезок дуги КК1 является интервалом разброса возможных перемещений по разлому с учетом данных о поясном распределении реализованных плоскостей и реконструкции палеонапряжений. В более древнем взбросовом поле № 1 разлом И1 при падении на ЮВ являлся сдвиго-взбросом с подчиненной право- или левосдвиговой компонентой перемещения в зависимости от величины ра.
Общепринято, что большинство борозд скольжения отражает следы последних перемещений на зеркалах скольжения. Соответственно, концентрация реализованных плоскостей по всему карьеру вблизи плоскостей действия максимальных касательных напряжений сдвигового общего поля напряжений подтверждает сделанные ранее выводы о его более молодом возрас-
Рис. 9. Поясное распределение полюсов плоскостей с бороздами скольжения и определение кинематического типа разлома на руднике «Железный», Ковдорский массив (сетка Вульфа, верхняя полусфера). 1 - пояс реализованных плоскостей РР1 и его ось (полюс) Р1; 2 - плоскость разлома и его полюс R1; 3 - векторы перемещений в плоскости разлома, определенные согласно: КО - поясному распределению трещиноватости, К1О - ориентировке осей главных нормальных напряжений общего поля II; 4-6 - оси максимальных (4), промежуточных (5) и минимальных (6) сжимающих главных нормальных напряжений II поля; 7 - полюса плоскостей действия максимальных касательных напряжений; 8 - изолинии плотности трещин с бороздами скольжения.
Fig. 9. Belt-type distribution pattern of poles of planes with striations, and determination of the kinematic type of the fault at the Zhelezny mining site of the Kovdor block. The Wulff net, upper hemisphere. 1 -belt of actual planes of РР1 and axis (pole) of Р1; 2 - fault plane and pole of R1; 3 - vectors of displacements at fault plane that are determined from: КО - belt-type fracturing pattern, К1О - orientation of axes of main normal stresses of General Stress Field II; 4-6 - axes of maximum (4), transitional (5) and minimum (6) compression of main normal stresses of General Stress Field II; 7 - poles of planes of maximum tangential stresses; 8 - isolines of density of fractures with striations.
те. В свою очередь, отсутствие максимумов реализованных плоскостей, близких по положению в пространстве к плоскостям максимальных касательных напряжений взбросового поля, также не противоречит тезису о его более древнем возрасте. Общая северовосточная ориентация субгоризонтальных осей сжатия в лежачем крыле прототектонической трещины (разлома), по которому произошло обрушение, а также крутые оси растяжения при условии, что т.н. 13 близка к СВ концу разлома, соответствуют правосдвиговому перемещению по этой плоскости в общем молодом поле напряжений. При этом положение плоскости обрушения по отношению к плоскостям действия максимальных касательных напряжений этого последнего поля значительно менее благоприятно для активизации по сравнению с положением разлома
Из сопоставления данных по трещиноватости с «реализованными плоскостями» очевидно, что в новейшем поле напряжений из множества систем разно-
возрастных трещин происходит их избирательная активизация - борозды скольжения образуются только на трещинах, благоприятно ориентированных относительно осей общего молодого поля напряжений.
Изучение разломной тектоники, трещиноватости и тектонических напряжений на руднике «Железный» дало возможность выделить этапы хрупкопластического деформирования исследованной части Ковдорского массива (рис. 10). Первый этап соответствует начальной стадии образования прототектониче-ских трещин; второй - ранний посткарбонатитовый (после завершения магматизма с ранними карбонати-тами). К этой стадии относится внедрение даек поздних карбонатитов и формирование разлома Иь крутое падение которого дает возможность предположить сдвиговое поле на этом этапе; по преобладающей ориентировке даек карбонатитов можно предположить субширотную ориентировку оси растяжения. Последние стадии С и Б восстановлены с помощью изучения тектонических напряжений.
Хибинский щелочной массив приурочен к тектоническому контакту Имандра-Варзугской зоны карелид (на юге и западе массива) и архейских гранитогнейсов Центрально-Кольского домена. Он локализован в трансрегиональной северо-восточной зоне, пересекающей весь Кольский п-ов и протягивающейся до грабена Осло. К пересечению долгоживущих глубинных разломов этой зоны с разломами других простираний тяготеют интрузии щелочно-ультраосновных пород и нефелиновых сиенитов. Хибинский массив является многофазным эллипсовидным в плане плутоном, длинная ось которого протяженностью 45 км вытянута по азимуту В 82; меридиональная ось равна 35 км (рис. 11). Морфологически он близок к асимметричному лополиту, у которого крутые контакты - восточный и северный, более пологие - южный и западный. На глубине он переходит в коническую структуру центрального типа. Внутреннее строение массива многофазное концентрически-зональное, от периферии к центру породы разного состава сменяют друг друга [Геология..., 2002].
Основные тектонофизические работы были проведены на месторождении Олений Ручей (рис. 11) в восточной части Хибинского массива нефелиновых сиенитов, на котором выделены четыре локальных стресс-состояния кинематическим методом и два - методом выделения сопряженных сколов (табл. 2).
Общим для этих стресс-состояний является пологая ось сжатия, ориентировка которой меняется от субши-ротной (точки 1 и 7а) до ССЗ (точки 6, 7 и ЭР) и суб-меридиональной до СВ растяжения, угол наклона которой меняется от пологого (т.н. ЭР ) через наклонное (т.н. 7) до крутого (т.н. 6). Исключением является стресс-состояние в т.н. 8, расположенной на ЮВ крутом склоне г. Суолуайв южнее рудника. Локальное стресс-состояние в этой точке характеризуется пологой СВ ориентировкой оси сжатия и СЗ пологой - рас-
Рис. 10. Этапы хрупко-пластического деформирования интрузии с апатит-магнетитовым месторождением (Ковдорский массив) [Жиров, Сим, 2010]: А - этап формирования прототектонических трещин; В - этап образования разлома Rb первых систем тектонических трещин, внедрение даек поздних карбонатитов; C - первый постмагматический этап, соответствующий установленному тектонофизи-ческими методами 1-му полю напряжений; D - неотектонический этап, соответствует 2-му сдвиговому полю напряжений.
Fig. 10. Stages of brittle-plastic deformation of the intrusion with apatite-magnetite deposit (Kovdor block) [Жиров, Сим, 2010]: А - phase of prototectonic fracturing; В - phase of formation of fault Rb the first systems of tectonic fractures, penetration of dykes of late carbonatites; C -the first postmagmatic phase corresponding to Stress Field I which is revealed by tectonophysical methods; D - neotectonic phase which corresponds to Stress Field II (strike-slip).
тяжения, практически соответствующей ориентировке оси сжатия в т.н. 7 и ЭР (табл. 2, рис. 12). Такая пере-индексация локальных осей главных нормальных напряжений характерна для разных крыльев разломов. Гора Суолуайв расположена в лежачем ЮЮВ крыле обновленного сложнопостроенного разлома, относящегося к краевым разломам Хибинского массива, частично вскрывающегося в т.н. 7. В последней, где вскрышные работы на карьере обнажили первые апатитовые тела, разлом представлен зоной дробления и интенсивной шпреуштейнизацией. Крутая ось растяжения в т.н. 6, возможно, также обусловлена влиянием небольшого разлома с азимутом простирания СЗ 305325°, который разработан левым притоком Оленьего ручья.
На рис. 12 (точка ЭР) представлен вариант определения локального стресс-состояния по трещинам в четвертичных водно-ледниковых отложениях, вскрытых в восточном борту дороги при въезде на карьер.
Реконструкция тектонических напряжений в т.н. 7 также была произведена с помощью выделения сколо-вых сопряженных трещин (объем 7а). Угол наклона осей сжатия и растяжения стресс-состояний, определенных разными методами, сохраняется пологий, но субмеридиональная ось растяжения разворачивается от меридиональной (в определении № 7а 13^20°) до ЮЗ (в определении № 7 на 226°Z30°). Азимут погружения оси сжатия меняется от ЗСЗ 283° до СЗ 320°. Сравнение результатов реконструкции локальных стресс-состояний на руднике «Олений Ручей», несмотря на малое количество определений, показывает в целом доминирование сдвигового поля напряжений с усредненной ЗСЗ ориентировкой оси сжатия и субме-ридиональной осью растяжения (пять определений из шести). Относительная непротиворечивость определений локальных стресс-состояний, восстановленных в коренных и четвертичных породах, позволяет говорить, по сути, о современном поле напряжений, вос-
Таблица 2. Ориентировка тектонических напряжений Кольского полуострова
Table 2. Orientation of tectonic stresses of the Kola Peninsula
№ Точка наблю- дений Район a1, азимут и угол падения, ° a2, азимут и угол падения, ° a3, азимут и угол падения ° Tmax. азимут и угол падения, °, тип Tmax. азимут и угол падения, °, тип
1 1 Олений ручей (карьер «СЗФК») 218Z20 338Z54 116Z3Q 346Z54 ps 258Z85 ls
2 3 г. Поачвумчорр (к северу от г. Кировска) 35QZ2Q 96Z3Q 232Z55 14Z6Q ls 134Z35 ps
3 4 г. Мурманск, карьер 118Z10 254Z76 27Z1Q 242Z88 lw 252Z76 p
4 5 Пулозеро, в выемке дороги М18 219Z6Q 61Z28 326Z1Q 115Z43 lw 35QZ59 pw
5 6 Олений ручей (карьер «СЗФК») 336Z88 242Z2 154Z2 332Z4Q w 152Z48 w
6 7 Олений ручей (карьер «СЗФК») 226Z3Q 67Z56 32QZ1Q 99Z6Q lw 359Z74 pw
7 7a Олений ручей (карьер «СЗФК») 13Z2Q 144Z6Q 283Z25 63Z87 ls 19QZ86 lw 144Z6Q
8 8 Олений ручей (естественное обнажение) 324Z13 234Z1Q 1Q7Z74 p 1QQZ72 p
9 9Q Олений ручей (карьер «СЗФК») 6QZ2 141Z72 32QZ2Q 184Z78 96Z76 l 1QQZ83 p
1Q 1Q Рудник «Центральный» (верхняя часть) 324Z2Q 174Z67 58Z1Q p 194Z7Q l
11 11 Рудник «Восточный» (Коашва) 29Z18 296Z6 197Z68 22QZ23 26Z66 ls 226Z43 ls
12 11а Рудник «Восточный» (Коашва) 16Z1Q 281Z28 124Z6Q s 353Z58 ps
13 11Q Рудник «Восточный» (Коашва) 282Z17 138Z70 16Z12 58Z86 pw 15QZ7Q L
14 13 Рудник «Центральный» (средняя часть, север) 57Z15 182Z64 321Z2Q 188Z63 p 98Z86 ls
15 13а Рудник «Центральный» (средняя часть, юг) Региональное поле напряжений 91Z15 356Z6Q 354Z14 254Z8 239Z6Q 159Z3Q 8QZ7Q ps 164Z74 pw 296Z23 ls 314Z17 lw
Примечание. Азимут и угол погружения осей главных нормальных напряжений: а! - ось растяжения, а2 - промежуточная ось, а3 -ось сжатия; Ттах - азимут и угол падения плоскости действия максимальных касательных напряжений; геолого-кинематические типы: w -взброс, 8 - сброс, 1 - левый сдвиг, р - правый сдвиг.
Note. Azimuth and dip angle of major normal stresses: 0^ - extension axis; o2 - transition axis; o3 - compression axis; Tmax - azimuth and dip angle of the plane of maximum tangential stresses; kinematic type: w - reserve fault; s - normal fault, l - sinistral strike-slip fault; p - dextral strike-slip fault.
становленном кинематическим методом.
Тектонофизические исследования на месторождении Расвумчорр (рудник «Центральный») проводились в его южной части на горизонтах +835 и +745 южного борта карьера в разных крыльях зоны субширотного разлома, оперяющего Главный разлом. Азимутальные характеристики изученного разлома (оперяющего по отношению к Главному разлому) варьируются от аз. пад. 19^35° на верхних уступах до 344^35-40° на нижнем, а его кинематический тип по бороздам скольжения определен как правый взбросо-сдвиг. Тип Главного разлома по геолого-структурным данным устанавливается как взбросо-надвиг с неоднократной циклической инверсией направления смещения. Всего на руднике выделено три локальных стресс-состояния (табл. 2, т.н. 10, 13 и 13а). Локальные стресс-состояния
в висячем (рис. 12, т.н. 10, горизонт +835) и лежачем (т.н. 13, горизонт +745) крыльях характеризуются пе-реиндексацией осей сжатия и растяжения. На рис. 13 показаны полюс разлома И, его плоскость И^-Иг и вектор перемещения ОА, зафиксированный в карьере на верхнем уступе. Оси локальных стресс-состояний в т.н. 10 и 13 лежат в разных крыльях, видна их переин-дексация. При этом они близки к плоскости разлома и стремятся занять положение, параллельное и перпендикулярное вектору перемещения разлома, что хорошо согласуется с данными математического и физического моделирования [Осокина, 1987, 1997] и полевыми исследованиями разломов [Сим, 1980, 1996].
На стереограмму (рис. 13) вынесены оси сжатия и растяжения общего поля напряжений (табл. 2), согласно которому перемещение по разлому должно быть
Рис. 11. Геологическая карта Хибинского массива. Составлена В.П. Павловым и др. (по: [Геология.., 2002]). 1 - нефелиновые фениты; 2 - роговики различного состава; 3 - щелочно-ультраосновные породы; 4 - трубки взрыва; 5 - фойяиты массивные; 6 - фойяиты тра-хитоидные; 7 - щелочные сиениты; 8 - лявочорриты; 9 - хибиниты массивные; 10 - хибиниты трахитоидные; 11 - хибиниты нечетко трахитоидные; 12 - апатит-нефелиновые руды; 13 - ювиты; 14 - рисчорриты массивные; 15 - рисчорриты гнейсовидные; 16 - уртиты массивные; 17 - ийолит-уртиты гнейсовидные; 18 - карбонатиты; 19 - месторождения и рудопроявления апатит-нефелиновых руд (1 -Валепахк-Намуайв, 2 - Лявочорр, 3 - Партомчорр, 4 - Куэльпор, 5 - Снежный цирк, 6 - Кукисвумчорр, 7 - Юкспор, 8 - Апатитовый цирк, Э - Расвумчорр, 10 - Эвеслогчорр, 11 - Коашва, 12 - Вуоннемйок, 13 - Ньоркпахк, 14 - Олений Ручей); 20 - Тульокская зона разломов; 21 - разрывные нарушения, установленные по геофизическим данным и аэрофотоснимкам; 22 - установленные разрывные нарушения (названия рудников с тектонофизическими данными выделены курсивом).
Fig. 11. Geological map of the Khibim block. The map is compiled by V.P. Pavlov et al. based on [Геология...., 2002]. 1 - nepheline fenite; 2 -hornblende varying in composition; 3 - alkaline-ultrabasic rocks; 4 - volcanic pipes; 5 - trachytoid foyaites; 7 - alkaline syenites; 8 - lyavor-chites; 9 - massive khibinites; 10 - trachytoid khibinites; 11 - indefinite trachytoid khibinites; 12 - apatite-nepheline ores; 13 - yuvites; 14 -gneissoid rischorrites; 15 - gneissoid rischorrites; 16 - massive urtites; 17 - gneissoid yolite-urtites; 18 - carbonatites; 19 - deposits and mineral occurrences of apatite-nepheline ores (1- Valepakh-Namuaiv, 2 - Lyavochorr, 3 - Partomchorr, 4 - Kuaelpor, 5 - snow cirque, 6 - Kukisvum-chorr, 7 - Yukspor, 8 - apatite cirque, 9 - Rasvumchorr, 10 - Eveslogchorr, 11 - Koashva, 12 - Vuonnemjok, 13 - Njorkpakh, 14 - Deer Brook); 20 - Tuljok fault zone; 21 - faults determined from geophysical data and aerial photos; 22 - determined fractures (names of mining sites with tectonophysical data are shown in italics).
взбросо-сдвигом по линии ОС при одноосном сжатии и по линии ОВ - при одноосном растяжении. Часть плоскости разлома ВС соответствует возможным выходам векторов перемещения по плоскости разлома в найденном общем поле напряжений при разных величинах коэффициента ра. Вектор ОА находится в створе угла ВОС. Из этого следует два вывода: 1 - зафиксированное смещение по разлому согласуется с общим полем напряжений; 2 - вид напряженного состояния
найденного общего поля напряжений близок к трехосному по положению реального векторы смещения примерно на равном расстоянии между векторами ВО и ОС.
На руднике «Коашва» выделены три однородно напряженных объема с определением их локальных стресс-состояний: два определения (т.н. 11, 11б) относятся к коренным расвумчорритам, одно (т. 1Щ) -к четвертичным флювиогляциальным отложениям,
Рис. 12. Стереограммы локальных стресс-состояний Хибинского массива (сетка Вульфа, верхняя полусфера). 1-4 - векторы перемещений на зеркалах скольжения: 1-3 - соответствующие - 1-2 - со знаком смещения, измеренные разными авторами, 3 - без знака смещения, 4 - противоречащие найденному локальному стресс-состоянию; 5-11 - оси главных нормальных напряжений и плоскости их действия: 5-7 - локальных стресс-состояний, восстановленных по измерениям в магматических, 8-9 - в четвертичных породах, 10-11 -общего поля напряжений; Ст! - ось минимальных (растяжения), <г2 - промежуточных, ст3 - максимальных (сжатия) сжимающих напряжений; 12 - вектор перемещения по плоскости действия максимальных касательных напряжений в точке полюса; 13-17 - простирания плоскостей действия максимальных касательных напряжений и их кинематические типы: 13 - взбросы, 14 - сбросы, 15 - сдвиги, 16 - взбросо-сдвиги, 17 - сбросо-сдвиги; двойными линиями выделена преобладающая составляющая перемещений, зубцы ориентированы по падению плоскостей; 18 - изолинии плотности трещин; 19 - сопряженные сколовые трещины.
Fig. 12. Stereograms of local states of stresses of the Khibin block. The Wulff net, upper hemisphere. 1-4 - vectors of displacements on slicken-sides: 1-3 corresponding to the local state of stresses (1-2 - with determined sign if displacement as per measurements by various authors, 3 -without displacement sign), 4 - contradicting with the local state of stresses; 5-11 - axes and planes of main normal stresses: 5-7 - local states of stresses determined from measurement in magmatic rocks, 8-9 - in the Quaternary rocks, 10-11 - general stress field; - axis of minimum
stresses (extension), ct2 - transitional stresses, ct3 - maximum stresses (compression); 12 - vector of displacement on plane of maximum tangential stresses at pole point; 13-17 - strike of planes of maximum tangential stresses and kinematic types of stresses: 13 - reverse faults, 14 - normal faults, 15 - strike-slip faults, 16 - strike-slip faults with reverse component, 17 - strike-slip faults with normal component. Dual lines show dominating components of displacement; «waves» are oriented by strikes of planes; 18 - isolines of density of fractures; 19 - conjugated shear fractures.
Рис. 13. Стереограмма, иллюстрирующая связь перемещений и осей главных нормальных напряжений вблизи разлома на руднике «Центральный». 1 - полюс (R) и плоскость разлома R: R2, на которой ВС - интервал возможного выхода векторов перемещений по разлому в общем поле напряжений Хибин при разном виде напряженного состояния («Центральный», горизонты +835 и +745); 2-3 - вектор перемещения по разлому: в точке полюса (2), на плоскости разлома (3); 4-5 - оси минимальных (4) и максимальных (5) сжимающих главных нормальных напряжений общего поля Хибин; 6-7 - оси минимальных (6) и максимальных (7) сжимающих главных нормальных напряжений локальных стресс-состояний; цифрами даны номера точек наблюдений; 8-10 - дуги больших кругов, соединяющих полюс разлома и ось растяжения (8), ось сжатия общего поля (9) и вектор перемещения по разлому (10).
Fig. 13. Stereogram showing the relationship between displacement and axes of main normal stresses near the fault on the Central mining site. 1 - pole (R) and plane of fault R: R2, where “ВС” is an range of potential outlet of displacement vectors on the fault in the general stress field of the Khibin with various states of stresses (Central, horizons +835 и +745); 2-3 - vector of displacement along the fault: at pole point (2), and on fault plane (3); 4-5 - axes of minimum (4) and maximum (5) compression main normal stresses of the general stress field of the Khi-bin; 6-7 - axes of minimum (6) and maximum (7) compression main normal stresses of local states of stresses. Numbers show points of measurements; 8-10 - arcs of large circles connecting the fault pole and extension axis (8), compression axis of the general field (9), and vector of displacement along the fault (10).
вскрытым в северо-западном борту карьера. Оси растяжения в т.н. 11 и 11а, разделенных разломом с азимутом падения 130-145°^70-80°, субгоризонтальные и субмеридиональные, а оси сжатия круто погружаются на юго-восток, обусловливая сдвиго-сбросовое поле напряжений в исследованной части рудника.
Сдвиговое поле напряжений в т.н. 1Щ характеризуется субмеридиональным близгоризонтальным сжатием (рис. 13), которое практически совпадает с осями растяжения, восстановленными в коренных породах, а пологая ось растяжения субширотна. Переиндексация осей главных нормальных напряжений в разнесенных в пространстве точках тектонофизических исследований может свидетельствовать о разном возрасте вос-
становленных напряжений. Но так как такая переин-дексация обычно характерна для разных крыльев разломов, можно предположительно говорить о связи этого явления с активизацией в молодом поле напряжений краевого разлома, в разных крыльях которого расположены т.н. 11 и 1Щ.
Трасса Мурманск - Санкт-Петербург. В архейских гранитогнейсах (кольско-беломорский нестратифици-рованный комплекс [Геологическая карта2001]) Мурманского блока реконструкция тектонических напряжений была произведена в карьере на западной окраине г. Мурманска (см. рис. 12, т.н. 4). В западной стенке карьера выделяется зона небольшого крутопадающего разлома с аз. пад. 205^85°. Ширина зоны дробления составляет 80-90 см, породы в ней обохре-ны и обводнены. В разных крыльях разлома измерения борозд скольжения произведены разными исследователями. На рис. 12 (т.н. 4) замеры измерения в северном и южном крыльях, произведенные разными исследователями, выделены в точках полюсов трещин отличающимися значками (стрелки с крестиком или с кругом). Все векторы перемещений идеально соответствуют сдвиговому полю напряжений с субгоризон-тальными осями ССВ сжатия и ЗСЗ растяжения. Это свидетельствует об отсутствии новейшей активизации разлома.
Другая реконструкция осуществлена по измерениям борозд скольжения в лопийских эндербитах [Геологическая карта., 2001]), вскрытых в восточном борту трассы Мурманск-Санкт-Петербург между пос. Пуло-зеро и г. Оленегорск (т.н. 5, табл. 2). Локальное стресс-состояние характеризуется взбросовым полем с крутой осью растяжения и осью сжатия, полого погружающейся на северо-запад.
6. ОБЩЕЕ ПОЛЕ НАПРЯЖЕНИЙ
Все полученные оси сжатия и растяжения локального уровня были вынесены на одну стереограмму (рис. 14), на которой виден значительный разброс в их ориентировках. Заметим, что в анализе участвуют стресс-состояния, восстановленные в породах архейского, девонского и четвертичного возраста. Тем не менее находится такое положение конусов сжатия и растяжения с телесным углом при вершине 45° и взаимно перпендикулярными осями, в которое не попадают противоположные локальные оси: в конусе растяжения нет осей сжатия и, наоборот, в конусе сжатия нет осей растяжения. Это свидетельствует, согласно основным принципам метода нахождения общего поля напряжений [Сим, 1980, 1996], о едином поле напряжений, в котором деформируются породы Хибин в неотектонический этап. Распространять это поле на весь Кольский п-ов представляется некорректным, хотя в его реконструкции и участвуют определения локальных стресс-состояний за пределами Хибин (из-за
T
Рис. 14. Стереограмма общего поля напряжений Хибинского массива. Условные обозначения приведены на рис. 11.
Fig. 14. Stereogram of the general stress field of the Khibin block. Legend same as in Figure 11.
малого числа последних). Пологая ось сжатия ССЗ простирания погружается на ЮЮВ по азимуту 159°^30°, крутая ось растяжения - 356°^60°, плоскости действия максимальных касательных напряжений являются взбросами с подчиненной горизонтальной компонентой перемещений, которая на плоскости 164°Z74° - правосдвиговая, а на плоскости 314°Z17° -левосдвиговая.
Локальные оси главных нормальных напряжений в т.н. 7 (хибиниты рудника «Олений Ручей») и в водноледниковых отложениях т.н. 9Р практически совпадают (см. рис. 13), а оси растяжения в этих точках имеют общее СВ простирание, но в хибинитах она более крутая. В целом принято, что расхождение в ориентации локальных осей напряжений в близлежащих точках наблюдения в разновозрастных породах незначительно. В таком случае в хибинитах восстановлено молодое поле напряжений (см. рис. 11, 13, табл. 2).
Крутое положение оси растяжения, определяющее в целом взбросовый тип общего поля напряжений, более уверенно относится только к Хибинскому массиву, в котором восстановлены 12 из 14 определений локального стресс-состояния. Возраст восстановленных по бороздам скольжения полей напряжений принимается за неотектонический и современный, так как в определении общего поля напряжений участвуют и локальные стресс-состояния в четвертичных породах. Воздымание Хибинского массива в неотектонический этап, сопровождающееся повышенной сейсмичностью
и фиксирующееся в современном рельефе, согласуется с найденным общим взбросовым полем напряжений.
Малая концентрация осей сжатия и растяжения в соответствующих конусах свидетельствует о трехосном напряженном состоянии общего поля напряжений, что согласуется с выводом, полученным ранее при анализе полей напряжений и векторы перемещения по разлому на руднике «Центральный» (см. рис. 12).
7. Обсуждение результатов
Приведенные результаты тектонофизических исследований позволяют говорить о двух общих полях напряжений на исследованной части Ковдорского массива. Более древнее поле № 1 взбросовое с субгори-зонтальной осью сжатия ЗСЗ простирания, более молодое поле № 2 - сдвиговое с тем же сжатием, но с переиндексацией оси растяжения и промежуточной оси. Поле напряжений № 2 хорошо коррелируется с региональным неотектоническим полем напряжений восточной части Балтийского (Фенноскандинавского) щита, установленным ранее [Сим, 199в, 2000; Sim et al., 2010]. Обрушение в юго-восточном борту карьера обусловлено в основном наличием благоприятно ориентированных плоскостей крупных трещин и, вероятно, возможным незначительным правосдвиговым перемещением с обстановкой локального растяжения на этом участке. Правосдвиговые перемещения также нашли подтверждение в южном борту карьера (l50-200 м от плоскости обрушения), где задокументирована серия прототектонических трещин и даек карбона-титов с правосдвиговой кинематикой.
Из-за переориентировки минералов на плоскостях сместителей субпараллельно линии перемещения, требующей особых РТ условий, можно предположить, что возраст взбросового поля напряжений синхронен возрасту карбонатитов последней стадии. Для определения этого возраста требуются дополнительные исследования.
В результате полевых тектонофизических исследований на рудниках Хибинского массива выделены l4 однородно нагруженных объемов с определением локальных стресс-состояний. Они характеризуются, так же как и на Ковдоре, существенной изменчивостью ориентировок осей главных нормальных напряжений. Наибольшее число определений относится к руднику «Олений Ручей», на котором локальные стресс-состояния восстановлены в коренных породах как по бороздам скольжения, так и по сколовым сопряженным трещинам в хибинитах, а также в четвертичных флювиогляциальных песках. По всем реконструкциям локальных стресс-состояний определено общее поле напряжений взбросового типа.
На карьере рудника «Железный» показана относительная активность крупных разломов, по-разному
ориентированных к современному полю напряжений. Разлом Кь близкий по положению к одной из плоскостей действия максимальных касательных напряжений последнего поля, наиболее активный.
На устойчивость бортов карьера влияет не только активность разломов, но и наличие крупных ослабленных плоскостей, неблагоприятно ориентированных относительно выработанного в процессе разработки карьера пространства, а также вид напряженного состояния в локальных объемах горных пород. Современная активность разломов на рудниках Хибинского массива фиксируется по резкой смене ориентировок осей локальных напряжений в разных крыльях.
Как известно из теоретического моделирования напряжений, связанных со смещением по разломам, вблизи зон сдвигов ориентировки осей напряжений внешнего поля (собственно, обусловливающего смещение) плавно изменяются по простиранию, достигая предельной переориентации на концах разлома. Во фронтальной части движущегося крыла ориентировка оси сжатия параллельна, а оси растяжения - перпендикулярна плоскости сдвига; в тыльной части противоположного крыла, наоборот, ось растяжения параллельна, а ось сжатия перпендикулярна плоскости сдвига [Осокина, 1987]. Данные о природных напряжениях полностью согласуются с теоретическими на разломе рудника «Центральный» в Хибинах и частично - на разломе (крупной трещине) в зоне обрушения на руднике «Железный». Поскольку перечисленные разломы не выходят за пределы магматических массивов, удалось зафиксировать возмущения локальных напряжений на концах этих разломов. Разлом К1 на руднике «Железный» продолжается по данным дешифрирования за границу карьера, в пределах которого вскрыта предположительно срединная часть разлома. По всей вероятности, разнообразие ориентировок в разных крыльях разлома Кь не согласующееся с теоретическими расчетами, обусловлено их значительно неоднородным геологическим строением. При сопоставлении данных моделирования сдвигов в однородной среде с данными активизации разломов в природных массивах необходимо иметь в виду значительно более сложное строение последних и более длительное развитие, сопровождавшееся, скорее всего, неоднократным обновлением разломов даже в едином для новейшего этапа общем поле напряжений. Косвенно это подтверждается тем, что вдоль сейсмогенных разрывов распределение амплитуд смещений может быть крайне неоднородным [Стром, Никонов, 1999]. В этом случае локальные напряжения могут быть так же неоднородны. Широкий спектр механизмов афтершоков после сильных землетрясений и разброс последних вдоль основного сейсмогенного шва тоже не противоречит «неупорядоченным» локальным напряжениям в крыльях разломов и нахождению по этим локальным стресс-состояниям общего поля. Имеющиеся данные по тектоническим напряжениям дают возможность
выделить только два ранга напряжений: локальные напряжения можно отнести ко II рангу, а общие - к I (или «региональному») рангу.
Неотектонический и современный возраст восстановленных по бороздам скольжения в архей-проте-розойских породах восточной части Балтийского щита напряжений установлен различными методами. В Центральной Карелии он определен по результатам специальных геохимических исследований. В пределах Кандалакшского залива он установлен в соответствии с новейшей структурой Кандалакшского грабена, ориентированного поперек «региональной» субго-ризонтальной оси растяжения. Последнее в сочетании с крутой осью «регионального» сжатия на этом участке исследований целиком согласуется со структурой грабена. Непротиворечивость локальных стресс-состояний, восстановленных в девонских магматических образованиях и голоценовых отложениях в пределах Хибинского массива, также дает основание относить возраст реконструированных полей напряжений к современному и неотектоническому.
На исследованной части Ковдорского массива выделяются четыре этапа деформирования: два более ранних - по структурным признакам, а два последних
- по ориентации осей главных нормальных напряжений общих полей. В отличие от Ковдора в Хибинах восстановлено единственное «региональное» поле. Отметим, что взбросовое поле напряжений 3-го этапа деформирования (и 1-е по реконструкции тектонических напряжений) на Ковдорском массиве похоже на взбро-совое поле на рудниках Хибинского массива. В обоих массивах ось сжатия субгоризонтальна. Предположительно крутая ось растяжения отражает этап хрупкого деформирования, проявившегося после остывания обеих интрузий на глубине, и последующее их возды-мание. Дальнейшая история развития привела массивы к разному современному напряженному состоянию. На неотектоническом и современном этапах механизм их деформирования может определяться двумя источниками субгоризонтального сжатия. Один из них -остаточные гравитационные напряжения, сохранившиеся в значительно эродированных магматических телах как «память» о пребывании в стесненных условиях на глубине [Rebetsky, 2008], а второй - возможное влияние процессов спрединга в Северной Атлантике. Современное сдвиговое поле на Ковдоре по ориентации оси сжатия не противоречит его обусловленности процессами раздвижения в Атлантике, хотя не исключает и наличия гравитационных напряжений. Хибинский массив, имеющий форму лополита, под воздействием горизонтального сжатия любого генезиса должен выдавливаться вверх, что фиксируется субвертикаль-ной осью растяжения «регионального» поля Хибин и подтверждается современным рельефом - наиболее высокие вершины Кольского п-ова находятся именно на этом массиве, он продолжает «всплывать».
8. Выводы
Проведенные на рудниках Кольского п-ова текто-нофизические исследования дали возможность сделать несколько выводов, имеющих теоретическое и практическое значение. Получены данные о современном напряженном состоянии, которое играет значительную роль в устойчивости горных выработок. На обоих массивах установлены «региональные» субгоризонталь-ные оси сжатия. При этом деформирование двух исследованных массивов на новейшем и современном этапах происходит по-разному: Ковдорский массив деформируется в сдвиговом поле напряжений с субго-ризонтальной запад-северо-западной ориентацией оси сжатия, не противоречащей влиянию рифтогенеза в Северной Атлантике, а Хибинский массив, имеющий лополитообразную форму, независимо от генезиса горизонтального сжатия, выдавливается наверх. Расчеты Ю.Л. Ребецкого величин остаточных гравитационных горизонтальных напряжений в Хибинах с использованием результатов полевых тектонофизических исследований, приведенных в данной работе, показали, что они хорошо коррелируются с данными инструментальных измерений напряжений в этом массиве и что роль дальнодействующих источников горизонтальных напряжений здесь практически исключена [Ребецкий и др., 2011].
9. Благодарности
Авторы выражают благодарность заместителю директора Горного института Кольского научного центра РАН профессору А.А. Козыреву и доценту этого института В.В. Рыбину за помощь в организации полевых работ, а также заведующему лаборатории текто-нофизики ИФЗ РАН Ю.Л. Ребецкому за обсуждение работы на всех стадиях и конструктивные замечания.
Работа выполнена при поддержке грантов РФФИ 09-05-12064 офи_м, 09-05-00687а, 09-05-1213а, 09-05-01022а.
10. Литература
Бабак В.И., Гарбар Д.И., Сим Л.А., Касаткин Ф.Г. Тектонические условия и тенденции развития Онежского озера // Изв. ВУЗов. Геология и разведка. 1979. № 2. С. 44-51.
Бабак В.И., Николаев Н.И. Пояснительная записка к карте геомор-фолого-неотектонического районирования Нечерноземной зоны РСФСР (за исключением горной части Урала и Калининградской области). М-б 1:500000. М.: МГУ, 1983.
Бабак В.И., Сим Л.А., Касаткин Ф.Г. Вариации тектонического поля напряжений восточной части Балтийского щита // Бюл. МОИП. Отд. геол. 1981. Т. 56, вып. 1. С. 150.
Геологическая карта Кольского региона. М-б 1:100000 / Ред. Ф.П. Митрофанов. Апатиты, 2001.
Геология рудных районов Мурманской области / В.И. Пожиленко, Б.В. Гавриленко, Д.В. Жиров, С.В. Жабин. Апатиты: КНЦ РАН, 2002. 359 с.
Гзовский М.В. Тектонические поля напряжений // Известия АН СССР. Серия геофизическая. 1954. № 5. С. 390-410.
Гзовский М.В. Основы тектонофизики. М.: Наука, 1975. 375 с.
Гущенко О.И. Анализ ориентировок сколовых тектонических смещений и их тектонофизическая интерпретация при реконструкции палеонапряжений // Доклады АН СССР. 1973. Т. 210, № 2. С. 331-334.
Гущенко О.И. Метод кинематического анализа структур разрушения при реконструкции полей тектонических напряжений // Поля напряжений в литосфере. М.: Наука, 1979. С. 7-25.
Гущенко О.И., Сим Л.А. Обоснование метода реконструкции напряженного состояния земной коры по ориентировкам сдвиговых тектонических смещений (по геологическим и сейсмологическим данным) // Механика литосферы: Тез. докл. Всесоюз. науч.-техн. совещ. Л.-М., 1974. С. 5-8.
Данилович В.Н. Метод поясов в исследовании трещиноватости, связанной с разрывными смещениями. Иркутск: Иркут. политехи. ин-т, 1961. 47 с.
Жиров Д.В., Сим Л.А. Мультикинематические дизъюнктивы как регуляторы транспозициональных полей напряжений на примере массива пород Ковдорского апатит-магнетитового месторождения // Триггерные эффекты в геосистемах: Тез. докл. се-минара-сов. (Москва, июнь 2010). М.: ИДГ РАН, 2010. С. 3132.
Карта новейшей тектоники СССР и сопредельных областей м-ба 1:1000000 / Ред. Н.И. Николаев. М.: ГУГК СМ, 1985.
Кольская сверхглубокая. М.: Недра, 1984. 490 с.
Кошечкин Б.И. Голоценовая тектоника восточной части Балтийского щита. Л.: Недра, 1979. 152 с.
Митрофанов Ф.П., Баянова Т.Б., Балабонин Н.Л., Сорохтин Н.О., Пожиленко В.И. Кольский глубинный раннедокембрийский коллизион: новые данные по геологии, геохронологии, геодинамике и металлогении // Вестник СПбГУ. 1997. Сер. 7, вып. 3. С. 5-18.
Мультимедийный справочник по минерально-сырьевым ресурсам и горно-промышленному комплексу Мурманской области: Цифровой информационный ресурс / Под ред. Ф.П. Митрофанова. Апатиты: ГИ КНЦ РАН, 2001. Ч. 1. (СБ-ШМ).
Николаев П.Н. Методика тектонодинамического анализа. М.: Недра, 1992. 295 с.
Осокина Д.Н. Об иерархических свойствах тектонического поля напряжений и деформаций в земной коре // Поля напряжений и деформаций в земной коре. / Под ред. А.С. Григорьева, Д.Н. Осокиной. М.: Наука, 1987. С. 136-151.
Осокина Д.Н. Парагенезы напряжений и вторичных структур в зонах активных разломов: Математическое моделирование. Структурные парагенезы и их ансамбли / Под ред. А.В. Лукьянова. М.: ГЕОС, 1997. С. 129-131.
Ребецкий Ю.Л., Сим Л.А., Козырев А.А., Рыбин В.В., Жиров Д.В. Первые оценки величин напряжений по геологическим данным // Современное состояние наук о Земле: М-лы междунар. конф., посв. памяти В.Е. Хаина. М.: Геол. факультет МГУ, 2011, С. 1553-1554. http://khain2011.web.ru.
Ронов А.Б., Хаин В.Е., Балуховский А.Н. Атлас литологопалеогеографических карт мира. Мезозой и кайнозой континентов и океанов / Отв. ред. В.Л. Барсуков, Н.П. Лаверов. Л.: ВСЕГЕИ, 1989. 79 с.
Рухлов А.С., Иванников В.В., Белл К. Рекуррентный дайково-диатремовый магматизм Кандалакшского грабена (Кольская щелочная провинция) // Рифтогенез, магматизм, металлогения докембрия. Корреляция геологических комплексов Фенноскан-дии. Тез. междунар. конф. Петрозаводск, 1999. С. 124-125.
Светов А.П., Свириденко Л.П. Магматизм шовных зон Балтийского щита. Л.: Наука. Ленинградское отд., 1991. 197 с.
Сим Л.А. Выявление новейших активизированных разломов в условиях слабообнаженных районов (на примере зоны сочленения Мезенской синеклизы и Среднего Тимана): Автореф. дис. ... канд. геол.-мин. наук. М.: МГУ, 1980. 24 с.
Сим Л.А. Некоторые особенности полей напряжений в зонах разломов (по геологическим и сейсмологическим данным) // Поля
напряжений и деформаций в земной коре. М.: Наука, 1987. С. 151-158.
Сим Л.А. Изучение тектонических напряжений по геологическим индикаторам (методы, результаты, рекомендации) // Изв. ВУЗов. Геология и разведка. 1991. № 10. С. 3-22.
Сим Л.А. Некоторые особенности взаимосвязи неотектонических движений Восточно-Европейской платформы с глобальным тектогенезом на основе изучения неотектонических напряжений) // Неотектоника и современная геодинамика контин. и океанов: Тез. докл. ХХ1Х тект. сов. М., 1996. С. 132-134.
Сим Л.А. Влияние глобального тектогенеза на новейшее напряженное состояние платформ Европы // М.В. Гзовский и развитие тектонофизики. М.: Наука, 2000. С. 326-350.
Сим Л.А. О возможной причине эрозионного сноса за мезозой -кайнозой и в начале новейшего этапа на Фенноскадинавском щите // Свойства, структура, динамика и минерагения литосферы Восточно-Европейской платформы: Материалы XVI международной конференции, 20-24 сентября 2010 г. Воронеж: Научная книга, 2010. С. 216-221.
Сим Л.А., Наумкин А.Н., Никулин В.И., Симонов Д.А. Тектонофизи-ческое изучение разломов Карелии в связи с выбором площадки для строительства Карельской АЭС // Тектонофиз. аспекты разломообразования в литосфере: Тез. Всес. сов. «Разломооб-разование в литосфере: тектонофиз. аспекты». Иркутск, 1990. С. 63.
Сим Л.А., Попова З.В., Шенкман Е.Я., Гущенко А.О., Серов А.В. Поля напряжений в зонах разломов по геологическим и экспериментальным данным // Эксперим. тект. в решении задач тео-рет. и практ. геологии: Тез. Второго Всес. симп. в Ялте. Киев, 1987. С. 150-151.
Сим Л.А., Свириденко Л.П., Брянцева Г.В. Об унаследованном развитии неотектонических разломов восточной части Балтийского щита // Разломообразование и сейсмичность в литосфере:
тектонофизические концепции и следствия: М-лы Всерос. сов. (18-21 августа 2009 г.). Т. 1. Иркутск: ИЗК СО РАН, 2009. С. 100-102.
Строение и динамика литосферы Восточной Европы. Результаты исследований по программе EUROPROBE. - М.: Геокарт: Геос, 2006. 736 с.
Стром А.Л., Никонов А.А. Распределение смещений вдоль сейсмо-генных разрывов и учет неравномерности подвижек при палео-сейсмологических исследованиях // Вулканология и сейсмология. 1999. № 6. С. 47-59
Шерман С.И. О потенциальной способности глубинных разломов к магмоконтролирующей деятельности // Вестник научной информации Забайкальского отд. Геогр. о-ва СССР. Чита, 1966. № 5. С. 16-24.
Garetsky R.G., R.E. Aizberg R.E., Karabanov A.K., Palienko V.P., Sli-aupa A.I. Neotectonics and Neogeodynamics of Central Europe // Geotectonics. V. 33. № 5. P. 343-352.
Rebetsky Yu.L. Possible Mechanism of Horizontal Compression Stress Generation in the Earth’s Crust // Doklady Earth Sciences. 2008. V. 423. № 2. P. 1448-1451. doi:10.1134/S1028334X08090274.
Sim L.A., Zhirov D.V., Smaglichenko T.A., Smaglichenko A.V. Results of the neotectonic stress state study in the eastern part of Baltic shield // Geoph. Research Abstr. V. 12, EGU2010-3599, 2010, EGU General Assembly 2010.
Zoback M.L., Zoback M.D., Adams J., Assumpgao M., Bell S., Bergman E.A., Blumling P., Brereton N.R., Denham D., Ding J., Fuchs K., Gay N., Gregersen S., Gupta H.K., Gvishiani A., Jacob K., Klein R., Knoll P., Magee M., Mercier J.L., Muller B.C., Paquin C., Rajenran K., Stephansson O., Suarez G., Suter M., Udias A., Xu Z.H. & Zhizhin M. Global patterns of tectonic stress // Nature. 1989. V. 341. № 6240. P. 291-298. doi:10.1038/341291a0.
Сим Лидия Андреевна, докт. геол.-мин. наук, в.н.с.
Институт физики Земли им. О.Ю. Шмидта РАН
123995, Москва, ул. Большая Грузинская, 10, стр. 1, Россия
И e-mail: [email protected]
Sim, Lidia A., Doctor of Geology and Mineralogy, Lead Researcher The Schmidt Institute of Physics of the Earth RAS 123995, Moscow, Bolshaya Gruzinskaya street, 10-1, Russia И e-mail: [email protected]
Жиров Дмитрий Вадимович, зав. отделом Геологический институт Кольского научного центра РАН 184209, Мурманская обл., Апатиты, ул. Ферсмана, д. 14, Россия e-mail: [email protected]
Zhirov, Dmitiy V., Head of Section Geological Institute The Kola Science Centre RAS 184209, Apatity, Fersman street, 14, Russia e-mail: [email protected]
Маринин Антон Витальевич, канд. геол.-мин. наук., сн.с..
Институт физики Земли им. О.Ю. Шмидта РАН 123995, Москва, Большая Грузинская ул., 10, стр. 1, Россия e-mail: [email protected]
Marinin, Anton V., Candidate of Geology and Mineralogy, Senior Researcher The Schmidt Institute of Physics of the Earth RAS 123995, Moscow, Bolshaya Gruzinskaya street, 10-1, Russia e-mail: [email protected]
