Generalización del cálculo de derivadas basado en β-VDR para una detección de bordes de fuente y estimación de profundidad robustas a partir de datos de campo potencial

Ver las estructuras ocultas bajo nuestros pies

Los minerales, las aguas subterráneas, los recursos geotérmicos y el petróleo suelen encontrarse en profundidad, ocultos a la vista directa. Los geofísicos usan cambios sutiles en los campos de gravedad y magnéticos de la Tierra para cartografiar estas estructuras enterradas sin perforar. Pero convertir esas mediciones delicadas en imágenes nítidas de fallas, diques y contactos es difícil, porque los métodos de procesamiento habituales tienden a amplificar el ruido junto con la señal. Este estudio presenta una forma de afinar las imágenes del subsuelo manteniendo el ruido bajo control, lo que hace que los mapas subterráneos sean más fiables para la ciencia y la exploración.

Por qué importan los bordes en el subsuelo

Cuando rocas de distintos tipos o densidades se encuentran, forman «bordes» en el subsuelo: fallas, contactos e intrusiones que suelen controlar dónde se acumulan fluidos, calor y minerales. Los levantamientos gravimétricos y magnéticos pueden detectar indirectamente estas fronteras como pequeñas variaciones, llamadas anomalías de campo potencial. Para localizar con precisión las posiciones y profundidades de las fuentes, los intérpretes calculan derivadas matemáticas de los datos, que destacan donde el campo cambia más bruscamente. Desafortunadamente, estas derivadas actúan como filtros que realzan el contenido de alta frecuencia, por lo que incluso pequeñas cantidades de ruido aleatorio pueden ocultar las características de interés. Las soluciones existentes o bien funcionan solo cuando los datos están muy limpios, o requieren cómputos pesados que son difíciles de aplicar a los levantamientos actuales de gran tamaño.

Una forma más inteligente de tomar diferencias

Un método anterior, conocido como β-VDR, ya ofrecía una manera más estable de calcular derivadas verticales combinando hábilmente versiones de los datos que habían sido continuadas hacia arriba—proyectadas matemáticamente a elevaciones mayores para suavizar el ruido. β-VDR proporciona derivadas verticales más limpias que los filtros estándar basados en Fourier, pero presentaba dos inconvenientes importantes. Primero, seguía apoyándose en fórmulas más frágiles de diferencias finitas para los componentes horizontales, lo que generaba un desequilibrio: las derivadas verticales eran robustas y las horizontales no. Segundo, la receta original requería cinco rondas separadas de pesados cálculos de Fourier, lo que la hacía lenta y costosa para mallas grandes.

Equilibrando la vista vertical y la lateral

Los autores reformularon la idea de β-VDR en un filtro compacto en el dominio de frecuencia que logra el mismo efecto con solo una transformada de Fourier directa y una inversa en vez de cinco. Este paso por sí solo reduce el tiempo de cómputo teórico en un factor cercano a cinco. Luego extendieron la misma lógica estabilizadora a las derivadas horizontales, creando una familia de filtros emparejados que denominan β-HDR. En conjunto, la β-VDR vertical y la β-HDR horizontal forman un esquema unificado, β-VDR-con-β-HDR, que trata todas las direcciones de derivada de forma coherente. En términos simples, el método suaviza el ruido lo justo en cada dirección al mismo tiempo que preserva las transiciones nítidas que marcan límites geológicos reales.

Poniendo el método a prueba

Para comprobar que el nuevo enfoque es correcto y útil, el equipo realizó extensos experimentos por ordenador. Comenzaron con modelos sintéticos—bloques subterráneos idealizados con formas, profundidades y propiedades físicas conocidas—y generaron sus respuestas gravimétricas y magnéticas. Al añadir distintos niveles de ruido aleatorio, recrearon el tipo de datos desordenados que se encuentran en levantamientos reales. Usando una técnica estándar de realce de bordes llamada gradiente total, que depende de las derivadas verticales y horizontales, compararon cuatro opciones: filtros tradicionales basados en Fourier, un método llamado ISVD, la β-VDR original combinada con diferencias horizontales convencionales, y la nueva β-VDR-con-β-HDR. El nuevo método reprodujo los resultados estándar cuando no se aplicó estabilización, confirmando la matemática. En condiciones ruidosas destacó claramente: los bordes permanecieron nítidos, los picos espurios fueron raros y las profundidades estimadas se mantuvieron cerca de los valores reales incluso cuando otros métodos fallaban.

De modelos de prueba a una cuenca sedimentaria real

Los autores aplicaron luego su técnica a datos aeromagnéticos de alta resolución del sector nigeriano de la cuenca del Chad, una región de sedimentos gruesos donde fallas e intrusiones influyen en el potencial geotérmico e hidrocarburífero. Sin recurrir al suavizado preliminar habitual, calcularon las derivadas estabilizadas y el gradiente total, y luego estimaron las posiciones y profundidades de las fuentes magnéticas tanto en perfiles como en vistas 3D completas. Las soluciones se alinearon con tendencias regionales conocidas y revelaron rasgos coherentes tipo falla e intrusiones, incluidas estructuras someras y lineamientos más profundos que pueden guiar el flujo de fluidos. Es importante destacar que las estimaciones de profundidad obtenidas desde perfiles 2D y mallas 3D concordaron de cerca, lo que sugiere que los resultados no son artefactos del método.

Imágenes subterráneas más claras para datos difíciles

Para un lector no especialista, la conclusión es que este trabajo ofrece un “filtro de enfoque” mejorado para mapas gravimétricos y magnéticos ruidosos del subsuelo. Al rediseñar cómo se calculan las diferencias verticales y horizontales, el método β-VDR-con-β-HDR extrae bordes y profundidades de estructuras enterradas con mayor fiabilidad, incluso cuando las mediciones están contaminadas por ruido significativo. Además, como es más eficiente de calcular, puede aplicarse a grandes conjuntos de datos modernos. Esto se traduce en imágenes más claras y de mayor confianza de lo que hay bajo nuestros pies—apoyando decisiones de perforación más seguras, mejores evaluaciones geotérmicas y una comprensión más profunda de la arquitectura oculta de la Tierra.

Cita: 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

Palabras clave: levantamientos gravimétricos y magnéticos, detección de bordes, estimación de profundidad, derivadas robustas frente al ruido, geología de la cuenca del Chad