Análisis gravimétrico y aeromagnético de las principales intrusiones en la carta El Batamote, Sonora Central, México: implicaciones en la exploración minera

Abstract

La carta El Batamote se ubica a 25 km en línea recta al noroeste de Hermosillo, Sonora, México, cubriendo un área de 899 km2. En la geología superficial se tienen aflorando rocas desde el Neoproterozoico al Holoceno, encontrándose aproximadamente la mitad de la carta cubierta por sedimentos de las rocas preexistentes. Con el fin de identificar las principales estructuras geológicas que se encuentran en el subsuelo y determinar su posible relación con las zonas mineralizadas presentes en la carta El Batamote, se realizó un estudio gravimétrico con datos del Instituto Nacional de Estadística y Geografía (INEGI) y de perfiles gravimétricos medidos en este trabajo, así como un análisis a partir de la aeromagnetometría del Servicio Geológico Mexicano (SGM). El procesamiento de los datos gravimétricos consistió en realizar la separación regional-residual por el método de Zeng (1989), posteriormente se aplicaron diferentes filtros a la anomalía de Bouguer con el fin de enfatizar las anomalías de interés. A partir de la anomalía residual se realizó la modelación bidimensional de perfiles perpendiculares a las anomalías, donde los modelos se ajustan a la morfología del Basin and Range. Las secuencias que se encuentran rellenando los valles presentan un espesor máximo de 1100 m. Las profundidades obtenidas mediante el método de deconvolución de Euler se correlacionaron con las obtenidas en la etapa de modelado. Del análisis de lineamientos aeromagnéticos se obtuvo que los principales lineamientos tienen una orientación principal NW-SE, donde el filtro que arrojó una mejor definición de éstos fue la derivada de Tilt (Miller y Singh, 1994). Para el análisis de estructuras profundas se emplearon continuaciones ascendentes, reflejando la presencia de dos altos magnéticos a profundidad en los laterales de la carta. Las profundidades promedio a la cima de los cuerpos obtenidas mediante el método de deconvolución de Euler se encuentran en el rango de 500 a 2000 m, las cuales se correlacionan con las obtenidas mediante el análisis espectral. A partir de este último se realizó la separación entre la anomalía magnética regional y residual. Una vez obtenida la anomalía residual se procedió a la modelación tridimensional (3D) de las anomalías de interés, donde el mejor ajuste se obtuvo con la geometría de diques verticales con profundidades a la cima de 500 a 1000 m y espesores promedio de 250 m.