CC BY
63
ПЕРСПЕКТИВЫ РАЗВИТИЯ ГЕОЛОГИЧЕСКИХ ИССЛЕДОВАНИЙ
УДК 550.8: 553.2: 552:535.33
КОМПЛЕКСНЫЕ ПРОГНОЗНО-ПОИСКОВЫЕ КРИТЕРИИ ГЛАВНОЙ МЕТАЛЛОГЕНИЧЕСКОЙ СПЕЦИАЛИЗАЦИИ (МАЛОСУЛЬФИДНОЙ Pt-Pd ИЛИ СУЛЬФИДНОЙ ^-N0 РАННЕПРОТЕРОЗОЙСКИХ БАЗИТ-ГИПЕРБАЗИТОВЫХ РАССЛОЕННЫХ МАССИВОВ КРИСТАЛЛИЧЕСКИХ ЩИТОВ
академик Ф.П. Митрофанов, Д.В. Жиров, Т.Б. Баянова
Г еологический институт КНЦ РАН
Аннотация
В результате многолетних фундаментальных и прикладных исследований разработаны комплексы индикаторов главной металлогенической специализации (малосульфидной Pt-Pd или сульфидной ^-М) раннепротерозойских базит-гипербазитовых расслоенных массивов кристаллических щитов, и на их основе сформулирована прогнозно-поисковая методика, позволяющая оперативно диагностировать главную металлогеническую специализацию расслоенных массивов кристаллических щитов на самых ранних стадиях геологических работ. Прогнозно-поисковая методика будет полезна для обоснования металлогенического прогноза, ранжирования очередности направлений и объектов поисков, а также для планирования мероприятий по воспроизводству минерально-сыревой базы платиноидов и цветных металлов.
Ключевые слова:
металлогения, докембрийский, базит-гипербазитовый, изотопный, петрология, геодинамика, стратегические полезные ископаемые, прогнозно-поисковые критерии.
Введение
Большинство достоверных запасов и ресурсов стратегических видов полезных ископаемых (СПИ), к которым относятся и такие "бренды" Кольского полуострова, как металлы платиновой группы (PGM), никель (№), медь (Си) и кобальт (Со) приурочены к основным и ультраосновным массивам [1-4]. Здесь необходимо особо подчеркнуть, что уровень добычи во всём мире определяют в первую очередь крупные и уникальные месторождения. При этом для преобладающего большинства эксплуатирующихся отечественных сырьевых объектов к началу - середине первого десятилетия XXI в. остро встала проблема истощения наиболее богатой и рентабельной
приповерхностной части месторождений. По оценкам ведущих экспертов [2-3, 5], запасов
разрабатывающихся в настоящее время уникальных месторождений, обеспечивающих основные потребности отечественной и всей мировой экономики, хватит всего на 10-30 лет. Проблема усугубляется тем, что промышленные возможности, а также степень хозяйственной емкости - нагрузки на окружающую среду традиционных горнорудных регионов близки к предельным и не могут обеспечить существенный рост производства на этих территориях. Кроме этого, в России негативно сказываются и последствия реформирования экономики и государственной промышленной политики 1990-2000 гг., вследствие которых появился и с каждым годом обостряется дефицит поисковых заделов. В этих условиях особую актуальность приобретают фундаментальные исследования в области разработки новых прогнознопоисковых методов и технологий, позволяющих уверенно диагностировать металлогению и обосновывать выбор направлений и объектов поиска крупных / уникальных месторождений СПИ.
Теоретические и фактографические предпосылки металлогенических прогнозов СПИ на древних кристаллических щитах
Древние кристаллические щиты играют роль одного из главных источников разнообразных видов стратегических полезных ископаемых. В базит-гипербазитах добывается большая часть товарной продукции платиноидов, цветных металлов, золота, хрома, ванадия, титана, алмазов и
других типов СПИ. Каждый из геолого-промышленных типов месторождений характеризуется набором специфичных геолого-генетических черт и параметров: геодинамическими условиями; типом, составом и режимом магматизма; внутренней структурой тела и др. В расслоенных комплексах находятся сульфидные платинометалльные и медно-никелевые, а также оксидные хромитовые и титан-ванадиевые объекты, образующие как объемные залежи, так и стратиформные структуры типа "риф" [3-4, 6]. Причем часто в рамках одного крупного массива могут сосуществовать несколько геолого-промышленных типов месторождений с различной металлогенической специализацией. В других случаях петрология, структура и вещественный состав интрузива определяют распространение только одного типа оруденения для множества геологических тел, однако промышленно значимые объекты обнаруживаются далеко не в каждом из них.
В последние годы результаты наших исследований, а также работы ведущих мировых специалистов [4, 7-8] подтверждают важность понимания закономерностей распределения крупных месторождений разных металлов не только в пространстве, но и во времени, что требует использования в геохронологических и петрологических исследованиях новейшей и весьма дорогостоящей прецизионной лабораторной техники.
Вероятно, наиболее благоприятные условия образования и сохранения крупных ортомагматических сульфидных Cu-Ni (+МПГ) и малосульфидных Pt-Pd (+Ni, Cu) месторождений прямо связаны с особо высокими температурами в мантии, приводящими к выплавлению высокомагнезиальных магм, обогащенных рудными элементами мантийного ряда, и с большой мощностью и плавучестью субконтинентальной литосферной мантии (SCLM). Такие условия в эволюционирующей Земле моделируются для коматиитовых провинций преимущественно позднего архея и для провинций основных-ультраосновных пород суперконтинентов - континентальных литосферных плит со зрелой земной корой докембрия и реже - фанерозоя. Это предположение иллюстрируется концентрацией основных мировых ресурсов Pt-Pd руд в расслоенных интрузивах в постархейских - палеопротерозойских месторождениях с возрастами 2.7-2.5 и 2.0-1.9 млрд лет, а Ni руд - в позднеархейских коматиитах, в мезопротерозойских и в позднепалеозойских месторождениях. Приведенные эпохи совпадают с эпохами существования в разных регионах наиболее мощных (250-150 км) континентальных литосфер суперконтинентов, обусловленных завершением коллизионных процессов и последующих интенсивных плюмовых процессов с продолжительностью развития до 200 и более млн лет. По нашему уточненному анализу структуры, вмещающие Pt-Pd месторождения, должны быть типично внутриплитные (пример расслоенных интрузий Восточно-Скандинавской изверженной провинции), а вмещающие богатые сульфидные Ni-Cu месторождения - характеризоваться развитием на активных окраинах континентальных плит (примеры Печенги, Джиньчуаня, Норильска), где, вероятно, в субдукционных условиях происходило обогащение мантийного вещества коровым (в первую очередь коровой серой).
Проведенные нами комплексные исследования базит-гипербазитов Кольского региона и аналитическое сопоставление их результатов с известными данными по крупнейшим металлогеническим провинциям мира показали хорошо выраженную стереотипность геологопромышленных типов PGE и Cu-Ni месторождений и проявлений раннепротерозойских расслоенных комплексов древних щитов, что является крайне важным обстоятельством, позволяющим обосновывать прогнозно-поисковые модели в отношении новых интрузивов [9-10]. Аналогии и/или подобие геодинамических, геологических и структурных условий образования, а также изотопногеохимических, петрологических и минералогических характеристик позволяют сделать вполне весомое постановочное обоснование для задания параметров и направлений геологических исследований малоизученных объектов - раннепротерозойских расслоенных базит-гипербазитовых массивов, залегающих в пределах докембрийских щитов.
Одной из наиболее перспективных на обнаружение новых крупных PGE и Cu-Ni месторождений является Восточно-Скандинавская базитовая обширная изверженная провинция Фенноскандинавского (Балтийского) щита (ВСкБОИП), в пределах которой выявлено порядка 20 крупных и сотни мелких перспективных базит-гипербазитовых тел (рис. 1). Многие из этих раннепротерозойских расслоенных интрузий хорошо изучены геофизически, геохимически и разбурены по детальной сети наблюдений, кроме того, в них обнаружены практически все геологопромышленные типы рудных объектов из когда-либо описанных в других провинциях мира.
В - раннепротерозойские расслоенные комплексы с платинметалльной минерализацией; КП - Кольский пояс; ФКП - Фенно-Карельский пояс;
1 - Федорово-Панский; 2 - Мончеплутон; 3 - Мончетундровский, Волчетундровский массивы, Габбро Главного хребта; 4 - гора Генеральская;
5- Кандалакшский и Колвицкий массивы; 6-Луккулайсваара; 7- Ковдорский массив; 8-Толстик; 9-Ондомозерский; 10-Песочный; 11 - Пялочный; 12-Кейвица; 13-Портимо комплекс (Контиярви, Сиика-Кямя; 14-Пеникат; 15-Кеми; 16-Торнио; 17 - Коиллисмаа комплекс; 18-Аканваара
Рис. 1. Схема геологического строения ВСкБОИП с выделением рудоконтролирующих рифтов и основных раннепротерозойских расслоенных комплексов с платинометалльной минерализацией
Разработка прогнозно-поисковой методики на малосульфидные Pt-Pd или сульфидные ^-№ руды
За отправную точку наших исследований было взято представление о том, что результаты проявления стереотипных геологических процессов от зарождения магмы до ее внедрения и кристаллизации, в близких и схожих условиях не всегда конвергентны, что делает невозможным нахождение одного универсального признака рудоносности, и требует комплексности и структурированности действий при построении прогнозно-поисковых моделей. Согласно этому представлению массивы основного - ультраосновного состава, образованные территориально близко друг к другу и имеющие похожие черты строения и вещественного состава пород, могут нести различную рудную минерализацию или вообще оказаться безрудными.
К числу основных факторов, влияющих на конечный результат, относятся:
• благоприятная предыстория магматогенного источника - промежуточного слоя «кора-мантия», способствующая накоплению в нем сидерофильных (в т.ч. платиноидов) и халькофильных элементов;
• особенности геодинамического режима (активность плюма, рифтогенез и т.п.) и локализации магматического очага;
• последовательность, длительность и направленность кристаллизации минералов, обуславливающие концентрацию и форму нахождения полезных компонентов;
• геохронологический (абсолютный) возраст и длительность формирования пород, что определяется по соотношению ряда маркирующих изотопов (U-Pb, Rb-Sr, Sm-Nd и др.);
• изменчивость металлогении магматических и гидротермальных процессов, обусловленная как закономерными эндогенными процессами - ликвацией, фракционированием и обеднением/обогащением исходных компонентов в результате кристаллизации и/или контаминации магмы, так и эволюцией и цикличностью "открытия/закрытия" системы по отношению к мантийным источникам и мн. др.
Все эти особенности истории развития массива фиксируются весьма многочисленными параметрами и характеристиками: степенью дифференциации строения, вариациями химического и петрологического состава различных слоев и тел, их халькофильностью, значением R-фактора,
изотопными соотношениями и мантийными трассерами £Nd (T), TDM, Sr87/Sr86, Os187/Os188, He3/He4 и др., а также содержанием и соотношениями редких и рассеянных элементов и многими другими индикаторами.
Таким образом, были разработаны комплексы индикаторов главной металлогенической специализации (малосульфидной Pt-Pd или сульфидной Cu-Ni) раннепротерозойских базит-гипербазитовых расслоенных массивов кристаллических щитов, и на их основе сформулирована прогнозно-поисковая методика. Эта методика позволяет оперативно диагностировать малосульфидную платино-палладиевую или сульфидную медно-никелевую главную металлогеническую специализацию базитовых расслоенных массивов кристаллических щитов путем выполнения комплекса аналитических и лабораторных исследований (рис. 2). Соответствующие комплексы разноранговых критериев металлогенической специализации систематизированы и представлены в таблицах 1 и 2.
Рис. 2. Алгоритм выполнения аналитических и лабораторных исследований с целью определения главной металлогенической специализации базит-ультрабазитовых массивов
Таблица 1
Прогнозно-поисковые признаки и индикаторы условий образования малосульфидной Pt-Pd (с попутными №, Си, Аи, Со, ЯЬ) комплексной промышленной минерализации
Наименование признака/индикатора Параметры признака/индикатора
Г еофизический контроль Наличие по данным глубинной геофизики в основании коры гранулит-базитового (анортозитового) слоя с характеристиками "кора-мантия" ^р = 7.7-7.1 км/с), сформированного в результате плюмового андерплейтинга (состав слоя определяется по глубинно-коровым ксенолитам в трубках взрыва)
Структурный контроль Региональный: распространение на обширной площади архейских доменов фундамента дискордантного ансамбля рифтогенных осадочно-вулканогенных прогибов, даек и многофазных расслоенных базитовых массивов. Локальный: рудные тела локализуются в базальных (нижних) контактах, протяженных рифовых залежах, в участках пегматоидных базитов, в жильных и офсетных проявлениях
Г еодинамическая обстановка Масштабное, долговременное и пульсационное проявление глубинных плюмовых или астеносферно апвеллинговых процессов, обуславливающих формирование обширной изверженной базитовой внутриплитной континентальной провинции (ЫРЪ) несубдукционного типа. Смена геодинамического режима архейского орогенеза на внутриконтинентальный рифтогенез (с формированием разноориентированных энсиалических поясов). Формирование рудоконтролирующих базит-гипербазитовых массивов происходит на инициальной (предрифтовой) стадии континентального рифтогенеза
Вещественный контроль Кремнистая высоко Mg (бонинитоподобная) и анортозитовая магмы. Цикличность (закономерная полистадийность) строения раслоенных комплексов и резкая изменчивость кумулусной стратиграфии (ассоциаций) и геохимической специализации расплава. Большинство палеопротерозойских расслоенных интрузивов имеют в своем строении от 2 до 5 и более мегациклических подразделений - закономерно расслоенных серий от ультрабазитовых разновидностей к менее основным (габброидам). Приуроченность оруденения к наиболее контрастным сериям чередования тонких прослоев пород, разных по составу - лейко- и мезократовых габбро, норитов, анортозитов, плагиопироксенитов, разнозернистых и неоднородных текстур (например, varitextured gabbro), лейкократовых разностей (лейкогаббро, анортозитов, «пятнистых габбро»), разнозернистых, грубозернистых и пегматоидных пород с эруптивными магматическими взаимоотношениями. Все известные стратиформные месторождения рифового типа приурочены к границам мегациклов, которые при этом, как правило, отражают смену высокохромистой на малохромистую магму. Характерно интенсивное проявление в породных комплексах глубинных восстановительных флюидов, обогащенных соединениями С, F, СІ, Н и др. Минералогические факторы: связь PGM с рассеянной (убогой) сульфидной минерализацией, аномально высокая концентрация платиноидов в сульфидах. Коэффициент распределения платиновых металлов между ликвирующими силикатными и сульфидными расплавами достигает 100 тыс. и более
Изотопно-геохимические характеристики Г лубинный мантийный источник магм, изначально обогащённый рудными компонентами (фертильный источник) и литофильными элементами, что проявляется в таких изотопных метках, как (Т) - малые отрицательные значения от -1 до -3, Кг = Sr87/ Sr86 (0.702-0.705), Не3/Не4 (п х 10-5 - п х 10-6). Источник магмы и руд, отличный от источников срединноокеанических хребтов и зон субдукции
Г еохронологическая характеристика Внутриплитные базитовые обширные изверженные провинции с малосульфидными платино-палладиевыми месторождениями (Восточно-Скандинавская на Фенноскандинавском (Балтийском) щите, Восточно-Саянская - на выступе фундамента Сибирской платформы, Гуронская - на Канадском щите) формируются в самом начале эпох разрушения суперконтинентов чаще всего на геохронологической границе "архей -палеопротерозой" 2600-2400 млн лет назад. Для Восточно-Скандинавской провинции -это эпоха сумия - начала сариолия: 2530-2420 млн лет. Рудно-магматические комплексы развиваются длительно и пульсационно (фазы 2515±10 млн лет; 2490±10 млн лет; 2470±10 млн лет; 2450±10 млн лет) со сменой бонинитовых магм на анортозитовые, а их рудной специализации - от Сг и Си+№ к Pt+Pd и Ti+V
Известные промышленные месторождения находятся в регионально неметаморфизованных породах, а в регионально метаморфизованных расслоенных базитовых комплексах известны только непромышленные И-Рс! рудопроявления.
Имеются данные, что превышение РТ параметров условий среднетемпературной амфиболитовой фации ведет к обеднению оруденения
Таблица 2
Прогнозно-поисковые признаки и индикаторы условий образования сульфидной Си-№
(с попутными Со, S, PGM, Se, Те и др.) комплексной промышленной минерализации
Наименование признака/индикатора Параметры признака/индикатора
Геофизический контроль Наличие по данным геофизики локальных гравитационных аномалий с концентрацией их в узких линейных зонах. Подъем границы Мохоровичича с 40-42 км в обрамлении до 39-38 км в рудоконтролирующих структурах
Структурный контроль Региональный: узкие протяженные пояса в общем композитном ансамбле палеопротерозойских орогенов кристаллических щитов (например, Печенга-Имандра-Варзугская структура). Рудоносные интрузивные тела инъецированы в верхнюю часть разреза вулканогенно-осадочных толщ раннего палеопротерозоя. Локальный: руда находится в базальных контактах инрузивов, в жильных телах переотложенных (эпигенетических) руд, в том числе в офсетном залегании
Геодинамическая обстановка Рудогенные процессы и магматизм локализованы в пространстве и во времени в период смены геодинамического режима от внутриконтинентального рифтогенеза (энсиалического) на режим раннего спрединга красноморского (энсиматического) типа. Формирование рудоконтролирующих базит-гипербазитовых массивов происходит на завершающей стадии континентального рифтогенеза
Вещественный контроль Первичная магма деплетирована и по резкоземельному спектру близка к БСОХ типу. Продукты ферропикритовой магмы, обогащенной Fe и Ti, образуют единую вулкано-плутоническую серию пород. Для интрузивных рудоносных тел характерен габбро-верлитовый состав, субвулканические и гипабиссальные условия кристаллизации, широкая дифференциация пород с образованием сингенетического ряда: верлит-клинопироксенит-габбро-ортоклазовое габбро
Изотопно- геохимические характеристики Верхнемантийный источник декретированной магмы с изотопными метками: £Nd (Т) -малые положительные значения (от +0.5 до +4), ISr = Sr87/Sr86 (0.703-0.704), 0s187/0s188=0.935±0.03 (единичное определение)
Геохронологическая характеристика Базитовый магматизм спредингового типа в кристаллических щитах проявлен в конце внутриконтинентального рифтогенеза, завершая "Переходный период" (Transiton period) и начиная эпохи типичной тектоники литосферных плит (22001980 млн лет). В пределах Фенноскандинавского щита эта эпоха начала формирования Свекофенского палеоокеана
Метаморфизм Коллизионный метаморфизм приводит к образованию переотложенных (ремобилизованных) рудных тел, как внутри самих рудоносных интрузивов, так и в офсетных обстановках
Обсуждение результатов
Необходимо отметить, что повсеместно применяющаяся практика опоискования и разведки основана на экстенсивном подходе: массовом бурении по всему разрезу и протяжению массивов с целью выявления рудных пересечений с высокими концентрациями полезных компонентов. Нередко для больших интрузивов объемы только поискового бурения составляют сотни скважин суммарной протяженностью десятки километров [11]. В свою очередь такой подход влечет за собой интенсивное негативное воздействие на окружающую среду и нерациональное использование финансовых, трудовых и материальных ресурсов. В отличие от этой практики, предложенный нами метод для обоснования и поисков месторождений малосульфидных или сульфидных руд не требует на ранних стадиях геологоразведочных работ валового опоискования всех массивов и их участков с проходкой многочисленных дорогостоящих и несущих экологические риски горных выработок и скважин.
Разработанная система критериев формирует базис для новой прогнозно-поисковой методики, которая отличается от традиционных иерархичностью и комплексностью подходов. Она предназначена для того, чтобы диагностировать на самых ранних стадиях геологических работ малосульфидную Pt-Pd или сульфидную Cu-Ni специализации базит-гипербазитовых расслоенных массивов кристаллических щитов и производить ранжирование объектов по потенциальной рудоносности. Соответствие характеристик конкретного исследуемого объекта какому-либо из перечней признаков/индикаторов (табл. 1, 2) отражает сочетание (подчас исключительное) целого комплекса благоприятных геологических, геодинамических и тектонических условий, способствующих формированию промышленно значимых по масштабам и концентрациям рудных тел малосульфидного (Pt-Pd специализации) или сульфидного (Cu-Ni-Co специализация) типов. Вместе с тем несоответствие или существенное отклонение от вышеупомянутых индикаторов принципиально не отрицает возможность обнаружения в исследуемом геологическом теле повышенных концентраций элементов и минералов Pt, Pd, Cu, Ni, Co, однако обусловливает весьма малую вероятность нахождения в рассматриваемом массиве промышленно значимых месторождений малосульфидного Pt-Pd или сульфидного Cu-Ni-Co типов руд. При этом не исключается возможность локализации в пределах этого геологического тела месторождений и проявлений других стратегических полезных ископаемых (Ti-V, Cr и др. типов руд).
Таким образом, разработанный метод носит избирательный характер и применяется, прежде всего, для обоснования металлогенического прогноза, ранжирования очередности и оптимизации ресурсных затрат при геологическом изучении базит-ультрабазитовых массивов и их частей, может быть использован с целью планирования воспроизводства минерально-сыревой базы платиноидов и цветных металлов.
Авторы благодарят большой коллектив сотрудников Геологического института КНЦ РАН, ОАО Пана и других геологических организаций Кольского региона, внесших на разных этапах ценный вклад в разрабатываемую методику своей квалифицированной работой, советами и обсуждением всех сторон проблематики.
1. Крупные и суперкрупные месторождения: закономерности размещения и условия образования / под ред. Д.В. Рундквист. М.: ИГЕМ РАН, 2004. 430 с. 2. Богатство недр России. Минерально-сырьевой и стоимостный анализ. Пояснительная записка к геолого-экономическим картам. СПб.: Изд-во ВСЕГЕИ, 2007. 550 с. 3. Налдретт А. Дж. Магматические сульфидные месторождения медно-никелевых и платинометалльных руд: перевод с англ. В.А. Федоренко. СПб.: СпбгУ, 2003. 487 с. 4. Robb Introduction to Ore Forming Processes. Blackwell Publishing. 2004. 368 p. 5. Государственный доклад о состоянии и использовании минерально-сырьевых ресурсов в Российской Федерации, выпуски за 2000-2010 гг. / центр «Минерал» ФГУНПП «Аэрогеология». 2002-2011 гг. Режим доступа: http://www.mineral.ru/ 6. Alapieti T., Lahtinen J. Platinum-Group Element Mineralization in Layered Intrusions of Northern Finland and the Kola Peninsula, Russia // The Geology, geochemistry, Mineralogy and Mineral Benefication of Platinum-Group Elements / Cabri, L. (Ed.). Canadian Institute of Mining, Metallurgy and Petroleum, 2002. Spec. vol. 54. 507-546. 7. Groves D.I, Ojala J. Lithospheric scale controls on mineral deposits (with Precambrian examples, including Fennoscandian shield) / Presentation on Fennoscandian Exploration and Mining (FEM). 2009, Rovaniemi, Finland, December 1, 2009. 30 p. 8. Bayanova T., Ludden J., Mitrofanov F. Timing and duration of Palaeoproterozoic events producing ore-bearing layered intrusions of the Baltic Shield : metallogenic, petrological and geodynamic implications // Geological Society Special Publications, 2009. № 323. P. 165-198. 9. Митрофанов Ф.П. Поисковые индикаторы новых промышленных месторождений Rh-Pt-Pd, Co-Cu-Ni и Cr руд на Кольском полуострове // Отечественная геология. 2006. № 4. С. 3-9. 10. Митрофанов Ф.П. Корреляция состава и рудообразующих этапов в раннепротерозойских базит-гипербазитовых расслоенных интрузиях Финляндии и Кольского полуострова (Россия) / Стратегические минеральные ресурсы Лапландии - основа устойчивого развития Севера: сб. материалов проекта Интеррег-Тасис. Апатиты: Изд. КНЦ рАн, 2008. С. 13-18. 11. Rasilainen K. et al. Quantitative mineral resource assessment of platinum, palladium, gold, nickel, and copper in undiscovered PGE deposits in mafic-ultramafic layered intrusions in Finland / K. Rasilainen , P. Eilu, T. Halkoaho, M. Iljina, T. Karinen // Geological Survey of Finland, Report of Investigation 180, 2010. 338 p.
Сведения об авторах
Митрофанов Феликс Петрович - академик, главный научный сотрудник, консультант-советник РАН; e-mail: fe1ix@geoksc.apatity.ru
Жиров Дмитрий Вадимович - и.о. руководителя отдела инноваций; e-mail: zhirov@geoksc.apatity.ru Баянова Тамара Борисовна - д.г.-м.н., зав. лабораторией геохронологии и изотопной геохимии; e-mail: tamara@geoksc.apatity.ru
Работа выполнена при поддержке и финансировании грантов Р И № 11-05-12012 офи-м и госконтракта Минобрнауки России № 16.515.11.5013.
ЛИТЕРАТУРА
