Teoría básica y aplicación de la tomografía sísmica
La tomografía sísmica es un método geofísico que se utiliza para "fotografiar" el interior de la Tierra mediante ondas sísmicas. Al igual que una tomografía computarizada en medicina, que crea una imagen tridimensional del cuerpo a partir de rayos X, la tomografía sísmica modela las estructuras del subsuelo utilizando datos sobre el tiempo de llegada, la amplitud o las formas de onda registradas por una red de sismómetros. Este método es crucial porque el interior de la Tierra no se puede observar directamente, mientras que fenómenos como la tectónica de placas, la actividad volcánica y los terremotos están fuertemente influenciados por las variaciones en las propiedades físicas de las rocas en profundidad.
1. Concepto básico de las ondas sísmicas
Las ondas sísmicas son ondas elásticas que se propagan a través de las rocas. Generalmente, las ondas sísmicas se dividen en:
1. Ondas corporales
– Ondas P (primarias/de compresión): se propagan más rápido, pueden atravesar sólidos y líquidos, y son sensibles a los cambios en la velocidad de compresión y la densidad.
– Ondas S (secundarias/de cizallamiento): más lentas que las ondas P, no pueden propagarse en fluidos, sensibles a la rigidez (módulo de cizallamiento) de la roca.
2. Ondas superficiales
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
La tomografía sísmica se basa en dos conceptos de cálculo:
a) Problema hacia adelante (problema hacia adelante)
A partir de un modelo de la Tierra (por ejemplo, una distribución de velocidad de las ondas), podemos calcular una predicción de los datos sísmicos: el tiempo de viaje de una onda desde una fuente (terremoto o fuente artificial) hasta una estación. En un enfoque sencillo basado en la teoría de rayos, se supone que la trayectoria de la onda es un rayo que sigue el principio de Fermat: la onda elige la trayectoria con el menor tiempo de viaje.
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
Lo contrario también es cierto: disponemos de datos de observación (tiempos de viaje, residuos de tiempo de viaje, dispersiones de ondas superficiales o formas de onda) y queremos estimar el modelo de velocidad que mejor se ajuste. El problema inverso suele ser no único y mal condicionado: varios modelos pueden ajustarse igualmente bien, y los datos contienen incertidumbre o ruido. Por lo tanto, se requiere regularización, como suavizado, amortiguación o restricciones geológicas.
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) Tomografía de forma de onda / inversión de forma de onda completa (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. Adquisición de datos
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. Identificación y selección de fases
Determinación del tiempo de llegada de ondas P, S o superficiales. La calidad de la selección determina la calidad del modelo.
3. Corrección inicial y modelado
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. Trazado de rayos / simulación de ondas
Calcula la trayectoria de los rayos o formas de onda sintéticas para construir la matriz de sensibilidad.
5. Inversión y regularización
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. Evaluación de la resolución y pruebas de fiabilidad
Por ejemplo, con la prueba del tablero de ajedrez, la prueba de picos o el análisis de la función de covarianza/dispersión de puntos para ver qué partes están realmente resueltas por los datos.
5. Interpretación del modelo de velocidad
El resultado de la tomografía suele ser un mapa de anomalías de velocidad relativa: las zonas de alta velocidad se interpretan a menudo como rocas más frías, densas o rígidas (por ejemplo, placas en subducción). Las zonas de baja velocidad pueden indicar altas temperaturas, rocas alteradas, zonas de fractura llenas de fluidos o fusión parcial, frecuentemente asociadas con volcanes activos o astenosfera caliente.
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
La tomografía sísmica es una herramienta clave para investigar la estructura interna de la Tierra a escalas que van desde lo local hasta lo global. Su base teórica reside en la relación entre las propiedades elásticas de las rocas y la propagación de ondas, así como en la resolución de problemas inversos que requieren regularización y pruebas de resolución. Sus aplicaciones son amplias: desde el mapeo de placas subducidas y sistemas de magma volcánico hasta la exploración energética y la evaluación del riesgo sísmico. Gracias a los avances en redes sísmicas, computación y métodos de inversión (incluida la inversión de forma de onda completa), la tomografía sísmica es cada vez más capaz de producir imágenes más nítidas y con mayor información, tanto para fines científicos como para la mitigación de desastres.
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).
