Технический справочник по спектральной алгебре (Band Ratioing) для локализации рудных узлов
1. Введение: Почему одного канала недостаточно?
В профессиональной среде дистанционного зондирования существует понятие «Проклятие альбедо». Если мы смотрим на одиночный канал спутникового снимка (например, SWIR Band 6 спутника ASTER), яркость пикселя зависит от двух факторов:
- Состав породы: Отражательная способность минералов.
- Топография: Угол падения солнечных лучей. Южный склон горы всегда ярче северного, даже если они сложены одной и той же породой.
Для геолога топография — это шум. Нам нужно знать химию, а не рельеф. Решением является метод Спектральных Отношений (Band Ratioing). Это арифметическая операция деления яркости одного канала на яркость другого.
Физический смысл деления: Если склон затенен, яркость падает и в Канале А, и в Канале Б пропорционально. При делении (А / Б) коэффициент затенения в числителе и знаменателе сокращается. Результат: Мы получаем «плоское» изображение, где яркость зависит только от спектральных свойств материала. Тени исчезают, химия остается.
В этом руководстве мы разберем промышленные стандарты спектральных формул, применяемых GeoJet Exploration для поиска золота, меди и полиметаллов.
2. Группа индексов оксидов железа (Iron Oxide Indices)
Цель: Поиск зон окисления сульфидов (Gossans), «железных шляп», зон гематитизации. Работаем в диапазоне VNIR (0.4–1.0 мкм).
2.1. Физика процесса
Ионы трехвалентного железа (Fe3+) в минералах гематит (Fe2O3) и гетит (FeOOH) имеют уникальную спектральную подпись:
- Сильное поглощение в ультрафиолете и синем диапазоне (< 0.5 мкм) из-за переноса заряда (Charge Transfer).
- Высокое отражение в красном диапазоне (0.63–0.69 мкм).
- Кристаллическое поле: Широкая полоса поглощения около 0.85–0.90 мкм (ближний ИК).
Мы ищем пиксели, где разница между «Красным» и «Синим» максимальна.
2.2. Формулы для Landsat 8/9 (OLI)
Сенсор OLI (Operational Land Imager) является рабочей лошадкой для картирования железа.
- Iron Oxide Ratio (Индекс оксидов железа):
Formula: Band 4 (Red) / Band 2 (Blue)- Интерпретация: Высокие значения (яркие пиксели) указывают на красноцветные породы, лимониты, гематиты.
- Риск: Этот индекс также подсвечивает полевые шпаты (розовые граниты), которые могут не нести оруденения.
- Ferrous Iron Index (Индекс двухвалентного железа): Ищем минералы типа хлорита или карбонатов с Fe2+. Они отражают в SWIR, но поглощают в NIR.
Formula: Band 6 (SWIR1) / Band 5 (NIR)- Применение: Полезно для выделения «свежих», неокисленных вулканитов.
2.3. Формулы для Sentinel-2 (MSI)
Европейский спутник Sentinel-2 имеет лучшее разрешение (10 м), но другие диапазоны.
- Iron Oxide (Sentinel adaptation):
Formula: Band 4 (Red) / Band 2 (Blue)- Нюанс: Band 2 у Sentinel очень чувствителен к атмосферной дымке. Требуется жесткая атмосферная коррекция (Sen2Cor) перед расчетом, иначе индекс покажет «облака» вместо железа.
[СОВЕТ ЭКСПЕРТА] Никогда не используйте индекс Fe3+ изолированно в пустынях. Красные пески (эоловые отложения) дадут мощный ложный сигнал. Всегда накладывайте маску рельефа: рудное тело обычно приурочено к выходам коренных пород, а не к дюнам.
3. Группа индексов глинистых изменений (Clay Alteration Indices)
Цель: Поиск гидротермальных метасоматитов — аргиллизитов, филлитов, пропилитов. Работаем в диапазоне SWIR (2.0–2.5 мкм). Платформа: Только ASTER (или WorldView-3). Landsat здесь слишком груб.
3.1. Физика гидроксилов
Гидротермальные минералы содержат группы OH, встроенные в решетку.
- Al-OH (Алюминиевые глины): Каолинит, Алунит, Пирофиллит. Поглощение при 2.16–2.20 мкм. (Это ASTER Band 5 и Band 6).
- Mg-OH / CO3 (Магниевые силикаты и Карбонаты): Хлорит, Эпидот, Кальцит. Поглощение при 2.30–2.35 мкм. (Это ASTER Band 8).
3.2. Базовые формулы ASTER (Standard GeoJet Stack)
- OHI (Al-OH Bearing Minerals Index) — Индекс Алунита/Каолинита: Мы берем каналы на «плечах» полосы поглощения (где отражение высокое) и делим на канал в центре поглощения.
Formula: (Band 5 + Band 7) / Band 6- Механика: Band 5 (2.16 мкм) и Band 7 (2.26 мкм) находятся по краям ямы поглощения каолинита (2.20 мкм, Band 6).
- Интерпретация: Чем глубже яма поглощения в B6, тем выше значение индекса. Яркие пиксели = зоны кислотного выщелачивания (центр системы).
- Kaolinite Index (Специфический индекс Каолинита): Для разделения каолинита и алунита (у алунита более широкое поглощение).
Formula: Band 6 / Band 5- Примечание: Часто используется инвертированная форма B5/B6, нужно смотреть на софт.
- MGI (Mg-OH Bearing Minerals Index) — Индекс Хлорита/Эпидота/Карбонатов:
Formula: (Band 6 + Band 9) / Band 8- Механика: Поглощение происходит в Band 8 (2.33 мкм). Каналы 6 и 9 — это фон.
- Интерпретация: Выделяет пропилитовые зоны (периферия порфировых систем) и скарны (карбонаты).
- Alunite Index (Специфический индекс Алунита): Алунит имеет дополнительные полосы поглощения.
Formula: (Band 5 + Band 8) / (Band 6 + Band 7)- Важность: Алунит — маркер верхней части эпитермальной системы (Lithocap).
4. Группа индексов кремнезема (Silica Indices)
Цель: Поиск окварцевания (Silicification), кварцевых жил, вторичных кварцитов. Проблема: Кварц прозрачен в VNIR и SWIR. У него нет линий поглощения. Решение: Использовать TIR (Тепловой диапазон) и эффект Reststrahlen.
4.1. Формула ASTER TIR
В диапазоне 8–9 мкм кварц имеет низкую излучательную способность (выглядит холодным), в то время как другие минералы излучают нормально.
- Quartz Index (QI):
Formula: Band 11 / Band 10- Или более сложный: (Band 11 * Band 11) / (Band 10 * Band 12)
- Интерпретация: Чем выше значение, тем «кислее» порода (больше SiO2). Позволяет отделить граниты от диоритов и найти кварцевые жилы (если они достаточно широкие).
5. Борьба с врагом: Индексы растительности
Растительность — главный враг спектрального геолога. Хлорофилл имеет поглощение в красном и мощное отражение в NIR. Глины (из-за воды в структуре) и растительность (из-за воды в клетках) могут путаться.
- NDVI (Normalized Difference Vegetation Index):
Formula: (NIR — Red) / (NIR + Red)- Применение: Мы не используем NDVI для поиска руды. Мы используем его для маскирования.
- Правило: Перед расчетом любого геологического индекса мы создаем маску: IF NDVI > 0.3 THEN NULL. Мы просто стираем лес с карты, чтобы не анализировать ложные аномалии.
6. Синтез: Создание композитов (RGB Composite)
Смотреть на черно-белую карту индекса «Каолинит» скучно и неинформативно. Геологу нужно видеть соотношение минералов. Мы создаем цветные композиты (False Color Composites), где каждому цветовому каналу (Red, Green, Blue) назначается свой минеральный индекс.
6.1. Классический «Alteration Composite» (GeoJet Standard)
Этот композит позволяет мгновенно оценить зональность месторождения.
- Red (Красный канал): OHI (Индекс глин / Al-OH). (Формула: (B5+B7)/B6).
- Смысл: Красным будут светиться зоны аргиллизитов (Каолинит/Алунит).
- Green (Зеленый канал): MGI (Индекс хлорита / Mg-OH). (Формула: (B6+B9)/B8).
- Смысл: Зеленым будут светиться пропилиты и карбонаты.
- Blue (Синий канал): Iron Oxide (Индекс железа). (Формула: B2/B1 для ASTER или аналог).
- Смысл: Синим будут светиться зоны окисления.
Интерпретация цветов на карте:
- 🟥 Чистый Красный: Интенсивная аргиллизация без окисления (глубокие части, не выветрелые).
- 🟨 Желтый (Красный + Зеленый): Смешанная зона, филлиты (серицит + хлорит).
- 🟪 Маджента (Красный + Синий): Аргиллизация + Окисление. Это Главная поисковая цель (Target Priority 1). Это окисленная «шляпа» над рудным телом с метасоматическим ореолом.
- 🟦 Синий: Просто оксиды железа (возможно, латеритная кора выветривания без гидротермалики).
- ⬛ Черный: Неизмененные породы или вода.
7. Продвинутые техники: Directed Band Ratioing
Простое деление каналов иногда дает сбой. Мы используем Нормализованные разностные индексы, которые работают стабильнее в условиях переменной освещенности.
- Relative Absorption Band Depth (RBD): Вместо простого отношения мы вычисляем «глубину» спектральной ямы относительно континуума (линии, соединяющей плечи спектра).
Formula: (Band_Left + Band_Right) / (2 * Band_Center) Это более надежный способ выделить именно минерал, отсекая влияние влажности почвы.
8. Практический кейс: Золоторудное поле (Чукотка)
Вводные: Поиск эпитермального золота в вулканическом поясе. Растительность скудная (тундра). Данные: ASTER L1T.
Алгоритм действий:
- Кросс-калибровка: Коррекция Crosstalk (известная проблема сенсора ASTER SWIR).
- Расчет маски: NDVI > 0.2 удален (кустарники в долинах).
- Расчет индекса Каолинита (B5+B7)/B6. Получена карта яркости.
- Расчет индекса Кремнезема (TIR B11/B10).
- Наложение: Мы ищем области, где High Kaolinite пересекается с High Silica.
- Геологический смысл: Это зоны вторичных кварцитов (Vuggy Silica) с каолинитовым оторочкой. Классическая модель High Sulfidation Epithermal Gold.
- Результат: Выделено 4 аномалии. Полевая проверка подтвердила 3 из них (одна оказалась осыпью кислых туфов). На одной из аномалий содержание Au в борозде составило 4.5 г/т.
9. Заключение: Алгебра как инструмент открытия
Использование мультиспектральных индексов — это самый дешевый способ конвертировать терабайты бесплатных спутниковых данных в векторные карты перспективных площадей.
Однако, Индустриальный Партнер должен помнить: Индекс — это не диагноз, это симптом. Высокий индекс оксидов железа может означать ржавую крышу ангара. Высокий индекс глин может означать распаханное поле. Индексы работают только в комплексе:
- В сочетании друг с другом (RGB).
- С учетом геологического контекста (структуры, разломы).
- С обязательной полевой заверкой (Ground Truthing).
В системе VEDART RS мы используем индексы как «первичное сито», отсеивающее 90% пустой породы, чтобы геолог мог сосредоточить свое внимание (и бюджет бурения) на оставшихся 10% аномалий.
FAQ: Ответы спектрального инженера
Вопрос 1: Можно ли рассчитать индексы в Excel? Ответ: Теоретически да, если вы экспортируете значения пикселей в таблицу. Но практически мы используем растровую алгебру в ГИС-пакетах (QGIS, ArcGIS, ENVI). Формула применяется ко всем 40 миллионам пикселей снимка одновременно.
Вопрос 2: Какой индекс лучше для поиска меди? Ответ: Для меди (порфировой) критически важен индекс серицита (Muscovite Index) и хлорита. Медь «сидит» в калиевой зоне, окруженной филлитовой (серицит) и пропилитовой (хлорит) зонами. Мы ищем «бычий глаз» (Bullseye): центр с низким индексом хлорита, окруженный кольцом высокого индекса хлорита.
Вопрос 3: Почему мои индексы показывают «шум»? Ответ: 99% проблем — в отсутствии атмосферной коррекции. Если вы делите сырые данные (DN values) без перевода в отражательную способность (Reflectance), вы делите «попугаев на попугаев». Атмосфера искажает синий канал и SWIR каналы по-разному. Сначала сделайте FLAASH/DOS коррекцию, потом считайте индексы.
