Моделирование рудообразующих систем как основа для прогнозирования крупных месторождений стратегических металлов

А.В. Волков — д.г.-м.н., заведующий лабораторией Геологии рудных месторождений, главный научный сотрудник ИГЕМ РАН.





Введение

Россия располагает минерально-сырьевым комплексом (МСК) мирового класса, представляющим собой фундамент экономики страны, а также является одним из ведущих экспортёров минерального сырья и его первичной продукции. Развитие экономики, национальная безопасность и уровень жизни населения Российской Федерации в значительной степени зависят от мощной горнодобывающей промышленности и надежной минерально-сырьевой базы (МСБ). Однако проблема обеспеченности высокотехнологичной промышленности России минеральным сырьем с каждым годом становится все более острой в связи с неуклонным сокращением богатых и относительно доступных в техническом и экономическом отношении месторождений и неустойчивостью его импорта из стран-продуцентов вследствие политической нестабильности в мире.

В статье обсуждаются результаты металлогенических исследований последних лет. В ходе исследований создана база данных месторождений и перспективных рудопроявлений стратегических и высокотехнологичных металлов России; выделены главные рудообразующие системы — источники высокотехнологичных критических металлов; выполнены металлогенические исследования на основе пространственно-статистического ГИС-анализа; составлены варианты прогнозно-металлогенических карт размещения рудообразующих систем стратегических и высокотехнологичных металлов, разработаны новые подходы к прогнозированию крупных месторождений стратегических металлов.

Концепция рудообразующей (минеральной) системы

Под рудообразующей системой мы, вслед за Д.В. Рундквистом [1], понимаем комплекс взаимосвязанных процессов генерации рудоносных флюидов, их миграции к поверхности и разгрузки на геохимических барьерах, определяющих формирование аномальных скоплений полезных ископаемых. Определение рудообразующей (минеральной) системы по мировым источниками [1–5] аналогично концепции нефтяной системы, но из-за особенностей природы рудных месторождений и вмещающих пород минеральная система гораздо более разнообразна и сложна (рис. 1).

Рис. 1. Концепция рудообразующей (минеральной) системы, модифицировано из [4] и [5]

В таких системах принято выделять корневую зону (область магмо- и флюидозарождения), зону транспорта (тепломассопереноса) и зону концентрированного рудоотложения. Для изучения, как правило, доступна лишь зона рудоотложения, которая реализуется в виде ареала развития оруденения, гидротермально измененных пород и ореолов рассеяния рудного вещества, выделяемых в качестве рудных узлов.

Для формирования рудного месторождения требуется источник металлов, агенты транспортировки (флюиды или магмы) и область осаждения или накопления руд. Объединяющая все эти процессы рудообразующая система включает все геологические и геодинамические факторы на всех масштабных уровнях, которые контролируют возникновение, развитие и сохранение рудных месторождений.

Выделению рудообразующих систем должны сопутствовать: (1) локальные исследования на известных месторождениях, включая изучение их геологического строения, типов горных пород, расположение потенциальных рудовмещающих структур, физико-химические процессы рудообразования; (2) исследования регионального масштаба, включающие изучение геодинамического контроля размещения рудных месторождений, физико-химических процессов в масштабе крупного геологического блока, определяющих формирование рудных месторождений; эволюции магм и других источников энергии и флюидов в масштабе рудообразующей системы.

Главные рудообразующие системы — источники стратегических металлов на территории России

К стратегическим видам в мире в основном относят минеральное сырье, которое лежит в основе высоких технологий и имеет важнейшее значение для поступательного развития многих отраслей промышленности [6].

Рис. 2. Основные месторождения стратегических металлов РФ. Крупными значками показаны крупные месторождения, мелкими – мелкие и средние. Крупные месторождения стратегических высокотехнологичных металлов РФ. Крупными значками показаны крупные месторождения, мелкими — мелкие и средние.

В таблице 1 приведена краткая характеристика 8 главных рудообразующих систем высокотехнологичных критических металлов России: расслоенных интрузивов основных пород (магматическая базитовая и ультрабазитовая); скарново-порфировая и порфирово-эпитермальная; орогенная (в том числе связанная с интрузивами гранитоидов); связанная с щелочными гранитами; ультраосновных-щелочных пород и карбонатитов; колчеданно-вулканогенная и осадочно-вулканогенная; осадочных бассейнов; россыпей и кор выветривания. Расположение основных месторождений высокотехнологичных критических металлов на территории России показано на обзорной карте (рис. 2).

Для каждой рудообразующей системы характерны свои промышленные и минеральные типы месторождений, отличающиеся набором основных и важнейших попутных компонентов (в том числе и критических) и по технологическим свойствам руд.

База данных месторождений и перспективных рудопроявлений стратегических металлов на территории России

В ходе исследований подготовлена база данных месторождений и перспективных рудопроявлений стратегических металлов, состоящая из связанных между собой информационных блоков, характеризующих примеры минеральных месторождений, геологическую среду их размещения, интегральные модели, отражающие условия рудообразования и геодинамические обстановки развития рудообразующих систем.

На территории Российской Федерации выявлены более 18 тыс. рудных объектов разного ранга, от рудопроявления до крупных месторождений, среди которых насчитывается около 8 тыс. золоторудных объектов, медных — более тысячи, свинцово-цинковых — около 2 тыс. и редкометалльных — более тысячи. Кроме того, известны Pb-Zn месторождения и рудопроявления скарнового, связанного с карбонатитами, и эпитермального в терригенных толщах типов — менее значимые для промышленности. При этом среди свинцово-цинковых месторождений миссиссипский тип (MVT) насчитывает более 80, тип SEDEX — более ста, а колчеданный тип (VMS) — около двухсот объектов. На основе этой базы данных были составлены обзорные карты (рис. 2).

Месторождений, в рудах которых стратегические высокотехнологичные металлы являются основными и попутными компонентами, насчитывают менее тысячи двухсот, из них около двухсот — крупные (рис. 2). В этой выборке выделяется до 200 мелких и средних месторождений, руды которых отличаются повышенной комплексностью: в них в виде основных и попутных насчитывается более пяти компонентов, входящих в группу стратегических высокотехнологичных металлов.

Анализ базы данных показал, что наиболее комплексные руды характерны для следующих минеральных типов: редкоземельно-редкометалльный апогранитовый, щелочных метасоматитов, апатит-редкоземельно-редкометалльный, редкометалльный пегматитовый, касситерит-вольфрамитовый грейзеновый, медно-колчеданный в вулканогенно-осадочных, медно-цинково-колчеданный, колчеданно-полиметаллический в терригенных и вулканогенных породах и базитовый титаномагнетитильменит-ванадиевый. В химический состав руд месторождений перечисленных минеральных типов входят Ag, Be, Bi, Cd, Ce, Co, Cs, Cu, Ga, Ge, Hf, Hg, In, La, Li, Nb, Rb, REE, Sc, Se, Sn, Sr, Ta, Te, Th, Ti, Tl, U, V, W, Y, Zr.

Металлогенические исследования на основе пространственно-статистического ГИС анализа

Рудообразующие системы представляют собой сложные иерархические объекты, состоящие из более мелких подсистем, меняющихся при взаимодействии, а также под воздействием внешних факторов. Единство состава и структурная общность определяют целостность системы. Границы системы, по которым она может быть выделена как единое целое, определяются контрастом связей между элементами системы и связей с элементами других систем. Фундаментальная особенность геологической среды заключается в том, что это иерархически структурированная среда. Наиболее полная характеристика минералообразующих систем раскрывается при системном подходе, основой чего служит построение их комплексных моделей.

В соответствии с полученными данными и составляются комплексные модели разномасштабных металлогенических объектов (провинций, зон, рудных районов и полей, месторождений, рудных тел), которые можно рассматривать как минеральные системы различных иерархических уровней. Региональные прогнозно-поисковые модели составляются с целью металлогенического районирования, что предполагает анализ региональной обстановки (карт), отражающей свойства минеральной системы, и выявление площадей с сочетанием ее признаков, перспективных для выявления месторождений. Для металлогенических зон и рудных районов сочетание признаков соответствует уровню масштабности (см. рис. 1), присущих этим таксонам; их набор характеризуется обобщенным сочетанием формационных и геолого-структурных факторов и геодинамических условий их проявления. Составление локальных прогнозно-поисковых моделей имеет целью использовать различные локальные проявления свойств рудообразующей системы, проявленной в эталонном месторождении, для поиска новых месторождений этого типа. В рудных полях характеристика модели более детальная и соответствует масштабу геологоразведочных работ.

Рис. 3. Этапы прогнозирования и оценки крупных месторождений и перспективных рудопроявлений стратегических металлов 

Выполнен сравнительный анализ современных отечественных и мировых тенденций прогнозирования крупных месторождений стратегических и высокотехнологичных металлов. Анализ показал, что в настоящее время наиболее перспективным направлением в прогнозировании — создание комплексных моделей рудообразующих систем. Это направление прогнозирования основывается на интеграции генетических моделей рудных месторождений со знаниями, полученными из других дисциплин, таких как геофизика, геохимия, дистанционное зондирование Земли, пространственный ГИС анализ и геолого-экономический анализ. Практика применения комплексных моделей минералообразующих систем, наряду с применением современной аппаратуры и компьютерных технологий при прогнозе и поисках новых месторождений показывает отличные результаты.

В связи с этим, пространственно-статистический ГИС-анализ включал:
1. Анализ закономерностей размещения месторождений стратегических металлов в различных геологических формациях.
2. Анализ возраста рудовмещающих формаций и геодинамических обстановок их образования.
3. Выявление пространственно-временных связей рудных месторождений с различными свойствами земной коры: строением, вариации мощности, неоднородности границ между слоями коры, коры и верхней мантии и др.
4. Выделение потенциально перспективных геодинамических обстановок формирования руд и вмещающих толщ.
5. Анализ позиции месторождений в различной по глубинному строению земной коре с выделением ее дислоцированных областей, потенциально перспективных на выявление новых рудных месторождений.

В результате анализа в составлении и применении комплексных моделей рудообразующих систем наметились несколько направлений. В первую очередь использование установленной закономерности — приуроченности рудных объектов к локальным утолщениям земной коры, особенно ее нижнего слоя. Во-вторых, выявление пространственным ГИС анализом крупных ареалов различных по возрасту и типу минеральных систем и их совмещенности в пространстве, которая позволяет прогнозировать металлогеническую специализацию еще не открытых потенциальных комплексных месторождений.

Общая методология ГИС-анализа заключается в изучении пространственно-временной роли различной геодинамической природы блоков земной коры в локализации и времени образования металлоносных геологических формаций и связанных с ними месторождений полезных ископаемых. Для проведения ГИС-анализа применялись методические приемы, заложенные в аналитический аппарат ARCMAP, MAPINFO и других картографических математико-аналитических систем. Наиболее эффективными показали себя устоявшиеся методы выявления и оценки пространственных связей (растровой алгебры, нечеткой логики, вероятностного анализа и др.).

Этапы прогнозирования и оценки крупных месторождений и перспективных рудопроявлений стратегических металлов показаны на рисунке 3.

Металлогенические исследования показали, что прогнозирование стратегических месторождений определяет комплекс признаков, включающий: (1) оценку потенциальной перспективности рудовмещающих геологических формаций на основе геодинамических обстановок их формирования в том числе; (2) исследование глубинных структурных преобразований земной коры, определяющих региональные критерии выявления новых крупных месторождений стратегических металлов, в т.ч. и не выходящих на поверхность.

Рис. 4. Размещение Pb-Zn-месторождений MVT, SEDEX и VMS на фоне сочетаний дислоцированных областей коры по границам Мохо, Конрада (нижнего), ее среднего и верхнего слоев. Составлена на основе модели CRUST1.0

Кроме того, необходимо учитывать, что для выявления закономерностей пространственных связей требуется использовать данные не только для изучаемой площади, но и для ее обрамления. Поэтому пространственно-статистический анализ осуществлялся с охватом значительно большего пространства, а также всего северного полушария Земли.

На территории РФ выделены три базовые геодинамические обстановки размещения стратегических месторождений высокотехнологичных металлов: выступы древнего основания, пассивная континентальная окраина и островодужные комплексы активной окраины. Дальнейшее развитие представлений о региональных рудообразующих системах предполагает выделение типовых областей — источников стратегических металлов, типовых геодинамических блоков земной коры, проницаемых для флюидов и магм.

Для ГИС-анализа пространственных соотношений геологической структуры и геодинамических обстановок формирования месторождений стратегических металлов России использованы результаты современных исследований литосферы на основе гравитационных данных, полученных космическим аппаратом GOCE (Gravity Field and SteadyState Ocean Circulation Explorer): глобальные карты глубины поверхности Мохо, мощности и строения осадочного чехла, модель CRUST1.0 [7]. Модель CRUST1.0 унаследована из модели CRUST2.0 [8] мощности и строения осадочного чехла [9] — содержит данные о глубине Мохо, плотности трехслойного осадочного чехла, а также верхнего, среднего и нижнего слоев консолидированной коры. В данной публикации результаты ГИС-анализа на основе этой модели продемонстрированы на примере Pb-Zn-месторождений MVT, SEDEX и VMS типов (рис. 4).

Общая пространственная статистика показывает, что в пределах локальных зон утолщения коры размещается более 65 % всех мировых месторождений и проявлений благородных, цветных, редких, черных и благородных металлов. В геофизической глобальной модели земной коры наиболее резкие локальные утолщения и утонения ее глубинных слоев проявлены в виде обособленных зон аномальной изменчивости (градиента) мощности средней и нижней коры (рис. 4).

Положение этих зон в целом соответствует коллизионным областям (Тетис, Урал, Байкало-Патом, Алданский щит, Таймыр, Приморье). По соотношению градиентных зон в размещения Pb-Zn-месторождений отмечается их отчетливая приуроченность к эпизонам над блоками наиболее изменчивой по мощности нижней корой (рис. 4). Причем колчеданное оруденение тяготеет при этом к областям с утоненной верхней корой. Таким образом, из пространственных соотношений размещения свинцовоцинковых месторождений и различия структуры вмещающих блоков коры вытекают основные глобальные закономерности размещения Pb-Zn месторождений.

Первое — это приуроченность MVTи SEDEX-месторождений к флангам крупных осадочных бассейнов, сложенных платформенными и субплатформенными комплексами с устойчивыми признаками нефтегазоносности. Такие территории отмечаются на Северо-Востоке России, где, в частности, на флангах Вилюйской нефтегазоносной провинции размещаются недоизученное месторождение Менгенилер на севере и объекты Сарданского узла на юге.

Второе — высокая изменчивость мощности нижней коры, вследствие интенсивного развития разрывных нарушений и диапиризма. При этом, месторождения MVT-типа, часто совместно с объектами SEDEX-типа, тяготеют к верхней и осадочной коре, а руды VMS-типа, часто вместе с эксгаляционными рудами (SEDEX) — над наиболее изменчивыми по мощности участками нижней и средней коры. Расположение этих потенциально перспективных областей показано на рисунке 5. Учитывая непрямой характер признака, следует рассматривать его лишь как косвенный инструмент для совместного использования с прямыми признаками (коренные проявления руд и геохимические аномалии).

Следующий этап в развитии представлений о перспективности и распространении региональных рудообразующих систем состоит в выделении типовых областей, являющихся источниками стратегических металлов, типовых геодинамических блоков земной коры, являвшихся проницаемыми каналами для продвижения флюидов и магм. Комплексное изучение месторождений позволит решить проблему сырьевых источников большинства стратегических металлов.

Рудообразующая система (тип месторождения)	Геологическая обстановка	Главные металлы	Попутные критические металлы	Примеры месторождений
Магматическая базитовая и ультрабазитовая (Ni–Cu коматиитовый, Ni–Cu, PGE и Fe–V–Ti в расслоенных базитовых и ультрабазитовых интрузиях)	Зеленокаменные докембрийские пояса, островодужные офиолиты, орогенные пояса	Ni, Cu, Со, PGE, Fe, V, Ti,	Au, Ag, Te, Se	Норильская и Печенгская группы, Кингаш, Чинейское, Медведковское
Скарново-порфировая и порфирово-эпитермальная (Cu–Mo–Au- и Cu-Moпорфировый, Au-Ag и PbAg эпитермальный, Cu-Au и Pb-Zn скарновый)	Островодужные, окраинно- и внутрикон- тинентальные вулканоплутонические пояса	Cu, Mo, Au, Ag	Re, Sb, Pt, Pd, Pb, Zn, Te, In, Se, Bi, Cd, W, REE	Песчанка, Малмыж, Быстринское, Михеевское, Синюхинское, Сорское, Бугдаинское, Дукат, Гольцовое,
Орогенная, связанная с интрузивами гранитоидов (Au–Bi, Sn–W, Mo-W- и Можильно-штокверковый, Mo-W скарновый), Ta-LiSn-Be пегматитовый	Магматические пояса, области и зоны тектономагматической активизации пассивных континентальных окраин	Au, Sn, W, U, Mo, B, Sb, Li, Та, Ве	Re, Bi, In, Cd, Ge, REE, Ag, Pb, Zn, Te, Co, U	Правоурмийское, Депутатское, Инкурское, Холтосонское, Вишняковское, Тигриное, Забытое, м-я Колымского пояса
Колчеданные вулканогенные и осадочно-вулканогенные (Cu–Zn и Cu-Zn–Pb)	Островодужные вулканические пояса, зоны рифтогенеза	Cu, Zn, Pb, Ag, Au	Bi, Cd, Te, In, Se	Уральская и Алтае-Саянская группы
Осадочных бассейнов (Миссисипи тип Zn–Pb–Ag, медистые песчаники (Cu– Co–Ag), несогласий (U), долинный (U) и др.)	Платформы, энсиалические бассейны, рифтогенные зоны, пассивные континентальные окраины	Cu, Co, Mn, Zn, Pb, U, Ag	Re, Bi, In, Cd, Ge, REE, Ag, Pb, Zn, Te, Sc, Sr,Ga, Ga	Сардана, Павловское, Удокан, Хиагдинская группа, Дулурское
Связанная с щелочными гранитами (Nb, Ta, Zr, Ве, REE-месторождения)	Внутриплитная, все типы бассейнов и фундамента	Nb, Ta, Ве, Zr, REE, криолит	U, Th, Au, Ag, Te, Co, Sr, Rb, Pb, Zn, Sn	Катугинское, Зашихинское, Улуг-Танзек
Связанная с массивами ультраосновных-щелочных пород и карбонатитов (Nb, Ta, апатитовые, железорудные, REE- и U-Th месторождения)	Внутриплитная, все типы бассейнов и фундамента	U, Th, Nb, REE, Ta, апатит, магнетит	Au, Ag, Te, Co, Sr,	Хибинская группа,  Томтор, Белозиминское, Чуктуконское и др.
Россыпей и кор выветривания (Au, Sn, W, PGE, Ti–Zr, палеороссыпи Au–U), вторичного обогащения (Cu, Zn и Pb безсульфидный, бокситы, Ni–Co–Cr латеритовый), соленых озер (Li–K–B)	Все	Au, Sn, Ti, W, Zr, Ni, Fe, Li, B, алмазы, PGE, K, бокситы	REE, Co, Cr, Ga, Sc, Mn, Re	Ичувеемская, Рывеемская, Берелехская, Чайюрьинская, Тирехтях, Одинокое, Млелювеем, Серовское, Красная шапочка, Баскунчак

Табл. 1. Краткая характеристика главных рудообразующих систем стратегических металлов России
Примечание: таблица составлена на основе [6]

Пространственно-статистический ГИС-анализ размещения Pb-Zn месторождений на территории Карело-Кольского региона Балтийского щита и сравнение геодинамических обстановок их формирования [10] позволил выделить ниже следующие закономерности. Позднеархейские зеленокаменные пояса, которые контролируют размещение медно-никелевых месторождений магматической ультрабазитовой рудообразующей системы, приурочены к утоненной верхней коре. Ареалы лопийских гранитоидных (порфировой и орогенной) рудообразующих систем с медно-порфировыми и медно-молибден-порфировыми, а также с серебро-полиметаллическими и скарновыми месторождениями и проявлениями, располагаются на наиболее утолщенной коре (рис. 6а). А палеопротерозойские стратифицированные комплексы (островодужные вулканические пояса и зоны рифтогенеза), вмещающие колчеданные и полиметаллические месторождения, располагаются над участками утолщения нижнего слоя и соответственного увеличения мощности земной коры до 42 и более км (рис. 6б). Полученные результаты имеют практическое значение для региональных прогнозно-металлогенических построений, поисков и оценки полиметаллических месторождений.

Заключение

Прогнозные карты, созданные на основе развития концепции рудообразующих систем, пространственно-статистического ГИС-анализа и глубинной геолого-геофизической модели земной коры — основной результат разработки новых подходов к прогнозированию металлогенических исследований последних лет — полезный инструмент для определения районов с высоким потенциалом размещения определенных рудообразующих систем, используют все знания и данные, доступные на территории проекта. Эти карты могут быть полезными геологоразведочным и горнодобывающим компаниям, занимающимся поисками и оценкой месторождений стратегических металлов, для выбора районов для новых проектов или уточнения их поисковых площадей, а также помогут ранжировать поисково-оценочные цели и выделения области с потенциалом различных типов минерализации. Кроме того, они могут использоваться правительственными структурами для предоставления предконкурсных материалов геологоразведочным компаниям, а также для принятия внутренних решений по недропользованию.

Рис. 5. Сводные признаки на выявление новых Pb-Zn-месторождений MVT, SEDEX и VMS 

Различные цели будут влиять на детализацию и масштаб проводимых работ, а также на разрешение данных, необходимых для выполнения задачи. Однако в каждом из этих сценариев крайне важно убедиться, что результаты — статистически достоверные, геологически значимые и практически полезные. Если конечный пользователь — геологоразведочная компания, то карты, которые сокращают область поиска до небольшой цели, могут считаться приоритетными. Напротив, если результаты используются для определения перспективных областей, для поисков новых месторождений и для целей планирования недропользования, то на картах выделяются большие площади.

Выполненные исследования показывают, что помимо выделения перспективных на открытие новых месторождений площадей в пределах существующих рудных районов, дополнительные рудные районы без известной экономической минерализации также могут быть выделены как перспективные. Эти районы предоставляют новые возможности для получения новых данных с более высоким разрешением для дальнейшего анализа с целью направления геологоразведочных работ регионального масштаба. Разработка более удобных для отображения на прогнозно-металлогенических картах репрезентативных исходных данных должно стать основным организующим принципом будущих металлогенических исследований.

Рис. 6. Поверхность Мохо, глубина в км (а). Соотношение мощности нижней и средней коры (б). По материалам [9]

Анализ показал, что успешное дальнейшее развитие отмеченных выше направлений связано с применением оригинальных компьютерных технологий: нейронных сетей, методов опорных векторов и случайного леса. Однако, несмотря на огромный потенциал в оценке ресурсов, современные модели на основе машинного обучения обладают недостатками, такими как переобучение, более длительное время вычислений и эффект сглаживания. Эти ограничения можно устранить, объединив их с другими методами машинного обучения и алгоритмами глубокого обучения.

1. Рундквист Д.В. Пульсационная гипотеза С.С. Смирнова в свете новых данных о процессах рудообразования // Проблемы региональной металлогении и эндогенного рудообразования. Л.: ВСЕГЕИ, 1968 г. C. 46–66.
2. Pirajno F. Hydrothermal processes and mineral systems. Springer, 2009, 1250 p.
3. Joly A., Porwal A., McCuaig T.C., Chudasama B., Dentith M.C., Aitken A.R.A. Mineral systems approach applied to GIS-based 2D-prospectivity modeling of geological regions: Insights from Western Australia // Ore Geology Reviews. 2015. Vol. 71. P. 673–702.
4. Hagemann S.G., Cassidy K.F. Archean orogenic lode gold deposits // Rev Econ Geol. 2000. Vol. 13. P. 9–68.
5. Hagemann S.G., Lisitsin V., Huston D.L. Mineral system analysis: Quo Vadis // Ore Geology Reviews. 2016. Vol. 76. P. 504–522.
6. Бортников Н.С., Волков А.В., Галямов А.Л., Викентьев И.В., Аристов В.В., Лаломов А.В., Мурашов К.Ю. Минеральные ресурсы высокотехнологичных металлов в России: состояние и перспективы развития // Геология рудных месторождений. 2016 г. T. 68. № 2. C. 97–119.
7. Волков А.В., Галямов А.Л., Савчук Ю.С. Применение моделей глубинного строения земной коры и верхней мантии, созданных на основе гравитационных данных спутника ГОСЕ, в металлогеническом анализе // Исследование Земли из космоса. 2020 г. № 4. С. 41–50.
8. Artemieva I.M., Meissner R., Crustal thickness controlled by plate tectonics: a review of crust– mantle interaction processes illustrated by European examples // Tectonophysics. 2012. Vol. 519. P. 3–34
9. Laske G., Masters G., Ma Z., Pasyanos M.E. Update on CRUST1.0 — A 1-degree global model of Earth’s crust. // Geophys. Res. Abstr. 2013. 15. EGU2013-2658.
10. Галямов А.Л., Волков А.В., Лобанов К.В. Размещение полиметаллических месторождений и глубинное строение земной коры Карело-Кольского региона // Труды Ферсмановской научной сессии ГИ КНЦ РАН. 2023. Т. 20. С. 304–311.

Опубликовано в журнале «Золото и технологии», № 3 (61)/сентябрь 2023 г.