Методология каскадного дистанционного зондирования на основе спутниковой рекогносцировки и гетерогенного роя БПЛА для геологоразведочных работ

Обоснованность выбора площадей для детализации наземных геофизических работ и поискового бурения во многом определяет эффективность геологоразведочных работ на поисковой стадии. Традиционные схемы поиска ограничены либо низкой производительностью наземных методов на обширных территориях, либо зависимостью оптической съемки от облачности, либо хаотичностью автономных облетов беспилотных летательных аппаратов (БПЛА). В работе предложена методология каскадного дистанционного зондирования, интегрирующая спутниковую рекогносцировку и целенаправленную верификацию аномалий гетерогенным роем БПЛА. На первом этапе по радиолокационным (SAR, InSAR) и мультиспектральным (NDVI, NDSI) данным формируется цифровая карта априорных вероятностей геологических аномалий; на втором этапе рой выполняет лидарную и мультиспектральную съемку в векторизованных полигонах приоритетного сканирования. Научная новизна состоит в архитектуре бортового контроля качества геологической информации на базе периферийного искусственного интеллекта и в алгоритме жесткого приоритета лидерства, обеспечивающем непрерывность миссии при отказах аппаратов. Имитационное моделирование показало, что интеграция данных InSAR сужает зону поискового бурения на 35 %, бортовой агент отсеивает до 38 % кадров с артефактами, а при отказе лидера сохраняется 100 % маршрутных точек в зонах аномалий при времени передачи лидерства около 5 с. Разработанная морфологическая матрица обосновывает применимость и риски подсистем комплекса. Практическая значимость результатов заключается в сокращении объемов «слепого» бурения и трудоемкости ручной постобработки материалов съемки.

На поисковой стадии геологоразведочных работ (ГРР) именно точность выбора площадей под детализацию наземной геофизикой и поисковым бурением определяет экономику всего последующего цикла. Обоснованному выбору препятствуют три взаимосвязанных ограничения: низкая производительность точных наземных методов на обширных труднодоступных территориях; зависимость оптической спутниковой съемки от облачности и сезонного снега; избыточное сканирование однородных участков беспилотными летательными аппаратами (БПЛА), способными обеспечивать сантиметровое разрешение, но применяемыми без априорной геологической привязки.
Отдельные компоненты будущего интегрированного решения уже апробированы: низковысотная аэромагниторазведка на базе БПЛА при структурном картировании и прогнозировании рудных объектов, роевая робототехника и планирование групповых траекторий, машинное обучение в прогнозно-поисковых построениях, литологическое картирование по гипер- и мультиспектральным данным, независимая от облачности интерферометрия радаров с синтезированной апертурой, или InSAR-интерферометрия. Однако спутниковый и беспилотный уровни сбора геопространственных данных рассматриваются, как правило, раздельно: маршруты БПЛА планируются равномерными сетками без учета априорной геологической информации, контроль качества материалов вынесен на камеральную постобработку, а замкнутый контур «спутниковая рекогносцировка — целенаправленная верификация роем» с бортовым принятием решений практически не исследован.
Чтобы закрыть обозначенный пробел, целью работы поставлена разработка методологии каскадного дистанционного зондирования, в которой карта априорных вероятностей аномалий, построенная по спутниковым данным, непосредственно определяет геометрию маршрутов гетерогенного роя БПЛА с бортовым контролем качества геологической информации.
Каскад открывается спутниковой ступенью: рекогносцировка выполнена по архивным и оперативным данным Sentinel-1 (режим IW, C-диапазон) за 2021–2025 гг., полученным в том числе через центр коллективного пользования «ИКИ-Мониторинг», и мультиспектральным данным Sentinel-2 (прибор MSI, разрешение 10 м). Радиолокационная обработка включает ко-регистрацию сцен, построение серий интерферограмм, оценивание скоростей смещений поверхности методами PS/SBAS и выделение линеаментов; мультиспектральная обработка — расчет индекса NDVI как биоиндикатора скрытых литологических и гидротермальных процессов и индекса NDSI для алгоритмического маскирования сезонного снега. Производные слои интегрируются взвешенным суммированием с весами, определенными по ретроспективному анализу эталонных участков, результат нормируется к вероятностной шкале от нуля до единицы, а пороговое отсечение и векторизация формируют полигоны приоритетного сканирования, загружаемые в бортовые вычислительные модули роя до начала миссии.

Верификацию аномалий выполняет гетерогенный рой трех ярусов: аппарат-ретранслятор с гибридной силовой установкой (эшелон 1000–1500 м), создающий воздушный mesh-хаб с дальностью до 100 км; носители лидара и мультиспектральной камеры (300–500 м); рабочие аппараты верификации и повторных пролетов (50–100 м). Узкополосный канал 900 МГц передает только отчеты контроля качества и статусы верификации, тогда как массивы лидарных и мультиспектральных данных записываются на борт. При подтверждении аномалии в целевом полигоне рой получает декларативную команду локальной концентрации ресурсов и увеличивает плотность съемки сектора без вмешательства наземного оператора.

Контроль качества геологической информации перенесен на борт: сверточная нейронная сеть семейства YOLO (англ. You Only Look Once), дообученная на отсев помех съемки (блики от солнца, влажных и заснеженных поверхностей, тени облачности, артефакты дымки), сегментирует каждый кадр и принимает решение «принять, отклонить, повторный пролет»; дополнительно в целевом полигоне контролируется плотность лидарного облака точек. При спектральной несогласованности либо недостаточной плотности облака аппарат автономно инициирует повторный пролет, не загружая канал связи нерелевантными данными. Непрерывность миссии обеспечивает алгоритм жесткого приоритета лидерства: при потере телеметрических «heartbeat-сигналов» ведущего аппарата его функции без участия наземного оператора передаются следующему по приоритету аппарату в нормативный интервал, а топология маршрутов адаптируется к ветровой нагрузке и рельефу.

Эффективность каскадной схемы оценена дискретно-событийным имитационным моделированием миссий на модельном полигоне, воспроизводящем залесенные ландшафты со сложным микрорельефом, развитой гидрографической сетью и сезонным снеговым покровом. Каждый сценарий воспроизведен в 30 независимых прогонах с варьированием ветровой нагрузки, облачности и отказов сенсоров; значимость различий с базовой схемой равномерного облета проверена критерием Уилкоксона при уровне значимости 0,05. Применимость и риски подсистем систематизированы морфологической матрицей (метод Цвикки).

Результатом спутниковой ступени стала цифровая карта априорных вероятностей аномалий, на которой устойчиво выделяются четыре класса объектов: карстово-суффозионные депрессии с риском просадки грунтов; компактные зоны спектральных аномалий, интерпретируемые как потенциальные рудные тела; линеаменты тектонических разломов по данным InSAR; протяженные зоны стресса растительности. На основе карты система генерирует не равномерную сетку облета, а векторизованные полигоны приоритетного сканирования, переводя рой в режим целенаправленной верификации.
В сериях прогонов бортовой агент корректно идентифицировал и отфильтровал в среднем 38 % кадров, содержащих блики и теневые артефакты, которые в стандартном камеральном конвейере потребовали бы трудоемкой ручной постобработки; ложных срабатываний на геологически значимые спектральные подписи не зафиксировано. Моделирование отказа ведущего аппарата подтвердило работоспособность алгоритма жесткого приоритета: время передачи лидерства не превысило 5 с, что меньше цикла актуализации данных, при этом сохранением маршрутных точек в зонах выявленных аномалий; деградация узкополосного канала компенсировалась перестроением mesh-топологии через ретранслятор без потери управляемости роя. Детализация процедуры: отказ лидирующего аппарата детектируется по пропаданию «heartbeat-сигналов» в течение двух циклов телеметрии, после чего его функции переходят к следующему по закрепленному приоритету аппарату, а маршрутные точки отказавшего аппарата перераспределяются между оставшимися без повторного планирования миссии.
Интеграция данных InSAR позволила сузить зону поискового бурения на 35 % за счет сосредоточения наземных работ исключительно в полигонах, верифицированных роем. Экономия летного ресурса относительно равномерной сетки облета достигала 40 %, объем временно хранимых на борту данных снизился на 30–40 %; различия схем по всем регистрируемым показателям статистически значимы при уровне значимости 0,05, без включения z-аппроксимации. К ограничениям исследования относятся имитационный характер оценки и идеализация поведения сенсоров; натурная апробация комплекса составляет предмет дальнейших работ. Использование воздушного пространства и эшелоны выше 150 м, требуемые для функционирования ретранслятора, подлежат согласованию в установленном порядке, что закладывается в план миссии при ее подготовке. Предложенная методология не привязана к конкретному типу ландшафта: для залесенных и открытых территорий изменяются лишь состав сенсоров и параметры полетных эшелонов, тогда как архитектура каскада и протоколы роевого взаимодействия сохраняются. Замкнутый контур «спутниковая рекогносцировка — карта априорных вероятностей — адресная верификация роем с бортовым контролем качества» обеспечивает статистически значимое сокращение площади «слепого» поиска и трудоемкости постобработки, создавая основу для перехода от экстенсивного сканирования территорий к адресным геологоразведочным работам.