Teoría básica e aplicación da tomografía sísmica
A tomografía sísmica é un método xeofísico empregado para "fotografar" o interior da Terra mediante ondas sísmicas. De xeito semellante a unha tomografía computarizada en medicina, que crea unha imaxe tridimensional do corpo a partir de raios X, a tomografía sísmica modela as estruturas do subsolo utilizando datos sobre o tempo de percorrido, a amplitude ou as formas de onda rexistradas por unha rede de sismómetros. Este método é crucial porque o interior da Terra non se pode observar directamente, mentres que as dinámicas como a tectónica de placas, a actividade volcánica e os terremotos están fortemente influenciadas polas variacións nas propiedades físicas das rochas en profundidade.
1. Concepto básico das ondas sísmicas
As ondas sísmicas son ondas elásticas que se propagan a través da rocha. As ondas sísmicas divídense xeralmente en:
1. Ondas de corpo
– Ondas P (primarias/de compresión): propáganse máis rápido, poden atravesar sólidos e líquidos, son sensibles aos cambios na velocidade e densidade de compresión.
– Ondas S (secundarias/de cizallamento): máis lentas que as P, non se poden propagar en fluídos, sensibles á rixidez (módulo de cizallamento) da rocha.
2. Ondas superficiais
Por exemplo, as ondas de Rayleigh e Love, que xeralmente son dominantes nos rexistros de terremotos distantes, teñen dispersión e son moi informativas para a estrutura da litosfera ata a astenosfera pouco profunda.
A velocidade das ondas sísmicas depende dos parámetros elásticos (módulo) e da densidade. As variacións de temperatura, composición mineral, presión, porosidade, fracturas e presenza de fluídos ou fusións parciais afectarán a velocidade. Esta é a base física para mapear as variacións de velocidade co fin de proporcionar pistas sobre as condicións xeolóxicas do subsolo.
2. Principios da tomografía: problema frontal e problema traseiro
A tomografía sísmica baséase en dous conceptos de cálculo:
a) Problema directo (problema directo)
Dado un modelo da Terra (por exemplo, unha distribución da velocidade das ondas), podemos calcular unha predición dos datos sísmicos: o tempo de percorrido dunha onda desde unha fonte (terremoto ou fonte artificial) ata unha estación. Nunha teoría simple dos raios, asúmese que a traxectoria da onda é un raio que segue o principio de Fermat: a onda escolle a traxectoria co tempo de percorrido mínimo.
Matematicamente, o tempo de viaxe T pódese escribir como:
\[
T = \int_{\text{ray}} \frac{ds}{v(\mathbf{x})}
\]
onde \(v(\mathbf{x})\) é a velocidade da onda na posición \(\mathbf{x}\) e \(ds\) é o elemento da lonxitude do percorrido.
b) Problema inverso
O contrario ocorre: temos datos observacionais (tempos de viaxe, residuos do tempo de viaxe, dispersións de ondas superficiais ou formas de onda) e queremos estimar o modelo de velocidade que mellor se axusta. O problema inverso xeralmente non é único e está mal formulado: varios modelos poden axustarse igual de ben e os datos conteñen incerteza/ruído. Polo tanto, requírese regularización, como suavización, amortiguación ou restricións xeolóxicas.
Na práctica, a inversión realízase a miúdo de xeito lineal: o modelo inicial perturbáse lixeiramente e, a continuación, a diferenza no tempo de percorrido calcúlase como:
\[
ΔT aprox. int_{text{ray}} Δs(\mathbf{x})\, ds
\]
onde Δs = Δ(1/v) é a perturbación de lentitude. Estas ecuacións están organizadas nun gran sistema lineal (d = Gm) e resólvense mediante o método dos mínimos cadrados amortecidos ou as súas variacións.
3. Tipos de tomografía sísmica
1) Tomografía en tempo de viaxe
O método máis común emprega chegadas en fase P e S de moitos terremotos. Os datos son residuos de tempo de viaxe contra un modelo de referencia (por exemplo, IASP91 ou ak135 a escala global). É axeitado para a modelización de velocidade en 3D na codia e no manto superior, especialmente se a rede de estacións é densa.
2) Tomografía de onda superficial
Usando a dispersión (cambio na velocidade de fase/grupo co período) das ondas de Rayleigh/Love. Sensible a estruturas superficiais e intermedias, esta tomografía é moi útil para mapear o grosor da litosfera, as zonas de baixa velocidade (LVZ) e as variacións de temperatura.
3) Tomografía telesísmica (tomografía telesísmica)
Usando terremotos distantes (telesísmicos) cuxas ondas viaxan a través do manto e son rexistradas por unha rede local. A vantaxe: teñen múltiples fontes e proceden de varias direccións, axudando así a "iluminar" o volume subterráneo da área de estudo, por exemplo debaixo de volcáns ou zonas de subdución.
4) Tomografía de forma de onda / inversión completa da forma de onda (FWI)
Usar a forma de onda completa, non só o tempo de chegada, teoricamente produce unha alta resolución porque utiliza plenamente a información de amplitude e fase, pero require grandes requisitos computacionais e un bo modelo inicial para evitar trampas mínimas locais.
4. Etapas xerais de procesamento e inversión
1. Adquisición de datos
Colección de rexistros de sismómetros permanentes ou temporais. A escala de exploración, a fonte pode ser unha explosión ou un vibrosismo; a escala rexional-global, a fonte primaria é un terremoto.
2. Identificación e selección de fases
Determinación do tempo de chegada das ondas P, S ou superficiais. A calidade da selección determina a calidade do modelo.
3. Corrección inicial e modelización
Corrección do tempo (desviación do reloxo), corrección da elevación e selección do modelo de referencia. O modelo inicial pode ser 1D ou 3D simple.
4. Trazado de raios / simulación de ondas
Calcula a traxectoria dos raios ou formas de onda sintéticas para construír a matriz de sensibilidade.
5. Inversión e regularización
Resolución do sistema para obter un modelo de velocidade. A regularización escóllese en función do obxectivo: se se queren destacar anomalías nítidas ou tendencias suaves.
6. Avaliación da resolución e probas de fiabilidade
Por exemplo, coa proba de damero, a proba de pico ou a análise da función de dispersión de covarianza/puntos para ver que partes están realmente resoltas polos datos.
5. Interpretación do modelo de velocidade
O resultado da tomografía adoita ser un mapa de anomalías de velocidade relativa: as zonas de alta velocidade adoitan interpretarse como rochas máis frías, máis densas ou máis ríxidas (por exemplo, placas en subdución). As zonas de baixa velocidade poden indicar altas temperaturas, rochas alteradas, zonas de fractura cheas de fluído ou fusión parcial, a miúdo asociada a volcáns activos ou astenosfera quente.
Non obstante, a interpretación debe ser cautelosa porque a velocidade está influenciada por moitos factores. Idealmente, a tomografía sísmica combínase con outros datos como a gravidade, a magnetotelúrica (MT), a xeodesia, a xeoloxía superficial e a petroloxía.
6. Aplicacións da tomografía sísmica
a) Estudo das zonas de subdución e da dinámica das placas
A tomografía global e rexional pode mapear a placa de subdución do manto, incluíndo a súa xeometría, profundidade e segmentación. Esta información é crucial para comprender as orixes dos grandes terremotos, os mecanismos de acoplamento de placas e a evolución tectónica dunha rexión.
b) Sistemas volcánicos e mitigación de catástrofes
En zonas volcánicas, a tomografía pode identificar zonas de baixa velocidade asociadas a cámaras magmáticas, vías de fluídos ascendentes ou rochas alteradas. Cunha monitorización regular, os cambios de velocidade (tomografía de lapso de tempo) poden proporcionar indicadores do movemento de fluídos/magma relevantes para a alerta temperá de erupcións.
c) Exploración de enerxía e recursos
A escala superficial da codia terrestre, a tomografía utilízase na exploración de hidrocarburos, xeotermia e minería. As variacións de velocidade axudan a mapear a litoloxía, as estruturas de falla, as zonas de fractura e os reservorios. Na xeotermia, a integración da tomografía coa MT adoita ser eficaz: a tomografía proporciona información elástica, a MT proporciona información sobre a condutividade dos fluídos.
d) Caracterización de fallas activas e riscos sísmicos
A tomografía local pode revelar zonas débiles, zonas de danos e heteroxeneidade arredor das fallas. Isto axuda a comprender a segmentación das fallas, o posible bloqueo e as variacións de velocidade que inflúen na amplificación das ondas (efectos do sitio).
e) Estrutura rexional da codia e da litosfera
Usando a tomografía de ondas superficiais, os investigadores poden mapear o grosor litosférico, o límite de Moho e as anomalías astenosféricas. Os resultados constitúen a base dos modelos xeodinámicos, incluíndo a formación de concas, a oroxénese e a evolución continental.
7. Limitacións e desafíos
A tomografía sísmica depende en gran medida da distribución das fontes e estacións. As zonas con poucos terremotos ou estacións dispersas terán unha iluminación deficiente, o que resultará nunha baixa resolución. Ademais, as suposicións da teoría dos raios poden ser menos precisas para unha forte heteroxeneidade ou certas frecuencias. O ruído, os erros de selección e a incerteza na localización do hipocentro tamén poden influír nos resultados da inversión. Polo tanto, a avaliación da resolución e a integración multimétodo son cruciais para evitar a sobreinterpretación.
Peche
A tomografía sísmica é unha ferramenta fundamental para investigar a estrutura interna da Terra desde escalas locais ata globais. A súa base teórica baséase na relación entre as propiedades elásticas das rochas e a propagación das ondas, así como na resolución de problemas inversos que requiren regularización e probas de resolución. As súas aplicacións son amplas: desde o mapeo de placas de subdución e sistemas de magma volcánico ata a exploración enerxética e a avaliación do risco de terremotos. Cos avances nas redes sísmicas, a computación e os métodos de inversión (incluída a inversión completa da forma de onda), a tomografía sísmica é cada vez máis capaz de producir imaxes máis nítidas e informativas tanto para fins científicos como de mitigación de desastres.
Se o desexa, podo adaptar este artigo a unha versión máis académica (con citas e bibliografía) ou centrarme nunha das aplicacións (por exemplo, tomografía para volcáns ou para enerxía xeotérmica).
