Tomografia seismică de refracție: ce este, de ce se face și cum se interpretează (cu Easy Refract + Tomo)

O prospecțiune seismică de refracție interpretată cu metodele clasice dă un model stratificat: doi sau trei refractori, fiecare cu viteza și adâncimea sa sub fiecare geofon. Funcționează bine când terenul este într-adevăr stratificat. Dar terenul real are depresiuni ale substratului, lentile de material diferit, zone de alterare care se îngroașă, fâșii fracturate de-a lungul unei falii. Tomografia seismică de refracție există tocmai pentru a vedea aceste lucruri: din aceleași prime sosiri obține o secțiune de viteză continuă, celulă cu celulă, în locul unui model stratificat. În acest articol explicăm ce este, de ce se face, cum se achiziționează datele, cum funcționează calculul și cum se citește rezultatul, cu exemplul noii ediții Easy Refract + Tomo de la GeoStru.

Ce este tomografia seismică de refracție

Tomografia seismică de refracție (în engleză seismic refraction tomography, SRT) este o metodă de investigare geofizică ce reconstruiește distribuția vitezei undelor seismice în teren de-a lungul unei secțiuni verticale. Datele sunt cele ale unei prospecțiuni de refracție obișnuite: o linie de geofoni (dispozitivul), o serie de lovituri (cu baros, greutate în cădere liberă sau alte surse) și, pentru fiecare lovitură și fiecare geofon, timpul primei sosiri a undei.

Diferența stă în interpretare. Terenul este împărțit într-o rețea de celule, fiecare cu viteza sa. Un algoritm caută distribuția de viteze care, propagând undele de la surse la geofoni, reproduce timpii măsurați în limita erorii de citire. Rezultatul este o imagine color a vitezei — de obicei a undelor P (Vp) — care urmărește topografia și arată variații atât pe verticală, cât și pe orizontală.

Termenul „tomografie” vine din medicină: așa cum tomograful (CT) reconstruiește interiorul corpului din multe măsurători cu raze X pe trasee diferite, tomografia seismică reconstruiește terenul din mulți timpi de parcurs de-a lungul unor raze seismice care se intersectează.

De ce se face: limitele modelelor stratificate

Metodele clasice de interpretare a refracției — metoda timpilor de intercepție, metoda reciprocă și metoda reciprocă generalizată (GRM) a lui Palmer — sunt rapide, robuste și rămân standardul pentru multe studii geotehnice. Se bazează însă pe câteva ipoteze:

• terenul este format din straturi separate de interfețe nete;
• viteza crește cu adâncimea de la un strat la următorul;
• în interiorul fiecărui strat viteza este constantă sau variază puțin lateral.

Când aceste ipoteze nu sunt valabile, modelul stratificat simplifică prea mult. Un strat mai lent decât cel de deasupra (inversiune de viteză) nu produce un refractor și rămâne invizibil; un strat subțire poate să nu apară niciodată ca primă sosire (strat ascuns); o creștere treptată a vitezei cu adâncimea, tipică pământurilor care se îndesează sau rocii care trece de la alterată la sănătoasă, nu are o interfață de identificat. Și mai ales, variațiile laterale — o depresiune a substratului, o lentilă, o zonă fracturată — sunt „întinse” pe un strat cu viteză medie.

Tomografia nu are nevoie de aceste ipoteze: nu întreabă câte straturi există și unde se află. De aceea se folosește când ne așteptăm la un teren eterogen sau pentru a verifica interpretarea stratificată cu o metodă independentă.

Aplicații tipice

• Alunecări de teren și versanți instabili: grosimea depozitelor deluviale sau a corpului alunecării, forma substratului stabil, zone decomprimate. Sunt date care intră direct în modelul geologic pentru analizele de stabilitate.
• Adâncimea și morfologia substratului pentru fundații, lucrări de susținere și excavații: depresiuni, trepte, albii vechi.
• Excavabilitatea pământurilor și rocilor: tabelele producătorilor de utilaje de terasamente corelează Vp cu posibilitatea de a scarifica (ripper) sau cu nevoia de ciocan hidraulic ori explozivi.
• Falii și zone fracturate, care apar ca zone cu viteză redusă în substrat.
• Goluri, doline și zone de dizolvare, ca anomalii de viteză mică (cu precauțiile privind rezoluția prezentate mai jos).
• Portaluri de tuneluri, baraje și diguri, depozite de deșeuri și situri contaminate, unde este nevoie de o secțiune continuă de-a lungul lucrării.
• Hidrogeologie: Vp crește la circa 1500 m/s în pământurile saturate, iar tomografia ajută la localizarea nivelului apei subterane în depozitele necoezive.

Cum se achiziționează datele pentru o tomografie

Echipamentul este cel al refracției: seismograf multicanal, geofoni verticali (pentru undele P) sau orizontali (pentru undele SH), sursă. Ce se schimbă cel mai mult este numărul de lovituri. Tomografia are nevoie de multe raze care traversează terenul din direcții diferite, deci:

• lovituri dese: pe lângă loviturile de capăt și centrală ale refracției clasice, o lovitură la câțiva geofoni de-a lungul dispozitivului (la fiecare 2–4 geofoni este o regulă bună);
• lovituri exterioare dispozitivului, dincolo de geofonii de capăt: fără ele, adâncimea substratului la marginile secțiunii este controlată de puține raze;
• cotele geofonilor măsurate cu atenție: topografia intră în calcul;
• prime sosiri citite cu precizie pe toate loviturile, cu o estimare a erorii de citire (picking), de regulă 0,5–2 ms.

Adâncimea de investigare depinde de lungimea dispozitivului și de contrastul de viteză; ca ordin de mărime este între o treime și o pătrime din lungime. Un dispozitiv de 24 de geofoni la 3 m (69 m) investighează de regulă primii 15–25 m. Rezoluția depinde de distanța dintre geofoni și de densitatea razelor: anomaliile mai mici decât distanța dintre geofoni nu se pot rezolva.

O veste bună: cine are o arhivă de prospecțiuni de refracție cu mai multe lovituri le poate adesea reprocesa ca tomografie fără să revină pe teren. Cu 5 lovituri și 12 geofoni rezultatul va fi neted și puțin adânc, dar arată deja variațiile laterale pe care modelul stratificat le ascunde.

Cum funcționează calculul

O tomografie de refracție alternează două probleme.

Problema directă: dat fiind un model de viteze, se calculează timpii de sosire la fiecare geofon. Timpul de parcurs T satisface ecuația eikonală |∇T| = 1/v, care spune că frontul de undă avansează peste tot cu viteza locală. Metoda cea mai folosită pentru a o rezolva pe o rețea este Fast Marching (Sethian, 1996): propagă frontul de la sursă celulă cu celulă, în ordinea timpului crescător, și găsește singură drumul cel mai rapid, fie că este o undă directă, refractată sau care ocolește o anomalie.

Problema inversă: se corectează modelul pentru a reduce diferența dintre timpii calculați și cei măsurați. Problema se liniarizează în jurul modelului curent (metoda Gauss-Newton) și se rezolvă un sistem mare, rar, de cele mai mici pătrate, cu algoritmi iterativi precum LSQR (Paige și Saunders, 1982). Deoarece datele sunt mai puține decât necunoscutele și au erori, se adaugă o regularizare: o constrângere de netezime care, dintre toate modelele compatibile cu datele, îl preferă pe cel mai simplu.

Cele două etape se repetă până când abaterea medie pătratică (RMS) dintre timpii măsurați și cei calculați coboară la nivelul erorii de picking. A merge mai departe nu are sens: ar însemna ajustarea modelului la zgomotul citirilor. Împreună cu secțiunea se obține acoperirea cu raze, care arată în ce celule viteza este cu adevărat constrânsă de date.

Cum se citește o secțiune tomografică

Trei reguli de citire evită cele mai frecvente greșeli.

1. Priviți razele. Unde razele sunt dese și se intersectează din direcții diferite, viteza este bine definită. La margini și în adâncime, unde razele sunt puține sau lipsesc, secțiunea arată mai ales modelul inițial: trebuie citită cu prudență sau ascunsă.
2. Tomografia netezește contrastele. Constrângerea de netezime transformă o trecere netă între două straturi într-o bandă de tranziție. Pentru a poziționa o interfață se folosește o izolinie de viteză aleasă cu discernământ sau comparația cu interpretarea stratificată.
3. Viteza nu este o litologie. Secțiunea spune cu ce viteză circulă undele, nu despre ce material este vorba. Trebuie calibrată cu foraje, aflorimente și geologia amplasamentului.

Pentru a citi secțiunea pe zone se folosesc clasele de viteză. Valori orientative pentru undele P:

Vp [m/s]	Interpretare orientativă
până la 800	Acoperiri, umpluturi, depozite necoezive uscate
800 – 1500	Depozite îndesate, rocă foarte alterată
1500 – 2000	Pământuri saturate, rocă alterată sau fracturată
peste 2000	Substrat stâncos

Atenție la apă: sub nivelul apei subterane Vp depășește circa 1500 m/s chiar și în pământurile necoezive, care pot intra în clasa rocii alterate. Este unul dintre motivele pentru care, în prezența apei, tomografia în unde P se completează cu o investigație în unde de forfecare.

Tomografia în unde P și clasa de teren

O precizare utilă pentru cei care redactează studii geologice și geotehnice: clasificarea terenului pentru acțiunea seismică (Eurocodul 8) se bazează pe viteza undelor de forfecare Vs, nu pe Vp. O tomografie în unde P nu o furnizează. Pentru Vs se folosesc metode ca MASW, raportul spectral HVSR sau refracția în unde SH (care poate fi și ea procesată ca tomografie, cu clase de viteză potrivite undelor S). În practică metodele se completează: tomografia în P descrie geometria terenului de-a lungul secțiunii, MASW și HVSR dau profilul Vs. Am tratat subiectul în articolul despre analiza combinată MASW și HVSR.

Tomografie sau GRM? O comparație

	Interpretare stratificată (GRM)	Tomografie seismică
Model	Straturi cu interfețe nete	Rețea de celule, viteză continuă
Variații laterale	Limitate, interfață cu interfață	Da, inclusiv treptate
Gradienți verticali	Nu (viteză constantă pe strat)	Da
Inversiuni de viteză	Nu	Parțial (zone slab acoperite de raze)
Contraste nete	Bine definite	Netezite într-o bandă de tranziție
Date necesare	Puține lovituri (chiar 2–5)	Cât mai multe lovituri, inclusiv exterioare
Rezultat tipic	Adâncimile și vitezele refractorilor	Secțiune Vp cu acoperirea cu raze

Cele două metode nu se exclud, se completează. Suprapunerea refractorilor GRM peste secțiunea tomografică este cel mai simplu mod de a vedea dacă cele două citiri concordă și unde terenul se abate de la modelul stratificat.

Easy Refract + Tomo: tomografia în programul de refracție

Easy Refract este software-ul GeoStru pentru interpretarea prospecțiunilor seismice de refracție: import de seismograme (SEG2, SEGY, SU și alte formate), citirea primelor sosiri, hodografe, metodele timpilor de intercepție, reciprocă și GRM, raport automat. Noua ediție Easy Refract + Tomo cuprinde tot Easy Refract și adaugă fila Tomografie seismică.

Punctul forte este că nu trebuie introdus nimic nou: geometria dispozitivului, cotele geofonilor și primele sosiri ale tuturor loviturilor sunt cele folosite deja pentru interpretare. Se deschide fila și se apasă Calculează.

Ce face

• Rețea care urmărește topografia, cu dimensiunea celulei și adâncimea stabilite automat (un sfert din distanța dintre geofoni, 35% din lungimea dispozitivului) sau manual.
• Timpi de parcurs cu metoda Fast Marching și inversie prin cele mai mici pătrate (Gauss-Newton, LSQR) cu o constrângere de netezime care se relaxează treptat pe parcursul calculului.
• Oprire la eroarea de picking (chi pătrat normalizat egal cu 1): programul nu ajustează modelul la zgomot.
• Trei modele inițiale: automat din date (recomandat), cu gradient sau din modelul stratificat GRM.
• Calcul în fundal cu iterația și RMS în timp real; programul rămâne utilizabil.
• Raze seismice, zone fără raze ascunse, interfețe GRM suprapuse, geofoni și surse.
• Clase de viteză cu limită, descriere și culoare personalizabile.
• Export al secțiunii ca imagine, grilă XYZ (x, cotă, adâncime, viteză, acoperire cu raze: se deschide cu Surfer sau cu o foaie de calcul) și DXF cu celule colorate, suprafață, interfețe și legendă pe layere separate.
• Capitol dedicat în raportul de calcul: metodă, parametri, calitatea ajustării, imaginea secțiunii.
• Tomografia se salvează în fișierul de proiect .rft, care rămâne lizibil și pentru versiunile anterioare ale Easy Refract.

Fluxul de lucru, în șase pași

1. Se importă seismogramele loviturilor și se citesc primele sosiri, ca la refracție.
2. Se introduc cotele geofonilor (topografia).
3. Se activează fila de tomografie seismică din meniul Proiect.
4. Se păstrează parametrii impliciți, potriviți pentru majoritatea cazurilor, și se apasă Calculează.
5. Se verifică RMS final (apropiat de eroarea de picking) și acoperirea cu raze; se compară secțiunea cu refractorii GRM.
6. Se activează clasele de viteză, se exportă în DXF sau XYZ și se tipărește raportul.

În dosarul de exemple există o prospecțiune reală cu 8 lovituri și 21 de geofoni, cu topografie, și tomografia deja calculată: RMS coboară de la circa 6 ms pentru modelul inițial la circa 1 ms. Calculul pe acest exemplu durează mai puțin de o secundă.

Un calcul validat și documentat

Înainte de lansare, motorul de calcul a fost supus unei validări cu toleranțe stabilite dinainte: compararea timpilor calculați cu soluții analitice (semispațiu omogen, două straturi, gradient liniar, refractor înclinat, suprafață înclinată), un model sintetic „treaptă a substratului cu falie” după testul orb SAGEEP 2011 (Zelt et al., 2013), cu și fără zgomot, și comparația cu un software comercial de tomografie pe aceleași date. Rezultatele, inclusiv cazurile în care toleranțele nu au fost respectate și motivul, sunt adunate într-un dosar public pe pagina Validarea codurilor de calcul GeoStru. Din validare vine și cel mai important sfat practic: faceți lovituri în afara dispozitivului.

Preț, upgrade și versiune de încercare

Easy Refract + Tomo este o ediție unică, cu licență perpetuă, la 990 € + TVA (licență Standard; pentru companii, instituții publice și revânzători sunt disponibile edițiile Professional, Enterprise și PA). Cine are deja Easy Refract poate trece la noua ediție printr-un upgrade în condiții speciale: pentru detalii scrieți biroului comercial GeoStru. Versiunea de încercare include tomografia, cu mențiunea „versiune de încercare” pe secțiune.

Manualul filei este online (în engleză): Seismic tomography – manualul Easy Refract. Pentru o privire de ansamblu asupra tuturor programelor GeoStru de geofizică și geologie, consultați ghidul celor mai bune programe software pentru geologi.

Descoperiți Easy Refract + Tomo pe pagina produsului →

Întrebări frecvente

Bibliografie

• Palmer D. (1980). The Generalized Reciprocal Method of Seismic Refraction Interpretation. Society of Exploration Geophysicists, Tulsa.
• Sethian J.A. (1996). A fast marching level set method for monotonically advancing fronts. Proceedings of the National Academy of Sciences, 93(4), 1591–1595.
• Paige C.C., Saunders M.A. (1982). LSQR: An algorithm for sparse linear equations and sparse least squares. ACM Transactions on Mathematical Software, 8(1), 43–71.
• Sheehan J.R., Doll W.E., Mandell W.A. (2005). An evaluation of methods and available software for seismic refraction tomography analysis. Journal of Environmental and Engineering Geophysics, 10(1), 21–34.
• Zelt C.A., Haines S., Powers M.H., et al. (2013). Blind test of methods for obtaining 2-D near-surface seismic velocity models from first-arrival traveltimes. Journal of Environmental and Engineering Geophysics, 18(3), 183–194.