MT duomenų apdorojimo ir interpretavimo metodai
Pendahuluanas
Magnetotelurinis metodas (MT) yra pasyvus geofizinis metodas, kuris naudoja natūralius Žemės elektrinių ir magnetinių laukų svyravimus, kad būtų galima nustatyti požeminio varžos pasiskirstymą. Varža yra glaudžiai susijusi su uolienų tipu, skysčių kiekiu, temperatūra ir geologinėmis struktūromis, tokiomis kaip lūžiai ar pakitimų zonos. Todėl MT plačiai naudojamas geoterminiams tyrimams, mineralizacijai, nuosėdinių baseinų (angliavandenilių) tyrimams ir tektoniniams bei plutos tyrimams.
Tačiau neapdoroti MT duomenys retai kada yra tiesiogiai panaudojami. Norint pagerinti signalo ir triukšmo santykį, reikia atlikti keletą apdorojimo veiksmų, o po to interpretuoti duomenis remiantis geologiškai suderintu varžos modeliu. Šiame straipsnyje aptariami MT duomenų apdorojimo ir interpretavimo metodai – nuo duomenų gavimo ir kokybės kontrolės iki impedanso apdorojimo, inversijos ir geologinės interpretacijos.
-
Duomenų bazė MT: E lauko, H lauko ir impedanso tenzorius
MT tyrimuose jutikliai registruoja elektrinio lauko dedamąsias (Ex, Ey) ir magnetinio lauko dedamąsias (Hx, Hy, kartais Hz) kaip laiko funkciją. Iš elektrinio ir magnetinio laukų santykio dažnio srityje gaunamas impedanso tenzorius:
\[
\begin{bmatrix} E_x \\ E_y \end{bmatrix}
=
\begin{bmatrix} Z_{xx} ir Z_{xy} \\ Z_{yx} ir Z_{yy} \end{bmatrix}
\begin{bmatrix} H_x \\ H_y \end{bmatrix}
\]
Šie tenzoriniai elementai tada išvedami į pagrindinius interpretavimo parametrus, tokius kaip tariamasis varžos koeficientas (ρa) ir fazė (φ). Paprastai geri duomenys rodo sklandžius ir nuoseklius ρa ir φ pokyčius skirtinguose dažniuose ir turi realistiškas paklaidų juostas.
-
1 etapas: kokybės kontrolė (KK) ir išankstinis apdorojimas
Geras MT apdorojimas prasideda nuo kokybės kontrolės (QC), nes lauko duomenysŠis etapas apima:
1. Laiko eilučių patikrinimas: ieškoma šuolių, dreifo, jutiklio prisotinimo ar periodinių trikdžių.
2. Kultūrinis triukšmas: dažni triukšmo šaltiniai yra 50/60 Hz elektros linijos, traukiniai, pramoninė veikla, elektrinės tvoros ir telekomunikacijos.
3. Elektrodų ir įžeminimo kontaktų būklė: didelė kontaktinė varža pablogins Ex/Ey kokybę.
4. Jutiklio orientacija ir padėtis: jutiklio azimuto paklaidos gali sukelti interpretavimo iškraipymus, ypač kryptinių struktūrų tyrimuose.
Šiame etape blogi duomenų segmentai paprastai yra nukerpami, atliekama poslinkio korekcija ir laiko sinchronizavimas, jei naudojama nuotolinė referencinė stotis.
-
2 etapas: Dažnių srities transformacija ir spektrinis įvertinimas
Kadangi MT analizuoja E ir H santykį dažnių srityje, laiko eilutė transformuojama naudojant spektrinius metodus, tokius kaip Furjė transformacija. Tada duomenys padalijami į langus (segmentus), siekiant užtikrinti statistinį stabilumą. Kai kurie dažniausiai naudojami metodai:
– Langų ir siaurėjančių diapazonų taikymas (pvz., Hanningo metodas), siekiant sumažinti spektrinį nuotėkį.
– Vidurkinimas tarp langų, siekiant gauti patikimesnį spektro įvertinimą.
– Kryžminės galios spektrai E ir H ryšiui nustatyti.
Tikslas yra gauti tikslius impedanso įvertinimus plačiame dažnių diapazone – nuo aukštų (negilių) iki žemų (gilių) dažnių.
-
3 etapas: patikimas apdorojimas ir nuotolinė nuoroda
Vienas iš mašininio apdirbimo (MT) iššūkių yra dažnai koreliuojantis triukšmas konkrečiame kanale. Todėl buvo sukurti patikimi apdorojimo metodai, skirti sumažinti išskirtinių reikšmių ir blogų duomenų segmentų įtaką. Patikimi metodai naudoja iteracinį svertinį vertinimą, kad sumažintų nenuoseklių duomenų įtaką.
Kita labai svarbi technika yra nuotolinis referencinis nustatymas (RR). Koncepcija yra ta, kad magnetinis laukas matuojamas stotyje, esančioje toli nuo vietinio triukšmo šaltinio. Koreliuojant E ir H pagrindinėje stotyje su H RR stotyje, galima sumažinti vietinio triukšmo poveikį H. RR ypač veiksmingas triukšmingose vietose, pavyzdžiui, šalia gyvenamųjų rajonų ar infrastruktūros.
Šio etapo rezultatas paprastai yra toks:
– ρa ir φ kreivės Zxy ir Zyx
– Paklaidos juosta (standartinis nuokrypis)
– Pirminių signalų suderinamumas
-
4 etapas: triukšmo šalinimas ir iškraipymų valdymas (statinis poslinkis)
Be triukšmo, MT dažnai veikia statinis poslinkis – vertikalus tariamosios varžos kreivės poslinkis dėl nedidelio nevienalytiškumo (pvz., plonų molio sluoksnių, žvyro ar elektrodų sąlygų). Statinis poslinkis reikšmingai nepakeičia fazės, bet daugybiškai paslenka ρa aukštyn arba žemyn.
Statinį poslinkį galima valdyti šiais būdais:
1. TDEM/CSAMT duomenimis pagrįsta korekcija kaip seklios varžos valdymas.
2. MT-TDEM jungtinė inversija, siekiant labiau apriboti sekliojo modelio ribas.
3. Patikimas inversijos metodas poslinkiui, pavyzdžiui, leidžiant poslinkio parametrus kiekvienai stočiai.
Be to, atliekamas ir dažnio nuokrypių aptikimas: tam tikri kreivės taškai, kurie smarkiai nukrypsta, paprastai atmetami arba jiems suteikiama didesnė paklaida.
-
5 etapas: Matmenų ir smūgių analizė
Prieš inversiją svarbu nustatyti, ar požeminė struktūra yra vienmatė, dvimatė ar trimatė. Tai turi įtakos inversijos metodo pasirinkimui ir interpretacijai. Įprastos analizės apima:
– Pasvirimo parametrai (pvz., Bahro pasvirimas), skirti 3D lygiui įvertinti.
– Fazės tenzorius, skirtas matyti dominuojančią struktūros kryptį, nepaveikiamą statinio poslinkio.
– Smūgio analizė, siekiant nustatyti dominuojančią 2D struktūros smūgio kryptį.
Jei duomenys pasižymi stipriomis 2D charakteristikomis, paprastai atliekama smūgio nukreipta impedanso tenzoriaus rotacija, siekiant išsiaiškinti pagrindinius komponentus (TE/TM). Jei 3D charakteristikos yra stiprios, pirmenybė teikiama 3D inversijai.
-
6 veiksmas: MT duomenų inversija (1D, 2D, 3D)
Kiekybinis MT interpretavimas paprastai atliekamas inversijos būdu, kurios metu randamas varžos modelis, geriausiai atitinkantis duomenų atsaką. MT inversija yra netiesinė ir netinkamai išplanuota, todėl reikia reguliarizuoti, kad modelis netaptų „laukinis“. Inversijos tipai:
1. Vienmatės inversijos metodas: tinka sluoksniuotoms sritims (pvz., horizontalioms nuosėdoms). Greitas, bet ribotas.
2. 2D inversija: tinka pailgoms struktūroms, tokioms kaip lūžiai, grabenai arba kryptingos geoterminės sistemos.
3. 3D inversija: realiausia sudėtingai geologijai, tačiau reikalauja tankių duomenų, didelių skaičiavimų ir griežtos kokybės kontrolės.
Atvirkštinė tikslo funkcija paprastai apima duomenų neatitikimą ir modelio šiurkštumą:
\[
\Phi = \Phi_d + \lambda \Phi_m
\]
kur λ yra reguliarizacijos parametras. λ pasirinkimas yra labai svarbus: per mažas modelis tampa pernelyg grubus (per didelis), per didelis – pernelyg glotnus ir praleidžia svarbius geologinius ypatumus.
-
7 veiksmas: inversijos rezultatų ir jautrumo įvertinimas
Po inversijos modelis turi būti įvertintas ne tik vizualiai. Šis įvertinimas apima:
– Neatitikimas (RMS): ar jis atitinka tikslinį rodiklį (pvz., RMS ~ 1–2, priklausomai nuo paklaidos apibrėžimo).
– Duomenų kreivių ir modelio atsako palyginimas kiekvienoje stotyje.
– Skyros testas: pavyzdžiui, šachmatų lentos testas arba jautrumo žemėlapis.
– Ankstesnių modelių įtaka: išbandykite kelis pradinius modelius, kad pamatytumėte rezultatų stabilumą.
Jei modelis smarkiai pasikeičia pakeitus nedidelį parametrą, tai reiškia, kad interpretacija turėtų būti atsargesnė ir gali prireikti papildomų duomenų.
-
Geologinis aiškinimas: varžos susiejimas su požeminėmis sistemomis
Varža nėra tiesioginis „uolienų tipas“, o sudėtinis fizikinis atsakas. Tačiau kai kurie bendri modeliai dažnai naudojami kaip gairės:
– Mažas varžos koeficientas (laidumas): molis, hidroterminis pakitimas (molio kepurė), druskos skysčio zona, grafitas arba sulfidų mineralai.
– Didelė varža (varžinė): masyvios magminės uolienos, sausos uolienos, silicuotos zonos arba kristalinis pamatas.
Pavyzdžiui, atliekant geoterminius tyrimus, klasikiniai modeliai dažnai rodo:
1. Laidus molinis dangtelis viršuje
2. Rezervuaro varža po apačia yra didesnė
3. Konstrukcijos kontroliuojama aukštyn srauto zona (gedimas)
4. Šilumos šaltiniai, kurie, priklausomai nuo litologijos ir temperatūros, kartais atrodo kaip varžinės arba sudėtingos anomalijos.
Geras interpretavimas visada susijęs su kitais duomenimis: paviršiaus geologija, terminėmis apraiškomis, geochemija, gravitacija, seisminiais ar gręžinių duomenimis.
-
Uždarymas
MT duomenų apdorojimo ir interpretavimo metodai yra tarpusavyje susiję: nuo laiko eilučių kokybės kontrolės (KK), spektrinio įvertinimo, patikimo/nuotolinio etaloninio apdorojimo, statinio poslinkio korekcijos, dimensijų analizės iki 2D/3D inversijos ir skiriamosios gebos įvertinimo. MT sėkmė priklauso ne tik nuo inversijos programinės įrangos, bet ir nuo matavimų kokybės, triukšmo supratimo ir integravimo į geologinį kontekstą.
Dėl disciplinuoto darbo eigos ir daugiaduomenės interpretavimo, MT yra galingas įrankis požeminių varžos struktūrų kartografavimui ir sprendimų priėmimo palengvinimui atliekant sudėtingus tyrinėjimus ir geomokslinius tyrimus.
