Обобщение вычисления производных на основе β-VDR для устойчивого обнаружения краев источников и оценки глубины по данным потенциального поля

Видеть скрытые структуры под нашими ногами

Минералы, грунтовые воды, геотермальные ресурсы и нефть часто залегают глубоко под землей, скрытые от прямого наблюдения. Геофизики используют тонкие изменения в гравитационном и магнитном полях Земли, чтобы картировать эти погребённые структуры без бурения. Но преобразовать деликатные измерения в чёткие изображения разломов, дайков и контактов непросто: обычные методы обработки склонны усиливать шум вместе с сигналом. В этом исследовании предложен способ, который позволяет «затачивать» подземные изображения, удерживая шум под контролем, что делает карты недр более надёжными для науки и разведки.

Почему важны подземные края

Когда породы разных типов или плотностей соприкасаются, на глубине образуются «края» — разломы, контакты и интрузии, которые часто контролируют, где аккумулируются жидкости, тепло и руды. Гравитационные и магнитные съёмки косвенно улавливают эти границы как небольшие вариации, называемые аномалиями потенциального поля. Чтобы точно определить положение и глубину источников, интерпретаторы вычисляют математические производные данных, которые выделяют участки с резкими изменениями поля. К сожалению, эти производные действуют как фильтры, усиливающие высокочастотную составляющую, поэтому даже небольшие уровни случайного шума могут заглушить интересующие особенности. Существующие исправления либо работают только при очень чистых данных, либо требуют тяжёлых вычислений, которые трудно применить к большим современным съёмкам.

Более умный способ брать разности

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

Уравновешивание вертикального и бокового обзора

Авторы переформулировали идею β-VDR в компактный фильтр в частотной области, который достигает того же эффекта с помощью всего одного прямого и одного обратного преобразования Фурье вместо пяти. Этот шаг сам по себе сокращает теоретическое время вычислений примерно в пять раз. Затем они распространили ту же стабилизирующую логику на горизонтальные производные, создав сопоставимую семью фильтров, которую назвали β-HDR. Вместе вертикальный β-VDR и горизонтальный β-HDR образуют единый подход — β-VDR-with-β-HDR, который последовательно обрабатывает производные во всех направлениях. Проще говоря, метод аккуратно сглаживает шум во всех направлениях, сохраняя при этом резкие переходы, обозначающие реальные геологические границы.

Испытание метода

Чтобы проверить корректность и полезность нового подхода, команда провела обширные компьютерные эксперименты. Они начали с синтетических моделей — идеализированных подземных блоков с известными формами, глубинами и физическими свойствами — и сгенерировали их гравитационные и магнитные отклики. Добавляя различные уровни случайного шума, они воссоздали типичные «грязные» данные полевых съёмок. Используя стандартную технику усиления краёв, называемую полной градиентной мерой (total gradient), зависящую от вертикальных и горизонтальных производных, они сравнили четыре варианта: традиционные фильтры на основе Фурье, метод ISVD, оригинальный β-VDR в сочетании с традиционными горизонтальными разностями и новый β-VDR-with-β-HDR. Новый метод воспроизводил стандартные результаты при отсутствии стабилизации, что подтвердило математическую корректность. В условиях шума он явно выделился: края оставались чёткими, ложные пики были редки, а оценённые глубины оставались близкими к истинным значениям даже там, где другие методы давали сбои.

От тестовых моделей к реальному осадочному бассейну

Затем авторы применили свою технику к высокоразрешающим аэромагнитным данным нигерийской части бассейна Чад — территории с толстыми отложениями, где разломы и интрузии влияют на геотермальный и углеводородный потенциал. Не прибегая к обычному предварительному сглаживанию, они рассчитали стабилизированные производные и полную градиентную меру, а затем оценили положения и глубины магнитных источников как по профилям, так и в полном 3D. Решения согласовались с известными региональными тенденциями и выявили когерентные разломоподобные структуры и интрузии, включая мелкие структуры и глубокие линейности, которые могут направлять поток флюидов. Существенно, оценки глубин по 2D-профилям и 3D-сеткам тесно совпали, что указывает на то, что результаты не являются артефактами метода.

Более ясные подземные изображения для сложных данных

Для неспециалиста основной вывод таков: данная работа предлагает улучшенный «фильтр повышения резкости» для шумных гравитационных и магнитных карт недр. Переработав способ вычисления вертикальных и горизонтальных разностей, метод β-VDR-with-β-HDR надёжнее извлекает края и глубины погребённых структур, даже когда измерения сильно зашумлены. Поскольку он также вычислительно более эффективен, его можно применять к большим современным наборам данных. Это означает более чёткие и более достоверные изображения того, что скрыто под нашими ногами — что поддерживает более безопасные решения при бурении, улучшенные оценки геотермального потенциала и более глубокое понимание подземной архитектуры Земли.

Цитирование: Falade, S.C., Falade, A.H. Generalizing β-VDR-based derivative computation for robust source edge detection and depth estimation from potential field data. Sci Rep 16, 5672 (2026). https://doi.org/10.1038/s41598-026-36635-7

Ключевые слова: гравиметрические и магнитные съемки, обнаружение краев, оценка глубины, устойчивые к шуму производные, геология бассейна Чад