Трехмерная геолого-геофизическая модель Сорского медно-молибденового комплекса по электромагнитным и другим геофизическим данным

Размещено на http://www.allbest.ru/

Центр геоэлектромагнитных исследований ИФЗ РАН

ТРЕХМЕРНАЯ ГЕОЛОГО-ГЕОФИЗИЧЕСКАЯ МОДЕЛЬ СОРСКОГО МЕДНО-МОЛИБДЕНОВОГО КОМПЛЕКСА ПО ЭЛЕКТРОМАГНИТНЫМ И ДРУГИМ ГЕОФИЗИЧЕСКИМ ДАННЫМ

Спичак В.В., Гойдина А.Г.

Аннотация

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

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

Основная часть

сейсмический скорость месторождение молибденовый

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

Рис. 1 Карта геологических формаций и рудных месторождений Кузнецкого Алатау [16]. Формации: 1 - гранитовая, 2 - тоналит-платогранит-гранодиоритовая, 3 - трахибазальтовая, 4 - формации кор выветривания (нерасчлененные), 5 - известковая, 6 - карбонатно- терригенно-пепловая (туффитовая), 7 - алеврито-пигановая, 8 - известняково- песчанико-алевритовая, 9 - карбонато-терригенно-пепловая, 10 - угленосная паралическая, 11 - гранитосиенит-гранитовая, 12 - трахибазальт-трахиатерит-трахилинарит; 13 - разломы; 14 - участок Сорского месторождения [1]

Трехмерная модель удельного сопротивления

В работе [2] по данным магнитотеллурического зондирования была построена блочная модель удельного сопротивления региона. На ее основе нами построена геоэлектрическая модель участка (см. его расположение на рис. 1), включающего Сорское месторождение (рис. 2а). Ее анализ позволяет сделать следующие выводы. Во-первых, фоновое удельное сопротивление среды составляет примерно 10⁴ Ом·м, т.е. область в целом является высокоомной. Во-вторых, на глубинах 10-40 км выделяется листрический разлом, простирающийся в северо-западном направлении и выполаживающийся к северо-востоку на глубинах Мохо. К юго-западу от него от поверхности до глубин больше 40 км прослеживается проводящая структура, а на северо-востоке выделяется высокоомный батолит с удельным сопротивлением 10⁵-10⁶ Ом·м. На глубинах 32-37 км расположена относительно проводящая зона диаметром примерно 5 км с удельным сопротивлением в центральной части от 10 до 40 Ом·м, а на поверхности - две области с удельным сопротивлением от 40 до 150 Ом·м, соединенные каналами с удельным сопротивлением 500-1000 Ом·м с нижнекоровой аномалией.

Рис. 2 Объемные модели удельного сопротивления (а), отношения скоростей сейсмических волн V_p/V_s (б), плотности (в) [1]

Трехмерные модели скоростей сейсмических волн

В работе [1] были также построены трехмерные модели скоростей продольных и поперечных сейсмических волн (V_p и V_s) по их двумерным разрезам [3]. На рис. 2б показана модель отношения V_p/V_s, которое изменяется от 1.74 до 1.88 (медианное значение составляет 1.81), что соответствует мощности консолидированной коры 43-45 км [3]. При этом в аномальных областях V_p/V_s составляет 1.77-1.78, т.е. уменьшение от медианного значения составляет примерно 2.3 %.

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

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

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

Трехмерная модель плотности

Применение уравнения регрессии, полученного на основе множественной корреляции параметров V_p, V_s и V_p/V_s, позволило построить модель плотности пород рассматриваемого участка (рис. 2в). Анализ этой модели показывает, что плотность во всей рассматриваемой области изменяется от 2.50 до 3.33 г/см³ (медианное значение - 2.91 г/см³), что покрывает практически весь диапазон типов пород - от кислых до ультраосновных. При этом общее распределение плотности достаточно однородно, за исключением областей, в которых наблюдаются отмеченные выше аномалии удельного сопротивления и сейсмических скоростей. В них значения плотности составляют от 2.67 до 2.75г/см³, что на 5.7-8.5 % меньше среднего значения.

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

Оценка литологического состава пород по содержанию кремнезема

Оценка типов пород может быть основана на общей классификации всех магматических пород земной коры по их химическому составу и, в первую очередь, по содержанию и соотношению в породах кремнезема и щелочей. Распределение содержания кремнезема (K, %) в коре рассматриваемого участка довольно неоднородно, что характерно для молодых складчатых структур, к которым он относится. При этом преобладают значения К в диапазоне 58-65 %, что соответствует представлению о «кислой - средней» коре. При этом в верхней и, в меньшей степени, средней коре заметно присутствие базальтовых пород с К < 54 %, в то время как в нижней коре заметно повышение концентрации более кислых пород с К > 65 %, связываемое с заложением геосинклиналей в складчатых областях континентальной коры. По содержанию кремнезема можно осуществить грубую оценку преобладающих типов пород в коре, не учитывающую содержание в них солей натрия и калия (рис. 4).

Рис. 4 Трехмерная модель литотипов пород [1]

Анализ построенной таким образом объемной модели литологических типов пород (рис. 4) показывает, что разрез представлен, главным образом, габбро-диоритами и диоритами. На этом фоне на западе участка выделяются области, сложенные, по-видимому, преимущественно кислыми (гранитными) магмами, которые прослеживаются от глубин примерно 32-37км, где можно предполагать расположение рудогенерирующего очага, и до поверхности.

Геолого-геофизическая модель месторождения

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

По геохимическим данным, в условиях литосферы континентальных окраин водосодержащий базальтоидный расплав поступал в кору с глубины около 70 км. Большой временной разрыв между лейкогранитами и габбро не позволяет допускать образование кислого расплава в процессе кристаллизационной дифференциации базальтоидной магмы. Можно предположить, что формирование кислой магмы происходило на нижнекоровом уровне при давлениях >10-12 кбар (что соответствует глубинам примерно 40 км) в связи с частичным плавлением предшествующих закристаллизовавшихся базитов. Обогащенные амфиболом породы в низах коры рассматриваются как потенциальный источник рудоносных флюидов, вовлекаемых в гидротермальный процесс на этапе развития лейкогранитов и рудоносных гранит-порфиров.

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

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

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

Список литературы

1. Спичак В.В., Гойдина А.Г. Геофизические индикаторы эндогенных рудных месторождений (на примере Сорского медно-молибденового комплекса). // Разведка и охрана земных недр. 2014. №10. C.47-51.

2. Белявский В.В., Гойдина А.Г. Трехмерная геоэлектрическая модель металлогенических зон Кузнецко-Алатауской складчатой области // Физика Земли. 2012. № 11-12. С. 97-117.

3. Кадурин И.Н., Белявский В.В. и др. Геофизические исследования сейсморазведочными и электроразведочными методами глубинного строения Алтае-Саянской складчатой области по профилям общей протяженностью 3300 км. // Отчет (Росгеолфонд № 487225). 2008.

Размещено на Allbest.ru