Teoria básica e aplicação da tomografia sísmica
Tomografi seismik adalah salah satu metode geofisika yang digunakan untuk “memotret” bagian dalam Bumi dengan memanfaatkan gelombang seismik. Seperti halnya CT-scan pada dunia medis yang menyusun citra tiga dimensi tubuh dari perambatan sinar-X, tomografi seismik menyusun model struktur bawah permukaan dari data waktu tempuh, amplitudo, atau bentuk gelombang (waveform) yang direkam oleh jaringan seismometer. Metode ini sangat penting karena interior Bumi tidak dapat diamati secara langsung, sementara dinamika seperti tectônica de placas, aktivitas vulkanik, dan gempa bumi sangat dipengaruhi oleh variasi sifat fisik batuan di kedalaman.
1. Conceito Básico Ondas sísmicas
Ondas sísmicas adalah gelombang elastik yang merambat melalui medium batuan. Secara umum gelombang seismik dibagi menjadi:
1. Gelombang badan (body waves)
– Gelombang P (primer/kompresi) : merambat paling cepat, dapat melalui padat maupun cair, sensitif terhadap perubahan kecepatan kompresi dan densitas.
– Gelombang S (sekunder/geser) : lebih lambat dari P, tidak dapat merambat dalam fluida, sensitif terhadap kekakuan (shear modulus) batuan.
2. Gelombang permukaan (surface waves)
Por exemplo, as ondas de Rayleigh e Love, que geralmente são dominantes em registros de terremotos distantes, apresentam dispersão e são muito informativas sobre a estrutura da litosfera até a astenosfera superficial.
A velocidade das ondas sísmicas depende de parâmetros elásticos (módulo) e da densidade. Variações na temperatura, composição mineral, pressão, porosidade, fraturas e a presença de fluidos ou fusões parciais afetam a velocidade. Essa é a base física para o mapeamento das variações de velocidade, que fornece indícios sobre as condições geológicas do subsolo.
2. Princípios da Tomografia: Problema Direto e Problema Inverso
Tomografia sísmica berdiri di atas dua konsep perhitungan:
a) Problema direto (problema direto)
Diberikan model Bumi (misalnya distribusi kecepatan gelombang), kita dapat menghitung prediksi dados sísmicos: waktu tempuh gelombang dari sumber (gempa atau sumber buatan) menuju stasiun. Dalam pendekatan sederhana ray theory, lintasan gelombang diasumsikan sebagai sinar (ray) yang mengikuti prinsip Fermat: gelombang memilih lintasan dengan waktu tempuh minimum.
Matematicamente, o tempo de viagem \(T\) pode ser escrito como:
\[
T = \int_{\text{ray}} \frac{ds}{v(\mathbf{x})}
\]
onde \(v(\mathbf{x})\) é a velocidade da onda na posição \(\mathbf{x}\), e \(ds\) é o elemento de comprimento do caminho.
b) Problema inverso
Kebalikannya: kita memiliki data pengamatan (waktu tempuh, residu waktu tempuh, dispersi gelombang permukaan, atau waveform) dan ingin menaksir model kecepatan yang paling cocok. Masalah balik umumnya tidak tunggal (non-unique) dan ill-posed : banyak model dapat sama-sama cocok, dan data memiliki ketidakpastian/noise. Karena itu diperlukan regularisasi, misalnya pemulusan (smoothing), damping, atau batasan geologi.
Na prática, a inversão é frequentemente realizada de forma linearizada: o modelo inicial é ligeiramente perturbado e, em seguida, a diferença no tempo de trânsito é calculada como:
\[
δT ≈ ∫_{\text{ray}} δs(\mathbf{x}) ds
\]
onde \(\delta s = \delta(1/v)\) é a perturbação de lentidão. Essas equações são organizadas em um grande sistema linear \( \mathbf{d} = \mathbf{Gm} \), então resolvido pelo método dos mínimos quadrados amortecidos ou suas variações.
3. Tipos de Tomografia Sísmica
1) Tomografia de tempo de deslocamento
O método mais comum utiliza as chegadas das fases P e S de vários terremotos. Os dados são resíduos de tempo de trânsito em relação a um modelo de referência (por exemplo, IASP91 ou ak135 em escala global). Adequado para modelagem de velocidade 3D na crosta e no manto superior, especialmente se a rede de estações for densa.
2) Tomografia de ondas de superfície
Utilizando a dispersão (variação na velocidade de fase/grupo com o período) das ondas de Rayleigh/Love, esta tomografia, sensível a estruturas rasas a intermediárias, é muito útil para mapear a espessura da litosfera, zonas de baixa velocidade (ZBV) e variações de temperatura.
3) Tomografia telesísmica (tomografia telesísmica)
Utilizando terremotos distantes (telessísmicos) cujas ondas viajam pelo manto e são então registradas por uma rede local. A vantagem: eles têm múltiplas fontes e vêm de várias direções, ajudando assim a "iluminar" o volume subterrâneo da área de estudo, por exemplo, sob vulcões ou zonas de subducção.
4) Tomografi waveform / full waveform inversion (FWI)
Utilizar a forma de onda completa, e não apenas o instante de chegada, teoricamente proporciona alta resolução, pois aproveita totalmente as informações de amplitude e fase, mas exige grande poder computacional e um bom modelo inicial para evitar armadilhas de mínimos locais.
4. Etapas gerais de processamento e inversão
1. Dados de Akuisisi
Coleta de registros de sismógrafos permanentes ou temporários. Em escala de exploração, a fonte pode ser uma explosão ou vibrossismo; em escala regional-global, a fonte primária é um terremoto.
2. Identifikasi fase dan picking
Determinar o tempo de chegada das ondas P, S ou de superfície. A qualidade da seleção determina a qualidade do modelo.
3. Koreksi dan pemodelan awal
Correção de tempo (desvio do relógio), correção de elevação e seleção do modelo de referência. O modelo inicial pode ser 1D ou 3D simples.
4. Ray tracing / simulasi gelombang
Calcule o percurso dos raios ou formas de onda sintéticos para construir a matriz de sensibilidade.
5. Inversi dan regularisasi
Resolvendo o sistema para obter um modelo de velocidade. A regularização é escolhida com base no objetivo: destacar anomalias acentuadas ou tendências suaves.
6. Evaluasi resolusi dan uji keandalan
Misalnya dengan checkerboard test , spike test, atau analisis covariance/point spread function untuk melihat bagian mana yang benar-benar teresolusi oleh data.
5. Interpretação do Modelo de Velocidade
Hasil tomografi biasanya berupa peta anomali kecepatan relatif: zona kecepatan tinggi sering diinterpretasikan sebagai batuan lebih dingin, lebih rapat, atau lebih kaku (misalnya slab subduksi). Zona kecepatan rendah dapat menunjukkan temperatur tinggi, batuan teralterasi, zona retakan dengan fluida, atau lelehan parsial—sering terkait gunung api aktif atau astenosfer panas.
No entanto, a interpretação deve ser cautelosa, pois a velocidade é influenciada por muitos fatores. Idealmente, a tomografia sísmica é combinada com outros dados, como gravimetria, magnetotelúrica (MT), geodesia, geologia de superfície e petrologia.
6. Aplicações da Tomografia Sísmica
a) Estudo das zonas de subducção e da dinâmica das placas tectônicas
A tomografia global e regional permite mapear a placa subductada do manto, incluindo sua geometria, profundidade e segmentação. Essas informações são cruciais para a compreensão das fontes de grandes terremotos, dos mecanismos de acoplamento das placas e da evolução tectônica de uma região.
b) Sistemas vulcânicos e mitigação de desastres
Em áreas vulcânicas, a tomografia pode identificar zonas de baixa velocidade associadas a câmaras magmáticas, trajetórias ascendentes de fluidos ou rochas alteradas. Com monitoramento regular, as mudanças de velocidade (tomografia time-lapse) podem fornecer indicadores do movimento de fluidos/magma relevantes para o alerta precoce de erupções.
c) Exploração de energia e recursos
Em escala crustal rasa, a tomografia é utilizada na exploração de hidrocarbonetos, geotérmica e de mineração. As variações de velocidade auxiliam no mapeamento da litologia, estruturas de falhas, zonas de fratura e reservatórios. Em geotérmica, a integração da tomografia com a MT (tomografia magneto-óptica) costuma ser eficaz: a tomografia fornece informações sobre a elasticidade do fluido, enquanto a MT fornece informações sobre a condutividade do fluido.
d) Caracterização de falhas ativas e riscos sísmicos
A tomografia local pode revelar zonas de fragilidade, zonas de dano e heterogeneidade em torno de falhas. Isso ajuda a compreender a segmentação de falhas, o potencial de bloqueio e as variações de velocidade que influenciam a amplificação de ondas (efeitos de sítio).
e) Estrutura regional da crosta e da litosfera
Utilizando a tomografia de ondas superficiais, os pesquisadores podem mapear a espessura da litosfera, o limite da descontinuidade de Mohorovičić (Moho) e anomalias astenosféricas. Os resultados servem de base para modelos geodinâmicos, incluindo a formação de bacias, a orogênese e a evolução continental.
7. Limitações e Desafios
A tomografia sísmica é altamente dependente da distribuição de fontes e estações. Áreas com poucos terremotos ou estações esparsas terão iluminação deficiente, resultando em baixa resolução. Além disso, as premissas da teoria de raios podem ser menos precisas em casos de forte heterogeneidade ou certas frequências. Ruído, erros de seleção e incerteza na localização do hipocentro também podem influenciar os resultados da inversão. Portanto, a avaliação da resolução e a integração de múltiplos métodos são cruciais para evitar interpretações excessivas.
Fechando
Tomografi seismik merupakan alat utama untuk menyelidiki struktur internal Bumi dari skala lokal hingga global. Dasar teorinya bertumpu pada hubungan antara sifat elastik batuan dan perambatan gelombang, serta pada penyelesaian masalah balik yang membutuhkan regularisasi dan uji resolusi. Aplikasinya luas: dari pemetaan slab subduksi, sistem magma gunung api, eksplorasi energi, hingga penilaian bahaya gempa. Dengan kemajuan jaringan seismik, komputasi, dan metode inversi (termasuk full waveform inversion), tomografi seismik semakin mampu menghasilkan citra yang lebih tajam dan informatif untuk kebutuhan sains maupun mitigasi bencana.
Se desejar, posso adaptar este artigo para uma versão mais acadêmica (com citações e bibliografia) ou focar em uma das aplicações (por exemplo, tomografia para vulcões ou para energia geotérmica).
