Исследование основано на данных геолого-­структурного анализа и изучения более чем 500 образцов пород, руд и метасоматитов, отобранных из поверхностных обнажений и более чем 3000 м задокументированных эксплуатационных, разведочных и поверхностных выработок месторождения золота Тунанца.

Петрографические и минералогические исследования включали макроскопическое изучение 315 образцов с месторождения Тунанца.

Оптическая микроскопия. Изучено 96 образцов, из которых 47 прозрачных шлифов и 10 аншлифов соответствовали фрагментам руды и пробам из обогащенной зоны; остальные 39 прозрачных шлифов относились к породам и метасоматитам. Исследования проводились на микроскопе Axioscope-5, оборудованном цифровой камерой и компьютером с программным обеспечением SIAMS на кафедре геологии месторождений полезных ископаемых геологоразведочного факультета Российского государственного гео­ло­го­раз­ве­доч­но­го университета имени Сер­го Орджо­ни­кид­зе (РГГРУ — МГРИ).

Физические параметры. Проанализировано 335 образцов, получены данные по плотности, магнитной восприимчивости и удельному электрическому сопротивлению. Анализы выполнены с использованием гидростатических весов, каппаметра Model KT-5 серия 7088 и Terrmоmmeter E6–3. Исследование проводилось на кафедре геофизики геофизического факультета Российского государственного гео­лого­раз­ведоч­ного университета имени Серго Орджони­кидзе (РГГРУ — МГРИ).

Анализ пространственного распределения физических параметров (плотности, магнитной восприимчивости и удельного электрического сопротивления) на исследуемой площади выполнен с использованием плагина Smartmap 1.5 для QGIS и метода интерполяции кригинга. Согласно [8][21], данный метод соответствует условиям настоящего исследования, поскольку данные характеризуются пространственной автокорреляцией, а сеть опробования является неоднородной: расстояние между пробами варьирует от 1 до 200 м. Анализируемая площадь соответствует полигону 1000 x 600 м, который плане включает крутопадающее тело брекчий, включает минерализованную зону месторождения и простирается до зон контакта брекчий с вмещающими породами, основана на данных 80 отобранных, проверенных и очищенных проб, со схожим размером, высокой целостностью и без признаков гипергенных изменений.

Золоторудное месторождение Тунанца, ранее известное как Пиунца, разрабатываемое по лицензии горного общества «Питука II», расположено на юге Эквадора, в провинции Самора-­Чинчипе, в 9 км к юго-востоку от города Самора (78°53’ в. д., 4°07’ ю. ш.), в Кордильере-дель-­Кондор. Район месторождения входит регион рифтогенеза пермо-­триасового возраста [17]. В этом же регионе находятся эпитермальное месторождение Фрута-дель-­Норте, с запасами 9,91 млн унций золота и 15,0 млн унций серебра [17], медно-­порфировый объект Мирадор с ресурсами 609 млн метрических тонн при содержании 0,58% Cu, 3,2 млн унций золота и 22 млн унций серебра [11], скарновое месторождение Намбиха и другие практически не изученные золотоносные зоны, к ним относится и месторождение Тунанца.

Месторождение Тунанца представляет собой мелко-среднезернистое эпитермальное место­рождение, расположенное внутри полимиктовой брекчиевой трубки фреатомагматического про­ис­хождения [2][4], которая прорывает гранитные породы батолита Самора (рис. 1). К югу от рудного поля расположены гнейсы и сланцы формации План-дель-­Осо предположительно неопротерозойско-­палеозойского возраста [15], слагающие останцы кровли над батолитом. На севере обнажаются порфировые андезиты, входящие в состав основания триасовой формации Пиунца [15]; они находятся в контакте с батолитом и также образуют останцы кровли над ним. На северо-­западе и юго-востоке наблюдаются афировые андезиты, залегающие несогласно на формации Пиунца, обрамляющие брекчиевый массив и входящие в состав формации Ла-Сакеа. 

Самора батолит простирается на 200 км в северо-северо-восточном направлении при ширине 50 км и является частью юрского известково-­щелочного вулканического пояса. Батолит состоит из двух различных интрузивных комплексов: 1 — ранний плутон преимущественно однородного состава, сложенный главным образом среднезернистыми гранодиоритами, которые обнажаются к юго-западу от рудного поля в контакте с брекчиями [4]; 2 — поздний комплекс в виде субвулканических интрузий, даек и штоками, образующими купольные структуры среднего или кислого состава, которые пересекают батолит в различных местах [11][17]. Именно они выделяются своей связью с зонами медно-­золотой минерализации [11]. Породы последнего типа относятся к риолитовым кварц-­полевошпатовым порфирам, присутствующим в составе брекчии в виде обломков и неправильных блоков размером более 5 м; они представляют собой реликты гипабиссального тела, внедрившегося в плутонические породы Самора батолита до образования брекчии [2][4].

Брекчии образуют неправильной формы массив с диаметром на поверхности примерно 1,4 км. Контактные поверхности круто падают к центру, что указывает на трубчатую, сужающуюся с глубиной форму. Они представлены преимущественно матриксными брекчиями, содержащими хаотично распределенные обломки разного размера. Их пересекают гидротермальные брекчии в виде субвертикальных прожилков и жил, беспорядочно распределенных в брекчиевом теле. На основании комплекса диагностических признаков эти брекчии относятся к фреатомагматическим [2][4].

Золоторудные зоны представлены двумя формами: первая — минерализованные брекчии с сульфидным цементом, образующие неправильные субвертикальные удлинённые блоки поперечником около 20 м; вторая — тонкие зигзагообразные кварц-­сульфидные жилы мощностью около 2 см, с почти вертикальным падением, протяжённостью до 30 м, расположенные радиально относительно первой формы. В целом минерализация локализована в центральной части брекчиевого тела. Рудные минералы включают пирит, сфалерит, галенит, халькопирит, арсенопирит и, в меньших количествах, другие сульфиды, которые в совокупности составляют более 60%. Нерудные минералы представлены кварцем и поздними карбонатами [2].

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

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

Физические свой­ства пород месторождения Ту­нанца, представленных обломками брекчий, имеют различные значения для каждого типа породы (табл. 2). 

Плотность изменяется от 2,44 до 2,75 г/см3 для пород и резко возрастает до 3,63 г/см3 для руд, что подчеркивает явный контраст между вмещающими породами и минерализованным телом. Магнитная восприимчивость сильно варьирует от 0,08×10-3 СИ для руды и роговиков до 5,05×10-3 СИ для базальта. Кварц-полевошпатовый риолитовый порфир имеет умеренные значения около 2,4×10-3 СИ, как и гранодиорит. Средние значения удельного электрического сопротивления в большинстве пород, как правило, превышают 1×106 Ом•м. Рудные образцы и сланцы имеют значения на порядок ниже 0,14×106 Ом•м.

Значения физических свой­ств брекчий варьируют в зависимости от степени их метасоматического изменения и присутствия рудных минералов (табл. 3). 

Средняя плотность пропилитизированных брекчий составила 2,46 г/см3, а охваченных кварц-серицитовым метасоматозом — 2,52 г/см3 и каолинит-алунитовым — 2,51 г/см3, что заметно выше. Надо отметить значительное увеличение плотности для брекчий с присутствием сульфидов со средними значениями, достигающими 2,64 г/см3, что обусловлено присутствием тяжелых рудных минералов пирита, галенита и сфалерита.

Магнитная восприимчивость отражает четкую нисходящую тенденцию: от значений, соответствующих неизмененным брекчиям, 2,15×10-3 СИ, которые прогрессивно снижаются при пропилитовом метасоматозе до 1,62×10-3 СИ, далее до 0,89×10-3 СИ при кварц-серицитовом изменении, достигая минимальных зна­че­ний при као­ли­нит-алу­ни­то­вом мета­со­ма­то­зе 0,25×10-3 СИ и в брекчиях с сульфидами, (0,18–0,25)×10-3 СИ. 

Удельное электрическое сопротивление брек­чий имеет различные значения. У неизме­нен­ного типа наблюдаются высокие средние значения 1,00×106 Ом•м, которые снижаются для пропилитов и кварц-серицитовых метасоматитов до 0,57×106 и 0,63×106 Ом•м соответственно. Брекчии, подверженные каолинит-алунитовым изменениям, демонстрируют увеличение сопротивления до 1,16×106 Ом•м. Брекчии с сульфидами обладают самым низким сопротивлением, 0,57×106 Ом•м, что обеспечивает уникальные характеристики, связанные с обогащенной золотом зоной. Рассчитаны коэффициенты корреляции Пирсона разных физических свой­ств брекчий (табл. 4).  

Корреляция между плотностью и магнитной восприимчивостью брекчий является средней или отрицательной, r = –0.63, что интерпретируется как общая тенденция к уменьшению магнитной восприимчивости при увеличении плотности. Корреляция между магнитной восприимчивостью и удельным сопротивлением является отрицательной, r = –0.48, вероятно, вследствие вариабельности значений удельного электрического сопротивления, как и корреляция между плотностью и удельным сопротивлением, которая является отрицательной, r = –0.33, демонстрируя слабую корреляционную зависимость.

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

Плотность варьирует в диапазоне от 2,42 до 2,52 г/см3. При пропилитовом метасоматозе сред­нее значение составляет 2,45 г/см3, которое незначительно снижается до 2,42 г/см3 при кварц-серицитовом изменении, а затем возрастает до 2,52 г/см3 для каолинит-алунитового мета­соматоза.

Магнитная восприимчивость отражает нисходящую тенденцию: от очень высоких значений, соответствующих неизмененным образцам 4,94×10-3 СИ, которые резко снижаются в пропилитах до 0,68×10-3 СИ, далее до 0,11×10-3 СИ в кварц-серицитовых метасоматитах, достигая минимальных значений в породах, охваченных каолинит- алунитовым изменениям 0,05×10-3 СИ. В образцах с сульфидами наблюдается незначительное увеличение до 0,12×10-3 СИ.

Удельное электрическое сопротивление демонстрирует сложную и контрастную картину в диапазоне от 0,01×106 до 2,00×106 Ом•м. В неизмененных образцах средние значения высокие, 1,08×106 Ом•м, которые снижаются при пропилитах до 0,25×106 Ом•м. В кварц-серицитовых мет­соматитах сопротивление возрастает до 2,00×106 Ом•м, а в образцах с каолинит-алунитовыми изменениями вновь снижается до 0,35×106 Ом•м. Образцы с сульфидами обладают самым низким сопротивлением, 0,01×106 Ом•м, создавая электрический контраст.

Коэффициенты корреляции Пирсона для значений физических свой­ств, различно измененных кварц-­полевошпатовых риолитовых порфиров представлены в таблице 6.

Корреляции между плотностью и магнитной восприимчивостью практически нет, r = –0,09. Практически то же надо отметить для корреляции между магнитной восприимчивостью и удельным электрическим сопротивлением, r = 0,20. Однако коэффициент корреляции между плотностью и удельным электрическим сопротивлением является сильной отрицательной, r = –0,71, что интерпретируется как связь между увеличением плотности и уменьшением сопротивления.

Распределение метасоматитов и руд с различными физическими свой­ствами проанализировано в плане и по рельефу участка месторождения Тунанца (рис. 2). 

Пространственный анализ выполнен с использованием метода кригинга в геоинформационной среде QGIS. Анализ пространственного распределения плотности показывает прогрессивное увеличение значений в рудной зоне, где величина достигает диапазона между 2,60 и 2,70 г/см3. На разрезе наблюдается пик значений, связанный с массовым присутствием сульфидов. В отличие от этого, в зонах, наиболее удалённых от рудного тела, соответствующих периферийным участкам и пропилитам, плотность значительно ниже в диапазоне 2,40–2,45 г/см3. Отсутствие прямой связи между вариациями значений плотности и границей между пропилитовой и кварц-серицитовой зонами позволяет предположить, что этот переход может быть размытым и постепенным. Данная особенность указывает на то, что изменения значений плотности в меньшей степени связаны с конкретным типом метасоматоза и в большей — с присутствием рудных минералов, которое не ограничивается разными метасоматитами.

Анализ пространственного распределения значений магнитной восприимчивости показывает значительное снижение в рудных зонах с величинами от 0,00 до 0,25×10-3 СИ, что проявляется как отрицательная аномалия на разрезе, эквивалентная магнитной аномалии, связанной с рудной зоной, лишённой ферромагнитных минералов. Значения прогрессивно увеличиваются в направлении к внешним частям кварц-серицитовой метасоматической зоны, достигая 2,00×10-3 СИ. Этот рост продолжается и сохраняется в диапазоне (2–4)×10-3 СИ в зоне пропилитов. Аномальные зоны со значениями выше 4×10-3 СИ связаны с присутствием базальтовых даек и блоков кварц-полевошпатового риолитового порфира. Выявленная закономерность демонстрирует связь между типом метасоматоза и магнитной восприимчивостью, указывая, что магнитная восприимчивость является хорошим маркером для поиска минерализованных зон и околорудных метасоматических ореолов.

Анализ пространственного распределения удельного электрического сопротивления показывает сложное пространственное зонирование. В области рудного тела значения находятся в диапазоне (0,50–0,75)×106 Ом•м. Они формируют поверхность, простирающуюся до определённых участков за пределами границы кварц-серицитовой зоны. Далее значения снижаются ниже <0,50×106 Ом•м в пропилитовой зоне и западной части кварц-серицитовой зоны, после чего на внешних участках постепенно возрастают, достигая диапазона (2,00–3,00)×106 Ом•м в областях, наиболее удаленных от рудного тела. Такое сочетание значений проявляется на разрезе, построенном по результатам анализа, в виде пологого поднятия, соответствующего рудной зоне, что предполагает наличие неконтрастной аномалии в ореоле пониженного сопротивления.

Определение физических свой­ств неизмененных пород, метасоматитов и руд, их пространственное распределение относительно рудного тела позволило установить определенные закономерности на месторождении золота Тунанца.

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

Установлено значительное влияние на изменения физических свой­ств метасоматических изменений и количества рудных минералов в рудовмещающих фреатомагматических брекчиях (см. табл. 3) и риолитового порфира (см. табл. 5). Плотность увеличивается, в то время как магнитная восприимчивость имеет тенденцию к снижению от пропилитов через кварц-серицитовые метасоматиты к рудному телу. Удельное электрическое сопротивление демонстрирует сложные вариации. Коэффициенты корреляции Пирсона (табл. 4, 6) показывают взаимосвязи: для брекчий плотность и магнитная восприимчивость имеют умеренную отрицательную корреляцию (r = –0,63), тогда как для порфира между плотностью и удельным электрическим сопротивлением наблюдается сильная отрицательная корреляция (r = –0,71). Это предполагает различные рудоконтролирующие признаки. В брекчиях присутствие рудных сульфидов ассоциировано с увеличением плотности, что совпадает со снижением магнитной восприимчивости, и, вероятно, является результатом метасоматического замещения ферромагнитных минералов, таких как магнетит, присутствующий в обломках брекчии, возможно, на пирит или сидерит [13][18]. Напротив, в порфире увеличение плотности, связанное с присутствием рудных сульфидов, коррелирует со снижением удельного сопротивления, а не с магнитной восприимчивостью, возможно, из-за того, что в этом типе породы не наблюдалось присутствия магнетита в метасоматитах. 

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

Определены физические свой­ства различных типов пород, метасоматитов и руд месторождения золота Тунанца. Установлено влияние метасоматических процессов и характера рудной минерализации на изменения этих параметров. Отчетливо проявлена зональность с ростом интенсивности процессов гидротермального изменения от периферийных участков с пропилитами к центральной части рудного тела. Соответственно установлен систематический рост плотности и резкое снижение магнитной восприимчивости. Это изменение связано с метасоматическим замещением первичных ферромагнитных минералов (магнетита в базальтах и обломках брекчий) сульфидами и карбонатами (пирит, сидерит) в сильно минерализованных и обогащенных золотом зонах.

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