CC BY
59
GEODYNAMICS & TECTONOPHYSICS
PUBLISHED BY THE INSTITUTE OF THE EARTH’S CRUST SIBERIAN BRANCH OF RUSSIAN ACADEMY OF SCIENCES
2014 VOLUME 5 ISSUE 4 PAGES 905-917
ISSN 2078-502X
http://dx.doi.org/10.5800/GT-2014-5-4-0162
Experimental study of decompression, permeability
AND HEALING OF SILICATE ROCKS IN FAULT ZONES
V. Ya. Medvedev1, L. А. Ivanova1, B. А. Lysov2, V. V. Ruzhich1, М. V. Marchuk1
1 Institute of the Earth's Crust of SB RAS, Irkutsk, Russia
2 Irkutsk State University, Irkutsk, Russia
Abstract: The article presents results of petrophysical laboratory experiments in studies of decompression phenomena associated with consequences of abrupt displacements in fault zones. Decompression was studied in cases of controlled pressure drop that caused sharp changes of porosity and permeability parameters, and impacts of such decompression were analyzed. Healing of fractured-porous medium by newly formed phases was studied. After experiments with decompression, healing of fractures and pores in silicate rock samples (3x2x2 cm, 500 °C, 100 MPa) took about 800-1000 hours, and strength of such rocks was restored to 0.6-0.7 of the original value. In nature, fracture healing is influenced by a variety of factors, such as size of discontinuities in rock masses, pressure and temperature conditions, pressure drop gradients, rock composition and saturation with fluid. Impacts of such factors are reviewed.
Key words: physical experiment, decompression, permeability, fracture healing, earthquake.
Citation: Medvedev V.Ya., Ivanova L.A., Lysov B.A., Ruzhich V.V., Marchuk M.V. 2014. Experimental study of decompression, permeability and healing of silicate rocks in fault zones. Geodynamics & Tectonophysics 5 (4), 905-917. doi:10.5800/GT-2014-5-4-0162.
Экспериментальное изучение декомпрессии, проницаемости
И ЗАЛЕЧИВАНИЯ СИЛИКАТНЫХ ПОРОД В ЗОНАХ РАЗЛОМОВ
1 12 1 1
В. Я. Медведев , Л. А. Иванова , Б. А. Лысов , В. В. Ружич , М. В. Марчук
1 Институт земной коры СО РАН, Иркутск, Россия
2 Иркутский государственный университет, Иркутск, Россия
Аннотация: Рассмотрены результаты петрофизических лабораторных экспериментов по изучению явлений декомпрессии, связанных с последствием скачкообразных смещений в зонах разломов. Измерялись и анализировались последствия декомпрессии при регулируемом спаде давлений, которая приводила к резкому изменению параметров пористости и проницаемости. Изучались процессы залечивания новообразованными фазами раскрывшегося тре-щинно-порового пространства. Показано, что после экспериментов с декомпрессией время заращивания трещин и
905
Tectonophysics
V.Ya. Medvedev et al.: Experimental study of decompression, permeability and healing of silicate rocks...
пор образцов силикатных пород размером 3x2x2 см при температуре 500 °С, давлении 100 МПа составляет порядка 800-1000 часов с восстановлением их прочности 0.6-0.7 по отношению к исходной. Рассматривается влияние основных факторов: размеров нарушений сплошности породных массивов, РТ-условий, скачков градиентов давления, состава пород и насыщенности флюидом на процессы залечивания в природных условиях.
Ключевые слова: физический эксперимент, декомпрессия, проницаемость, залечивание разрывов, землетрясение.
1. Введение
Физико-химические процессы, происходящие в глубинных сегментах зон разломов литосферы, где формируются рудные месторождения, а также очаги землетрясений, остаются недостаточно изученными до настоящего времени. Отсюда высоковостребованными можно считать экспериментальные исследования, с помощью которых воспроизводятся и оцениваются разнообразные параметры декомпрессии, проницаемости и залечивания горных пород с учетом тех термодинамических и трибохимических условий, которые предполагаются на разных глубинных уровнях и в различных геодинамических обстановках, сопровождающих деформации в зонах разломов. С позиций текто-нофизики и мезомеханики деструкция в литосфере является многоуровневым процессом, т.е. имеет место последовательное объединение мелких зародышевых трещин в более крупные, которые на следующем масштабном уровне также объединяются до стадий формирования внутриплитных и межплитных разломов в литосфере. Важно акцентировать внимание на следующих особенностях деструкции в литосфере. Морфология бортов магистральных разрывов при слиянии разномасштабных разрывов низшего иерархического уровня всегда характеризуется шероховатостью, то есть идеально гладких разрывных нарушений в геологической обстановке не бывает. Следовательно, режимы смещений берегов разрывов при наличии исходной шероховатости, как правило, являются прерывистыми, о чем свидетельствуют многолетние высокоточные измерения смещений в разломах [Ruzhich et al., 1999]. Этому также способствуют разнообразные режимы деформирования литосферы Земли под влиянием эндогенных и внеземных энергетических источников, которые имеют выраженную волновую природу. Строение зон разломов представлено сочетанием выступов с пятнами контактного скольжения и участков относительного растяжения. Таким образом, в глубинных сегментах зон разломов при неравномерности смещений шероховатых поверхностей и взаимодействии неровностей всегда возможны быстрые перепады давлений и явления декомпрессии, приводящие к скачкообразному снижению литостатического и флюидного давления и интенсивному восходящему тепломассопереносу. Именно такие бароградиентные геологические
системы моделировались лабораторными физическими экспериментами. В качестве объектов исследования авторами использовались две группы силикатных пород: амфибол-биотитовые гранодиориты (как равномерно-зернистые) и граниты-рапакиви (как неравномерно-зернистые породы).
2. Методы проведения эксперимента
И ИЗМЕРЕНИЙ
Эксперименты по декомпрессии силикатных пород проводились с целью петрохимического моделирования эндогенных процессов, происходящих в зонах разломов в условиях быстрых изменений скорости смещений. Использовалась аппаратура высокого давления, позволяющая работать в диапазоне температур 200-800 °С и давлений 20-200 МПа при одновременном воздействии флюида системы С-О-Н с контролируемым соотношением компонентов c амфибол-биоти-товыми гранодиоритами и гранитами рапакиви. Сброс давления после необходимой изобарно-изотермической выдержки осуществлялся через специальный клапан, время сброса давления от 100 МПа до атмосферного 2-1800 с. Моделировались различные режимы диссипации (от постепенного приведения системы к стандартным условиям до шоковой декомпрессии). До и после экспериментов на образцах определялась проницаемость; кроме того, исследовались прочностные характеристики пород (удельная прочность на сжатие), и просматривались шлифы.
При изучении восстановления сплошности образцов путем залечивания их в гидротермальных условиях применялись стандартные автоклавы и экзоклавы с рабочим объемом 200 см3 по следующей схеме (рис. 1). Образцы, ранее подвергнутые шоковой декомпрессии в условиях: давление 100 МПа, скорость сброса 50 МПа/с, температура 500 °С, помещались в автоклав, где и происходило их залечивание в течение необходимого времени.
Изучение петрофизических свойств горных пород проводилось в лаборатории физики нефтяного и газового пласта (ФНП) Иркутского государственного университета на стандартной аппаратуре типа АКМ, а также на созданных в лаборатории приборах. Исследования были выполнены на образцах, изготовленных в
906
Geodynamics & Tectonophysics 2014 Volume 5 Issue 4 Pages 905-917
| Рис. 1. Схема рабочей установки. | Fig. 1. Installation scheme.
виде цилиндров диаметром d=(2.44-2.51)-10-2 м, высотой h=(1.0-2.5)-102 м и объемом V=(4.9-10.0)-106 м3. Всего проанализировано порядка ста образцов горных пород. Определялось шестнадцать параметров, представленных в таблицах 1, 2.
Емкостная характеристика порового пространства пород изучалась по газу (N2) с помощью коэффициента открытой пористости, который определялся по формуле:
1-^ .100%,
Рт.ф./
где Коп - коэффициент открытой пористости (%); робр -объемная плотность породы; рт.ф. - плотность твердой фазы (минеральной части) образца.
Определение объемной плотности образцов проводилось тремя независимыми методами: гидростатическим взвешиванием; денсимером; по массе и объему, найденному простым геометрическим обмером образцов. Плотность твердой фазы образцов находилась по массе и объему минеральной части, который определялся на приборе КОФСП (комплексный определитель физических свойств породы), созданном в лаборатории ФНП. Принцип действия прибора основан на использовании закона Менделеева-Клапейрона, рабочим агентом служит газ. Следует отметить, что использо-
вание газа является единственным надежным способом определения малых значений емкости порового пространства представленных образцов горных пород. Точность прибора определялась эталонированием.
Динамические (фильтрационные) характеристики порового пространства образцов исследовались с помощью коэффициента проницаемости (Кпр), определяемого в основном на приборе КОФСП и в единичных случаях - на УИПК (универсальный исследователь проницаемости кернов) по газу (N2). Максимальное рабочее давление равнялось 0.6 МПа, давление обжима образцов - 1.11 МПа при температуре 293 K, т.е. условия фильтрации были близки к стандартным (лабораторным).
Для нахождения Кпр использовалась двучленная формула:
I = aQ + bQ2,
где I - градиент давления; a, b - фильтрационные коэффициенты; Q - расход газа.
Данный способ определения Кпр отличается от общепринятого, основанного на линейном сопоставлении среднего градиента давлений и среднего расхода газа, проходящего через образец. Коэффициент проницаемости рассчитывался по фильтрационному коэффициенту а, определенному графическим или аналитическим способом в комбинированных координатах:
(Р!+Р2+2Рб) (Р1-Р2) =
где Р1 - давление на входе в образец; Р2 - давление на выходе из образца; Рб - барометрическое давление. Полученные таким образом значения фильтрационного коэффициента а меньше определяемых по стандартной методике, поэтому коэффициент проницаемости, найденный следующим образом:
7Т 2Рб^-^
кпр = ~^г,
где ц - вязкость газа (МПа-с); h - длина образца; F -площадь фильтрации, оказывается выше рассчитанного традиционным путем.
Фильтрационный коэффициент b характеризует структуру порового пространства пород и степень турбулентности потока газа.
Точность определения проницаемости зависит от герметичности изоляции образца и чувствительности приборов, определяющих перепад давления на концах образца и расход газа. Герметичность в экспериментах проверялась по фильтрации газа при максимально возможных перепадах давления через непроницаемый образец-эталон. Минимально возможное определяемое значение проницаемости - 1-10-3 фм2. В коллекции представленных образцов наименьшее значение коэф-
907
908
Таблица 1. Физические свойства гранодиоритов Table 1. Physical properties of granodiorites
N обр h, CM d, CM F, cm2 V06p.> CM3 m, г Робр.) г/см3 Vx*, CM Рг.ф.5 г/см3 Vnop.5 см3 Кол., % a b 4 о t—1 « 2 mkm n-106 S-104 CM2 SyV104 cm-1 SyflT04 CM-1 bt, MKM 1, CM P-102 cm-1 S-103 CM2 syVio3 cm-1 44°3 cm-1
11 1.71 2.48 4.828 8.23 22.59 2.74 7.85 2.88 0.379 4.6 <1
12 1.75 8.45 23.05 2.73 8.15 2.86 0.397 4.7 - - <1 - - - - - - - - - - -
13 1.76 8.47 23.11 2.73 8.10 2.85 0.373 4.4 - - <1 - - - - - - - - - - -
14 1.77 8.55 23.44 2.74 8.15 2.88 0.393 4.6 - - <1 - - - - - - - - - - -
15 1.76 8.50 23.27 2.74 8.10 2.87 0.400 4.7 - - <1 - - - - - - - - - - -
16 1.76 - " - - " - 8.50 23.26 2.74 8.12 2.87 0.374 4.4 - - <1 - - - - - - - - - - -
1 1.65 2.48 4.828 7.97 21.89 2.75 7.71 2.84 0.266 3.2 1050 8730 10.0 0.050 408.0 102.00 1.280 383.50 _ _ _ _ _ _
2 1.68 8.11 22.22 2.74 7.75 2.87 0.360 4.4 1570 8800 7.1 0.040 1081.0 19.87 2.450 55.19 - - - - - -
3 1.71 2.49 4.870 8.32 22.71 2.73 7.92 2.87 0.400 4.8 2130 6202 5.3 0.030 1699.0 26.88 3.230 67.20 - - - - - -
4 1.62 2.48 4.828 7.82 21.38 2.73 7.50 2.85 0.321 4.1 1660 8800 6.8 0.040 1008.0 17.60 2.250 54.83 - - - - - -
5 1.71 8.26 22.54 2.73 7.84 2.87 0.416 5.0 830 8700 14.0 0.050 721.0 17.43 2.110 41.90 - - - - - -
6 1.56 - " - - " - 7.53 20.63 2.74 7.20 2.86 0.332 4.4 2280 8900 4.6 0.030 1666.0 22.89 3.040 68.95 - - - - - -
41 1.44 2.47 4.790 6.84 19.23 2.81 6.70 2.87 0.192 2.8 86 27 110 0.177 28.47 2.16 0.316 11.26 0.217 6176 12.90 17.64 2.58 91.90
42 1.39 4.770 6.59 18.51 2.80 6.47 2.86 0.119 1.8 109 47 86 0.196 14.92 1.19 0.184 9.97 0.239 3586 7.52 9.91 14.92 83.30
43 1.40 4.790 6.61 18.55 2.81 6.45 2.88 0.159 2.4 120 47 79 0.162 29.12 1.93 0.295 12.14 0.199 5783 12.08 15.96 29.12 100.40
44 1.44 6.86 19.28 2.80 6.70 2.87 0.165 2.4 100 13 97 0.180 23.59 1.83 0.267 11.10 0.220 5219 10.90 14.95 23.59 90.60
45 1.47 4.770 6.94 19.49 2.81 6.72 2.90 0.215 3.1 109 35 91 0.153 42.17 2.81 0.405 13.08 0.188 7879 16.52 22.93 42.17 106.60
46 1.32 2.46 4.750 6.27 17.61 2.81 6.00 2.94 0.276 4.4 77 34 120 0.148 63.98 3.74 0.596 13.53 0.181 11553 24.32 30.50 63.98 110.50
31 1.55 2.47 4.770 7.37 20.78 2.82 7.99 2.85 0.081 1.1 116 48 90 0.256 5.35 0.63 0.086 7.83 0.313 1675 3.51 5.18 0.70 63.90
32 1.45 2.47 4.790 6.94 19.53 2.81 6.80 2.87 0.146 2.1 107 143 91 0.186 19.33 1.56 0.225 10.70 0.228 4411 9.51 12.78 1.84 87.60
33 1.28 6.11 17.21 2.82 6.00 2.87 0.104 1.7 95 23 91 0.207 12.64 1.00 0.164 9.64 0.253 3212 6.71 8.20 1.34 78.80
34 1.25 5.96 16.73 2.81 5.83 2.87 0.125 2.1 53 50 160 0.247 10.96 1.01 0.170 8.11 0.302 3326 6.95 8.28 1.39 66.20
35 1.40 6.64 18.60 2.80 6.50 2.87 0.232 3.5 56 59 170 0.197 28.72 2.36 0.355 10.17 0.241 6942 14.50 19.25 2.90 83.00
36 1.31 4.770 6.19 17.34 2.80 6.00 2.87 0.217 3.5 109 44 81 0.136 60.27 3.18 0.514 14.67 0.167 10018 21.00 26.00 4.20 119.80
Примечание. 11-16 - исходные образцы гранодиоритов; 1-6 - контрольный эксперимент (снятие давления за 4 часа); 41-46 - сброс давления за 3.25 с; 31-36 -сброс давления за 2.5 с.
Note: 11-16 - original granodiorite samples; 1-6 - test experiment (pressure relief within four hours); 41-46 - pressure relief within 3.25 sec; 31-36 - pressure relief within 2.5
sec.
V.Ya. Medvedev et al.: Experimental study of decompression, permeability and healing of silicate rocks..
909
Таблица 2. Физические свойства гранитов-рапакиви Table 2. Physical Properties of Rapakivi Granites
N обр h, см d, см F, см2 Vo6p.> см3 ш, г Робр.) г/см3 СМ Рт.ф .> г/см3 VnOp.5 см3 Кол., % а Ъ кпр40-3, мД мкм п-106 S-104 см2 Syfl404 см-1 Ч40* см-1 К мкм 1, СМ Р-102 см-1 S-103 см2 Syflj'lO3 см-1 ч4°3 см-1
Г 1 1.61 2.48 4.828 7.77 20.46 2.63 7.65 2.67 0.124 1.6 1340.0 10209.00 8.1 0.064 125.8 3.90 0.50 31.25
Г 2 1.53 7.39 19.08 2.58 7.20 2.64 0.163 2.2 2550.0 6971.00 4.0 0.038 481.7 8.60 1.16 52.63 - - - - - -
Г 3 1.64 7.92 20.84 2.63 7.80 2.67 0.119 1.5 450.0 7220.00 25.0 0.115 35.8 2.10 0.26 17.40 - - - - - -
Г 4 1.46 7.05 18.76 2.66 7.00 2.68 0.071 1.0 4000.0 18182.00 2.4 0.022 663.7 6.41 0.91 90.91 - - - - - -
Г 5 1.67 8.06 20.86 2.59 7.90 2.64 0.161 2.0 1120.0 9333.00 10.0 0.032 637.0 10.10 1.25 62.50 - - - - - -
Г 6 1.80 8.69 22.66 2.61 8.45 2.68 0.243 2.8 113.0 14.80 110.0 0.089 113.5 5.47 0.63 22.47 - - - - - -
Г 21 1.67 2.49 4.870 8.13 20.81 2.56 7.70 2.70 0.428 5.3 6.8 3.5 1600 _ _ _ _ _ 0.602 4288 8.81 14.318 1.761 33.453
Г 22 1.69 2.48 4.828 8.16 21.15 2.59 8.06 2.62 0.100 1.2 21.5 3.8 520 - - - - - 0.721 804 1.66 2.716 0.333 27.158
Г 23 1.57 7.58 19.53 2.58 7.24 2.70 0.340 4.5 18.5 5.1 570 - - - - - 0.390 5575 11.54 17.498 2.308 51.465
Г 24 1.74 8.40 21.96 2.61 8.26 2.66 0.141 1.7 11.0 0.82 1060 - - - - - 0.865 949 1.97 3.301 0.393 23.416
Г 25 1.44 6.95 17.87 2.57 6.69 2.67 0.292 3.8 8.7 1.0 1100 - - - - - 0.589 3114 6.45 8.962 1.289 30.692
Г 26 1.37 6.61 17.29 2.61 6.60 2.62 0.014 0.2 8.2 0.5 1100 - - - - - 2.569 37.6 0.078 0.103 0.0156 7.350
Г 31 1.46 7.05 18.09 2.57 6.69 2.70 0.359 5.1 5.6 0.5 1700 _ _ _ _ _ 0.632 3895 8.06 11.37 1.613 31.671
Г 32 0.98 4.73 12.40 2.62 4.70 2.64 0.032 0.6 5.7 1.3 1100 - - - - - 1.483 195 0.40 0.383 0.809 11.958
Г 33 1.40 6.76 17.48 2.59 6.55 2.67 0.209 3.1 10.0 1.5 930 - - - - - 0.600 2496 5.17 6.985 1.033 33.422
Г 34 1.60 7.73 19.95 2.58 7.53 2.65 0.200 2.5 12.5 1.5 850 - - - - - 0.639 1890 3.91 6.051 0.783 30.255
Г 35 1.69 8.16 21.25 2.60 7.72 2.75 0.439 5.5 15.4 2.0 730 - - - - - 0.399 6656 13.78 22.491 2.756 61.232
Г 36 1.30 6.28 16.20 2.58 5.88 2.75 0.400 6.3 4.0 1.3 2200 - - - - - 0.647 4700 9.73 12.223 1.946 30.559
Г 11 1.66 8.02 20.57 2.57 7.80 2.64 0.216 2.7 7.8 0.59 1420.0 _ _ _ _ _ 0.790 1650 3.42 5.48 0.68 25.3
Г 12 1.49 7.19 18.56 2.58 6.90 2.69 0.294 4.1 8.0 2.60 1250.0 - - - - - 0.605 3272 6.78 9.76 1.36 33.1
Г 13 1.70 8.21 21.30 2.60 7.90 2.70 0.312 3.8 4.6 1.60 2500.0 - - - - - 0.888 2065 4.28 7.02 0.86 22.5
Г 14 1.80 8.69 22.25 2.56 8.45 2.63 0.243 2.8 5.4 0.70 2200.0 - - - - - 0.971 1392 2.88 5.01 0.58 20.6
Г 15 1.74 8.40 21.76 2.59 8.20 2.65 0.202 2.4 3.5 20.10 3300.0 - - - - - 1.285 902 1.87 1.57 0.18 7.8
Г 16 1.65 7.97 20.68 2.60 7.80 2.65 0.167 2.1 9.4 1.60 1200.0 — — — — — 0.828 1224 2.54 4.04 0.51 24.2
Примечание. Г-1-Г-6- исходные образцы гранитов-рапакиви; Г-21 - Г-26 - контрольный эксперимент (снятие давления за 4 часа); Г-31 - Г-36 - сброс давления за 13.5 с; Г-11 - Г-16 - сброс давления за 1.5 с.
N о t е: Г-1 - Г-6 - virgin rapakivi granite samples; Г-21 - Г-26 - test experiment (pressure relief within four hours); Г-31 - Г-36 - pressure relief within 13.5 sec; Г-11-Г-16-pressure relief within 1.5 sec. I
I
I
Geodynamics & Tectonophysics 2014 Volume 5 Issue 4 Pages 905-917
V.Ya. Medvedev et al.: Experimental study of decompression, permeability and healing of silicate rocks...
фициента проницаемости определено по образцу Г-4 и составило Кпр=2.4-10-3 фм2, т.е. в два раза превышает минимально возможное значение.
Колебание проницаемости исследованных образцов значительное: разброс значений достигает 3000 фм2. Для образцов с минимальной проницаемостью определялись средний радиус капилляров (пор) RK, число пор на единицу поверхности образца п, суммарная поверхность поровых каналов S, удельная поверхность на единицу объема образца и на единицу объема по-
рового пространства 5у^.
Средний радиус пор рассчитывался по формуле:
_ Isk^
где T - извилистость поровых каналов. Для нахождения параметра T необходимо проводить дополнительные исследования образца, которые пока не были осуществлены в силу технических трудностей, поэтому во всех расчетах T принимается равным единице.
Размеры пор малы и изменяются от 0.02 до
0.70 мкм. Число пор на единицу поверхности образца определялось по формуле:
п _
K(IJ
nRl
1012.
При минимальных значениях число пор достигало 6-109 п/м2, что говорит об однородной структуре порового пространства исследуемых горных пород, представленной огромной массой чрезвычайно малых пор.
Удельная поверхность в единице объема образца рассчитывалась через гидравлический радиус по формуле:
^ _ 2Уп
^ “ «к .
Так как поверхность электростатически заряжена, очевидно то, что исследуемая порода может стать хорошим катализатором поверхностных явлений на границах сред.
Удельная поверхность на единицу объема порового пространства определялась по формуле:
S _ —
ЭУД2 Rk .
Естественно, что она значительно больше Sw но имеет более плавный характер колебаний, что позволяет с большой достоверностью разбивать образцы на группы с близкими адгезионными свойствами.
Для образцов, у которых проницаемость резко увеличивалась, вместо пяти последних параметров определялись шесть дополнительных параметров, характеризующих трещиноватость. При незначительных изменениях Кпр рассчитывались все характеристики.
Определение раскрытия трещин bt и длины системы трещин l осуществлялось по закону Буссинеска:
Q _
Ь2 АР 12ц ' ДТ '
, Др
где bt - раскрытость трещин; — - градиент давления;
Q - расход газа; l - длина системы трещин.
Поверхностная плотность трещин (р) определялась по формуле:
где F - площадь фильтрации образца.
Для единой трещинной системы поверхностная плотность трещин совпадает с объемной, поэтому площадь поверхности трещин можно определить следующим образом:
S=2рV0бр,
где S - площадь поверхности всех трещин в образце; Уобр - объем образца.
Удельная поверхность единицы объема образца рассчитывалась по формуле:
-ТД1
_ 2р.
Удельная поверхность в единице объема порового пространства определялась по формуле:
где Упор - объем пор в образце.
Ошибки измерения всех физических свойств в абсолютном выражении в 2-4 раза меньше минимальных определенных значений соответствующих параметров, т.е. все определения являются значимыми.
3. Результаты исследований и их анализ
В более ранних экспериментах по декомпрессии силикатных расплавов было установлено, что при определенных режимах (в зависимости от скорости сброса давления) за счет вскипания расплава и резкого увеличения его объема (иногда почти десятикратно по отношению к первичному) происходит структурирование порового пространства формирующегося стекла [Letnikov et al., 1990]. При шоковой декомпрессии возникал столбик расплава, пронизанный системой тонких капилляров, вытянутых вдоль длинной оси (рис. 2). Таким путем формировалась высокоэффективная система флюидного массопереноса, реализация которой в подвергшейся декомпрессии магматической колонне весьма вероятна.
910
Geodynamics & Tectonophysics 2014 Volume 5 Issue 4 Pages 905-917
I Рис. 2. Проводящие структуры в гранитном расплаве после сброса давления. Т=850 °С, Р=2500 атм., скорость сброса давления 2500 атм/с. I
I Fig- 2. Conductive structures in granitic melt after pressure drop. T=850 °C; P=2500 atm; pressure relief rate = 2500 atm/sec.
С этой точки зрения интерес представляет проблема физического состояния горной породы при снятии литостатического и флюидного давления. Предполагалось, что декомпрессия образцов силикатных пород, помещаемых в условия высоких температур и давлений, соответствующих глубоким горизонтам земной коры, также приведет к определенным изменениям порового пространства [Medvedev et al., 1995; Balyshev et al., 1996].
При микроскопических исследованиях амфибол-биотитовых гранодиоритов, подвергнутых шоковой декомпрессии, сравнительно с исходными породами обнаруживаются следующие изменения. В биотите отчетливее проявляются трещины спайности, и повсеместно по этим трещинам выделяется мелкий рудный минерал черного цвета (магнетит или гематит). Чешуйки биотита становятся расщепленными, он как бы «разбухает», однако цвет его не меняется, а остается коричневым или светло-коричневым, как и в исходных породах. В плагиоклазах почти в каждом зерне появляются довольно широкие (до 0.02 мм) неминерализованные трещины. В крупных зернах наблюдается до 67 трещин, рассекающих зерно полностью или почти полностью. Нередко крупные трещины сопровожда-
ются более мелкими «оперяющими». Зерна амфибола разбиты грубыми субпараллельными трещинами, в которых заметна (хотя и не повсеместно) хлоритиза-ция, возникшая за время проведения эксперимента, что представляет несомненный интерес. Часто наблюдаются не открытые трещины, а зоны, где амфибол деформирован. Ширина таких зон порядка 0.012 мм. Многие крупные трещины не только рассекают отдельные зерна, но и проходят сквозь все минералы шлифа. Заметно увеличивается пористость породы за счет и трещин, и довольно крупных межзерновых пор овальной формы, причем размер последних достигает
0.8 мм.
В подвергнутых шоковой декомпрессии микроклиновых гранитах наблюдаются аналогичные явления: обилие мелких внутризерновых трещин, распространенных по всем слагающим породу минералам; развитие более крупных (до 0.12-0.15 мм) межзерновых трещин, переходящих в поры, и сквозных субпараллельных трещин, секущих все минералы шлифа.
Из изученных физических свойств наиболее закономерно и однозначно зависит от режима снятия давления проницаемость пород (рис. 3). При этом в случае максимально быстрого приведения системы к
911
V.Ya. Medvedev et al.: Experimental study of decompression, permeability and healing of silicate rocks...
I Рис. 3. Изменение проницаемости и открытой пористости в гранодиорите в зависимости от режима сброса давления.
IFig. 3. Changes of permeability and open porosity in granodio-rite samples depending on pressure release mode.
стандартным условиям проницаемость гранодиоритов возрастает от 0.00797 мД в контрольном эксперименте до 0.11383 мД в образцах, полученных при шоковой (снятие давления за 2 секунды) декомпрессии, т.е. более чем на порядок (см. табл. 1): для гранитов рапа-киви, соответственно, от 0.0233 мД и до 3.300 мД (табл. 2). Использованы средние значения по ряду экспериментов.
Средний радиус капилляров Rk (см. табл. 1) стабильно растет от 0.040 до 0.169 и 0.204 мкм с уменьшением времени декомпрессии и одновременным повышением проницаемости, при этом расчетное число пор резко уменьшается, что связано с методикой расчета, поскольку вклад в проницаемость капилляров большого сечения на порядки выше «слепых» или пор минимального размера.
У гранодиоритов коэффициент открытой пористости практически одинаков в исходных образцах и составляет 4.4-4.7 %. В контрольных экспериментах (см. табл. 1, образцы 1-6) коэффициент практически не меняется. С увеличением скорости декомпрессии до 500 атм/с резко увеличивается разброс значений коэффициента открытой пористости, уменьшаясь в образцах с максимальными значениями проницаемости. Вероятнее всего, это связано с тем обстоятельством, что при измерении проницаемости происходило значительное обжатие образца, что и обусловливало уменьшение измеряемого значения открытой пористости.
Отметим, что шоковая декомпрессия приводит к структурированию, выраженному не столько в количественном увеличении порового пространства, сколько в качественном изменении его конфигурации. Так, при малых скоростях декомпрессии проницаемость значительно увеличивается (см. табл. 1, образцы 1-6), но отсутствует система трещин по образцу. С увеличением скорости декомпрессии до 30 МПа/с начинает формироваться система пор или микро- и макротрещин, имеющая довольно отчетливую проводящую ориентировку (см. табл. 1, образцы 41-46). При дальнейшем уменьшении времени сброса давления и соответствующем увеличении скорости декомпрессии происходит закономерное увеличение проницаемости и величины раскрытия трещин bt до 0.3 мкм (см. табл. 1, образцы 31-36.)
Аналогичные явления наблюдаются для гранитов-рапакиви (табл. 2), хотя влияние декомпрессии на свойства порового пространства отмечается при меньших ее скоростях. Вероятнее всего, это происходит ввиду значительного различия в размерах кристаллов, слагающих образцы гранита-рапакиви.
Следует также отметить, что при шоковой декомпрессии и изменении структуры пространства закономерно уменьшается прочность образцов. Для гранодиорита эта зависимость имеет нелинейный характер (рис. 4), что подчеркивает значимость скорости декомпрессии.
В последующем было проведено залечивание образцов, подвергнутых шоковой декомпрессии со скоростью 50 МПа/с при начальном давлении 100 МПа и температуре 500 °С с целью определения возможности восстановления прочностных свойств породы в гидротермальных условиях.
В природе залечивание трещин может происходить в результате действия различных процессов. Основным фактором, определяющим закрытие микротрещин, является неоднородность поля напряжений, вызывающая растворение вещества в зонах концентрации напряжений, его перенос и отложение на бортах трещин. При подъеме гидротермального раствора в результате охлаждения растворимость кремнезема и других компонентов уменьшается, что приводит к выпадению осадка, осаждению его на стенках трещин и зарастанию или залечиванию трещин. Кроме того, упругие и другие физические свойства пород, возможности массопереноса связаны c геометрией порового пространства и общей пористостью и проницаемостью. Изученные процессы, вызывающие изменение морфологии трещин, могут сказываться и на величине этих характеристик. Экспериментально показано [Kotelnikova, Kotelnikov, 2008; Geguzin et al., 1988], что залечивание зависит от температуры, состава флюида и ориентации трещины. Этот процесс в силикатной среде может происходить путем диффузии через решетку и транспорта материала матрицы вдоль поверхности разграничения твердое - флюид или непосредственно через флюид и
912
Geodynamics & Tectonophysics 2014 Volume 5 Issue 4 Pages 905-917
| Рис. 4. Изменение прочности гранодиорита в зависимости от скорости сброса давления. | Fig. 4. Changes of granodiorite strength depending on pressure release rates.
не требует вовлечения нового матричного материала. При этом морфология трещины изменяется. Существует разность химического потенциала между передвигающимся концом трещины и ее стенками вследствие неровностей поверхности трещин:
Ap=y(Ce-Cm)Q,
где у - поверхностная энергия, Q - молекулярный объем, Ce, Cm - кривизна поверхности у перегибов и в средней части соответственно. Анализ уравнения показывает, что перенос вещества происходит от выпуклых или вогнутых участков с большой кривизной поверхности к выпуклым частям с низкой кривизной. Разная плотность дислокаций на стенках также является причиной разности химического потенциала [Geguzin et al., 1988]. В этой работе показано, что размер сужающейся поры в термоградиентном поле уменьшается до критического значения:
Y-D ~ (VT)-1/2.
В случае залечивания эта зависимость является проявлением пороговости по отношению к растворению боковых стенок трещины. Практически из изложенного следует, что чем больше неровностей в трещине, тем мельче образованные включения. Если плотность дислокаций низка (105-106 см2), они практически не влияют на перемещение ступени растворения. В случае сильнодислоцированного кристалла можно считать, что дислокаций так много, что они равномерно распределены вдоль поверхности. Ступени растворения (роста) тогда образуются выходами дислокаций, и перемещаться эти ступени могут, если расстояние
между дислокациями не меньше диаметра критического зародыша. Показано [Nichols, Mullins, 1965], что при залечивании трещин с образованием включений по отношению расстояния между включениями к радиусу A/R можно судить о преобладании формы массопереноса: при A/R=12.96 имеет место диффузия через кристаллическую решетку, если отношение равно 9.02 -диффузия идет через поровый раствор, при A/R=8.89 преобладает обмен поверхность-флюид. З.А. Котельникова [Kotelnikova, Kotelnikov, 2008], работая с призмами, выпиленными из монокристалла кварца и подвергнутыми термоудару для получения сети трещин, показала, что в случае залечивания флюидообразующей смесью разного состава средняя величина A/R для 500-600 °С составила 9.0, что отвечает массообмену через флюид.
Поскольку в нашем случае предполагается перенос вещества через флюид, заполняющий поры, для описания его можно воспользоваться видоизмененной схемой, предложенной Я.Е. Гегузиным с соавторами [Geguzin et al., 1988] для оценки скорости движения включений в термоградиентном поле (тот же тип массопереноса), и дополненной уравнениями [Smith, Evans, 1984] для определения времени образования первых включений при залечивании. Длину залеченного участка L можно определить:
L=A[(t exp(-Q/kT))/T],
где Л и и - константы, зависящие от начальной геометрии трещины; Q - энергия активации, для кварца по оценкам [Smith, Evans, 1984] равна 50 кДж/моль; T - абсолютная температура; k - константа Больцмана; t - время.
913
V.Ya. Medvedev et al.: Experimental study of decompression, permeability and healing of silicate rocks...
Тогда скорость перемещения замыкания трещины
(v):
v= A[(Q/kT)/T] = pw(Y-Y*)(dc/dT)AT,
где в - коэффициент, определяемый экспериментально по углу наклона на графике v-R. Анализ этих зависимостей приводит к выводу, что скорость перемещения замыкания трещины (фактически скорость залечивания) линейно зависит от ширины трещины. В работе [Nichols, Mullins, 1965] показано, что коэффициент в прямо зависит от плотности дислокаций в кристалле. Таким образом, залечивание в сильнодислоцированных кристаллах происходит быстрее. Скорость залечивания является функцией температуры, времени и концентрации раствора.
З.А. Котельниковой с А.Р. Котельниковым [Kotelni-kova, Kotelnikov, 2008] была проведена серия экспериментов с целью разработки способов образования синтетических флюидных включений путем залечивания трещин в минерале. В качестве затравки использовались кристаллы кварца (пьезокварц, р. Алдан). Шихта приготовлялась из оксидов Si и Al, карбоната Na и гидрокарбоната K. В систему добавлялся фтор в виде NaF и AlF3. Эксперименты проводились в гидротермальных условиях при давлении 2 кбар (200 МПа). Изначально система была нагрета до 675 °С и выдержана при одной температуре в течение 100 часов для уравновешивания расплава с материалом затравки. Затем проводилось ступенчатое (по 7-8 °C) снижение температуры, и при каждой новой температуре система выдерживалась в течение 24-30 часов. По достижении 650 °C температура выдерживалась еще в течение 100 часов, и затем производилась закалка опытов. В результате в образце из пересыщенной водой системы все трещины и поры были полностью залечены.
Анализ приведенных выше данных позволил предположить, что в нашем случае следует ожидать интенсивного залечивания образцов за время порядка 1000 часов при создании оптимальных условий массопереноса.
Эксперимент осуществлялся по методике изобарноизотермической выдержки в контролируемой атмосфере в термоградиентных условиях.
Образцы, предварительно подвергнутые шоковой декомпрессии, в виде прямоугольников размещались вдоль автоклава (см. рис. 1). Залечивание гранодиоритов производилось в области гидротермального роста кварца и слюды при температуре 500 °С, давлении 100 МПа в присутствии щелочного раствора. В качестве шихты как исходного материала для залечивания применялась смесь кварца, слюды и амфибола, которая помещалась на дно автоклава. Длительность эксперимента составляла 240 и 800 часов. После цикла залечивания на образцах отмечается рост новых фаз (преимущественно кварца) и замыкание в некоторых случаях устьев трещин. С увеличением времени гидро-
термальной обработки прочность образцов увеличивается до 60 % от исходной прочности гранодиорита. Ввиду того, что с увеличением длительности эксперимента происходит выполаживание кривой, дальнейшего увеличения прочности образцов не следовало ожидать и эксперименты длительностью более 800 часов не проводились.
4. Обсуждение результатов
Проведенные петрофизические эксперименты по декомпрессии показали, что в образцах, подвергнутых шоковой декомпрессии со скоростью до 50 МПа/с, наблюдается обилие мелких внутризерновых трещин, распространенных по всем слагающим породу минералам; развитие более крупных (до 0.12-0.15 мм) межзерновых трещин, переходящих в поры, и сквозных субпараллельных трещин, секущих все минералы образца. Иными словами, такая декомпрессия приводит к структурированию, выраженному не столько в количественном увеличении порового пространства, сколько в качественном изменении его конфигурации. Формируется система пор или микро- и макротрещин, имеющая довольно отчетливую проводящую ориентировку, что закономерно отражается на их физических свойствах. Резко меняются проницаемость, открытая пористость и другие параметры. Соответственно меняется и прочность образцов. Можно отметить, что эти эффекты нелинейно зависят от скорости сброса давления и, вероятно, при больших скоростях декомпрессии, существующих в природе, должны проявляться более значительно.
Как показывает опыт изучения условий формирования разнообразных рудных месторождений гидротермального типа, приуроченных к зонам разломов в интервале глубин 1-5 км, их становление происходило многостадийно и на протяжении длительного времени. [Kushnarev, 1982]. В природных условиях процесс залечивания сегментов разломов во многом сходен с условиями формирования рудных месторождений. Таким образом, небольшие трещины и поры при благоприятных условиях могут залечиваться в течение недель-месяцев в прерывистом режиме, тогда как крупные сегменты разломов, вскрывшиеся, например, при сильных землетрясениях, могут находиться в стадии залечивания в условиях зеленосланцевой фации метаморфизма в течение многих десятков, сотен или тысяч лет. Влияние скорости тектонического крипа на длительность залечивания разрывных нарушений в зонах разломов отмечалось в работе [Kocharyan, 2010]. Ее автором было показано, что высокие скорости крипа, порядка V > 25 мм/год, способны значительно снизить эффективность залечивания.
Таким образом, ответы на вопросы о длительности залечивания разломов всегда требуют уточнения различных параметров с учетом термо- и барометрических градиентов, из-за разнообразных сочетаний кото-
914
Geodynamics & Tectonophysics 2014 Volume 5 Issue 4 Pages 905-917
рых может быть дан только приблизительный ответ [Ruzhich et al., 1990, 2014; Kocharyan, 2009]. Важно учитывать, что при гидротермально-метасоматическом залечивании на глубинах порядка 5-10 км прочность пород на сжатие может восстанавливаться примерно до уровня 0.6-0.7 от исходной, что указывает на возможность соответствующего уменьшения сдвиговой прочности в тех сегментах разломов, где проявлялась декомпрессия.
Следует заметить, что в контактных пятнах скольжения на неровностях плоскостей разломов при высоких давлениях и темпах деформаций, а также при быстрых смещениях маловероятно проявление декомпрессионного залечивания трещинно-поровых образований. В подобных условиях реализуется иной трибохимический процесс. Имеются в виду процессы ускоренного преобразования вещества в условиях сильного сжатия и трения в плоскостях зеркал скольжения при высокоскоростных динамических подвижках. В таких случаях в глубинных сегментах сейсмоактивных разломов на поверхности зеркал скольжения при сверхтонком измельчении минералов до размеров, измеряемых первыми нанометрами, фиксируется наличие тончайших обводненных слоев из преобразованных минеральных фаз, обладающих пониженным сопротивлением сдвигу [Sobolev et al., 2011, 2014].
5. Заключение
На основе полученных результатов о шоковой декомпрессии, запускающей возникновение активных
транспортных путей в земной коре и проявленной в зонах внутрикоровых разломов и разломных пересечений, необходимо подчеркнуть,что подобные процессы могут быть кратковременными, но за период своего существования они способны за короткое время, например при сейсмической активизации, переместить значительные количества как инертного, так и химически активного вещества в разрывных нарушениях разного масштаба.
В рассмотренных выше экспериментах были созданы достаточно идеальные условия для быстрого, в течение нескольких недель, залечивания небольших по объемам образцов. Этот параметр может использоваться в качестве отправной точки при расчетах времени восстановления сплошности земной коры с учетом сложной комбинации различных факторов. Значимость каждого из перечисленных факторов еще предстоит изучить при проведении дальнейших экспериментов.
Интерпретируя полученные результаты, следует заметить, что в зонах разломов явления декомпрессионного залечивания возникают за пределами пятен контактного скольжения на неровностях. Однако именно в них в условиях сильнейшего деформирования и разогрева происходят сложные трибохимические процесс-сы, которые определяют уровень сдвигового сопротивления и режим концентрации сдвиговых напряжений. Процессы декомпрессионного залечивания в разломах земной коры, по-видимому, не являются определяющими при подготовке масштабного динамического разрушения высоконапряженных массивов горных пород.
6. Литература / References
Balyshev S.O., Medvedev V.Ya, Savel’eva V.B., 1996. On the problem of shock decompression of silicate rocks. Doklady AN 347 (1), 66-68 (in Russian) [Балышев С.О., Медведев В.Я., Савельева В.Б. К проблеме декомпрессии силикатных пород // Доклады АН. 1996. T. 347. № 1. С. 66-68].
Geguzin Ya.E., Kaganovsky Yu.S., Kruzhanov V.S., 1988. About the temperature relationship for NaCl crystals of the equilibrium shape. Kristallografiya 33 (6), 1495-1498 (in Russian) [Гегузин Я.Е., Кагановский Ю.С., Кружанов В.С. О температурной зависимости равновесной формы кристаллов NaCl // Кристаллография. 1988. Т. 33. Вып. 6. С. 14951498].
Kocharyan G.G., 2009. Mechanisms of strength recovery of seismogenic faults. In: Problems of interacting geospheres. GEOS, Moscow, p. 77-86 (in Russian) [Кочарян Г.Г. Механизмы восстановления прочности сейсмогенных разломов // Проблемы взамодействующих геосфер. М.: ГЕОС, 2009. C. 77-86].
Kocharyan G.G., 2010. Fault zone as a nonlinear mechanical system. Fizicheskaya Mezomekhanika 13 (Special Issue), 5-17 (in Russian) [Кочарян Г.Г. Разломная зона как нелинейная механическая система // Физическая мезомеханика. 2010. Т. 13. Специальный выпуск. С. 5-17].
Kotelnikova Z.A., Kotelnikov A.R., 2008. NaF-bearing fluids: Experimental investigation at 500-800 °C and P = 2000 bar using synthetic fluid inclusions in quartz. Geochemistry International 46 (1), 48-61. http://dx.doi.org/10.1134/S00167 02908010047.
Kushnarev I.P., 1982. Depths of Formation of Endogenous Ore Deposits. Nedra, Moscow, 166 p. (in Russian) [Кушнарев И.П. Глубины образования эндогенных рудных месторождений. М.: Недра, 1982. 166 с.].
Letnikov F.A., Kuznetsov K.E., Medvedev V.Ya., 1990. On the problem of shock decompression of fluidized silicate melts. Doklady AN SSSR 313 (3), 682-684 (in Russian) [Летников Ф.А., Кузнецов К.Е., Медведев В.Я. К проблеме шоковой декомпрессии флюидизированных силикатных расплавов // Доклады АН СССР. 1990. Т. 313. № 3. С. 682-684].
915
V.Ya. Medvedev et al.: Experimental study of decompression, permeability and healing of silicate rocks...
Medvedev V.Ya., Balyshev S.O., Savelieva V.B., 1995. Change of physical properties of rocks under shock decompression. In: The 13th Russian Meeting on Experimental Mineralogy. IEM RAS, Chernogolovka, p. 68 (in Russian) [Медведев В.Я., Балышев С.О., Савельева В.Б. Изменение физических свойств горных пород при шоковой декомпрессии // XIII Российское совещание по экспериментальной минералогии. Черноголовка: ИЭМ РАН, 1995. С. 68].
Nichols F.A., Mullins W.M., 1965. Morphological changes of a surface of revolution due to capillarity-induced surface diffusion. Journal of Applied Physics 36 (6), 1826-1835. http://dx.doi.org/10.1063/1.1714360.
Ruzhich V.V., Ivanova L.A., Medvedev V.Ya., 2014. On the study of tribological processes in slickensides from fault zones. In: Faulting in the lithosphere and related processes. Tectonophysical analysis. Abstracts of the All-Russia meeting with visiting researchers from other countries (11-16 August 2014, Irkutsk). IEC SB RAS, Irkutsk, p. 65 (in Russian) [Ружич В.В., Иванова Л.А., Медведев В.Я. Об изучении трибологических процессов в зеркалах скольжения из зон разломов // Разломообразование в литосфере и сопутствующие процессы: тектонофизический анализ: Тезисы докладов Всероссийского совещания с участием приглашенных исследователей из других стран (11-16 августа 2014 г., г. Иркутск). Иркутск: ИЗК СО РАН, 2014. С. 65].
Ruzhich V.V., Medvedev V.Ya., Ivanova L.A., 1990. Healing of seismogenic fractures and repeatability of earthquakes. In: Seismicity of the Baikal rift: prognostic aspects. Nauka, Novosibirsk, p. 44-50 (in Russian) [Ружич В.В., Медведев В.Я., Иванова Л.А. Залечивание сейсмогенных разрывов и повторяемость землетрясений // Сейсмичность Байкальского рифта: прогностические аспекты. Новосибирск: Наука, 1990. С. 44-50].
Ruzhich V.V., Truskov V.A., Chernykh E.N., Smekalin O.P., 1999. Recent movements in the fault zones of Pribaikalie and mechanisms of their initiation. Geologiya i Geofizika (Russian Geology and Geophysics) 40 (3), 360-372.
Smith D., Evans B. 1984. Diffusional crack healing in quartz. Journal of Geophysical Research 89 (В6), 4125-4135. http://dx.doi.org/10.1029/JB089iB06p04125.
Sobolev G.A., Genshaft Yu.S., Kireenkova S.M., Morozov Yu.A., Smul’skaya A.I., Vittegren’ V.I., Kulik V.B., 2011. Effects of high pressure and temperature on properties of nanocrystals in rocks: Evidences from Raman spectroscopy. Izvestiya, Physics of the Solid Earth 47 (6), 465-474. http://dx.doi.org/10.1134/S1069351311050053.
Sobolev G.A., Vettegren V.I., Ruzhich V.V., Ivanova L.A., Mamalimov R.I., Shcherbakov I.P., 2014. The study of nanocrystals in slickensides. In: The 15th International Conference "Physical, chemical and petrophysical studies in Earth Sciences". Moscow, September 29 - October 01, 2014. Conference Proceedings. IGEM RAS, Moscow, p. 211-213 (in Russian) [Соболев Г.А., Веттегрень В.И., Ружич В.В., Иванова Л.А., Мамалимов Р.И., Щербаков И.П. Исследование нанокристаллов в зеркале скольжения // Пятнадцатая международная конференция «Физико-химические и петрофизические исследования в науках о Земле». Москва, 29 сентября - 1 октября 2014 г. Материалы конференции. М.: ИГЕМ РАН, 2014. С. 211-213].
Медведев Владимир Яковлевич, канд. геол.-мин. наук, с.н.с.
Институт земной коры СО РАН
664033, Иркутск, ул. Лермонтова, 128, Россия
И e-mail: med@crust.irk.ru
Medvedev, Vladimir Ya., PhD in Geology and Mineralogy, Senior Researcher Institute of the Earth's Crust SB RAS 128 Lermontov street, Irkutsk 664033, Russia И e-mail: med@crust.irk.ru
Иванова Лариса Александровна, канд. геол.-мин. наук, с.н.с., Институт земной коры СО РАН 664033, Иркутск, ул. Лермонтова, 128, Россия e-mail: liva@crust.irk.ru
Ivanova Larisa A., PhD in Geology and Mineralogy, Senior Researcher Institute of the Earth’s Crust SB RAS 128 Lermontov street, Irkutsk 664033, Russia e-mail: liva@crust.irk.ru
916
Geodynamics & Tectonophysics 2014 Volume 5 Issue 4 Pages 905-917
Лысов Борис Антонович, канд. геол.-мин. наук, доцент кафедры геологии нефти и газа ИГУ Иркутский государственный университет, геологический факультет 664003, Иркутск, ул. Ленина, 3, Россия Тел.: 8(3952)243280
Lysov, Boris A., PhD in Geology and Mineralogy, associate professor Irkutsk State University, Geological Faculty 3 Lenin street, Irkutsk 664003, Russia Tel.: 8(3952)243280
Ружич Валерий Васильевич, докт. геол.-мин. наук, г.н.с.
Институт земной коры СО РАН
664033, Иркутск, ул. Лермонтова, 128, Россия
Тел. (3952)422776; e-mail: ruzhich@crust.irk.ru
Ruzhich, Valery V., Doctor of Geology and Mineralogy, Chief Researcher
Institute of the Earth's Crust SB RAS
128 Lermontov street, Irkutsk 664033, Russia
Tel. +7(3952)422776; e-mail: ruzhich@crust.irk.ru
Марчук Марина Валерьевна, канд. геол.-мин. наук, м.н.с.
Институт земной коры СО РАН
664033, Иркутск, ул. Лермонтова, 128, Россия
Marchuk, Marina V., PhD in Geology and Mineralogy, Junior Researcher
Institute of the Earth’s Crust SB RAS
128 Lermontov street, Irkutsk 664033, Russia
917
