Gravimetria

Gravimetria
P R O F. D R . V I N I C I U S L O U R O
Gravimetria
Investigação de contrastes de densidade
◦ Baseado em medidas das variações do campo gravimétrico da Terra
◦ Alcance
◦ Até o núcleo da Terra
Aplicações
◦ Exploração de petróleo
http://www.geologyforinvestors.com/wp-
content/uploads/GravitySurveys.jpg
◦ e.g. Domos salinos
◦ Exploração mineral
◦ e.g. Caracterização de intrusões
◦ Tectônica
◦ e.g. Delimitação de blocos tectônicos
GEOFÍSICA APLICADA - GRAVIMETRIA 2

Princípios Físicos
3a lei de Newton
◦ Atração Universal
http://petgeo.weebly.com/uploads/9/9/1/0/99107128/10-1_orig.jpg
𝐺𝐺𝑚𝑚1 𝑚𝑚2
𝐹𝐹⃗ =
𝑟𝑟⃗ 2
◦ Em que
◦ F ≡ força de atração que a massa m1 exerce em m2
◦ r ≡ distância entre o centro de massa desses corpos

Adaptando ao sistema Terra simplificadamente
◦ Assumimos
◦ Forma da Terra ≡ esférica
◦ Movimentação Estática ≡ Sem efeitos de rotação
◦ Raio da Terra ≡ R m
◦ Massa da Terra ≡ Homogênea de valor M
http://www.bbc.co.uk/schools/gcsebitesize/geography/images/tec_001.g
F
R
◦ Assim
𝐺𝐺𝐺𝐺𝐺𝐺
𝐹𝐹⃗ = 𝑅𝑅2 = 𝑚𝑚𝑎𝑎⃗𝐺𝐺
if
◦ Em que
◦ F ≡ força da gravidade Massa da Terra = M
◦ M ≡ massa da Terra
◦ R ≡ raio da Terra
◦ G ≡ constante de gravitação universal (6.67408 . 10-11 m3.kg-1.s-2)
◦ aG ≡ GM/R² ≡ aceleração da gravidade

Porém a Terra não é tão simples
◦ Assumimos
m
◦ Forma da Terra ≡ esférica
http://www.bbc.co.uk/schools/gcsebitesize/geography/images/t
◦ Quase elipsoidal F
R
◦ Movimentação Estática ≡ Sem efeitos de rotação
◦ Rotacional
◦ Raio da Terra ≡ R
ec_001.gif
◦ Topografia irregular
Massa da Terra = M
◦ Massa da Terra ≡ Homogênea de valor M
◦ Heterogênea
◦ Tudo isso gera variações no campo gravimétrico

http://www.spacedaily.com/images/earth-gravitymaps-desk-1024.jpg

Método Potencial
◦ Potencial Gravimétrico
◦ Energia potencial (EP)
◦ Força (F) exercida a uma determinada distância (δr)

𝐸𝐸𝑃𝑃 = 𝐹𝐹. 𝛿𝛿𝛿𝛿 = 𝑚𝑚. 𝑎𝑎. 𝛿𝛿𝛿𝛿
◦ Quando falamos de potencial gravitacional (UG)

𝐸𝐸𝑃𝑃 = 𝑚𝑚. 𝑈𝑈𝐺𝐺
◦ Então
𝐸𝐸𝑃𝑃 = 𝑚𝑚. 𝑈𝑈𝐺𝐺 = 𝑚𝑚. 𝑎𝑎𝐺𝐺 . 𝛿𝛿𝑟𝑟
◦ Lembrando que
𝐺𝐺𝐺𝐺𝐺𝐺
𝐹𝐹⃗ = = 𝑚𝑚𝑎𝑎⃗ 𝐺𝐺
𝑅𝑅 2
◦ Temos
𝑀𝑀
𝑈𝑈𝐺𝐺 = 𝐺𝐺
𝑅𝑅

Propriedade física
◦ Dado um corpo de densidade ρ(x,y,z)
◦ Potencial gravitacional (UG) no ponto r(x,y,z)
ρ ( x, y , z )
UG = G ∫ ∫ ∫ dxdydz
XY Z
r ( x, y , z )
◦ Para um ponto fora da esfera
◦ UG e aG serão os mesmos considerando a massa do corpo concentrada em seu centro
𝐸𝐸 𝐸𝐸
𝑈𝑈𝐺𝐺 = 𝐺𝐺 𝑟𝑟 𝑎𝑎𝐺𝐺 = 𝐺𝐺 𝑟𝑟𝑟

Superfícies equipotenciais
◦ Potencial gravitacional
◦ Campo vetorial conservativo
◦ Superfícies constantes ≡ superfícies equipotenciais.
◦ Para uma esfera com massa m

◦ Potencial gravimétrico varia apenas em função da distância r do seu centro
◦ Assim temos superfícies equipotenciais concêntricas à esfera, tal que:

◦ Equipotencial un dista rn do centro
◦ u0 coincide com o limite da esfera
◦ Energia do trabalho ao longo dessa superfície é zero
𝑑𝑑𝜏𝜏 = −𝐹𝐹. 𝑑𝑑𝑑𝑑 = 𝐹𝐹 cos 𝜃𝜃 𝑑𝑑𝑑𝑑

Potencial centrífugo
◦ Atração gravitacional e do efeito de rotação terrestre
◦ A aceleração cetrífuga (ac) é definida em função do eixo de rotação
◦ Considere um ponto de rotação na Terra com raio r
◦ O ângulo θ entre o raio e o eixo de rotação é chamado de colatitude (complementar a latitude λ)
◦ O potencial centrífugo é dado por
𝜕𝜕𝜕𝜕𝜕𝜕
𝑎𝑎⃗ 𝑐𝑐 = − 𝑥𝑥⃗ = 𝑤𝑤 2 𝑥𝑥 𝑥𝑥⃗
𝜕𝜕𝜕𝜕

Forma da Terra
◦ O formato da Terra é bem irregular
◦ Com uma distribuição de massa heterogênea
◦ Influenciando diretamente na gravidade
◦ Dois modelos terrestres

◦ Elipsóidal
◦ Geóidal
Por que não usar um único modelo para a superfície terrestre?

◦ A forma real da Terra leva em consideração infinitos parâmetros
◦ Matematicamente inviável
◦ Então utiliza-se dois modelos principais

◦ Modelo elipsoidal
◦ Geometria para uma distribuição de massas homogênea
◦ Modelo geoidal
◦ Incorpora as deformações geométricas decorrentes da distribuição heterogênea da massa

Modelo elipsoidal ou Terra Normal:
◦ Constituído por um elipsóide de revolução
◦ Distribuição homogênea de massa
◦ Forma Geométrica
◦ Que melhor se aproxima da forma real da Terra
◦ Uso
◦ Cálculos de posicionamento geodésico
◦ Entretanto
◦ Por não ter a mesma distribuição das massas da Terra real
◦ Inadequado para representar a massa terrestre em equilíbrio
◦ Calculado com base em observações de satélites artificiais

◦ Elipsóide que melhor ajusta o formato da Terra
◦ International Reference Ellipsoid

O modelo geoidal
◦ É o modelo terrestre que mais se aproxima da forma real da Terra
◦ Definição de superfície geoidal
“Superfície quipotencial do campo da gravidade que coincide com o nível médio dos mares
não perturbados e, imaginariamente se prolonga pelos continentes”
◦ Incorpora as deformações geométricas
◦ Associadas à distribuição das massas da Terra
◦ Inapropriado para o posicionamento geodésico de pontos na superfície física

Modelo geoidal
◦ Sobre os continentes o geóide é afetado pela massa sobre o nível do mar
◦ Ex: montanhas, planaltos, etc.
◦ Esse excesso de massa gera um aumento na atração da gravidade
◦ Gera uma elevação local do geóide sobre o elipsóide
◦ Denominado ondulação geoidal (N)

Convenções
◦ Altura geoidal (A)
◦ Distância sobre a normal entre um ponto do geóide e sua projeção no elipsóide
◦ Convenção
◦ Positivas no exterior e negativas no interior do geoide
H
◦ Caracterizam as ondulações geoidais (N)
◦ Altitude ortométrica (H) A
◦ Distância na vertical de um ponto na superfície terrestre e o geóide

Gravidade Normal
◦ Se o interior da Terra fosse uniforme
◦ Valor da gravidade no elipsóide internacional de referência variaria em função da latitude geodésica (λ)
◦ Denominada de gravidade normal (gN).
◦ Fórmula da Gravidade Teórica
◦ Geodetic Reference System 1967 – em mGal = 10-5 m/s2
𝑔𝑔𝑁𝑁 = 978031.8 1 + 5.28 . 10−3 . sin2 𝜆𝜆 − 2.34 . 10−5 sin4 𝜆𝜆

Valor da Gravidade Normal
◦ Importante pois fornece o valor da variação teórica da gravidade
◦ Em função da latitude na superfície elipsoidal
◦ Definida por 3 fatores

◦ Distância até o centro de massa da Terra
◦ Nos polos < No equador
◦ Concentração de massa
◦ Aumenta ligeiramente a aceleração da gravidade
◦ Força centrífuga
◦ Portanto...
◦ A gravidade é maior nos Polos do que no Equador
Superfície gravimétrica de Marte
http://geodesy.curtin.edu.au/research/models/mgm2011/

Gravimetria
Na prática
◦ Investigação de variações na densidade da litosfera (ρ)
◦ Mais denso = Anomalia positiva
◦ Dique de diabásio (ρ ≈ 3.1 g/cm³)

vs.
◦ Embasamento metassedimentar (ρ ≈ 2.6 g/cm³)
◦ Menos denso = Anomalia negativa
◦ Espessas cadeias de montanhas

◦ Composição andesítica (ρ ≈ 2.7 g/cm³)
vs.
◦ Sobre o manto listosféricos (ρ ≈ 3.3 g/cm³)

Gravimetria
Na prática
◦ Excelente método para delimitação lateral de estruturas

◦ Limitado pela cobertura e resolução do levantamento
◦ Para investigações de profundidade

◦ Da superfície ao topo do alvo em profundidade
◦ É possível estimar
◦ Não é o ponto forte do método
◦ Ambiguidade inerente ao método

◦ Dois corpos de densidades diferentes
◦ Se o mais denso estiver mais profundo
◦ Se o menos denso estiver mais raso

Gravimetria
Na prática
◦ O sinal gravimétrico pode ser o mesmo
◦ Dependendo
◦ Contrastes de densidade
◦ Profundidades
◦ Formas dos corpos

Gravimetria
Aquisições de Dados
◦ Mede-se a aceleração da gravidade
https://www.geomatrix.co.uk/cms/resources/gall
content/uploads/2017/03/A10-product.jpg
http://microglacoste.com/wp-
◦ Absoluta
ery/cg5-1-w1200h1200.jpg
◦ Gravímetro absoluto
◦ Satélite/modelos
◦ Relativa a uma medida absoluta
PVAwkNFjeys5q_mqSBcQmh9HVrURVSUaJrSM9
tbn0.gstatic.com/images?q=tbn:ANd9GcQPFrnA
◦ Gravímetro diferencial
https://encrypted-
U1MpLNJ
◦ Mede-se a altitude ortométrica local
◦ Distância entre a medida e o geóide
https://media.springernature.com/original/springer-
static/image/art%3A10.1007%2Fs12210-015-0382-
9/MediaObjects/12210_2015_382_Fig9_HTML.jpg

Gravimetria
http://microglacoste.com/wp-content/uploads/2017/03/A10-product.jpg
◦ Gravímetros absolutos
◦ Sistema de queda-livre
◦ Para um objeto em queda a partir de z0 com velocidade inicial u
◦ E equação do movimento para o ponto z no tempo t
1
𝑧𝑧 = 𝑧𝑧0 + 𝑢𝑢𝑢𝑢 − 𝑔𝑔𝑔𝑔𝑔
2
◦ Gravidade absoluta
◦ Ajuste de uma função quadrática para vários tn e zn
◦ O tempo de queda é medido através de um laser.
◦ Precisão
◦ 0,005 a 0,010 mGal

Gravimetria
repair.com/portals/0/Images/upgrade.jpg
http://www.gravitymeter-
http://www.geoequipment.info/image.php?id=1082
◦ Gravímetros Diferenciais
◦ Sistema Massa-Mola
◦ Os primeiros gravímetros se baseavam na lei de Hooke
Gravímetro LaCoste &
Gravímetro Scintrex
◦ A massa m suspensa por uma mola de comprimento s0 Romberg, Modelo G
CG-5
◦ Causa uma distensão s
◦ A extensão é então proporcional a g
𝐹𝐹 = 𝑚𝑚𝑚𝑚 = −𝑘𝑘 𝑠𝑠 − 𝑠𝑠0
◦ onde k é a constante da mola
◦ O gravímetro é calibrado para um ponto com g conhecido
◦ Possível medir a diferença de g para dois pontos diferentes
◦ Precisão: 0,010 mGal

Gravimetria
www2.unb.ca/gge/News/2006/Gravimeter/1.jpeg
content/uploads/2016/09/CG5_TREPIED-1.jpg
www.georeva.eu/wp-
◦ Nivelamento
Gravímetro LaCoste & Romberg Gravímetro

Modelo G Scintrex CG-5
◦ Gravímetro sempre nivelado
◦ Obtenção da aceleração vertical da gravidade
◦ Efeitos laterais de corpos com contraste de densidade
◦ Praticamente desprezíveis

Gravimetria
◦ Gravímetro LaCoste & Romberg
◦ Gravímetro de aquisição analógica
http://www.gravitymeter-repair.com/portals/0/Images/upgrade.jpg
◦ Extremamente sensível
◦ Nunca deve ser deixado sem uma fonte de energia

◦ Tomada
◦ Bateria
◦ Medidas tomadas pelo operador

◦ Normalmente 3 medidas
◦ Feita a média

Gravimetria
◦ Gravímetro Scintrex CG-5
◦ Gravímetro de aquisição automática

◦ Digital
◦ Extremamente sensível
◦ Nunca deve ser deixado sem uma fonte de energia

◦ Tomada
◦ Bateria
◦ Número de medidas por minuto definida pelo usuário

◦ Realizada a media dos valores adquiridos e o desvio padrão

Gravimetria
https://geotop.com.pe/wp-content/uploads/2017/03/497b_gs08-gs14uhf.jpg
◦ Medição da altitude ortométrica
◦ GPS diferencial
◦ Nivelado
◦ Com ao menos 4 satélites na visada
◦ Tempo de aquisição variável

◦ Depende
https://c1.staticflickr.com/3/2504/3890982064_eae33d446c_z.jpg?zz=1
◦ Tipo de processamento
◦ Número de GPS’s disponíveis
◦ Distância da base
◦ Real Time Kinematic (RTK)

◦ 3 minutos
◦ Posicionamento por Ponto Preciso (PPP)
◦ 10 a 15 minutos

Gravimetria
◦ Calibração
◦ LaCoste & Romberg
◦ Linha de pontos de referência
◦ Grande variação topográfica
◦ “g” é conhecido
◦ Ex: Linha do Observatório Nacional até o topo da Serra de Itatiaia
◦ CG-5
◦ Autocalibração no próprio gravímetro

Campo de Gravimetria
1) Checar a região
◦ Geologia
◦ Geomorfologia
◦ Pontos de Acesso
2) Checar Referências de Nível (RN; IBGE)
3) Ir para a região
4) 1a medida
◦ RN (Porta da Igreja, marcos em estradas, etc.)
◦ Não precisa medir altitude com GPS
5) Fazer medições (instalar estações gravimétricas)

◦ Precisa medir altitude com GPS
◦ Marcar a hora
6) Última medida
◦ RN (Porta da Igreja, marcos em estradas, etc.)
◦ Não precisa medir altitude com GPS
Importante!!!
◦ Marquem sempre as medidas em mais de um local

Gravimetria
Correções
◦ 1) Deriva instrumental
◦ 2) Correção de efeitos de latitude
◦ 3) Efeitos de mare sólida

◦ Atração luni-solar
https://pubs.usgs.gov/of/2000/of00-304/htmldocs/chap02/images/gravity6.gif
◦ 4) Efeito da altitude ortométrica
◦ Anomalia Ar-Livre
◦ 5) Efeito da massa referente à altitude ortométrica

◦ Anomalia Bouguer Simples
◦ 6) Efeitos de topografia
◦ Anomalia Bouguer Completa
◦ 7) Efeitos de deslocamento
◦ Anomalia Eötvos
◦ Levantamentos transportados
◦ 8) Anomalia isostática

Gravimetria
Deriva instrumental
◦ Variações nas propriedades da mola

◦ Temperatura
◦ Tempo de uso do equipamento
◦ Correção
◦ Medidas iniciais e finais em campo no mesmo local

◦ Marcos do IBGE
◦ Altitude
◦ Aceleração da Gravidade
◦ Anotação das horas de todas as medidas

◦ Curva de deriva

Gravimetria
Correção de efeitos de latitude
◦ A Terra não é perfeitamente esférica
◦ Força centrífuga exercida em um corpo sobre a superfície varia em função da latitude
◦ Gravidade varia em função da latitude
◦ Pode ser feita em função a partir da fórmula internacional da gravidade (gN)
𝑔𝑔𝑁𝑁 = 978031.8 1 + 5.28 . 10−3 . sin2 𝜆𝜆 − 2.34 . 10−5 sin4 𝜆𝜆

F
F aC

Gravimetria
Efeitos de mare sólida (atração luni-solar)
◦ Durante o dia
◦ Medições sujeitas a atração de maré
◦ Esse efeito varia com o tempo
◦ Pode ser calculado com precisão
◦ Importante saber a hora das medições
Representação da variação da amplitude da maré em função da sobreposição

do efeito lunar e solar (Lowrie, 2007).

Gravimetria
Efeito da altitude ortométrica (Anomalia Ar-Livre)
◦ Anomalia ar-livre (CA):
◦ Corrige a diferença entre...

◦ Altitude do ponto amostrado (na superfície terrestre) e o geóide de referência
◦ Determinada por...
◦ Gradiente vertical da aceleração da gravidade (∂g/∂h)
◦ Altitude ortométrica do ponto medido

◦ h = distância vertical
◦ Ponto de medida na superfície da Terra e o geóide
𝜕𝜕𝜕𝜕
𝐶𝐶𝐹𝐹𝐹𝐹 = − 𝜕𝜕𝜕 ℎ ≈ 0.3086ℎ
◦ Assim...
∆𝑔𝑔𝐹𝐹𝐹𝐹 = 𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔 + 0.3086ℎ − 𝑔𝑔𝑔𝑔

Gravimetria
Efeito da massa referente à altitude ortométrica
◦ Anomalia Bouguer Simples (CB):
◦ Diferença de massa topográfica entre o geóide e a superfície física da Terra
◦ Tanto ausência (ex. vales) quanto excesso de massa (alto topográfico)
◦ A aceleração gravitacional exercida pelas massas topográficas é expressa por:
𝐶𝐶𝐵𝐵 = 2𝜋𝜋𝜋𝜋𝜋𝜋𝜋 ≈ 0.041909. 10−8
◦ onde
◦ G = 2.670 x 10-8 cm3g-1s-2
◦ Constante de Gravitação Universal
◦ ρ = 2.67 g.cm-3
◦ Densidade média da crosta

Gravimetria
Efeitos de topografia (Anomalia Bouguer Completa)
◦ Correção de Terreno (ΔgT)
◦ Se a medida for feita próximo a uma montanha

◦ Influência dessa deve ser corrigida
◦ Cálculo
◦ Dividindo do terreno em diferentes prismas
◦ Considerando a contribuição de cada um
∆gT = Gρφ0 [( r 2
)(r
+ h 2 − r1 − 2
+ h 2 − r2 )]
◦ Geralmente necessárias quando

◦ Diferença topográfica dentro de um setor
◦ Superior a cerca de 5% da sua distância entre estações

Gravimetria
Efeitos de deslocamento (Anomalia Eötvos)
◦ Correção de Eötvös (gE)
◦ Deve ser aplicada quando o gravímetro encontra-se em movimento

◦ Ex.: Navio
◦ Depende da direção do movimento.
◦ Essa correção é dada por
g E = 7.503v cos λsenα + 0.004154v 2 em nós
g E = 14.585v cos λsenα + 0.015696v 2 em m/s
◦ Onde
◦ v = velocidade do navio
◦ α = azimute
◦ λ = latitude.

Gravimetria
A anomalia Bouguer é dada pela aplicação de todas as correções (necessárias)
∆g B = g OBS + (∆g FA − ∆g B + ∆gT + ∆gTide ) − g N
Onde:
◦ gOBS e gN são os valores da gravidade observada e teórica,
◦ ΔgFA é a correção Ar-Livre ou Faye,
◦ ΔgB é a correção Bouguer,
◦ ΔgT écorreção de terreno e
◦ ΔgTide é a correção de maré.

Gravimetria
Isostasia
◦ Processo de compensação à imposição de forças sobre a crosta

◦ Como a crosta reage ao acúmulo de massa/carga sobre ela
◦ e.g.: regiões de orógenos
◦ Modelos
◦ Modelo de Airy-Heiskanen
◦ Modelo de Pratt-Hayford
◦ Modelo de placa elástica de Vening-Meinesz

Gravimetria
Modelo de Airy-Heiskanen
◦ Crosta terrestre fina e sólida pouco resistente

◦ Apoiada sobre um substrato de densidade maior e de consistência plástica
◦ Manto
◦ Inserção de carga
◦ E.g.: orógeno extenso
◦ Rompimento da crosta e afundamento do orógeno

◦ Até equilíbrio por empuxo
◦ Compensação da massa do relevo

◦ Falta de massa criada pelo deslocamento de parte do material do substrato
◦ Substituição pela crosta menos densa

Gravimetria
Modelo de Pratt-Hayford
◦ Espessura da crosta constante
◦ Divisão em colunas
◦ Variação lateral da densidade
◦ Para uma dada profundidade
◦ Peso das várias colunas de diferentes densidades

◦ Equilibrado pela força de empuxo
◦ Quanto mais elevada for a coluna acima de uma base
◦ Menor a densidade

Gravimetria
Modelo de placa elástica de Vening-Meinesz
◦ Camada superior mais leve
◦ Flutua sobre um substrato fluido mais denso
◦ Age como uma placa elástica

◦ Sobre um fluido pouco rígido
◦ Compensação se estende lateralmente
◦ Carga topográfica inclina a placa para baixo afundando no substrato

◦ Substrato é deslocado lateralmente
◦ O empuxo empurra a placa para cima

◦ Suportando a inclinação da placa em locais afastados da depressão central

Gravimetria
Produtos Finais
◦ Mapa de Anomalia Ar-Livre

◦ Levantamentos aéreos e marinhos
◦ Mapa de Anomalia Bouguer Simples

◦ Levantamentos Terrestres em terrenos pouco acidentados
◦ Mapa de Anomalia Bouguer Completa

◦ Levantamentos terrestres em terrenos de topografia complexa
◦ Mapa de Anomalias Gradiométricas e de Anomalia Eötvos

◦ Levantamentos aéreos e marinhos
◦ Mapa de Anomalia Isostática

Gravimetria
Interpretação
◦ Anomalia Ar-Livre (Free-Air)
◦ Correção Ar-Livre
◦ Remove variações gravimétricas causadas por diferenças de elevação em relação ao geóide
◦ Geóide
◦ Superfície equipotencial equivalente ao nível médio dos mares não perturbados

Gravimetria
Interpretação
◦ Anomalia Bouguer
◦ Remove o efeito de excesso de massa

◦ Referente à elevação corrigida na correção Ar-Livre
◦ Massa referente à elevação em relação ao nível do geóide
◦ Simples
◦ Assume um terreno aproximadamente plano
◦ Calota Bouguer
◦ Completa
◦ Assume uma topografia complexa
◦ Algoritmos subdividem a topografia em prismas
◦ Cálculo individual das contribuições de massa

Gravimetria
Interpretação
◦ Anomalia gravimétrica
◦ Gerada pela distribuição heterogênea da densidade em profundidade
◦ Contraste de densidade
Shimabara-Beppu graben
Δρ = ρ – ρ0 Kyushu, Japão
Complexo Alcalino de Tapira

◦ Dois tipos
◦ Anomalia positiva
◦ Anomalia negativa

Gravimetria
Interpretação
◦ A anomalia deriva da presença de diferentes densidades na Terra
◦ A “aparência” da anomalia depende de
◦ Dimensões
◦ Contraste de densidade
◦ Profundidade do corpo anômalo
◦ Corpos mais rasos

◦ Tendem a produzir anomalias estreitas e de amplitude mais definida
◦ Pequenos comprimentos de onda
◦ Anomalias de grande comprimento de onda

◦ Tendem a pertencer a grandes estruturas ou corpos em altas profundidades

Gravimetria
Pós-Processamento
◦ Filtragem regional-residual
◦ Objetivo: remover a influencia da componente regional dos dados
◦ Diversas técnicas como: espectro de potência e continuação para cima
◦ Técnicas de realce da anomalia

◦ Objetivo:
◦ Realçar uma informação (lineamentos e contatos)
◦ Derivadas do campo gravimétrico
◦ Estimativa de propriedade
◦ Profundidade de topo
◦ Deconvolução de Euler
◦ EHD-Depth
◦ Densidade
◦ Modelagem
◦ Inversão

Gravimetria
Pós-Processamento
◦ Filtragem regional-residual
◦ Objetivo
◦ Isolar a contribuição da anomalia de interesse
◦ Procedimento
◦ Ajusta uma superfície
◦ Melhor aproximação do comportamento regional
◦ Remoção do campo original
◦ Atenção!
◦ Não falsear a anomalia

Gravimetria
Interpretação
◦ Upward continuation
◦ Estima o comportamento do sinal de uma região para alturas maiores
◦ Objetivo
◦ Atenuar as anomalias de pequeno comprimento de onda (superficiais)
◦ A filtragem é feita ponto a ponto
http://www.ga.gov.au/ausgeonews/ausgeonews201006/geology.jsp

Gravimetria
Interpretação
◦ Após a filtragem regional-residual

◦ Técnicas de realce
◦ Ex: Mapa da
Anomalia Bouguer
◦ Derivadas horizontais
◦ Filtros passa-alta, passa-baixa ou passa-banda
◦ Sinal Analítico (2D e 3D)
◦ Etc.
◦ Importante
◦ Qual o significado geológico do sinal observado?

◦ Como o filtro aplicado influencia o sinal?
◦ Gera falseamento?
◦ É estável?
◦ Aumenta a influência de ruídos de alta frequência?
Derivada Horizontal Derivada Vertical

Gravimetria
Interpretação
◦ Meia Largura
◦ Definida como o ponto ao longo de um perfil

◦ Onde a amplitude da anomalia cai pela metade do valor máximo
g max =
4
3
πG ( ) = GM
R 3∆ρ
z2
z 2
◦ Para uma esfera centrada em (0, 0, z)

◦ Amplitude máxima ocorre sobre o centro de massa
◦ Em x=x1/2
◦ gZ=gmax/2
◦ Sendo z máximo dado por:
x1
z< 2
1
 4 13 − 1 2
 

Gravimetria
Interpretação
◦ Thompson (1982) e Reid et al. (1990)
◦ Dada pelo desenvolvimento da equação diferencial homogênea de Euler
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕

𝑥𝑥 − 𝑥𝑥0 + 𝑦𝑦 − 𝑦𝑦0 + 𝑧𝑧 − 𝑧𝑧0 = 𝜂𝜂(𝑅𝑅 − 𝑀𝑀)
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
◦ Onde:
◦ x0, y0 e z0 = coordenadas da fonte

◦ M = Anomalia gravimétrica
◦ R = Campo gravimétrico residual
◦ Η = Índice Estrutural
◦ Η = Índice Estrutural
◦ Medida da taxa de decaimento da anomalia gravimétrica entre a fonte e o ponto de medida e a fonte
◦ Depende diretamente da geometria
◦ Esfera – 2
◦ Cilindro – 1
◦ Dique – 0
◦ Soleira (Sill) – 0

Gravimetria
Interpretação
◦ Modelo 3D
◦ Objetivo:
◦ caracterizar o comportamento em profundidade
◦ Estimar a região com maior contraste de densidade
◦ Importante:
◦ Ambiguidade dos métodos potenciais!
◦ Sinal associado a uma fonte rasa com pouco contraste X fonte profunda com altas densidades
◦ Como evitar:
◦ Máximo de informações possíveis sobre o alvo de estudo
◦ Ex:
◦ Informações sobre a geologia
◦ Local e/ou regional
◦ Se aflorante
◦ Amostras para estimar densidade em laboratório
◦ Estimar os limites laterais
◦ Sinal analítico, derivadas horizontais, etc.

Gravimetria
Interpretação
Modelo 3D
◦ Modelagem direta
◦ 2D (ex: algoritmo de Talwani et al. 1959)
◦ 3D (ex: Barnett 1976)
◦ Inversão 3D
◦ 3D (ex. Li e Oldenbrug, 2002)

Gravimetria
Interpretação
◦ Platô do Tibet
◦ Escala de placas tectônicas
◦ Anomalia de gravidade Bouguer completa
◦ Ordem de -1000 mGal
◦ O planalto fica a altitudes muito elevadas
◦ Devido ao espessamento da crosta associado à colisão dos continentes Índia e Eurásia

◦ A crosta continental espessada
◦ ρc = 2670 kg.m−3
◦ Manto em horizontes mais profundos
◦ ρm = 3300 kg.m− 3
(GRACE gravity data, Mishra et al., 2012, Journal of Asian Earth Sciences, 48:93–110)

Gravimetria
Interpretação
◦ Costa do Pacífico na América Central
◦ Trincheira oceânica e o limite da placa transformante (placas do Caribe e da América do Norte)

◦ Dados GRACE
◦ Dados regionais usados para modelar variações de densidade na zona de subducção
◦ Associadas a
◦ Arco vulcânico
◦ Desenvolvimento da crosta
◦ Mudanças na placa subdutada
(modelo do site GFZ)

Gravimetria
Interpretação
◦ Bacia de Nordkapp (Noruega)
◦ Modelo de diapiros de sal

◦ Domo dalino 'Uranus' na Bacia Nordkapp (Noruega)
◦ Dados de sísmica + gradiometria gravitacional
◦ Perfil sísmico
◦ Ilustra a falta de refletores na base do sal
◦ Apresenta modelos de extensão de sal
◦ Topo e base da interpretação do pré-poço
https://www.geoexpro.com/articles/2016/04/gravity-for-hydrocarbon-exploration
◦ Dados do poço
◦ Provaram que o modelo mínimo estava errado
◦ Aquisição de novos dados de gravidade e gradiometria gravimétrica

◦ Modelo de densidade
◦ Desenvolvido usando limites sedimentares e todo do sal da sísmica e densidades de poços

Gravimetria
Interpretação
◦ Bacia de Nordkapp (Noruega)
◦ O sal
◦ Contraste de densidade negativa com a encaixante
◦ Gerando uma anomalia de gravidade negativa
◦ Gradiente vertical de gravidade (Gzz)

◦ Detecta a parte mais rasa do sal
◦ Gravidade Gz
◦ Sensível às partes mais profundas do corpo de sal
◦ Planejamento de levantamento
◦ Modelos gravitacionais iniciais
https://www.geoexpro.com/articles/2016/04/gravity-for-hydrocarbon-exploration
◦ Orientação ideal e processamento sísmico para QC das velocidades sísmicas

◦ Velocidades podem ser convertidas em um modelo de densidade
◦ Resultando em anomalias de gravidade correspondentes
◦ Interpretação sísmica
◦ Pode ser melhorada em um loop interativo com modelagem e inversão gravitacional

Gravimetria
Interpretação na Tectônica de Placas
◦ Comportamentos específicos em limites de placas
◦ De acordo com o tipo de rochas envolvidas

◦ Limites Compressivos
◦ Limites Distensivos
◦ Limites Transformantes

Gravimetria
Extra
◦ Passo a passo de como configurar e executar medições
◦ Gravímetro Diferencial Lacoste & Romberg
◦ Gravímetro Diferencial Scintrex CG-5

Gravimetria
Lacoste & Romberg
◦ Configuração em campo
◦ Início da medição
http://www2.unb.ca/gge/News/2006/Gravimeter/1.jpeg
◦ LaCoste e Romberg
◦ 1) Ligar o equipamento
◦ 2) Destravamento da mola
◦ 3) Giro do dial várias vezes no sentido anti-horário

◦ 3.1) Soltando a massa
◦ 4) Giro do dial no sentido horário

◦ 4.1) Até o capilar atingir o valor definido pela companhia produtora
◦ 4.1.1) Informado na aquisição do equipamento
◦ 5) Travamento da mola
◦ 6) Desligamento do equipamento

Gravimetria
Scintrex CG-5
◦ Configuração anterior ao campo
◦ Definição de data e hora
◦ Definição de tempo de medição

◦ Read Time
◦ Definição da frequência de medidas (medidas por minuto)

◦ # of cycles
◦ Em caso de Levantamento regular

◦ Separação de linhas
◦ Separação de estações
◦ Mudança automárica de estações
◦ Modo de visualização dos dados

◦ Numérica
◦ Gráfica

Gravimetria
Scintrex CG-5
◦ Configuração em campo
◦ Inserção de parâmetros básicos
◦ Código do Levantamento
Botão
◦ Cliente
◦ Operador
◦ GPS do gravímetro ou manual (opcional)
◦ Localização
◦ Altitude aproximada

Gravimetria
Scintrex CG-5
◦ Início da medição
◦ 1) Aperte o botão MEASURE CLR

Botão
◦ 2) Nivele o gravímetro
◦ 3) Clique em READ GRAV
◦ 4) Aguarde o tempo configurado
◦ 5) Clique em FINAL DATA
◦ 6) Cheque os valores
◦ 6.1) Se ruins (S.D. > 0.800)

◦ 6.1.1) Aperte CANCEL
◦ 6.2) Se adequados (S.D. ≤ 0.800)

◦ 6.2.1) Aperte RECORD

Gravimetria - Estudos

Enviado por

Direitos autorais:

Formatos disponíveis

Gravimetria - Estudos

Enviado por

Dados do documento

Título original

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Gravimetria - Estudos

Enviado por

Direitos autorais:

Formatos disponíveis

◦ Baseado em medidas das variações do campo gravimétrico da Terra

GEOFÍSICA APLICADA - GRAVIMETRIA 2

◦ F ≡ força de atração que a massa m1 exerce em m2

◦ r ≡ distância entre o centro de massa desses corpos

GEOFÍSICA APLICADA - GRAVIMETRIA 3

GEOFÍSICA APLICADA - GRAVIMETRIA 4

◦ Tudo isso gera variações no campo gravimétrico

GEOFÍSICA APLICADA - GRAVIMETRIA 5

◦ Energia potencial (EP)

◦ Força (F) exercida a uma determinada distância (δr)

◦ Quando falamos de potencial gravitacional (UG)

GEOFÍSICA APLICADA - GRAVIMETRIA 6

◦ Dado um corpo de densidade ρ(x,y,z)

◦ Potencial gravitacional (UG) no ponto r(x,y,z)

◦ Para um ponto fora da esfera

◦ UG e aG serão os mesmos considerando a massa do corpo concentrada em seu centro

GEOFÍSICA APLICADA - GRAVIMETRIA 7

◦ Para uma esfera com massa m

◦ Assim temos superfícies equipotenciais concêntricas à esfera, tal que:

◦ Energia do trabalho ao longo dessa superfície é zero

𝑑𝑑𝜏𝜏 = −𝐹𝐹. 𝑑𝑑𝑑𝑑 = 𝐹𝐹 cos 𝜃𝜃 𝑑𝑑𝑑𝑑

GEOFÍSICA APLICADA - GRAVIMETRIA 8

◦ Atração gravitacional e do efeito de rotação terrestre

◦ A aceleração cetrífuga (ac) é definida em função do eixo de rotação

◦ Considere um ponto de rotação na Terra com raio r

◦ O ângulo θ entre o raio e o eixo de rotação é chamado de colatitude (complementar a latitude λ)

◦ O potencial centrífugo é dado por

GEOFÍSICA APLICADA - GRAVIMETRIA 9

◦ Dois modelos terrestres

Por que não usar um único modelo para a superfície terrestre?

◦ Então utiliza-se dois modelos principais

GEOFÍSICA APLICADA - GRAVIMETRIA 10

◦ Calculado com base em observações de satélites artificiais

GEOFÍSICA APLICADA - GRAVIMETRIA 11

◦ É o modelo terrestre que mais se aproxima da forma real da Terra

◦ Definição de superfície geoidal

◦ Incorpora as deformações geométricas

◦ Associadas à distribuição das massas da Terra

◦ Inapropriado para o posicionamento geodésico de pontos na superfície física

GEOFÍSICA APLICADA - GRAVIMETRIA 12

◦ Sobre os continentes o geóide é afetado pela massa sobre o nível do mar

◦ Ex: montanhas, planaltos, etc.

◦ Esse excesso de massa gera um aumento na atração da gravidade

◦ Gera uma elevação local do geóide sobre o elipsóide

◦ Denominado ondulação geoidal (N)

GEOFÍSICA APLICADA - GRAVIMETRIA 13

◦ Altura geoidal (A)

◦ Distância sobre a normal entre um ponto do geóide e sua projeção no elipsóide

◦ Altitude ortométrica (H) A

◦ Distância na vertical de um ponto na superfície terrestre e o geóide

GEOFÍSICA APLICADA - GRAVIMETRIA 14

◦ Se o interior da Terra fosse uniforme

◦ Denominada de gravidade normal (gN).

◦ Fórmula da Gravidade Teórica

◦ Geodetic Reference System 1967 – em mGal = 10-5 m/s2

𝑔𝑔𝑁𝑁 = 978031.8 1 + 5.28 . 10−3 . sin2 𝜆𝜆 − 2.34 . 10−5 sin4 𝜆𝜆

GEOFÍSICA APLICADA - GRAVIMETRIA 15

◦ Definida por 3 fatores