UNIVERSIDADE DE BRASÍLIA /INSTITUTO DE GEOCIÊNCIAS

TESE DE DOUTORADO No 30

MARCELO DE LAWRENCE BASSAY BLUM

PROCESSAMENTO E INTERPRETAÇÃO DE DADOS DE GEOFÍSICA AÉREA NO BRASIL CENTRAL E SUA APLICAÇÃO À GEOLOGIA REGIONAL E À PROSPECÇÃO MINERAL

Acessar TEXTO INTEGRAL

Palavras-chave:Geofísica Aplicada, Gravimetria, Magnetometria Aérea, Gamaespectrometria Aérea, Geologia Regional, Prospecção Mineral

DATA DA DEFESA:03/09/99
ÁREA DE CONCENTRAÇÃO: GEOLOGIA ECONÔMICA E PROSPECÇÃO
ORIENTADOR: PROF. ROBERTO ALEXANDRE V. DE MORAES (UnB)
CO-ORIENTADOR: PROF. AUGUSTO CESAR BITTENCOURT PIRES (UnB)
EXAMINADORES:PROF. REINHARDT A FUCK (UnB)
PROF. VASILE MARZA (UnB-SISMOLOGIA)
PROF. LUIZ FERNANDO S. BRAGA (GEOMAG)
PROF. CARLOS ALBERTO MENDONÇA (USP)

RESUMO
Com o fito de melhorar o conhecimento da geologia regional, foram utilizados os dados de geofísica aérea do Projeto Geofísico Brasil – Canadá (PGBC). Esses dados foram avaliados, tratados e interpretados à luz do conhecimento geológico da área estudada localizada no Estado de Goiás e no sudoeste do Estado do Tocantins, abrangendo terrenos da Província Estrutural Tocantins.
Na avaliação dos dados do PGBC foi proposta uma metodologia para a determinação de “picos” (spikes) na gravação original dos dados magnéticos, que revelou ser uma técnica eficiente. A partir daí, os dados foram interpolados em malha regular utilizando o método de splines bi-cúbicos. Os dados foram micronivelados utilizando um algoritmo desenvolvido, que mostrou eficiência na retirada dos desnivelamentos entre linhas de vôo.
Grande parte das filtragens aplicadas sobre os dados aeromagnéticos foi feita no domínio do número de onda. O filtro de redução ao polo somente foi eficiente para aquelas anomalias que provavelmente não possuíam magnetização remanescente. Mesmo assim, o resultado da aplicação do filtro de pseudo-gravidade, baseado nessa redução ao polo, pôde ser comparado ao mapa Bouguer da região.
A aplicação de técnicas de realce de fonte magnética (gradiente horizontal total, amplitudes e fases dos sinais analíticos de primeira e segunda ordens e terraceamento), mostrou eficiência na determinação de descontinuidades físicas que podem ser traduzidas em unidades geológicas. A amplitude do sinal analítico de segunda ordem mostrou melhores resultados do que de primeira ordem por revelar melhor os limites dos corpos causativos e discriminar as unidades de rochas máficas e ultramáficas.
A aplicação de terraceamento sobre a amplitude do sinal analítico de segunda ordem permitiu situar a resposta anômala sobre o traçado dos afloramentos em mapa, definindo melhor as unidades magnéticas em relação à geologia conhecida.
A fase do sinal analítico exibiu toda a complexidade das estruturas locais traduzidas em um emaranhado de estruturas lineares (lineamentos e limites de propriedade magnética) que auxiliaram na interpretação e correlação às conhecidas da geologia.
Na estimativa da profundidade da superfície Curie da área de estudo foram desenvolvidos programas específicos, acrescidos de rotinas que o tornaram mais eficiente e que possibilitou cobrir uma área equivalente à do trabalho sem comprometer o resultado. Nessas rotinas incluiu-se a retirada automática de superfície de tendência de grau 1 e a expansão das bordas por dobramento das janelas de dados no domínio de Fourier, garantindo maior estabilidade nas estimativas espectrais. As estimativas em foco, juntamente com alguns dados geotérmicos pontuais, permitiram que se extrapolasse e interpolasse esses últimos para toda a área de estudo. Isto mostra a importância das estimativas das profundidades Curie.
Da análise das imagens magnéticas do campo anômalo e de suas transformadas, foi possível definir quatro tipos principais de relevo magnético: relevo magnético suave, relevo magnético acidentado e incorrelato, relevo magnético acidentado e correlato e relevo magnético muito acidentado.
A análise das principais estruturas magnéticas interpretadas a partir, principalmente, das imagens da fase do sinal analítico, da amplitude do sinal analítico de segunda ordem e do campo magnético anômalo, aliado ao conhecimento geológico da região, sugere que essas estruturas estejam ligadas a um evento com tensão principal, ?1, na direção ~E-W.
O uso de composição do tipo CMY invertido nos mapas gamaespectrométricos ternários surtiu bom resultado. A melhor combinação foi a que usou o canal do K associado à cor ciano (C), o Th a magenta (M) e o U a amarelo (Y), que permitiu a identificação dos limites e subdivisões dos complexos ortognáissicos de Caiçara, da Anta e Uvá.
Para estudo dos greenstone belts locais, foi realizado nas região de Crixás um levantamento gravimétrico em semi-detalhe (cerca de 400 estações sobre 100.000 km2). A malha gravimétrica revelou os três greenstone belts com maior nitidez que a obtida dos dados regionalizados disponíveis. Metodologias de determinação indireta de densidades e de correção de terreno foram desenvolvidas e aplicadas na área estudada.
Modelagem 2-D (modelos de partida) e inversões 2,5-D e 3-D sobre os greenstone belts com dados gravimétricos e modelos baseados na geologia local, permitiram a elaboração de seções verticais com a geologia provável destas. As profundidades obtidas em dois métodos de inversão diferentes indicam que a base das anomalias pode chegar a 2 km de profundidade.
A análise integrada das imagens ternárias de K, Th e U, em conjunto com os mapas magnéticos transformados, os mapas Bouguer, o mapa de potássio anômalo e o mapa de integração magnético - gamaespectrométrica, à luz do conhecimento geológico da região, sugere que a região tenha sido afetada por uma seqüência de eventos, iniciando-se pela formação dos terrenos granitóide-greenstone e passando pela colisão brasiliana, indo até eventos mais recentes.
Usando os dados gamaespectrométricos, foi possível estimar o chamado potássio anômalo, que revelou alinhamentos coincidentes com algumas ocorrências de ouro primário. Essas mesmas ocorrências foram adicionadas aos mapas interpretativos. Observou-se que parte delas se associa a estruturas magnéticas, incluindo lineamentos na fase do sinal analítico, e a alguns contatos entre as unidades. Esse fato, aliado aos alinhamentos encontrados no potássio anômalo, podem indicar zonas de cisalhamento, principalmente na direção N55-65W, na qual haveria possibilidade da ocorrência de minerais hidrotermais.
Mapas de interpretação magnética, gamaespectrométrica e integrada da área de estudo e da área do levantamento gravimétrico representam esboços regionais de unidades geológicas e elementos estruturais. A validade dessas interpretações deverá ser verificada em campo.
 


  
 
UNIVERSITY OF BRASILIA / INSTITUTE OF GEOSCIENCES

PhD THESIS No 30

MARCELO DE LAWRENCE BASSAY BLUM

PROCESSING AND INTERPRETATION OF AIRBORNE GEOPHYSICAL DATA IN CENTRAL BRASIL AND ITS APPLICATION TO REGIONAL GEOLOGY STUDIES AND TO MINERAL EXPLORATION

KeyWords: Applied Geophysics, Gravity, Airborne Magnetometry, Airborne gamma spectrometry, Regional Geology, Mineral Prospection

DATE OF ORAL PRESENTATION:03/09/99
TOPIC OF THE THESIS: EXPLORATION AND ECONOMIC GEOLOGY
SUPERVISOR:PROF. ROBERTO ALEXANDRE V. DE MORAES (UnB)
CO-SUPERVISOR: PROF. AUGUSTO CESAR BITTENCOURT PIRES (UnB)
COMMITTEE MEMBERS:PROF. REINHARDT A FUCK (UnB)
PROF. VASILE MARZA (UnB-SISMOLOGIA)
PROF. LUIZ FERNANDO S. BRAGA (GEOMAG)
PROF. CARLOS ALBERTO MENDONÇA (USP)

ABSTRACT
In an attempt to shed some light on the regional geology, data of airborne geophysics of the Brazil - Canada Geophysical Project (PGBC) were used. Those data were evaluated, processed and interpreted in the light of the geological knowledge of the studied area which is located in the central and northwestern region of Goiás and in the southwestern region of Tocantins, Brazil, involving lands of the Tocantins Structural Province.
In the evaluation of the data of PGBC, a method was proposed for the determination of spikes in the original recording of the magnetic data that revealed to be an efficient technique. Following this, the data were interpolated in regular grid using the method of bi-cubic splines. The data were microleveled using a developed algorithm, which showed efficiency in treating bad-leveled data among flight lines.
Large part of the filters applied on the aeromagnetic data were in the wave number domain. The reduction to the pole filter was used to obtain the pseudo-gravity map, allowing the comparison to the Bouguer map of the area.
The application of enhancement techniques of magnetic sources (horizontal gradient, amplitude of the analytic signals of first and second orders, phase of the analytic signal and terracing), proved to be efficient in the determination of physical discontinuities, translated in geological units. The amplitude of the analytic signal of second order showed better results than the amplitude of the analytic signal of first order in revealing the limits of the causative magnetic bodies and in discriminating mafic and ultramafic rocks.
The application of terracing on the amplitude of the analytic signal of second order allowed placing the anomalous answer on the plan of the outcrops in map, defining the magnetic units better in relation to the well-known geology.
The phase of the analytic signal exhibited an entangle of linear structures (lineaments and limits of magnetic property) showing the complexity of the structural framework and in its interpretation and correlation to the geology.
For estimation of the Curie surface depths in the studied area, specific programs were developed, with routines that turned it more efficient and to recover the data on the borders. A routine was included for subtraction of a first-degree trend surface. The focused estimates, together with some punctual geothermal data, allowed the extrapolation and interpolation for the whole studied area. This shows the importance of the Curie depth estimates.
The analysis of the magnetic images for the anomalous field and its transformed ones made it possible to define four main types of magnetic relief: smooth, rough and uncorrelated, rough and correlated, and very rough.
The analysis of the main magnetic structures interpreted in the images of the phase of the analytic signal, of the amplitude of the analytic signal of second order and of the anomalous magnetic field, together with the geological information available, suggest that these structures must be tied up to an event with main tension, ?1, in the direction approximately E-W.
The use of color composition of the type CMY inverted in the gamma-ray spectrometric ternary maps showed good results. The best combination is K in the cyan (C), Th in the magenta (M) and U in the yellow (Y). This allowed the identification of the limits and subdivisions of the orthogneissic complexes of Caiçara, Anta and Uvá.
A detailed gravity survey (some 400 stations over 100,000 km2) was accomplished in the area of Crixás to study the local greenstone belts. The gravity grid showed better each of the three greenstone belts than available gravity data. Algorithms for terrain correction together with another for indirect determination of densities were developed.
A 2-D modeling (starting models) and 2,5-D and 3-D inversions were made on the greenstone belts using the gravity data and models based on the local geology. These inversions allowed the construction of vertical sections with the probable geology for them. The depths to the base obtained in two inversion methods are compatible and gave indication of about 2 km.
The integrated analysis of the ternary images of K, Th and U, together with the transformed magnetic maps, the Bouguer map, the map of anomalous potassium and the magnetic and gamma-ray spectrometric integration map, using the geological knowledge for the area, suggests that the area has been affected by a sequence of events, beginning with the formation of the granitoid-greenstone belt terrain, the Brasiliano collisional Cycle and finishing with more recent events.
Using the gamma-ray spectrometric data it was possible to estimate what is called anomalous potassium, which revealed lineaments that coincide with some occurrences of primary gold. These same occurrences were added to the interpretative maps. They are associated to magnetic structures, including lineaments picked up in the phase of the analytic signal, and along some contacts between different lithologic units. This fact, allied to the lineaments found in the anomalous potassium, can indicate shear zones, mainly in the directions N55W and N65W, in which there would be possibility of the occurrence of hidrotermal minerals.
Maps of magnetic interpretation, gamma-ray spectrometric interpretation and of integration of the study area represent regional sketches of geological units and structural elements. The validity of these interpretations should be verified in field.