Teoría básica y aplicación de la tomografía 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. Konsep Dasar 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 ejemplo, las ondas de Rayleigh y Love, que generalmente predominan en los registros de terremotos lejanos, presentan dispersión y son muy informativas sobre la estructura de la litosfera hasta la astenosfera superficial.
La velocidad de las ondas sísmicas depende de los parámetros elásticos (módulo) y la densidad. Las variaciones de temperatura, composición mineral, presión, porosidad, fracturas y la presencia de fluidos o fundidos parciales afectan la velocidad. Esta es la base física para mapear las variaciones de velocidad y obtener información sobre las condiciones geológicas del subsuelo.
2. Principios de la tomografía: Problema directo y problema inverso
Tomografía sísmica berdiri di atas dua konsep perhitungan:
a) Problema hacia adelante (problema hacia adelante)
Diberikan model Bumi (misalnya distribusi kecepatan gelombang), kita dapat menghitung prediksi datos 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.
Matemáticamente, el tiempo de viaje \(T\) se puede escribir como:
\[
T = \int_{\text{ray}} \frac{ds}{v(\mathbf{x})}
\]
donde \(v(\mathbf{x})\) es la velocidad de la onda en la posición \(\mathbf{x}\), y \(ds\) es el elemento de longitud de trayectoria.
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.
En la práctica, la inversión se realiza a menudo de forma linealizada: el modelo inicial se perturba ligeramente y, a continuación, la diferencia de tiempo de viaje se calcula como:
\[
\delta T \approx \int_{\text{ray}} \delta s(\mathbf{x})\, ds
\]
donde \(\delta s = \delta(1/v)\) es la perturbación de lentitud. Estas ecuaciones se organizan en un gran sistema lineal \( \mathbf{d} = \mathbf{Gm} \), que luego se resuelve mediante el método de mínimos cuadrados amortiguados o sus variaciones.
3. Tipos de tomografía sísmica
1) Tomografía de tiempo de viaje
El método más común utiliza las llegadas de las fases P y S de numerosos terremotos. Los datos corresponden a los residuos del tiempo de viaje con respecto a un modelo de referencia (por ejemplo, IASP91 o ak135 a escala global). Es adecuado para el modelado de velocidad 3D en la corteza y el manto superior, especialmente si la red de estaciones es densa.
2) Tomografía de ondas superficiales
Utilizando la dispersión (cambio en la velocidad de fase/grupo con el período) de las ondas de Rayleigh/Love, esta tomografía, sensible a estructuras superficiales e intermedias, resulta muy útil para cartografiar el espesor litosférico, las zonas de baja velocidad (ZBV) y las variaciones de temperatura.
3) Tomografía telesísmica (tomografía telesísmica)
Se utilizan terremotos lejanos (telesísmicos) cuyas ondas viajan a través del manto y son registradas por una red local. La ventaja es que tienen múltiples fuentes y provienen de diversas direcciones, lo que ayuda a "iluminar" el volumen del subsuelo del área de estudio, por ejemplo, bajo volcanes o zonas de subducción.
4) Tomografi waveform / full waveform inversion (FWI)
En teoría, utilizar la forma de onda completa, y no solo el tiempo de llegada, proporciona una alta resolución porque aprovecha al máximo la información de amplitud y fase, pero requiere grandes recursos computacionales y un buen modelo inicial para evitar trampas de mínimos locales.
4. Etapas generales del procesamiento e inversión
1. Datos de Akuisisi
Recopilación de registros de sismómetros permanentes o temporales. A escala de exploración, la fuente puede ser una explosión o un vibrosismo; a escala regional-global, la fuente principal es un terremoto.
2. Identifikasi fase dan picking
Determinación del tiempo de llegada de ondas P, S o superficiales. La calidad de la selección determina la calidad del modelo.
3. Koreksi dan pemodelan awal
Corrección horaria (desviación del reloj), corrección de elevación y selección del modelo de referencia. El modelo inicial puede ser unidimensional o tridimensional simple.
4. Ray tracing / simulasi gelombang
Calcula la trayectoria de los rayos o formas de onda sintéticas para construir la matriz de sensibilidad.
5. Inversi dan regularisasi
Resolver el sistema para obtener un modelo de velocidad. La regularización se elige en función del objetivo: resaltar anomalías pronunciadas o suavizar las tendencias.
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. Interpretación del modelo de velocidad
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.
Sin embargo, la interpretación debe ser cautelosa, ya que la velocidad está influenciada por muchos factores. Idealmente, la tomografía sísmica se combina con otros datos como gravedad, magnetotelúrica (MT), geología, geología superficial y petrología.
6. Aplicaciones de la tomografía sísmica
a) Estudio de las zonas de subducción y la dinámica de las placas.
La tomografía global y regional permite cartografiar la placa subducida del manto, incluyendo su geometría, profundidad y segmentación. Esta información es crucial para comprender el origen de los grandes terremotos, los mecanismos de acoplamiento de placas y la evolución tectónica de una región.
b) Sistemas volcánicos y mitigación de desastres
En zonas volcánicas, la tomografía permite identificar zonas de baja velocidad asociadas a cámaras magmáticas, conductos de fluidos ascendentes o rocas alteradas. Mediante un monitoreo regular, los cambios de velocidad (tomografía de lapso de tiempo) pueden proporcionar indicadores del movimiento de fluidos y magma relevantes para la alerta temprana de erupciones.
c) Exploración de energía y recursos
A escala de la corteza superficial, la tomografía se utiliza en la exploración de hidrocarburos, energía geotérmica y minería. Las variaciones de velocidad ayudan a cartografiar la litología, las estructuras de falla, las zonas de fractura y los yacimientos. En el ámbito geotérmico, la integración de la tomografía con la magnetotelúrica suele ser eficaz: la tomografía proporciona información elástica y la magnetotelúrica, información sobre la conductividad de los fluidos.
d) Caracterización de fallas activas y riesgos sísmicos
La tomografía local puede revelar zonas débiles, zonas dañadas y heterogeneidad alrededor de las fallas. Esto ayuda a comprender la segmentación de las fallas, el posible bloqueo y las variaciones de velocidad que influyen en la amplificación de las ondas (efectos de sitio).
e) Estructura regional de la corteza y la litosfera
Mediante la tomografía de ondas superficiales, los investigadores pueden cartografiar el espesor de la litosfera, el límite de la discontinuidad de Moho y las anomalías astenosféricas. Los resultados constituyen la base de modelos geodinámicos, incluyendo la formación de cuencas, la orogénesis y la evolución continental.
7. Limitaciones y desafíos
La tomografía sísmica depende en gran medida de la distribución de las fuentes y las estaciones. Las áreas con pocos terremotos o estaciones dispersas presentan una iluminación deficiente, lo que resulta en una baja resolución. Además, las suposiciones de la teoría de rayos pueden ser menos precisas en presencia de una fuerte heterogeneidad o ciertas frecuencias. El ruido, los errores de selección y la incertidumbre en la ubicación del hipocentro también pueden influir en los resultados de la inversión. Por lo tanto, la evaluación de la resolución y la integración de múltiples métodos son cruciales para evitar la sobreinterpretación.
Clausura
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.
Si lo desea, puedo adaptar este artículo a una versión más académica (con citas y bibliografía) o centrarme en una de las aplicaciones (por ejemplo, tomografía para volcanes o para energía geotérmica).
