MT Datuen Prozesatzeko eta Interpretatzeko Teknikak
Pendahuluan
Magnetoteluriko Metodoa (MT) teknika geofisiko pasibo bat da, Lurraren eremu elektriko eta magnetikoen aldaketa naturalak erabiltzen dituena lurpeko erresistentziaren banaketa mapatzeko. Erresistentzia oso lotuta dago arroka motarekin, fluido edukiarekin, tenperaturarekin eta egitura geologikoekin, hala nola failak edo alterazio eremuak. Hori dela eta, MT oso erabilia da esplorazio geotermikorako, mineralizaziorako, arro sedimentarioetarako (hidrokarburoak) eta ikerketa tektoniko eta lurrazalerako.
Hala ere, MT datu gordinak gutxitan erabil daitezke zuzenean. Seinale-zarata erlazioa hobetzeko hainbat prozesatzeko urrats behar dira, eta ondoren, geologikoki koherentea den erresistentzia-eredu batean oinarritutako interpretazioa. Artikulu honek MT datuak prozesatzeko eta interpretatzeko teknikak aztertzen ditu, eskuratzetik eta kalitate-kontroletik hasi eta inpedantzia-prozesamenduraino, inbertsio eta interpretazio geologikoraino.
-
MT Datu-basea: E eremua, H eremua eta inpedantzia-tentsorea
MT inkestetan, sentsoreek eremu elektrikoaren osagaiak (Ex, Ey) eta eremu magnetikoaren osagaiak (Hx, Hy, batzuetan Hz) denboraren funtzio gisa erregistratzen dituzte. Maiztasun-domeinuko eremu elektriko eta magnetikoen arteko erlaziotik, inpedantzia-tentsorea lortzen da:
\[
\begin{bmatrix} E_x \\ E_y \end{bmatrix}
=
\begin{bmatrix} Z_{xx} eta Z_{xy} \\ Z_{yx} eta Z_{yy} \end{bmatrix}
\begin{bmatrix} H_x \\ H_y \end{bmatrix}
\]
Tentsore-elementu hauek interpretazio-parametro gakoetan eratortzen dira, hala nola, itxurazko erresistentzia (ρa) eta fasea (φ). Oro har, datu onek joera leun eta koherenteak erakusten dituzte ρa eta φ-n maiztasunetan zehar, eta errore-barra errealistak dituzte.
-
1. etapa: Kalitate Kontrola (KK) eta Aurre-prozesamendua
MT prozesamendu ona QC-rekin hasten da, landa-datuekin hasita. Etapa honek honako hauek ditu barne:
1. Denbora-serieen ikuskapena: puntak, desbideratzea, sentsoreen saturazioa edo aldizkako asaldurak bilatzea.
2. Zarata kulturala: zarata-iturri ohikoenak 50/60 Hz-ko linea elektrikoak, trenak, industria-jarduerak, hesi elektrikoak eta telekomunikazioak dira.
3. Elektrodoen eta lurrerako kontaktuen egoera: kontaktu-erresistentzia handiak Ex/Ey-ren kalitatea okerrera egingo du.
4. Sentsorearen orientazioa eta posizioa: sentsorearen azimut erroreek interpretazio distortsioak sor ditzakete, batez ere norabide-egituren azterketetan.
Fase honetan, datu-segmentu txarrak normalean moztu egiten dira, desplazamendu-zuzenketa egiten da eta denbora-sinkronizazioa egiten da urruneko erreferentzia-estazio bat erabiltzen bada.
-
2. etapa: Maiztasun-domeinuaren eraldaketa eta espektro-estimazioa
MT-k E eta H arteko erlazioa maiztasun-domeinuan aztertzen duenez, denbora-seriea Fourier Transformatua bezalako teknika espektralak erabiliz eraldatzen da. Ondoren, datuak leihoetan (segmentuetan) banatzen dira egonkortasun estatistikoa bermatzeko. Ohiko teknika batzuk hauek dira:
– Leihoen murrizketa eta konikotasuna (adibidez, Hanning) espektro-ihesa murrizteko.
– Leihoen arteko batez bestekoa egitea espektroaren estimazio sendoagoa lortzeko.
– E eta H-ren arteko erlazioa ezartzeko gurutzadura-potentzia espektroak.
Helburua inpedantzia-estimazio zehatzak lortzea da maiztasun-tarte zabal batean, maiztasun altuetatik (azalekoetatik) baxuetara (sakonetara).
-
3. etapa: Prozesamendu sendoa eta urruneko erreferentzia
MT-aren erronketako bat kanal jakin batean askotan korrelazionatzen den zarata da. Hori dela eta, prozesatzeko teknika sendoak garatu dira outlierren eta datu-segmentu txarren eragina murrizteko. Metodo sendoek pisu iteratiboa erabiltzen dute datu inkoherenteen ekarpena kentzeko.
Beste teknika oso garrantzitsu bat urruneko erreferentzia (RR) da. Kontzeptua da eremu magnetikoa tokiko zarata-iturritik urrun dagoen estazio batean neurtzen dela. Estazio nagusiko E eta H RR estazioko H-rekin korrelazionatuz, tokiko zaratak H-n duen eragina murriztu daiteke. RR bereziki eraginkorra da gune zaratatsuetan, adibidez, bizitegi-eremuen edo azpiegituren ondoan.
Etapa honetako emaitza normalean hau da:
– ρa eta φ kurbak Zxy eta Zyx-erako
– Errore-barra (desbideratze estandarra)
– Lehen mailako seinaleen arteko koherentzia
-
4. etapa: Zarata kentzea eta distortsioaren kudeaketa (desplazamendu estatikoa)
Zarataz gain, MT askotan desplazamendu estatikoak eragiten du, hau da, erresistentzia itxurazko kurban desplazamendu bertikala heterogeneotasun txikiengatik (adibidez, buztin geruza meheak, legarra edo elektrodoen baldintzak). Desplazamendu estatikoak ez du fasea nabarmen aldatzen, baina ρa gora edo behera biderkatzen du.
Txanda estatikoen kudeaketa honela egin daiteke:
1. TDEM/CSAMT datuetan oinarritutako zuzenketa erresistentzia txikiko kontrol gisa.
2. MT-TDEM-ren inbertsio bateratua sakonera txikiko eredua lotuagoa izan dadin.
3. Aldaketa-inbertsio bidezko ikuspegi sendoa, adibidez, estazio bakoitzeko aldaketa-parametroak ahalbidetuz.
Horrez gain, maiztasun-abstentzioen detekzioa ere egiten da: kurban nabarmen desbideratzen diren puntu batzuk baztertu edo errore handiagoa ematen zaie normalean.
-
5. etapa: Dimentsioen eta kolpeen azterketa
Inbertsioa egin aurretik, garrantzitsua da lurpeko egitura 1D, 2D edo 3D den zehaztea. Horrek inbertsio-metodoaren aukeraketa eta interpretazioa eragiten du. Ohiko analisiak hauek dira:
– 3D maila ebaluatzeko okertze-parametroak (adibidez, Bahr-en okertzea).
– Fase-tentsorea egituraren norabide nagusia ikusteko, desplazamendu estatikoak eraginik izan gabe.
– 2D egituraren kolpe-norabide nagusia zehazteko kolpe-analisia.
Datuek 2D ezaugarri sendoak erakusten badituzte, normalean inpedantzia-tentsorearen errotazio zuzendua egiten da osagai nagusiak (TE/TM) argitzeko. 3D ezaugarriak sendoak badira, 3D inbertsioa hobesten da.
-
6. urratsa: MT datuen alderantzikatzea (1D, 2D, 3D)
MT-ren interpretazio kuantitatiboa, oro har, inbertsioaren bidez egiten da, eta horrek datuen erantzunari hobekien egokitzen zaion erresistentzia-eredua aurkitzea dakar. MT inbertsioa ez-lineala eta gaizki planteatuta dago, eta erregularizazioa behar du eredua "basati" bihurtzea saihesteko. Inbertsio motak:
1. 1D inbertsioa: geruzadun eremuetarako egokia (adibidez, sedimentu horizontaletarako). Azkarra, baina mugatua.
2. 2D inbertsioa: faila, graben edo norabideko sistema geotermikoen antzeko egitura luzangetarako egokia.
3. 3D inbertsioa: geologia konplexurako errealistenetakoa, baina datu trinkoak, kalkulu handiak eta QC zorrotza behar ditu.
Alderantzizko helburu-funtzioak normalean datuen desegokitzapena eta ereduaren zimurtasuna barne hartzen ditu:
\[
\Phi = \Phi_d + \lambda \Phi_m
\]
non λ erregularizazio-parametroa den. λ-ren aukera funtsezkoa da: txikiegia bada, eredua zakarregia (gehiegi egokitzea) bihurtzen du, eta handiegia bada, eredua leunegia eta ezaugarri geologiko garrantzitsuak galtzen ditu.
-
7. urratsa: Inbertsio-emaitzen eta sentikortasunaren ebaluazioa
Alderantzikatzearen ondoren, eredua ebaluatu behar da, ez bakarrik bisualki. Ebaluazio honek honako hauek barne hartzen ditu:
– Desegokitzapena (RMS): helburuarekin bat datorren ala ez (adibidez, RMS ~ 1–2, errorearen definizioaren arabera).
– Datuen kurben eta modeloaren erantzunaren konparaketa estazio bakoitzean.
– Bereizmen-proba: adibidez, dama-taularen proba edo sentikortasun-mapa.
– Aurreko ereduen eragina: saiatu hasierako hainbat eredu emaitzen egonkortasuna ikusteko.
Parametro txiki bat aldatzen denean eredua izugarri aldatzen bada, interpretazioa zainduagoa izan behar dela esan nahi du eta datu gehigarriak behar izan ditzakeela.
-
Interpretazio geologikoa: Erresistentzia lurpeko sistemekin erlazionatzea
Erresistentzia ez da zuzeneko "arroka mota" bat, baizik eta erantzun fisiko konposatu bat. Hala ere, eredu orokor batzuk jarraibide gisa balio ohi dute:
– Erresistentzia baxua (eroalea): buztina, alterazio hidrotermalak (buztin-geruza), gatz-fluidoaren eremua, grafitoa edo sulfuro mineralak.
– Erresistentzia handia (erresistiboa): arroka igneo masiboak, arroka lehorrak, silizeztatutako eremuak edo oinarri kristalinoa.
Geotermiaren esplorazioan, adibidez, eredu klasikoek sarritan erakusten dute:
1. Buztinezko estalki eroalea goialdean
2. Urtegiak erresistentzia handiagoa du azpian
3. Egiturak kontrolatutako goranzko fluxu-eremua (faila)
4. Litologiaren eta tenperaturaren arabera batzuetan erresistentziazko edo anomalia konplexu gisa agertzen diren bero-iturriak.
Interpretazio ona beti lotuta dago beste datu batzuekin: gainazaleko geologia, agerpen termikoak, geokimika, grabitatea, sismika edo putzuen datuak.
-
Itxiera
MT datuak prozesatzeko eta interpretatzeko teknikak elkarri lotuta daude: denbora-serieen kalitate-kontroletik (QC), espektro-estimaziotik, erreferentzia-prozesamendu sendo/urrunetik, desplazamendu estatikoaren zuzenketatik, dimentsio-analisitik, 2D/3D inbertsio eta bereizmen-ebaluazioraino. MTren arrakasta ez da soilik inbertsio-softwarearen mende, baita neurketen kalitatearen, zarataren ulermenaren eta testuinguru geologikoarekin integrazioaren mende ere.
Bere lan-fluxu diziplinatuari eta datu anitzetan oinarritutako interpretazioari esker, MT tresna indartsua da lurpeko erresistentzia-egiturak mapatzeko eta esplorazio eta geozientzia ikerketa konplexuetan erabakiak hartzen laguntzeko.
