Grav Mag Prism: A Matlab /Octave Program to Generate Gravity and Magnetic Anomalies Due to Rectangular Prismatic Bodies

Abstract

ABSTRACT. In this paper we present the GRAV MAG PRISM code to generate synthetic gravity and magnetic anomalies from rectangular prismatic bodies with arbitrary dimensions, densities and magnetizations. The code has been developed in a MATLAB/OCTAVE environment and provides a graphical and a numerical output, as well. The data are written in xyz (ASCII) format and can be contaminated with noise. We also summarize the theory and show some examples of the program’s applications. The code can be used for research and educational purposes and is freely-distributed by the authors on request.Keywords: gravity and magnetic anomalies, potential-field modeling, program. RESUMO. Neste trabalho é apresentado o programa GRAV MAG PRISM, o qual gera anomalias de gravidade e magnéticas sintéticas a partir de corpos prismáticos retangulares com dimensões, densidades e magnetizações arbitrárias. O código foi desenvolvido em ambiente MATLAB/OCTAVE, com saídas gráficas e numéricas. Os dados são escritos no formato xyz (ASCII) e podem ser contaminados por ruídos. Neste trabalho também é apresentado um resumo da teoria e alguns exemplos de aplicação do programa. O código acompanha um tutorial e pode ser utilizado para fins educacionais e de pesquisa, sendo de distribuição gratuita através de solicitação aos autores.Palavras-chave: anomalias gravimétricas e magnéticas, modelagem de campos potenciais, programa.

GRAV-MAG SUITE - An open source MATLAB-based program for processing potential field data

“…The magnetic anomaly shown in Fig. 3 was generated using GRAV MAG PRISM, a MATLAB-based forward modeling program for both magnetic and gravity fields generated by prismatic bodies (Bongiolo et al 2013).…”

Optimizing depth estimates from magnetic anomalies using spatial analysis tools

“…To illustrate our methodology, we calculated the magnetic anomalies over two hypothetical magnetic sources using Matlab program developed by Bongiolo et al (2013). Model A used geomagnetic inclination values of 90° and included two prisms (M1 and M2) with equal dimensions, equal magnetization (0.13 A/m), different depths, and sides oriented north-south or east-west ( Figure 1).…”

O lineamento transbrasiliano na Bacia do Paraná : compartimentação crustal do embasamento e reativações fanerozóicas

“…Figura 3.15 Aproximação para um corpo bidimensional por um polígono de N lados Figura 4.A1 Airborne geophysical surveys compiled for this study Figura 4.1 (a) Synthetic Model A, consisting of two magnetic prisms, M1 and M2, with same dimensions, infinite depth extent, and tops at depths of 3 and 5 km, respectively. Anomalies calculated with the method of Bongiolo et al (2013) Estimated depth values of magnetic anomalies observed from the Bacia do Paraná airborne survey, using tilt-depth technique (Salem et al, 2007), distance analysis tool and spline interpolation. Mapped Cretaceous units of Paraná Basin are represented in green , including volcanic sills………… 85 Basin generalized from CPRM (Serviço Geológico do Brasil) 1:1,000,000 scale mapping and classified according to geological periods.…”

“…A maior parte das técnicas de realce foi analisada e avaliada por meio de testes efetuados no programa desenvolvido em ambiente Matlab, GRAV_MAG_PRISM (Bongiolo et al, 2013), no qual dados magnéticos e gravimétricos sintéticos foram gerados em malha de dados de 250 x 250 m. O objetivo dos testes é entender a geometria e intensidade das anomalias em diferentes situações do campo potencial terrestre e avaliar empiricamente o desempenho das técnicas nessas diferentes condições. As fontes causativas correspondem a prismas com variáveis…”

“…Many methods for edge detection and depth-to-source estimation rely on horizontal and vertical derivatives of magnetic anomalies Blakely and Simpson, 1986;Thurston and Smith, 1997;Fedi and Fiori, 2000;Phillips, 2000;Ferreira et al, 2013). The total-gradient method, for example, calculates the analytic signal by combining horizontal and vertical derivatives of magnetic anomalies, calculated either from profile or gridded data.…”