Необходимость наращивания ресурсной базы благороднометалльного сырья требует развития технологий анализа спутниковых изображений, что наиболее важно для удаленных и труднодоступных регионов. Островодужные системы привлекают исследовательский интерес в связи со своим контрастным геологическим строением, а также наличием как традиционных, так и необычных типов минерального сырья. Рудные месторождения редко- и благороднометалльного типов сырья, связанные с измененными вулканогенно-осадочными и магматическими породами [1], могут быть выражены в ореолах вторичных изменений, фиксируемых на космических изображениях, следовательно, прогноз и освоение минерального сырья островодужных систем возможны с привлечением методов дистанционного зондирования Земли. Природно-климатические факторы в данном случае выступают как осложняющие прогноз и оценку, однако их влияние может быть редуцировано тщательным подбором сезонных изображений для анализа, а также использованием техник машинного обучения. Целью настоящего исследования является составление и апробация карты минерагенического потенциала дистанционных изображений на примере относительно хорошо изученной в минерагеническом отношении территории острова Кунашир.
Материалы и методы исследования
Курильская островодужная система (КОС) простирается в северо-восточном направлении от острова Хоккайдо до полуострова Камчатка на 1175 км. КОС включает две субпараллельные островные дуги Большой и Малой Курильских гряд. Большая Курильская гряда (БКГ) характеризуется наличием большого количества вулканов и расположена в тылу КОС, со стороны задугового Охотоморского бассейна. Малая Курильская гряда (МКГ) является фронтальной, ограничена Курило-Камчатским желобом Тихого океана, и вулканизм для нее менее характерен. Область исследования настоящей работы включает в себя территорию острова Кунашир (рис. 1). КОС была сформирована в ходе надсубдукционного вулканического процесса, включает более чем 200 вулканов (105 являются субаэральными), большая часть локализована в пределах Кунашир-Урупского сегмента. Мощность континентальной коры в основании БКГ достигает 35–40 км. В структуре КОС выделяют фундамент (поздний мел – ранний палеоген), несогласно перекрытый олигоцен-неоген-четвертичными островодужными комплексами [1; 2]. Среди разрывных нарушений преобладают правые сдвиги северо-восточного простирания, сформированные косым взаимодействием Тихоокеанской плиты и окраины континента. Структуры второго порядка характеризуются левосторонним смещением и северо-западным простиранием [1].
Острова БКГ сложены олигоцен-неоген-четвертичными породами и представляют собой горстово-антиклинальную структуру с углами падения на крыльях до 30–40 °, осложненную изолированными вулканическими постройками. В пределах БКГ обосабливаются комплексы «зеленых туфов», вулканогенно-кремнисто-диатомитовый, базальтовый и андезитовый [1; 3; 4]. Вулканогенные породы острова Кунашир (рис. 1) обычно описываются как супрасубдукционные вулканиты Тихоокеанского типа [1; 4; 5]. Мелководно-морские высококремнеземистые осадочные породы плиоцена, включающие различные обломочные фракции и диатомиты, насыщены пирокластикой. Они переслаиваются пемзами и тефрами различного возраста и прорываются многочисленными экструзиями и субвулканическими телами (до 63-64 % SiO2) и обладающие выраженной столбчатой отдельностью. До 60 % площади островов сложено лава-пирокластическим материалом плейстоцена-голоцена состава с существенным преобладанием двупироксеновых базальтовых андезитов и андезитов [1; 3; 4]. Южные острова Большой Курильской гряды вмещают скопления рассеянных элементов (Mo, Bi, Cd, In, Ge, Re), благородных (Au-Ag) и цветных (Zn, Cu, Sn, Zn, Pb) металлов, а также полуметаллов (Se, Te, Sb, As). Распространены малоглубинные (≤1 км) золото-серебряные, золото-полуметаллические, золото-алунитовые, золото-адуляровые, золото-сульфидные рудопроявления. Интерес представляют также сольфатарные поля и термальные источники центров современного вулканизма (рис. 1). Территория имеет потенциал для обнаружения порфировых Cu-Mo-Au месторождений (глубина ≤ 4 км), связанных интрузивными массивами (глубина внедрения 5–10 км) [1; 6]. Вулканогенно-пирокластические породы комплекса «зеленых туфов» вмещают стратиформные эксгаляционно-осадочные и жильные гидротермально-метасоматические тела сульфидных руд Zn, Pb, Cu типа Куроко. Прорывающие их габброплагиограниты вмещают золото-серебряные руды с Se и Te. Экструзивные, субвулканические и интрузивные тела фельзитового состава, а также ореолы их вторичных изменений (аргиллизации, окварцевания, серицитизации, ожелезнения и другие) связаны с многочисленными золото-серебряными ( ± Sn) проявлениями [1]. Крупным известным объектом – примером эндогенной минерализации юга БКГ считается Прасоловский рудный узел, расположенный на Охотоморском побережье острова Кунашир (рис. 1, Б) [7]. Прасоловский рудный узел приурочен к интрузивно-купольному поднятию, образованному габбродиорит-плагиогранитным массивом (его возраст 33 ± 8 – 30 ± 2 млн лет (K-Ar) и 31 ± 1 млн лет по цирконам (U-Pb) [1]. На площади Прасоловского рудного узла присутствуют сульфидная медно-цинковая минерализация, а также благороднометалльная минерализация различных типов, связанная с жильно-метасоматическими телами и ореолами вторичных изменений, окружающих субвулканические тела [7].
Площадное распределение спектральных особенностей дистанционного изображения позволяет получить картину наличия определенного вещественного признака, включая развитие вторичных изменений и минерализации определенного типа. Для оконтуривания интересующих площадей возможно использовать оценку главных компонентов изображения и спектральных индексов, т.н. метод Кроста [8], сопоставляемых с данными рудоперспективности территории. Анализ производится с помощью авторской прикладной программы, написанной на языке Python. Для проведения исследования подбирались сцены Landsat 8 OLI без облаков, прошедшие геометрическую и спектральную коррекцию. Ввиду пересеченного рельефа территории исследований и невысоких углов возвышения светила изображения требуют топографической коррекции с применением алгоритма SCS + C [9].
Главные компоненты, рассчитанные для изображения, отражают вклад в интегральную картину отдельных каналов, их анализ достаточно известен для картирования вторичных изменений при металлогенических исследованиях [10]. Глинистые минералы характеризуются сильным поглощением в 7 канале Landsat 8 OLI (2,11–2,29 µm) и значительным отражением в канале 6 (1,57–1,65 µm), отношение этих спектральных каналов (7/6) представляют собой индекс CM (clay minerals – глинистые минералы). Отношение каналов 7/2 представляет собой индекс-индикатор развития минералов-гидроксидов (HA, hydroxyl alteration), а 4/2 – индикатор присутствия минералов – оксидов железа (IO, Iron Oxide) [10; 11]. Характеристика 7 спектральных каналов Landsat 8 OLI (105/029 path/row, дата съемки 27.10.2015), использованных для анализа участка интереса «Кунаширский», приведена в табл. 1.
Расчет главных компонент изображения осуществлялся с помощью функции pca модуля decomposition библиотеки scikit-learn Python. Сведения о ковариации каналов изображения и его главных компонентов приведены в табл. 2.
Рис. 1. A – геологическая карта Курильской островодужной системы (с использованием [3]); Б – геологическая карта и металлогения острова Кунашир
Таблица 1
Статистическая характеристика спектральной яркости изображения
|
|
band1 |
band2 |
band3 |
band4 |
band5 |
band6 |
band7 |
|
Minimun |
7660,00 |
6952,00 |
5963,00 |
5161,00 |
4993,00 |
4674,00 |
4628,00 |
|
Maximum |
33884,00 |
36697,00 |
38807,00 |
42869,00 |
42474,00 |
17273,00 |
13330,00 |
|
Mean |
8632,79 |
7911,12 |
6784,60 |
6153,18 |
6150,50 |
5405,09 |
5191,28 |
|
Median |
8545,00 |
7801,00 |
6618,00 |
5976,00 |
5444,00 |
5067,00 |
5017,00 |
|
Std.dev. |
671,39 |
776,81 |
904,92 |
1016,75 |
2241,78 |
1062,90 |
532,93 |
Таблица 2
Ковариация спектральных каналов изображения и их главных компонентов (PC)
|
PC1 |
PC2 |
PC3 |
PC4 |
PC5 |
PC6 |
PC7 |
|
|
band1 |
0,021 |
0,082 |
0,158 |
-0,044 |
-0,011 |
0,019 |
0,000 |
|
band2 |
0,023 |
0,098 |
0,182 |
-0,056 |
-0,003 |
0,026 |
0,000 |
|
band3 |
0,013 |
0,146 |
0,193 |
-0,073 |
0,008 |
0,009 |
0,000 |
|
band4 |
0,011 |
0,170 |
0,213 |
-0,079 |
-0,028 |
0,009 |
0,000 |
|
band5 |
-0,136 |
0,609 |
0,072 |
-0,065 |
-0,012 |
0,009 |
0,000 |
|
band6 |
-0,091 |
0,245 |
-0,096 |
-0,106 |
-0,016 |
0,004 |
0,000 |
|
band7 |
-0,044 |
0,119 |
-0,044 |
-0,057 |
-0,016 |
0,003 |
0,000 |
Отношения спектральных каналов Landsat 8 OLI и главные компоненты сопоставлялись с матрицей обратных евклидовых расстояний до эпитермальных рудопроявлений и месторождений, пересчитанные в пиксели и ограниченные радиусом 50 пикселей для каждого объекта, которые рассматривались как выражение рудоносности территории (IOD). Для близко расположенных рудных объектов перекрывающиеся радиусы суммировались. Значения IOD применялись для обучения и валидации прогнозной модели логистической регрессии (также известной как «максимальной энтропии», MaxEnt).
Территория исследования покрывалась сеткой квадратных окон осреднения (ОО) площадью 1 км2, для каждого полигона рассчитывалось среднее значение признака, хранящегося в растровых слоях (CM, IOA, PC1-7, HA, NDVI), включая значения отношений спектральных каналов, главных компонентов космического изображения и евклидовых расстояний до эпитермальных рудопроявлений и месторождений. Данные по каждому ОО помещались в атрибутивную таблицу, которая затем рассматривалась как таблица экспериментальных данных. Алгоритм оценки экспериментальных данных включает использование методов анализа особенностей отражательной способности поверхности, а также обучение и модели логистической регрессии, перспективной для прогнозных минерагенических исследований [12]. Используется класс OneClassSVM, входящий в библиотеку scikit-learn (рис. 2). Этот метод позволяет использовать неполные данные и оценивает вероятность принятия значения целевой переменной с максимальной энтропией и контролируется набором ограничений, представленных неполнотой информации о целевом распределении [13; 14].
Результаты исследования и их обсуждение
Для визуальной оценки и дешифрирования территории применяются композиции каналов RGB432 («естественные цвета»), а также псевдоцветные композиции: RGB642, RGB742, RGB765, RGB764, позволяющие визуально различить растительность, водные объекты, гидротермальные изменения субстрата (RGB765), развитие минералов гидроокислов и минералов – оксидов железа (HA, IOA, (HA+IOA)/2, «зона изменений») (рис. 3).
Подготовленная модель логистической регрессии была апробирована на территории острова Кунашир для данных, не использовавшихся при ее обучении. Результат построения моделью карты прогнозной перспективности территории острова Кунашир, показан на рис. 4, А. Оценка качества осуществлялась оценкой площади под ROC-кривой. Полученное нами значение 0,6540 характеризует приемлемо обученную модель (рис. 4, Б).
Рис. 2. Алгоритм прогнозного исследования с использованием спектральных данных и логистической регрессии, реализованный на языке Python
Рис. 3. Материалы для визуальной оценки территории обучения, Landsat 8 OLI (105/029 path/row, дата съемки 27.10.2015)
Рис. 4. А – результат, полученный при работе модели логистической регрессии, обученной по 60 % данных, на острове анализа продуктов Landsat 8 OLI; Б – ROC-кривая модели, обученной по 40 % исходных данных
Видно, что распределение рудных объектов в целом показано верно, показанные в оверлее рудопроявления и месторождения находятся в осевых частях выявленных «рудоперспективных» кластеров. Выделение этих кластеров производилось с помощью выделением пяти классов окон осреднения с помощью классификации естественных границ Дженкса. Три класса наименьших значений рассматривались как фоновые, два класса наибольших значений – условно перспективные. Качество применяемой модели также подтверждалось наличием внутренней зональности кластеров – значение IOD уменьшалось от центра кластера к периферии.
Выводы
Северо-восточное простирание наиболее перспективных кластеров (рис. 4А) соответствует разрывным нарушениям первого порядка и позволяет предположить их взаимосвязь с проницаемыми правосдвиговыми структурами. Это позволяет рассматривать простирание выделенных зон в качестве дополнительного поискового признака. Карта прогнозной перспективности территории острова Кунашир показывает соответствие перспективных кластеров известным рудным объектам и вулканическим структурам. Кластер на юго-западе площади, по периферии вулкана Менделеева, может рассматриваться в качестве объекта для дальнейшего геологического опоискования. Используемая модель эффективно работает по исходным данным, включающим спектральные характеристики и только известные объекты в условиях частичной обнаженности, позволяя воссоздать распределение эпитермального оруденения по дистанционной основе.