P. Núñez Demarco, C. Prezzi, L. Sánchez Bettucci
2 017
Latinmag Letters, Volume 7, Special Issue (2017), GEP04, 1-6. Proceedings Juriquilla, Qro, Mexico
Se desarrolló un programa basado en MatLab que computa automáticamente el análisis espectral en 2-D de mapas magnéticos y estima la profundidad del punto de Curie, el gradiente geotérmico y el flujo térmico, según los métodos y ecuaciones establecidos por Spector y Grant (1970), Bhattacharyya y Leu (1977), Okubo et al. (1985) y Tanaka et al. (1999). El programa carga el mapa magnético, interpola los datos faltantes si es necesario, y divide el mapa en subregiones de acuerdo con el tamaño establecido por el usuario. Además, diferentes niveles de solapamiento entre las subregiones (de entre 0% y 99.9%) pueden ser evaluados simultáneamente, según las preferencias del usuario. La potencia espectral de la transformada rápida de Fourier en dos dimensiones (FFT 2-D) es calculada para cada subregión. La profundidad del tope y del centroide de la fuente magnética es derivada de la media radial del espectro de potencia escalado por la frecuencia y de la media radial del espectro de potencia, respectivamente. La profundidad de la isoterma de Curie es obtenida a partir de los valores de profundidad del tope y centroide y el gradiente geotérmico calculado. Todos los datos son graficados. Asimismo, se visualiza información sobre el procedimiento que permite su validación. Esta información incluye parámetros estadísticos, coeficientes de correlación y una visualización rápida de los ajustes lineales. Este programa fue testeado y comparado con otros software como ser Oasis Montaj, obteniendo un aumento en resolución y una reducción substancial en tiempo de procesamiento.
Palabras Clave: Análisis Espectral, Espectro de Potencia, Profundidad de la Isoterma de Curie,
Gradiente Geotérmico.
We developed a new MatLab based program that automatically computes the 2-D spectral analysis of magnetic maps, estimates the Curie point depth, the geothermal gradient and heat flux according to the methods established by Spector and Grant (1970), Bhattacharyya and Leu (1977), Okubo et al. (1985) and Tanaka et al. (1999). The program loads the magnetic data, interpolates the missing data if required, and divides the map in subregions according to the size established by the user. Also, different levels of overlapping between subregions (ranging from 0% to 99.9%) are tested simultaneously, according to the user customization. The 2-D Fast Fourier Transform (FFT) power spectrum is calculated for each subregion. The centroid and top depth of magnetic sources are derived from the slopes of the longest and second-longest wavelength of the radial averaged frequency-scaled power spectra and the radial averaged power spectra, respectively. The Curie point depth is derived from top and centroid depths (Okubo et al., 1985) and the geothermal gradient calculated. All the data are plotted in different visualization forms. In addition, information about the procedure is visualized in order to validate it. This information includes statistical parameters, correlation coefficients and fast visualization of the linear adjustments. This program was tested and compared with the results of other software like Oasis Montaj, obtaining an improvement in resolution and substantial reduction in processing time.
Keywords: Spectral Analysis, Power Spectrum, Curie Point Depth, Geothermal Gradient

http://www.geofisica.unam.mx/LatinmagLetters/LL17-01-SP/GEP/GEP04.pdf