Геология рудных месторождений, 2020, T. 62, № 5, стр. 462-474

Эпитермальное Cu–Au месторождение Дудица (Республика Северная Македония)

Г. Тасев a, Д. Серафимовский a, А. В. Волков b*, Б. Боев a, Т. Джорджевич a, У. Колич cd, Т. Серафимовский a

a Университет “Госе Делчева”
Штип, Р. Македония

b Институт геологии рудных месторождений, петрографии, минералогии и геохимии РАН
119017 Москва, Старомонетный пер., 35, Россия

c Институт минералогии и кристаллографии, Венский университет
1090 Вена, Альтанстр, 14, Австрия

d Отдел минералогии и петрографии Музея Естественной истории
1010 Вена, Бургринг, 7, Австрия

* E-mail: alexander@igem.ru

Поступила в редакцию 27.08.2019
После доработки 08.12.2019
Принята к публикации 14.03.2020

Полный текст (PDF)

Аннотация

Месторождение Дудица расположено в южной части Республики Северная Македония (РСМ), недалеко от границы с Грецией в Козуф-Аридском вулканическом районе. Генетически месторождение связано с неогеновым вулканизмом (≥6 Ma, K/Ar), проявившимся на пересечении северо-восточных разрывных нарушений с Вардарской зоной разломов северо-западного простирания. РЗЭ-спектр вулканических пород характеризуется слабой отрицательной Eu-аномалией и обогащением легкими РЗЭ относительно тяжелых РЗЭ. Минерализация встречается в жилах, прожилках и вкрапленности, концентрация Cu – около 0.5%, Au – до 1 г/т. Доминирующие рудные минералы – халькопирит, борнит, энаргит, ковеллин, халькозин и дигенит, реже встречаются пирит, галенит, сфалерит, марказит и самородное золото. Изотопный состав δ34S находится в узком интервале значений (от +1.00 до +2.50‰), что указывает на ее магматическое происхождение. Полученные результаты свидетельствуют о принадлежности месторождения Дудица к высокосульфидизированному классу эпитермальных месторождений.

Ключевые слова: Республика Северная Македония, Вардарская зона, Au–Cu месторождение Дудица, РЗЭ, минеральная ассоциация, изотопия серы, порфирово-эпитермальная система

ВВЕДЕНИЕ

Месторождение Дудица расположено в 2 км от границы РСМ с Грецией (фиг. 1). Месторождение периодически изучалось (1916–1917, 1940–1944, 1946–1948 гг.), но не разрабатывалось, за исключением добычи самородной серы в незначительных объемах. На месторождении выполнен достаточно большой объем геологоразведочных работ: пройдены несколько штолен, одна из которых длиной более 400 м, и пробурено более 4 км колонковых скважин, глубина которых варьировала от 50 до 690 м. Работы проводились исходя из концепции, что месторождение Дудица принадлежит к колчеданно-полиметаллическому классу (аналог месторождения Бор в Сербии). Краткая информация о месторождении Дудица была опубликована в статьях: Hiessleintner, 1945; Ivanov, 1965; Jankovic, 1967.

Фиг. 1.

Позиция месторождения Дудица в Вардарской зоне. 1 потенциальные рудные районы; 2 третичные вулканические породы; 3 Cu–Au-месторождение; 4 Cu-месторождение; 5 населенный пункт; 6 перспективные рудные месторождения; 7 Козуфская вулканическая зона; 8 региональные линеаменты; 9 участки опробования. N – неогеновые и палеогеновые осадочные толщи; a, q – вулканиты (третичные); C – меловые осадочные толщи; ν, ββ – габбро-диабазы (юрские); Gr – гранитоиды (юрские); Pz – палеозойские метаморфические породы; G – гнейсы (докембрийские).

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

МЕТОДИКА ИССЛЕДОВАНИЙ

Каменный материал для литогеохимии и изотопного датирования был отобран в 2015 г. из рудных отвалов разведочных штолен и вмещающих пород. Микроэлементный состав отобранных образцов был определен с использованием комбинации методов индуктивно связанной плазмы (ICP-AES), масс-спектрометрии с индуктивно связанной плазмой (ICP-MS) и инструментального нейтронно-активационного анализа (INAA). Породообразующие оксиды определены методом ICP-AES. Точность, определяемая воспроизводимостью лабораторных стандартов и дубликатов, в пределах 5% для основных породообразующих оксидов и 10% – для второстепенных и микроэлементов.

K–Ar-возраст горных пород определен в лаборатории Института ядерных исследований ATOMK, Дебрецен, Венгрия. Калий определен методом пламенной фотометрии с буфером Na и внутренним стандартом Li с использованием пламенного фотометра (Sherwood M420). Точность и воспроизводимость этого метода находятся в пределах 2%. Аргон извлечен в условиях сверхвысокого вакуума путем радиочастотного индукционного нагрева и плавления образцов горных пород в молибденовых тиглях. Отношения изотопов аргона измерялись масс-спектрометрическим методом разбавления изотопов 38Ar, предварительно откалиброванным с использованием атмосферного аргона и международных стандартов горных пород.

Для получения значений δ34S для одного минерала сульфиды отделяли от образцов руды путем ручного отбора под бинокулярным микроскопом. Образцы пирита, халькопирита, галенита и энаргита были измельчены до 200 меш с использованием агатовой ступки. Анализ изотопов серы выполнен в лаборатории ACME, департамента геологических наук, Королевского университета г. Кингстон, Онтарио, Канада.

РЕГИОНАЛЬНОЕ ГЕОЛОГИЧЕСКОЕ СТРОЕНИЕ И МЕТАЛЛОГЕНИЯ

Рудное поле Дудица расположено на западном фланге Вардарской рифтогенной зоны, на ее контакте с Пелагонским кристаллическим массивом (фиг. 1) и находится в пределах плиоценового вулканоплутонического андезит-кварц-латитового комплекса в Аридии (Греция) и Козуф в РСМ. Детально изучена только северная Козуфская часть этого района (Ivanov, 1965; Jankovic et al., 1997; Bogojevski, 1998; Serafimovski, Tasev, 20131; и др.).

Козуф-Аридская металлогеническая зона плиоценовой минерализации занимает особое положение на крайнем юге РСМ, располагаясь вдоль границы с Грецией. Вулканические структуры и субвулканические интрузии известково-щелочного состава приурочены к неотектоническим разломам северо-восточного направления (Кочнева и др., 2006).

Докембрийские альбититовые гнейсы со спорадическими линзами амфиболитов относятся к самым древним породам в Козуф-Аридском районе. Изредка в гнейсах встречаются блоки мраморов. Локально развиты палеозойские сланцы, филлиты, метапесчаники, глинистые сланцы и кварциты. Триасовые отложения занимают большую часть района. Они представлены двумя основными фациями: 1) мраморизованными известняками и доломитами и 2) аргиллитами и песчаниками со спорадическими включениями диабазов и изредка зеленых сланцев. Позднемеловые отложения представлены баремиан-альбскими конгломератами и туронскими известняками. Верхнеэоценовые толщи состоят из базальных конгломератов, перекрытых флишевыми отложениями (алевролитами, глинами, песчаниками и пластами известняка). Плиоценовые отложения состоят из конгломератов и глинистых песчаников с пропластками известковых глин. Вулканогенно-осадочные пирокластические породы и отложения глинистых песчаников развиты в отдельных плиоценовых бассейнах. Большая часть бывшего вулканического кратера Аридия перекрыта плиоценовыми осадками. Четвертичные отложения слагают террасы.

Третичный магматизм проявился после закрытия мезозойского океанического бассейна (Karamata, 1984). Это закрытие связано со сближением Динаридской плиты и Карпато-Балканского террейна с Сербо-Македонским массивом и последовавшей коллизией этих континентальных сегментов (Dimitrijevic, 1974). Завершающий этап третичного магматизма в зоне Вардар (включая Козуфский вулканизм) в позднем плиоцене имеет типичные субдукционные характеристики (Boev, Yanev, 2001).

Характеристика Козуфского вулканического комплекса подробно рассмотрена в статье Б. Боева (Boev, 1988). Его состав представлен различными типами андезитов и кварцевых латитов. Доля интрузивных обнажений на площади развития этого комплекса незначительна, вследствие низкого уровня эрозии. Вулканоплутоническая деятельность в районе продолжалась от 7 до 1.8 млн лет.

Верхнеюрская субдукция сопровождалась известково-щелочным магматизмом во время среднего и верхнего мела (Boev, Jelenkovic, 2012). Послесубдукционная экстенсиональная тектоника генерировала внутриконтинентальные синорогенные разломы, рифтовые бассейны и обширные орогенные третичные магматические пояса плутонических и субвулканических пород I типа, которые в течение миоцена частично дифференцировались до задугового вулканизма (Marchev et al., 2005).

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

Система кольцевых разломов нашла отражение в нескольких морфологически отрицательных формах рельефа (что видно на космоснимках) и положительной структуре в районе Дудицы (Кочнева и др., 2006).

Система разломов от СВ-ЮЗ до СЮ ориентировки – более молодая, чем Вардарская система, и характеризуется недавней сейсмической активностью.

ГЕОЛОГИЧЕСКОЕ СТРОЕНИЕ МЕСТОРОЖДЕНИЯ

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

Фиг. 2.

Геологическая карта месторождения Дудица (Ivanov, 1965; модифицирована: Serafimovski, Tasev, 20131). 1 – аллювий; 2 – гумусовый и склоновый материал; 3 – андезит; 4 – гидротермально измененный андезит; 5 – вторичный кварцит; 6 – дацит, кварцевый порфир; 7 – мраморизованный известняк; 8 – известняк; 9 – песчаник; 10 – палеозойский сланец; 11 – разломы; 12 – положение минерализации.

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

В южной части рудного поля палеозойская толща перекрыта верхнемеловыми (сеннонскими) местами битуминозными известняками. К западу от рудного поля Дудица эта толща известняков (мощностью 400–600 м) вмещает Fe–Ni-месторождение Ржаново (Серафимовский и др., 2013).

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

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

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

Рудоносная зона реки Мирчевица составляет около 2 км в длину и 500–700 м в ширину. Линзовидные вытянутые залежи медных руд с небольшими раздувами расположены вдоль реки Мирчевица (фиг. 2, 3) и локализованы в интенсивно гидротермально измененных андезитах и меловых мраморированных известняках. Мощность минерализованных приконтактовых зон варьировала от 0.2 до 20 м. Глубина распространения детально не исследована.

Фиг. 3.

Геологический разрез через месторождение Дудица (Ivanov, 1965; модифицирован). 1 – аллювий; 2 – андезит; 3 – гидротермально измененный андезит; 4 – известняк; 5 – палеозойский сланец; 6 – зона смятия и дробления; 7 – рудное тело; 8 – разлом; 9 – буровая скважина.

На рудном поле Дудица установлены следующие морфологические типы медного оруденения:

1. Гидротермальные жилы образовались вдоль разломов, в основном в измененных вулканитах (фиг. 3). Преобладают жилы мощностью от 0.1 до 0.3 м, в раздувах до 1.0 м, обычно длиной от нескольких до 30 м, иногда более. В штольне № 6 обнаружено несколько энаргит-пирит-кварцевых жил, мощность которых колеблется от 0.1 до 1.0 м, которые в ряде случаев залегают субпараллельно друг к другу. Жилы содержат: от 3 до 5% Cu, до 18 г/т Ag и Au (в основном около 0.1 г/т, местами 1.1 г/т).

2. Система прожилков (местами в виде линейных штокверков) обычно развивается вдоль разломов, ширина которых может достигать нескольких метров. Содержание Cu обычно 0.5–1.5%. Содержание Ag и Au не определялось.

3. Вкрапленная минерализация встречается преимущественно в сланцах и/или вдоль их контакта с андезитами. Штокверково-вкрапленный тип оруденения (Cu-порфировое) установлен на глубоких горизонтах колонковыми скважинами. Халькопирит и пирит – основные минералы в этом типе минерализации. Содержание Cu, по данным бурения, колеблется от 0.1 до 0.5%. Содержание Ag и Au не определялось.

Пространственное распределение ряда сульфидных минералов на месторождении Дудица представляет особый интерес. Оно изучено с поверхности (отметка 1800 м) и до горизонта 700 м (Jancovic et al., 1997). В интервале 1800–850 м установлены иллит, каолинит, серицит, ангидрит и пирит. В интервале 1650–1230 м развит энаргит. Медьсодержащая группа минералов (борнит, ковеллин, халькозин и тетраэдрит) отлагались в относительно узком диапазоне от 1500 до 1220 м. Халькопирит – основной Cu-содержащий минерал – отлагался в интервале 1400–850 м.

В южной части рудного поля установлены джаспероиды, образовавшиеся в результате окварцевания известняков и мраморов, которые могут содержать золотую минерализацию Карлинского типа, такую же, как на месторождении Алшар (Волков и др., 2006). До сих пор эти джаспероиды не были исследованы на золото. Однако, по данным двух штуфных проб (одна проба с обнажения), эти породы содержат 7.2 и 13.4 г/т золота (Jankovic et al., 1997).

Линзовидные залежи самородной серы – продукт сольфатарной активности, довольно часто встречается вдоль разломов на рудном поле Дудицы. В течение 1917 года здесь было добыто около 20 т самородной серы (Jankovic et al., 1997). Проявление самородной серы на рудном поле представляет скорее генетический, нежели экономический интерес.

ГЕОХИМИЯ ТРЕТИЧНЫХ ВУЛКАНИЧЕСКИХ ПОРОД

Химический состав (основные, микро- и редкоземельные элементы) восьми образцов вулканических пород месторождения Дудица представлен в табл. 1 и согласуется с опубликованными ранее данными по соседним участкам Козуфского района (Boev, Yanev, 2001; Boev, Jelenkovic, 2012). Фотографии типичных образцов приведены на фиг. 4.

Таблица 1.  

Химический состав вмещающих пород месторождения Дудица

Элементы D1S1 D1S2 D2S1 D2S2 Lat YB3 Lat YB4 Trdac YB5 Trdac YB6
SiO2 (%) 63.02 60.82 62.12 59.69 58.67 60.86 62.72 65.08
Al2O3 (%) 16.89 16.89 17.14 17.26 17.81 18.20 17.84 17.04
Fe2O3(T) (%) 4.34 5.1 4.54 6.02 5.51 4.64 4.12 3.39
MnO (%) 0.073 0.061 0.092 0.061 0.11 0.11 0.08 0.08
MgO (%) 1.66 2.45 1.76 2.4 1.50 1.11 0.79 0.47
CaO (%) 3.09 4.83 2.95 4.67 5.48 4.10 3.64 5.04
Na2O (%) 3.49 3.59 3.39 3.58 4.05 4.35 4.09 4.34
K2O (%) 3.62 3.58 3.61 3.55 4.71 4.75 4.77 3.84
TiO2 (%) 0.473 0.648 0.487 0.651 0.71 0.52 0.50 0.43
P2O5 (%) 0.27 0.29 0.27 0.3 0.68 0.56 0.54 0.54
LOI (%) 2.58 2.21 2.67 2.06 0.78 0.80 0.90 0.47
Total (%) 99.51 100.5 99.03 100.2 100.01 100.00 99.99 100.72
Sc (г/т) 6 11 7 11 15 10 15 10
Be (г/т) 4 3 4 3
V (г/т) 78 108 79 106
Cr (г/т) 100 50 100 50 26 25 20 20
Co (г/т) 10 14 11 12 20 20 10 10
Ni (г/т) < 20 20 < 20 20 30 20 <20 20
Cu (г/т) 20 20 20 50        
Zn (г/т) 60 60 50 50 80 100 20 20
Ga (г/т) 20 21 20 20
As (г/т) < 5 < 5 < 5 9
Rb (г/т) 132 117 133 116 174 180 190 200
Sr (г/т) 994 1056 987 1070 1100 1170 1250 1250
Y (г/т) 10 15 11 15 34 34 26 23
Zr (г/т) 208 207 199 197 200 210 210 220
Nb (г/т) 9 9 8 10 7 6 4 1
Sn (г/т) 1 1 1 2
Sb (г/т) < 0.5 1.1 0.6 1.6
Cs (г/т) 7.3 4.8 7.2 4.8 42 41 40 39
Ba (г/т) 1520 1539 1522 1544 1800 1760 1950 2100
La (г/т) 66.9 59.8 56.9 66.1 85 85 63 66
Ce (г/т) 130 114 115 122 145 157 125 115
Pr (г/т) 12.4 12.1 11.3 13 13.6 15.1 12.4 12.9
Nd (г/т) 42.7 44.5 39.8 47.6 54.8 56.9 41.3 51.2
Sm (г/т) 5.7 6.4 5.9 7 8.13 9.1 7.2 6.8
Eu (г/т) 1.43 1.68 1.41 1.78 2.0 1.9 1.42 1.38
Gd (г/т) 3.7 4.6 3.6 4.7 4.3 4.9 3.9 4.6
Tb (г/т) 0.4 0.6 0.5 0.6 0.75 0.78 0.7 0.7
Dy (г/т) 2.3 3 2.4 3.2 3.3 3.8 2.9 3.3
Ho (г/т) 0.4 0.6 0.4 0.6 0.8 0.9 0.9 0.8
Er (г/т) 1.3 1.6 1.3 1.7 2.0 2.3 2.2 1.6
Tm (г/т) 0.2 0.23 0.21 0.26 0.33 0.36 0.34 0.28
Yb (г/т) 1.3 1.5 1.4 1.6 2.01 1.85 1.8 1.7
Lu (г/т) 0.22 0.25 0.21 0.25 0.30 0.28 0.34 0.38
Hf (г/т) 4.5 4.4 4.4 4.3 6 5 5 4
Ta (г/т) 0.9 0.8 0.9 0.8 0.8 0.8 0.7 0
W (г/т) 3 3 3 4 3 5 6 4
Tl (г/т) 0.8 0.7 0.8 0.7 1.1 1.5 4.1 3.3
Pb (г/т) 61 52 59 54 55 63 55 62
Th (г/т) 28.3 22.2 26.6 22.9 28 31 28 28
U (г/т) 8.2 7 8.1 6.9 8 9 6 8
Sr/Y 99.4 70.4 89.73 71.33 32.35 34.41 48.08 54.35
La/Yb 51.46 39.87 40.64 41.31 42.29 45.95 35.00 38.82

Примечание. D1S1, D1S2, D2S1 и D2S2 – последние пробы; Lat YB3, Lat YB4, Trdac YB5 и Trdac YB6 (Boev, Yanev, 2001).

Фиг. 4.

Фотографии последних образцов из района Дудицы, отправленных на химический анализ: а – образец D1S1; б – образец D2S1.

Все образцы могут быть охарактеризованы как высококалиевые щелочно-шошонитовые породы: SiO2 58–65 мас. %; К2О 3.55–4.77 мас. % (табл. 1, фиг. 5; LeMaitre et al., 1989; Müller et al., 1992; Pearce, 1982; и др.). В породах отмечены незначительные вариации содержаний SiO2, а концентрация некоторых основных элементов стабильна (мас. %): Na2O – 3.39–4.35; CaO – 2.95–5.48 и MgO – 0.47–2.45.

Фиг. 5.

Классификация горных пород района Дудица (наши анализы) и опубликованные данные по Козуфскому району (Lat YB3, Lat YB4, Trdac YB5, Trdac YB6; Boev, Yanev, 2001).

На основной классификационной диаграмме ЅіО2 – (K2O+Na2O), (Middlemost, 1994) вулканические породы Дудицы и прилегающих магматических зон попадают в область трахиандезитового и трахидацитового состава (фиг. 5а). График зависимости K2O от SiO2 (фиг. 5б) показывает, что большинство образцов относятся к области высококалиевых пород и шошонитов; однако мы не можем исключить, что сдвиг в шошонитовое поле не был вызван метасоматическими изменениями анализируемых пород.

Редкоземельные элементы (РЗЭ) вмещающих пород (табл. 1) были нормированы к стандартным значениям РЗЭ хондрита (Boynton, 1984). Значения и характер распределения РЗЭ в новых пробах из района Дудицы и пробах Козуфской области (Boev, Yanev, 2001) практически не различаются (фиг. 6). Обогащение образцов легкими РЗЭ относительно тяжелых РЗЭ указывает на их происхождение из одного источника базальтовой магмы в пределах субконтинентальной области. Почти все образцы показывают слегка отрицательную Eu-аномалию.

Фиг. 6.

РЗЭ-нормализованные графики проб различных горных пород района Дудицы (наши анализы) и опубликованные данные по Козуфскому району (Boev, Yanev, 2001).

Геохронологические результаты K/Ar-датирования для четырех образцов горных пород месторождения Дудица представлены в табл. 2. Измеренные значения радиометрического возраста (5.0–6.5 Ма) перекрывают друг друга в пределах аналитической погрешности.

Таблица 2.  

K/Ar геохронологические анализы вмещающих пород месторождения Дудица

Проба K, % 40Ar rad (ccSTP/g ) 40Ar rad (%) K–Ar возраст ± 1 sigma
1 D1S1 3.484 8.1973 × 10–7 18.9 6.04 ± 0.44
2 D2S1 3.392 7.7663 × 10–7 19.8 5.88 ± 0.41
3 BBD1 L 4.360 8.0230 × 10–7 30.3 5.00 ± 0.40
4 BBD2 LQL 2.550 6.4080 × 10–7 25.1 6.50 ± 0.41

Примечание. № 1–2 наши пробы; № 3–4 пробы Б. Боева (Boev, 1988).

ИЗОТОПНЫЙ СОСТАВ СЕРЫ

Для изучения δ34S были приготовлены и отправлены на анализ четыре монофракции пирита, три – халькопирита, две – галенита и две – энаргита. Результаты δ34S приведены в табл. 3. Значения δ34S для пирита варьируют от 1.0 до 2.3‰, среднее 1.80‰, для халькопирита – от 1.5 до 2.1‰, среднее 1.70‰, для галенита – от 2.1 до 2.5‰, среднее 2.3‰, а для энаргита – от 1.8 до 2.4‰, среднее 2.1‰. Все анализируемые сульфидные минералы показали значения δ34S от 1.0 до 2.5‰, среднее 1.9‰, с узким диапазоном всего 1.5‰.

Таблица 3.  

Значения δ34S в сульфидах месторождения Дудица

№ пробы Минерал δ34S ‰
D1 Пирит 1.0
D2 Пирит 2.3
D4 Пирит 2.0
D5 Халькопирит 1.6
D6 Халькопирит 2.1
D7 Халькопирит 1.5
D8 Галенит 2.5
D9 Галенит 2.1
D10 Энаргит 1.8
D11 Энаргит 2.4

ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ И ЗАКЛЮЧЕНИЕ

Как уже упоминалось выше, нормированные на хондрит РЗЭ-значения показали нормальную дифференциацию с очень слабой отрицательной аномалией Eu (фиг. 6). Отсутствие значительных аномалий Eu, вероятно, связано с высокой степенью магматического окисления (Hanson, 1980; Carmichael, Ghiorso 1990; Moore, Carmichael 1998). Незначительное увеличение концентрации Sr в этих образцах с увеличением содержания SiO2 исключает возможность значительного фракционирования плагиоклаза в процессе эволюции пород района Дудицы. Расчеты по формулам (см. ниже) также показали отрицательные аномалии Eu и Ce. За исключением одного значения Eu* = 1.034, слегка положительного (образец 5-Lat YB3).

${\text{Eu}}* = \frac{{{\text{E}}{{{\text{u}}}_{N}}}}{{\sqrt {{\text{S}}{{{\text{m}}}_{N}} \cdot {\text{C}}{{{\text{d}}}_{N}}} }},$
${\text{Ce}}* = \frac{{{\text{3C}}{{{\text{e}}}_{N}}}}{{2{\text{L}}{{{\text{a}}}_{N}} + {\text{N}}{{{\text{d}}}_{N}}}}.$

В изученных образцах Eu* = 0.7547–1.0343, в среднем 0.9077, Ce* = 0.8354–0.9851, в среднем 0.9151, что указывает на слабо-окислительный характер рудообразующей среды (Bau, 1991). Полученные результаты аналогичны данным по другим участкам Козуфского района (Jankovic et al., 1997; Boev, Yanev, 2001).

Диаграммы Y к Rb, Ba к Rb и Sr к Rb для этих пород (фиг. 7a–в) показывают увеличение содержания Rb и Y с более или менее постоянными концентрациями Sr и Ba. Эти особенности вместе с очень слабой аномалией Eu свидетельствуют об ограниченном участии плагиоклаза в процессе дифференциации при формировании латитов и трахидацитов. Клинопироксен, вероятно, основной фракционировавший минерал при формировании этих пород (Boev, Yanev, 2001). На основе дискриминационной РЗЭ-диаграммы (фиг. 7д), козуфские породы показывают характеристики, типичные для связанных с субдукцией вулканических дуг.

Фиг. 7.

Диаграммы Козуфских трахидацит-трахиандезитов: a – Y к Rb; б – Ba к Rb; в – Sr к Rb, и д – дискриминационная диаграмма (Pearce et al., 1984). WPG – внутриплитные граниты, syn-COLG – синколлизионные граниты, VAG – граниты вулканических дуг и орогенов.

Образцы вмещающих пород (латит, трахидацит, андезит и доломит) из других участков Козуфского района характеризуются близкими, субпараллельными спектрами распределения РЗЭ (Jankovic et al., 1997; Boev, Yanev, 2001). Профили РЗЭ всех проб (фиг. 6) показали, что они обогащены LREE и имеют относительно плоские участки HREE со слабо отрицательными аномалиями Eu. Это согласуется с гипотезой о том, что HREE менее подвижны, чем LREE, и поэтому остаются преимущественно во вмещающих породах (McLennan, 1989). Обогащение легких РЗЭ (LREE) относительно тяжелых РЗЭ (HREE) – общая особенность всех образцов месторождения Дудица и Козуфского района (Jankovic et al., 1997; Boev, Yanev, 2001).

Высокие значения Sr/Y (32.35–99.40) и La/Yb (35.0–51.46) (фиг. 8) могут быть индикатором высокого давления и/или водосодержащего фракционирования (Kolb et al., 2013). Обогащение легкими РЗЭ (LREE) относительно тяжелых (HREE) в изученных образцах (фиг. 6), по сравнению с океаническими толеитовыми базальтами, указывает на участие во фракционировании корового материала (Najime et al., 2012). Такой процесс фракционирования подтверждается графиками (La/Yb)N к LaN и (La/Sm)N к SmN (фиг. 8a и б соответственно), как показано в (Chen et al., 2017).

Фиг. 8.

Диаграммы: a – La/Sm к Sm; б – La/Yb к La, месторождение Дудица.

ΣРЗЭ в анализируемых породах относительно высока (240.33–340.17 г/т), причем ΣLREE имеет более высокие значения (230.31–325.00 г/т), чем ΣHREE (9.82–15.17 г/т), что отвечает высоким значениям ΣLREE/ΣHREE = ~18.958–26.388.

Как уже упоминалось ранее, измеренные K–Ar радиометрические возрасты всех образцов (D 1; D 2; BBD1 L и BBD2 LQL) перекрываются друг с другом в пределах аналитической погрешности и предполагают возраст магматизма в районе Дудицы от миоцена до плиоцена (табл. 2). Все образцы очень схожи с петрографической точки зрения (трахидациты).

Таким образом, K–Ar-возраст можно считать минимальным возрастом (реальный геологический возраст может быть немного старше, чем аналитический возраст).

Полученные значения возраста показали высокую сходимость с предыдущими данными для пород в Козуфском районе 5.0–6.5 Ма (Boev, 1988), а также для отдельных минералов (биотит, полевой шпат) от 3.9 до 5.1 Ма (Lipolt, Fuhrman, 1986; Boev, Jelenkovic, 2012).

Что касается данных по изотопии серы, то мы сравнили значения δ34S сульфидных минералов с аналогичными данными в пределах Сербо-Македонской металлогенической провинции (Тулар-Киселяк, Плавица, Бучим, Боров-Дол и Скуриес; фиг. 9). Наиболее близкие значения δ34S с таковыми из Дудицы от 0.00 до +2.50‰ известны в Бучимском месторождении (Serafimovski, 1993; Serafimovski, Tasev, 20132; Serafimovski et al., 2016) и от +0.41 до +1.24‰ в месторождении Плавица (Serafimovski, 1993). В то время как для месторождения Скуриес значения δ34S находятся в интервале от –3.20 до +1.28‰ δ34S (Serafimovski, 1993). Для месторождений Тулар-Киселяк и Боров-Дол характерно незначительное расширение диапазонов в сторону отрицательных значений δ34S (Serafimovski, 1993), от –3.36 до +2.55‰ и от –4.05 до +0.72‰ соответственно (фиг. 9). Таким образом, полученный спектр значений δ34S близок к таковому для магматической серы – 0 ± 5‰ (Ohmoto, Rye, 1979).

Фиг. 9.

Значения δ34S сульфидов из Тулар-Киселякского, Плавицкого, Бучимского, Боров-Дола, Дудицы и Скуриеского месторождений.

Определенный диапазон значений δ34S в сульфидных минералах из месторождений Дудицы – репрезентативен для рудной минерализации в породах гранитного состава и порфировых Cu(±Mo)-месторождениях (Douglas et al., 2003). Узкий диапазон значений δ34S для пирита, халькопирита, галенита и энаргита от +1.0 до 2.5‰ предполагает относительно однородный источник серы (Ohmoto, Rye, 1979; Hedenquist et al., 2017).

Результаты, полученные в конце прошлого века, как и наши исследования, указывают на сложную структуру месторождения Дудица, которую можно наиболее адекватно сравнить с таковой на месторождении Плавица (Серафимовский и др., 2017), где самые высокие концентрации золота связаны с телами вторичных кварцитов. Подобные рудные тела широко развиты на высокосульфидизированных эпитермальных месторождениях по всему миру (Sillitoe, Hedenquist, 2003; Richards, 2009). Аналогичная Дудице минерализация обнаружена в Тулярском месторождении Сербии (Janković, 1997), которая на глубине трансформируется в Cu-порфировые руды. Кроме того, аналогичные проявления найдены на рудном поле Петрошница вблизи Куманово, Республика Северная Македония (M’nkov et al., 2015).

Наличие алунитизации, наряду со вторичными кварцитами, пористым кварцем, залежами самородной серы, а также аргиллитовый тип изменений свидетельствуют о развитии на рудном поле Дудица высокосульфидизированной эпитермальной минерализации (энаргит, пирит, золото) (Sillitoe, Hedenquist, 2003).

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

  1. Волков А.В., Серафимовский Т., Кочнева Н.Т., Томсон И.Н., Тасев Г. Эпитермальное Au–As–Sb–Tl месторождение Алшар (Южная Македония) // Геология руд. месторождений. 2006. Т. 48. № 3. С. 205–224.

  2. Кочнева Н.Т., Волков А.В., Серафимовский Т., Тасев Г., Томсон И.Н. Тектоническая позиция Au–As–Sb–Tl месторождения Алшар (Македония) // ДАН. 2006. Т. 407. № 2. С. 175–178.

  3. Серафимовский Т., Волков А.В., Боев Б., Тасев Г. Метаморфизованное латеритное Fe–Ni-месторождение Ржаново (Республика Македония) // Геология руд. месторождений. 2013. Т. 55. № 5. С. 443–460.

  4. Серафимовский Т., Волков А.В., Серафимовский Д., Тасев Г, Ивановский И. Мурашов К.Ю. Эпитермальное Аu–Аg–Сu месторождение Плавица (Восточная Македония): геологическое строение и 3D-модель распределения полезных компонентов в рудах // Геология руд. месторождений. 2017. Т. 59. №. 4. С. 296–304.

  5. Bau M. Rare-earth element mobility during hydrothermal and metamorphic fluid-rock interaction and the significance of the oxidation state of europium // Chemical Geol. 1991. V. 93. № 3–4. P. 219–230.

  6. Boev B., Yanev Y. Tertiary magmatism within the Republic of Macedonia: a review // Acta Vulcanologica. 2001. V. 13. № 1–2. P. 57–71.

  7. Boev B., Jelenkovic R. Allchar Deposit in Republic of Macedonia. Petrology and Age Determination // Petrology – New Perspectives and Applications. 2012. V. 64. P. 132–168.

  8. Boev B. Petrological, geochemical and volcanic features of volcanic rocks of the Kozuf Mountain. Doctoral Thesis, Štip: University “St. Cyril and Methodius”, 1988.

  9. Bogojevski K. Gold in Macedonia: Geology, Metallogenic features, Ore occurrences, Gold Ore deposits and their evaluation. Skopje: University “St. Cyril and Methodius”, 1998.

  10. Boynton W.V. Geochemistry of the rare earth elements: meteorite studies // Henderson, P. (ed.), Rare Earth Element Geochemistry. Oxford: Elsevier, 1984. P. 63–114.

  11. Carmichael I.S.E., Ghiorso M.S. The effect of oxygen fugacity on the redox state of natural liquids and their crystallizing phases // Rev Mineral, Min Soc America. 1990. V. 24. P. 191–212.

  12. Chen B., Long X., Wilde A.S. et al. Delamination of lithospheric mantle evidenced by Cenozoic potassic rocks in Yunnan, SW China: A contribution to uplift of the Eastern Tibetan Plateau // Lithos. 2017. V. 284–285. P. 709–729.

  13. Dimitrijević M. The Serbo-Macedonian Massif // Tectonic of the Carpathian-Balkan Regions. Bratislava: Geological Institute of Dyoniz Stura, 1974. P. 291–296.

  14. Douglas T.A., Chamberlain C.P., Poage M.A. et al. Fluid flow and the Heart Mountain fault: A stable isotopic, fluid inclusion, and geochronologic study // Geofluids. 2003. V. 3. P. 13–32.

  15. Hanson G.N. Rare earth elements in petrogenetic studies of igneous systems // Annu. Rev. Earth Planet. Sci. 1980. V. 8. P. 371–406.

  16. Hedenquist W.J., Arribas R.A., Aoki M. Zonation of sulfate and sulfide minerals and isotopic composition in the Far Southeast porphyry and Lepanto Epithermal Cu–Au Deposits, Philippines // Resource Geology. 2017. V. 67. P. 174–196.

  17. Hiessleitner G. Das Enargitvorkommen Dudice in Mazedonien // Geol. Bundesanst. 1945. V. 1. P. 53–93.

  18. Ivanov T. Metallogeny of the southern part of the Vardar zone. Doctoral thesis, Belgrade: Faculty of Mining and Geology, 1965.

  19. Jankovic S. Deposits of metalic mineral resources. Belgrade: Faculty of Mining and Geology, 1967.

  20. Jankovic S. The Carpatho-Balkanides and adjacent area: a sector of the Tethyan Eurasian metallogenic belt // Mineralium Deposita. 1997. V. 32. P. 426–433.

  21. Janković S., Boev B., Serafimovski T. Magmatism and Tertiary mineralization of the Kozuf metallogenetic district, the Republic of Macedonia with particular reference to the Alsar deposit. Skopje: University “Sts. Cyril and Methodius”. Special Issue, № 5. 1997.

  22. Karamata S. Plate tectonic phenomena in the regions of the Tethys type // Geotectonics. 1984. V. 17. P. 52–66.

  23. Kolb M., Von Quadt A., Peytcheva I. et al. Adakite-like and normal arc magmas: distinct fractionation paths in the East Serbian Segment of the Balkan-Carpathian Arc // J. Petrol. 2013. V. 54. № 3. P. 421–451.

  24. LeMaitre R.W., Bateman P., Dudek A. et al. A classification of igneous rocks and glossary of terms. Oxford: Blackwell Scientific Publishers, 1989.

  25. Lippolt H.J., Fuhrmann U. K–Ar age determination on volcanics of Alsar mine/Yugoslavia. Workshop on the feasibility of Solar Neutrino Detection with 206Pb by geochemical and mass spectroscopical measurements. Nolte E., ed. Report GSI-86-9. Munchen: Technische Univer, 1986.

  26. Marchev P., Kaiser-Rohrmeier M., Heinrich C., et al. Hydrothermal ore deposits related to post-orogenic extensional magmatism and core complex formation: The Rhodope Massif of Bulgaria and Greece // Ore Geol. Rev. 2005. V. 27. P. 53–89.

  27. McLennan S.M. Rare earth elements in sedimentary rocks: influences of provenance and sedimentary processes // Rev. in Mineral. and Geoch. 1989. V. 21. P. 169–200.

  28. Middlemost E.A.K. Naming materials in the magma/igneous rock system // Earth-Science Rev. 1994. V. 37. № 3–4. P. 215–224.

  29. M’nkov Sl., Antonov M., Siroshtan D. et al. Report for the geological mapping and geophysical studies of the territory of the concession area Petroshnitsa, Republic of Macedonia-Project 2015. Sofia: Geomac-Ltd, 2015.

  30. Moore G., Carmichael I.S.E. The hydrous phase equilibria (to 3 kbar) of an andesite and basaltic andesite from western Mexico: constraints on water content and conditions of phenocryst growth // Contrib Mineral Petrol. 1998. V. 130. P. 304–319.

  31. Müller D., Rock N.M.S., Grooves D.I. Geochemical discrimination between shoshonitic and potassic volcanic rocks in different tectonic settings: a pilot study // Mineralogy and Petrology. 1992. V. 46. P. 259–289.

  32. Najime T., Abaa S.I., Magaji S., Petrology and Rare Earth Elements (REE) distribution patterns of magmatic rocks in Gboko area, Lower Benue trough Nigeria: implication for tectonic evolution // Global Journal of Geological Sciences. 2012. V. 10. P. 47–58.

  33. Ohmoto H., Rye R.O. Isotopes of sulfur and carbon / Barnes, H.L. (ed.), Geochemistry of Hydrothermal Ore Deposits, 2 nd ed. New York: John Wiley, 1979. P. 509–567.

  34. Pearce J.A., Harris N.B.W., Tindle A.G. Trace element discrimination diagrams for the tectonic interpretation of granitic rocks // J. of Petrology. 1984. V. 25. P. 956–983.

  35. Richards J.P. Post-subduction porphyry Cu-Au and epithermal Au deposits: Products of remelting of subduction-modified lithosphere // Geology. 2009. V. 37. P. 247–250.

  36. Serafimovski T., Tasev G. Metallogeny of the Kozuf ore district, R. Macedonia. Eds. Serafimovski and Boev. 1st International Workshop on The Project Enigma, Proceedings Book, 20131. P. 11–22.

  37. Serafimovski T., Tasev G. Sulfur isotope compositions from different type of deposits in the Buchim-Damjan-Borov Dol ore district, Eastern Macedonia: 10th Applied Isotope Geochemistry Conference, Hungarian Academy of Sciences, 22–27th September 2013, Budapest, Hungary, 20132. P. 8–13.

  38. Serafimovski T. Structural metallogenetic characteristics of the Lece-Chakidiki zone: Types of Mineral Deposits and Distribution. Stip: Special Edition, 1993.

  39. Serafimovski T., Tasev G., Srmić-Palinkaš, S. et al. Porphyry Cu mineralization related to the small Tertiary volcanic intrusions in the Bučim ore deposit, Eastern Macedonia // Geologica Croatica. 2016. V. 69. P. 103–121.

  40. Sillitoe H.R., Hedenquist W.J. Linkages between volcanotectonic settings, ore-fluid compositions, and epithermal precious metal deposits // Soc. Econ. Geol. Sp. Publ. 2003. V. 10. P. 315–343.

Дополнительные материалы отсутствуют.