Teoria de bază și aplicațiile tomografiei seismice
Tomografia seismică este o metodă geofizică utilizată pentru a „fotografia” interiorul Pământului folosind unde seismice. La fel ca o tomografie computerizată în medicină, care creează o imagine tridimensională a corpului din raze X, tomografia seismică modelează structurile subterane folosind date despre timpul de deplasare, amplitudine sau formele de undă înregistrate de o rețea de seismometre. Această metodă este crucială deoarece interiorul Pământului nu poate fi observat direct, în timp ce dinamica, cum ar fi tectonica plăcilor, activitatea vulcanică și cutremurele, este puternic influențată de variațiile proprietăților fizice ale rocilor la adâncime.
1. Conceptul de bază al undelor seismice
Undele seismice sunt unde elastice care se propagă prin rocă. Undele seismice sunt în general împărțite în:
1. Undele corpului
– Undele P (primare/de compresie): se propagă cel mai rapid, pot traversa solide și lichide, sunt sensibile la modificările vitezei de compresie și ale densității.
– Unde S (secundare/de forfecare): mai lente decât P, nu se pot propaga în fluide, sensibile la rigiditatea (modulul de forfecare) a rocii.
2. Unde de suprafață
De exemplu, undele Rayleigh și Love, care sunt în general dominante în înregistrările cutremurelor îndepărtate, au dispersie și sunt foarte informative pentru structura litosferei până la astenosfera superficială.
Viteza undelor seismice depinde de parametrii elastici (modulul) și de densitate. Variațiile de temperatură, compoziția mineralelor, presiune, porozitate, fracturi și prezența fluidelor sau a topirilor parțiale vor afecta viteza. Aceasta este baza fizică pentru cartografierea variațiilor de viteză pentru a oferi indicii despre condițiile geologice ale subsolului.
2. Principiile tomografiei: Problema anterioară și problema posterioară
Tomografia seismică se bazează pe două concepte de calcul:
a) Problemă directă (problemă directă)
Având în vedere un model al Pământului (de exemplu, o distribuție a vitezei undei), putem calcula o predicție a datelor seismice: timpul de deplasare al unei unde de la o sursă (cutremur sau sursă artificială) la o stație. Într-o abordare simplă a teoriei razelor, se presupune că traiectoria undei este o rază care urmează principiul lui Fermat: unda alege traiectoria cu timpul de deplasare minim.
Matematic, timpul de călătorie ∫(T) poate fi scris ca:
\[
T = \int_{\text{rază}} \frac{ds}{v(\mathbf{x})}
\]
unde \(v(\mathbf{x})\) este viteza undei la poziția \(\mathbf{x}\), iar \(ds\) este elementul lungimii traiectoriei.
b) Problema inversă
Reversul este adevărat: avem date observaționale (timpi de deplasare, reziduuri ale timpului de deplasare, dispersii ale undelor de suprafață sau forme de undă) și dorim să estimăm modelul de viteză cel mai potrivit. Problema inversă este, în general, neunică și greșit pusă: mai multe modele se pot potrivi la fel de bine, iar datele conțin incertitudine/zgomot. Prin urmare, este necesară regularizarea, cum ar fi netezirea, amortizarea sau constrângerile geologice.
În practică, inversiunea este adesea efectuată într-o manieră liniarizată: modelul inițial este ușor perturbat, iar apoi diferența de timp de deplasare este calculată ca:
\[
ΔT ∫⁻¹ ...
\]
unde Δs = Δ(1/v) este perturbația de lentoare. Aceste ecuații sunt aranjate într-un sistem liniar mare (d = Gm), apoi rezolvate prin metoda celor mai mici pătrate amortizate sau prin variantele acesteia.
3. Tipuri de tomografie seismică
1) Tomografie în timp de călătorie
Cea mai comună metodă folosește sosirile în faza P și S de la numeroase cutremure. Datele sunt reziduuri ale timpului de călătorie față de un model de referință (de exemplu, IASP91 sau ak135 la scară globală). Potrivită pentru modelarea 3D a vitezei în scoarță și mantaua superioară, în special dacă rețeaua de stații este densă.
2) Tomografie cu unde de suprafață
Folosind dispersia (schimbarea vitezei de fază/grup în funcție de perioadă) undelor Rayleigh/Love. Sensibilă la structuri de la superficiale la intermediare, această tomografie este foarte utilă pentru cartografierea grosimii litosferei, a zonelor de viteză redusă (LVZ) și a variațiilor de temperatură.
3) Tomografie teleseismică (tomografie teleseismică)
Folosind cutremure îndepărtate (teleseismice) ale căror unde călătoresc prin manta și sunt apoi înregistrate de o rețea locală. Avantajul: au surse multiple și provin din diverse direcții, contribuind astfel la „iluminarea” volumului subteran al zonei studiate, de exemplu sub vulcani sau zone de subducție.
4) Tomografie cu formă de undă / inversiune completă a formei de undă (FWI)
Utilizarea formei de undă complete, nu doar a timpului de sosire, produce teoretic o rezoluție ridicată, deoarece utilizează pe deplin informațiile despre amplitudine și fază, dar necesită cerințe computaționale mari și un model inițial bun pentru a evita capcanele minime locale.
4. Etapele generale de procesare și inversare
1. Achiziția de date
Colecție de înregistrări de la seismometre permanente sau temporare. La scară de explorare, sursa poate fi o explozie sau un vibroseism; la scară regională-globală, sursa principală este un cutremur.
2. Identificarea și selectarea fazelor
Determinarea timpului de sosire al undelor P, S sau de suprafață. Calitatea selecției determină calitatea modelului.
3. Corecție inițială și modelare
Corecția timpului (abaterea ceasului), corecția elevației și selectarea modelului de referință. Modelul inițial poate fi 1D sau 3D simplu.
4. Trasarea razelor / simularea undelor
Calculați traiectoria razelor sau formelor de undă sintetice pentru a construi matricea de sensibilitate.
5. Inversiune și regularizare
Rezolvarea sistemului pentru a obține un model de viteză. Regularizarea este aleasă în funcție de obiectiv: dacă se doresc a evidenția anomalii bruște sau tendințe line.
6. Evaluarea rezoluției și testarea fiabilității
De exemplu, cu testul de tip tablă de șah, testul spike sau analiza funcției de covarianță/împrăștiere punctuală pentru a vedea care părți sunt cu adevărat rezolvate de date.
5. Interpretarea modelului de viteză
Rezultatul tomografiei este de obicei o hartă a anomaliilor de viteză relativă: zonele de mare viteză sunt adesea interpretate ca roci mai reci, mai dense sau mai rigide (de exemplu, plăci în subducție). Zonele de mică viteză pot indica temperaturi ridicate, roci alterate, zone de fractură umplute cu fluid sau topire parțială - adesea asociate cu vulcani activi sau astenosferă fierbinte.
Totuși, interpretarea trebuie să fie precaută, deoarece viteza este influențată de mulți factori. În mod ideal, tomografia seismică este combinată cu alte date, cum ar fi gravitația, magnetotelurica (MT), geodezia, geologia suprafeței și petrolologia.
6. Aplicații ale tomografiei seismice
a) Studiul zonelor de subducție și al dinamicii plăcilor
Tomografia globală și regională poate cartografia placa de subducție a mantalei, inclusiv geometria, adâncimea și segmentarea acesteia. Aceste informații sunt cruciale pentru înțelegerea surselor cutremurelor mari, a mecanismelor de cuplare a plăcilor și a evoluției tectonice a unei regiuni.
b) Sisteme vulcanice și atenuarea dezastrelor
În zonele vulcanice, tomografia poate identifica zone cu viteză redusă asociate cu camere magmatice, căi ascendente ale fluidelor sau roci alterate. Prin monitorizare regulată, schimbările de viteză (tomografie time-lapse) pot oferi indicatori ai mișcării fluidelor/magmei relevanți pentru avertizarea timpurie a erupției.
c) Explorarea energiei și a resurselor
La scară superficială a crustei terestre, tomografia este utilizată în explorarea hidrocarburilor, geotermală și minieră. Variațiile de viteză ajută la cartografierea litologiei, a structurilor faliilor, a zonelor de fractură și a rezervoarelor. În geotermală, integrarea tomografiei cu MT este adesea eficientă: tomografia oferă informații elastice, MT oferă informații despre conductivitatea fluidelor.
d) Caracterizarea faliilor active și a hazardelor cutremure
Tomografia locală poate dezvălui zone slabe, zone de deteriorare și eterogenitate în jurul faliilor. Acest lucru ajută la înțelegerea segmentării faliilor, a blocării potențiale și a variațiilor de viteză care influențează amplificarea undelor (efecte de amplasament).
e) Structura regională a crustei și litosferei
Folosind tomografia cu unde de suprafață, cercetătorii pot cartografia grosimea litosferei, limita Moho și anomaliile astenosferice. Rezultatele stau la baza modelelor geodinamice, inclusiv formarea bazinelor, orogeneza și evoluția continentală.
7. Limitări și provocări
Tomografia seismică depinde în mare măsură de distribuția surselor și a stațiilor. Zonele cu puține cutremure sau stații rare vor avea o iluminare slabă, rezultând o rezoluție scăzută. În plus, ipotezele teoriei razelor pot fi mai puțin precise pentru o eterogenitate puternică sau anumite frecvențe. Zgomotul, erorile de selectare și incertitudinea locației hipocentrului pot influența, de asemenea, rezultatele inversiunii. Prin urmare, evaluarea rezoluției și integrarea multi-metodă sunt cruciale pentru a evita suprainterpretarea.
Închidere
Tomografia seismică este un instrument cheie pentru investigarea structurii interne a Pământului de la scară locală la scară globală. Baza sa teoretică se bazează pe relația dintre proprietățile elastice ale rocilor și propagarea undelor, precum și pe rezolvarea problemelor inverse care necesită regularizare și testare a rezoluției. Aplicațiile sale sunt largi: de la cartografierea plăcilor de subducție și a sistemelor de magmă vulcanică până la explorarea energiei și evaluarea riscului de cutremure. Odată cu progresele înregistrate în rețelele seismice, calcul și metodele de inversie (inclusiv inversiunea completă a formei de undă), tomografia seismică este din ce în ce mai capabilă să producă imagini mai clare și mai informative, atât în scopuri științifice, cât și în scopuri de atenuare a dezastrelor.
Dacă doriți, pot adapta acest articol într-o versiune mai academică (cu citări și bibliografie) sau mă pot concentra pe una dintre aplicații (de exemplu, tomografia pentru vulcani sau pentru energia geotermală).
