DECONVOLUÇÃO DE EULER: NOVA PERSPECTIVA PARA A INTERPRETAÇÃO 3D DE

ESTRUTURAS GEOLÓGICAS
Luciano Costa Gonçalves, Roberto Alexandre Vitória de Moraes, Elton Luiz Dantas, Augusto César Bittencourt Pires –
Universidade de Brasília e Adalene Moreira Silva – Universidade Estadual de Campinas

Copyright 2005, SBGf - Sociedade Brasileira de Geofísica

colagem de pequenos fragmentos crustais, envolvendo
th
núcleos arqueanos, blocos gnáissicos migmatíticos
This paper was prepared for presentation at the 9 International Congress of the
Brazilian Geophysical Society held in Salvador, Brazil, 11-14 September 2005. paleoproterozóicos e seqüências supracrustais
neoproterozóicas. Todo este conjunto foi afetado por
Contents of this paper were reviewed by the Technical Committee of the 9th
International Congress of the Brazilian Geophysical Society. Ideas and concepts of the intrusões graníticas brasilianas e, geralmente, estão
text are authors’ responsibility and do not necessarily represent any position of the associadas a grandes zonas de cisalhamento
SBGf, its officers or members. Electronic reproduction or storage of any part of this
paper for commercial purposes without the written consent of the Brazilian transcorrentes (Figura 1). Os modelos de evolução
Geophysical Society is prohibited. geodinâmica para esta região sugerem a acreção e o
____________________________________________________________________
amalgamento dos diferentes terrenos tectono-
ABSTRACT estratigráficos durante a Orogênese Brasiliana – Pan
Africana (Jardim de Sá 1994, Dantas et al. 2003, Van
This paper deals with a 3D magnetic data (Seridó Project,
Schmus et al. 1995, 2003).
1973) interpretation of regional geologic structures (Patos
Lineament and João Câmara-Picuí Shear Zone) in Grandes zonas de cisalhamento, como o Lineamento
Northeast Brazil based on horizontal position and depth Patos, Zona de Cisalhamento João Câmara-Picuí e a
for magnetic sources given by Euler deconvolution. Best Zona de Cisalhamento Porto Alegre, destacam-se como
parameters for Euler deconvolutions were found from limites dos diferentes terrenos que constituem o Domínio
different combinations of structural index, window sizes, Rio Grande do Norte e apresentam trends dextrais com
and error tolerances. Best product yelded solutions for the direção E-W e NE-SW (Vauchez et al. 1995, Jardim de
structural trends with depth through 5 km. The Patos Sá 1994). Neste estudo de caso foram selecionadas duas
Lineament (E-W) and João Câmara-Picuí Shear Zone pequenas áreas dentro deste sistema de cisalhamento
(NE-SW) showed maximum depth in the order of 3 km para aplicar a técnica de deconvolução de Euler no
and 2 km, respectively. Deeper magnetic sources (> 5 entendimento do comportamento das zonas de
km) were found under supracrustal sequence of the cisalhamento em profundidade, em diferentes situações
Seridó Belt, in the portions Southwest of Currais Novos tectono-estrutrais. Na Área I, a Zona de Cisalhamento
city. 3D plots proved to be useful to get a glimpse on the João Câmara-Picuí ocorre no contato entre as unidades
geometry and distribution of magnetic material within supracrustais da Faixa Seridó, de idade neoproterozóica,
these structures. com os terrenos paleoproterozóicos do Terreno Cálcio-
alcalino Santa Cruz (Dantas 1997, CPRM 2001, Figura 1
INTRODUÇÃO
e 2). Na Área II, a Zona de Cisalhamento de Patos,
O presente trabalho visou a interpretação tridimensional ocorre no contato entre o terreno paleoproterozóico do
de estruturas geológicas regionais (Lineamento Patos e Domínio Rio Grande do Norte e os terrenos
Zona de Cisalhamento João Câmara-Picuí) na região mesoproterozóicos do Domínio da Zona Transversal
Nordeste do Brasil, baseada em resultados obtidos com a (Brito Neves 2000, Figura 1, 2).
deconvolução de Euler. Este estudo de caso foi efetuado DECONVOLUÇÃO DE EULER
utilizando-se os dados aeromagnéticos do Projeto Seridó
(1973). A melhor parametrização obtida veio das A deconvolução de Euler é uma técnica de mapeamento
diferentes combinações de índices estruturais, tamanho espacial das fontes anômalas de campos potenciais,
de janelas e tolerâncias de erros, com foco naquelas que baseada na equação de homogeneidade de Euler. Esta
mostraram agrupamentos de soluções mais coerentes técnica é usada como ferramenta quantitativa para
para as indicações sobre as estruturas geológicas da estimar a posição horizontal e profundidades das fontes
área de estudo. magnéticas ou gravimétricas (Reid et al. 1990).
CONTEXTO GEOLÓGICO REGIONAL As soluções da deconvolução de Euler são estimadas
para grupos de índices estruturais, que expressam a
No contexto geotectônico da Plataforma Sul-Americana, relação de homogeneidade da equação. Os diferentes
a área em estudo está localizada na porção Nordeste da índices estruturais aliados às posições estimatidas
Província Borborema (Almeida et al., 1981), limitada podem identificar e caracterizar uma grande variedade de
pelos meridianos 35°45’ W e 37°45’ W e pelos paralelos estruturas geológicas, como contatos, foliações, falhas,
5°33’ S e 7°10’ S (Figura 1, 2). diques, soleiras, pipes, intrusões e outras feições
O arcabouço tectônico da porção Nordeste da Província magnéticas (Thompson 1982, Reid et al. 1990).
Borborema é caracterizado por um arranjo complexo de
domínios tectono-estratigráficos limitados por zonas de
cisalhamento regionais. Estas zonas representam uma

Ninth International Congress of the Brazilian Geophysical Society

 DECONVOLUÇÃO DE EULER E INTERPRETAÇÃO 3D DA FAIXA SERIDÓ 2
________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________

A equação de Euler é descrita como nominal de 135 m e orientados na direção E-W, com
espaçamento de 1 km ([Link]).
(x − x0 ) ∂T + ( y − y0 ) ∂T + (z − z 0 ) ∂T = N (B − T ) (1)
PROCESSAMENTO DOS DADOS
∂x ∂y ∂z
Foram realizadas quatro etapas de processamento. Na
onde ( x0 , y0 , z0 ) representa a posição da fonte magnética primeira etapa, os dados foram avaliados quanto à
cujo campo total T foi detectado em ( x, y, z ) sobre um presença de inconsistências (“spikes”) e à distribuição
campo regional B , com N interpretado como o índice de espacial das linhas de vôo. A segunda etapa consistiu na
homogeneidade da equação ou o índice estrutural (Reid definição do algoritmo (Mínima Curvatura) e do tamanho
et al. 1990, GEOSOFT 2000). da célula (250 m) de interpolação. Na terceira etapa
foram corrigidos os desnivelamentos presentes nos
A tabela 1 correlaciona os índices estruturais de dados, visando minimizar estas imperfeições e,
homogeneização da equação de Euler com os principais conseqüentemente, homogeneizar a distribuição espacial
modelos físicos e estruturas geológicas associadas. dos dados. Para tanto, foi utilizada a rotina de
micronivelamento desenvolvida por Blum (1999) com
Tabela 1 – Índices estruturais usados na deconvolução
base no método proposto por Minty (1991). Na quarta
de Euler (Reid et al. 1990).
etapa, foi gerado o campo magnético anômalo usado nas
deconvoluções de Euler (GEOSOFT 2000, Figura 3).
Modelo Físico Estrutura Geológica Índice Estrutural (N) As deconvoluções de Euler foram calculadas para o
Linha de pólos Contatos, falhas e foliações 0 (<0,5)
modelo físico de linha de pólos, com os conjuntos de
parâmetros variando o índice estrutural (N) de 0 e 0,5, o
Linha de pólos Diques e soleiras 1 tamanho da janela de 6 a 15 células e a tolerância de
Ponto polar Corpos 2D e pipes verticais 2 erro 3 a 10%. Os melhores resultados para o
agrupamento das soluções Euler sobre o Lineamento
Ponto dipolar Corpos 3D e intrusões 3
Patos e a Zona de Cisalhamento João Câmara-Picuí
foram N=0, 10 células e 7% de erro (Figura 4).

BANCO DE DADOS 37°30'W 37°00'W 36°30'W 36°00'W

Y
#
João Câmara

Os dados aeromagnetométricos utilizados neste trabalho Lajes

Y
#

foram adquiridos no Projeto Seridó (1973), que recobriu a nT

porção central dos estados do Rio Grande do Norte e da Baixo Alto

Paraíba (Figura 2). Área I

Bom Jesus
6°00'S#Y

39°W 38°W 37°W 36°W 35°W

Currais São José
Y
#
Novos do Campestre
Y
#

75°W 60°W 45°W N

Fortaleza
Caicó
Y
# N
# Equador 6°30'S

4°S
Brasil Y
# Parelhas
15°S
25 0 25 km

30°S Área II
7°00'S
5°S CE
Y
# Patos
Área em
Estudo
RN
Figura 3 – Imagem do campo magnético anômalo.
# Natal
6°S Área I INTERPRETAÇÃO
Projeto As soluções calculadas na deconvolução de Euler
Seridó
apresentaram bons agrupamentos espaciais com as
7°S
estruturas geológicas da área (CPRM 2001),
João #
Pessoa principalmente com os sistemas de cisalhamentos
Área II
regionais, indicando o mapeamento de fontes magnéticas
PB
infracrustais com profundidades variando da superfície
8°S
até 5 km.
Recife#
PE 50 0 50 km As feições mais expressivas correspondem à zona de
cisalhamento com direção E-W associada ao Lineamento
Patos, marcada pelo agrupamento das soluções Euler,
Figura 2 – Localização das áreas de estudo I e II que indicaram profundidades máximas para as fontes
inseridas no Projeto Seridó. magnéticas entre 2 e 3 km (LP, Figura 4). A Zona de
Cisalhamento João Câmara-Picuí apresentou fontes
O Projeto Esperança recobriu 25.000 km2 de área em magnéticas com profundidades máximas de 2,0 km
28.000 km de perfis de produção, adquiridos à altura (ZCJP, Figura 4).

Ninth International Congress of the Brazilian Geophysical Society

 GONÇALVES, MORAES, DANTAS, PIRES E SILVA 3
________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________

Alinhamentos secundários de menor expressão foram A interpretação, em escala regional, mostrou que a
determinados pelas anomalias realçadas nos produtos deconvolução de Euler forneceu bons agrupamentos de
aerogeofísicos. Um forte alinhamento trend NE-SW, soluções correlacionadas aos trends geológicos
localizado na região do Seridó, apresentou profundidades conhecidos na região.
variando de 0 a 5 km. As anomalias que evidenciam as
maiores profundidades foram observadas a SW de
Currais Novos, onde fontes magnéticas com 5 km de
profundidades foram encontradas sob as rochas
metassedimentares da Faixa Seridó (A, Figura 4).
Um forte alinhamento das soluções Euler com direção
N65°E, é visível desde a região de Caicó, passando por
Currais Novos e seguindo em direção a borda do núcleo
arqueano em Bom Jesus (B, Figura 4).
Na região a Oeste de Lajes foi observado um
alinhamento NE-SW de soluções com profundidades
máximas da ordem de 2 km (C, Figura 4). Na região N-
NW de Currais Novos foram observadas alinhamentos
concêntricos de soluções Euler com profundidades
máximas entre 0 e 2 km (D, Figura 4). As estruturas com
direção NW-SE foram marcadas por alinhamentos de
soluções Euler na região das cidades de Currais Novos e
Lajes (E, Figura 4).
Os resultados da deconvolução de Euler para o Projeto
Seridó e, mais especificamente, para a Zona de
Cisalhamento João Câmara-Picuí (Área I) e o
Lineamento Patos (Área II) são observados em detalhe,
com figuras tridimensionais que auxiliam a visualização
das estruturas em profundidade (Figura 5, 6).
As visualizações tridimensionais mostram a geometria
das estruturas em profundidade e ajudam a entender
como o material magnético se distribuiu ao longo das
zonas de cisalhamento. Pode-se observar que a fonte
magnética é representada por minerais magnéticos
cristalizados nos planos da foliação cisalhante
verticalizada, por veios com minerais magnéticos
discordantes aos planos da foliação e estruturas
sigmoidais desenvolvidas nos planos da foliação vertical
(Figura 5, 6). São notáveis o mergulho verticalizado da
Zona de Cisalhamento João Câmara e os mergulhos
para N e S divergentes na Zona de Cisalhamento de
Patos.
CONCLUSÕES
A definição dos parâmetros utilizados na deconvolução
de Euler é o grande problema desta técnica. Os diversos
conjuntos de parâmetros podem fornecer diferentes
grupos de soluções, com diferentes profundidades. Para
uma definição apropriada dos parâmetros utilizados
deve-se gerar um modelo sintético que simule os alvos a
serem investigados ou o arcabouço estrutural da área de
estudo, em seguida, aplicar a deconvolução de Euler e
analisar qual o melhor resultado. Entretanto, em áreas
extremamente complexas do ponto de vista geológico e
estrutural, a definição de um modelo sintético que se
ajuste à geologia é um trabalho extremamente complexo.
Testes empíricos com diferentes combinações de
parâmetros podem ser realizados e, conseqüentemente,
analisados para determinar os melhores agrupamentos
de soluções. Neste trabalho, a definição dos parâmetros Figura 5 – Visualizações 3D da deconvolução de Euler
da deconvolução de Euler foi estritamente empírica. (N=0, 10 células e tolerância 7%) para a Área I. A – Visa
em planta. B e C – Visada inclinada da superfície.

Ninth International Congress of the Brazilian Geophysical Society

 DECONVOLUÇÃO DE EULER E INTERPRETAÇÃO 3D DA FAIXA SERIDÓ 4
________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________

O Lineamento Patos e a Zona de Cisalhamento João Evolution of South American, Rio de Janeiro, pp.:
Câmara-Picuí foram caracterizadas por fontes 151-182.
magnéticas com profundidades máximas de 3 km e 2 km,
respectivamente. A Sudoeste de Currais Novos, CPRM-Companhia de Pesquisa de Recursos Minerais
alinhamentos de soluções Euler com direção NE-SW e 2001. Mapa Geológico do Brasil. Rio de Janeiro.
profundidades de 5 km sugerem que estas fontes (Escala 1:2.500.000). CD-ROM.
magnéticas possam representar o embasamento da CPRM-Companhia de Pesquisa e Recursos Minerais,
Faixa Seridó nesta região. Na região Oeste de Lajes, o acessada em 20/03/2004 no endereço eletrônico
trend NE-SW com profundidades aproximadas de 2 km [Link].
caracterizou o contato tectônico entre as unidades
neoproterozóicas da Faixa Seridó e as rochas gnáissicas Dantas E. L. 1997. Geocronologia U-Pb e Sm-Nd de
migmatíticas do embasamento paleoproterozóico. terrenos arqueanos e paleoproterozóicos do
Agrupamentos concêntricos de soluções Euler na região Maciço Caldas Brandão, NE do Brasil. Instituto de
a N-NW de Currais Novos indicaram profundidades Geociências e Ciências Exatas, UNESP, Tese de
variando entre 0 e 2 km para os granitóides brasilianos e Doutoramento, 208 p.
sugerem a presença em profundidade de intrusões nos Dantas E. L., Van Schmus W. R., Hackspacher P. C.,
gnaisses migmatíticos paleoproterozóicos. Alinhamentos Fetter A. H., Brito Neves B. B., Cordani U., Nutman
de soluções com direção NW-SE indicam a reativação A. P., Williams I. S. 2003. The 3.4-3.5 Ga São José
pós-brasiliana de estruturas antigas. do Campestre Massif, NE Brazil: remnants of the
As visualizações 3D nas Áreas I e II mostraram a oldest crust in South America. Precambiran
geometria em profundidade dos sistemas transcorrentes Research, 2439: 1-25.
do Lineamento Patos e da Zona de Cisalhamento João TM
GEOSOFT 2000. OASIS Montaj , v. 5.1.7. GEOSOFT,
Câmara-Picuí, com os mergulhos divergentes, planos Inc., Toronto.
verticalizados bem evidentes e interação de feições
secundárias, como veios discordantes com a foliação Jardim de Sá E. F. 1994. A Faixa Seridó (Província
principal da zona de cisalhamento. Borborema, NE do Brasil) e o seu significado
geodinâmico na cadeia Brasiliana/Pan-Africana.
A deconvolução de Euler por ser um método quantitativo Instituto de Geociências, Universidade de Brasília,
de localização horizontal e profundidade das fontes Tese de Doutoramento, 804 p.
magnéticas ou gravimétricas pode trazer importantes
contribuições a respeito da geometria de corpos em Minty B. R. S. 1991. Simple micro-levelling for
profundidade ou do arcabouço tectônico regional. aeromagnetic data. Exploration Geophysics, 22:
591-592.
A aplicação da deconvolução de Euler em dados
aeromagnéticos e gravimétricos, que recobrem a porção Reid A. B., Allsop J. M., Granser H., Millet A. J.,
Nordeste da Província Borborema, fornecerá novos Somerton I. W. 1990. Magnetic interpretation in
argumentos para futuras interpretações e discussões three dimensions using Euler deconvolution.
sobre o arcabouço tectônico e sua evolução geodinâmica Geophysics, 55 (1): 80-91.
nesta região.
Thompson D. T. 1982. EULDPH – A new technique for
AGRADECIMENTOS making computer-assisted depth estimates form
magnetic data. Geophysics, 47: 31-37.
À Companhia de Pesquisa e Recursos Minerais (CPRM)
pela cessão e permissão do uso dos bancos de dados Van Schmus W. R., Brito Neves B. B., Hackspacher P.C.,
aerogeofísicos. Ao Laboratório de Geofísica Aplicada Babinski M. 1995. U/Pb and Sm/Nd geocronological
(LGA) do Instituto de Geociências (IG) da Universidade studies of the eastern Borborema Province,
de Brasília (UnB) pela utilização da infra-estrutural. A northeastern Brazil: initial conclusions. Journal of
CAPES pela bolsa de mestrado do primeiro autor. Ao South America Earth Sciences, 8 (3/4): 267-288.
CNPq pelo apoio financeiro projeto 475717-01-5 e
Van Schmus W. R., Brito Neves B. B., Williams I. S.,
471144/2003-7.
Hackspacher P. C., Fetter A. H., Dantas E. L.,
REFERÊNCIAS BIBLIOGRÁFICAS Babinski M. 2003. The Seridó Group of NE Brazil, a
late Neoproterozoic pre-to syn-collisional basin in
Almeida F. F. M., Hasui Y., Brito Neves B. B., Fuck H. A. West Gondwana: insights form SHRIMP U-Pb
1981. Brazilian structural provinces: an introduction. detrital zircon ages and Sm-Nd crustal residence
Earth Sciences Reviews, 17: 291-317. (TDM) ages. Precambrian Research, 127: 287-327.
Blum M. L. B. 1999. Processamento e interpretação de Vauches A., Neves S., Caby R., Corsini M., Egydio-Silva
dados de geofísica aérea no Brasil Central e sua M., Arthaud M. Amaro V. 1995. The Borborema
aplicação à geologia regional e à prospecção shear zone system, NE Brazil. Journal of South
mineral. Instituto de Geociências, Universidade de America Earth Sciences, 8 (3/4): 247-266.
Brasília, Tese de Doutoramento, 229 p.
Brito Neves B. B. 2000. Tectonic history of the Borborema
Province, NW Brazil. In: Cordani U. G., Milani E. J.,
Thomaz Filho A., Campos D. A. (eds) Tectonic

Ninth International Congress of the Brazilian Geophysical Society

 GONÇALVES, MORAES, DANTAS, PIRES E SILVA 5
________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________

Figura 1 – Mapa geológico da área em estudo localizada na porção Nordeste da Província Borborema, que mostra os
principais blocos crustais e as principais zonas de cisalhamento da região (Lineamento Patos, ZCJCP – Zona de
Cisalhamento João Câmara-Picuí, CPRM 2001).

37°30'W 37°00'W 36°30'W 36°00'W

Y
#
João Câmara

P
Lajes

JC
Y
#

ZC
m C

262 435 546 668 816 1034 1593

C
Bom Jesus
Y
#
Área I 6°00'S

E
P

D
JC

D
ZC

B
E Y
#
Currais São José
Novos do Campestre
Y
#
B
E
Caicó
Y
# N
6°30'S
B
A
B Y Parelhas
#
P

25 0 25 km
JC
ZC

Área II
7°00'S
LP YPatos
#
LP LP

Figura 4 – Tema da deconvolução de Euler, para o índice estrutural N=0, janela de 10 células e tolerância 7%.

Ninth International Congress of the Brazilian Geophysical Society

 DECONVOLUÇÃO DE EULER E INTERPRETAÇÃO 3D DA FAIXA SERIDÓ 6
________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________

Figura 6 – Visualizações 3D da deconvolução de Euler (N=0, 10 células e tolerância 7%) para a Área II. A – Vista em planta.
B, C e D – Visada inclinada da superfície, detalhe em C da estrutura sigmoidal presente no plano verticalizado da foliação
principal. E – Visada inclinada em profundidade.

Ninth International Congress of the Brazilian Geophysical Society