MT-gegevensverwerkings- en interpretatietechnieken
Pendahuluan
De magnetotellurische methode (MT) is een passieve geofysische techniek die gebruikmaakt van natuurlijke variaties in de elektrische en magnetische velden van de aarde om de verspreiding van de weerstand in de ondergrond in kaart te brengen. De weerstand is nauw gerelateerd aan het gesteentetype, het vloeistofgehalte, de temperatuur en geologische structuren zoals breuken of alteratiezones. Daarom wordt MT veelvuldig gebruikt voor geothermische exploratie, mineralisatie, sedimentaire bekkens (koolwaterstoffen) en tektonische en korststudies.
Ruwe MT-gegevens zijn echter zelden direct bruikbaar. Ze vereisen een reeks verwerkingsstappen om de signaal-ruisverhouding te verbeteren, gevolgd door interpretatie op basis van een geologisch consistent resistiviteitsmodel. Dit artikel bespreekt technieken voor de verwerking en interpretatie van MT-gegevens, van acquisitie en kwaliteitscontrole tot impedantieverwerking, inversie en geologische interpretatie.
-
Dasar Data MT: Medan E, Medan H, dan Tensor Impedansi
Pada survei MT, sensor merekam komponen medan listrik (Ex, Ey) dan medan magnet (Hx, Hy, kadang Hz) sebagai fungsi waktu. Dari hubungan antara medan listrik dan medan magnet di domain frekuensi, diperoleh tensor impedansi :
\[
\begin{bmatrix} E_x \\ E_y \end{bmatrix}
=
\begin{bmatrix} Z_{xx} & Z_{xy} \\ Z_{yx} & Z_{yy} \end{bmatrix}
\begin{bmatrix} H_x \\ H_y \end{bmatrix}
\]
Elemen tensor inilah yang kemudian diturunkan menjadi parameter interpretasi utama seperti apparent resistivity (ρa) dan phase (φ). Secara umum, data yang baik menunjukkan tren ρa dan φ yang mulus dan konsisten antar frekuensi, serta memiliki error bar yang realistis.
-
Fase 1: Kwaliteitscontrole (QC) en voorbewerking
Pengolahan MT yang baik dimulai dari QC sejak data lapangan. Tahap ini mencakup:
1. Pemeriksaan time series : melihat adanya lonjakan (spike), drift, saturasi sensor, atau gangguan periodik.
2. Kebisingan budaya : sumber noise umum adalah jaringan listrik 50/60 Hz, kereta, aktivitas industri, pagar listrik, dan telekomunikasi.
3. Kondisi elektroda dan kontak tanah : resistansi kontak tinggi akan memperburuk kualitas Ex/Ey.
4. Orientasi dan posisi sensor : kesalahan azimut sensor dapat menyebabkan distorsi interpretasi, terutama pada studi struktur berarah.
In deze fase worden doorgaans slechte datasegmenten verwijderd, wordt een offsetcorrectie uitgevoerd en wordt tijdsynchronisatie uitgevoerd indien gebruik wordt gemaakt van een extern referentiestation.
-
Fase 2: Frequentiedomeintransformatie en spectrale schatting
Karena MT menganalisis hubungan E dan H dalam domain frekuensi, time series diubah dengan teknik spektral seperti Fourier Transform . Data kemudian dibagi menjadi window (segmen) agar stabil secara statistik. Beberapa teknik yang sering dipakai:
– Windowing dan tapering (mis. Hanning) untuk mengurangi spectral leakage.
– Averaging antar window untuk memperoleh estimasi spektrum yang lebih robust.
– Cross-power spectra untuk membangun hubungan E terhadap H.
Het doel is om nauwkeurige impedantieschattingen te verkrijgen over een breed frequentiebereik, van hoge (lage) tot lage (diepe) frequenties.
-
Fase 3: Robuuste verwerking en referentie op afstand
Salah satu tantangan MT adalah noise yang seringkali berkorelasi pada kanal tertentu. Karena itu berkembang teknik robust processing , yang mengurangi pengaruh outlier dan segmen data buruk. Metode robust menggunakan pembobotan iteratif untuk menekan kontribusi data yang tidak konsisten.
Teknik lain yang sangat penting adalah remote reference (RR) . Konsepnya: medan magnet diukur juga pada stasiun yang jauh dari sumber noise lokal. Dengan mengkorelasikan E dan H di stasiun utama terhadap H di stasiun RR, efek noise lokal pada H dapat dikurangi. RR sangat efektif di area berisik, misalnya dekat permukiman atau infrastruktur.
De output van deze fase is doorgaans:
– ρa- en φ-curven voor Zxy en Zyx
– Foutbalk (standaarddeviatie)
– Coherentie tussen de primaire signalen
-
Fase 4: Ruisonderdrukking en vervormingscorrectie (statische verschuiving)
Selain noise, MT sering dipengaruhi static shift , yaitu pergeseran vertikal pada kurva apparent resistivity akibat heterogenitas dangkal (mis. lapisan lempung tipis, kerikil, atau kondisi elektroda). Static shift tidak banyak mengubah fase, tetapi menggeser ρa naik atau turun secara multiplikatif.
Statische shift-afhandeling kan op de volgende manieren worden gedaan:
1. Koreksi berbasis data TDEM/CSAMT sebagai kontrol resistivitas dangkal.
2. Inversi bersama (joint inversion) MT-TDEM agar model dangkal lebih terikat.
3. Pendekatan inversi yang robust terhadap shift , misalnya memperkenankan parameter shift per stasiun.
Selain itu, dilakukan pula deteksi outlier frekuensi : titik-titik tertentu pada kurva yang menyimpang tajam biasanya dibuang atau diberi error yang lebih besar.
-
Fase 5: Dimensionaliteit en stakingsanalyse
Voordat een inversie wordt uitgevoerd, is het belangrijk om te bepalen of de ondergrondse structuur 1D, 2D of 3D is. Dit beïnvloedt de keuze van de inversiemethode en de interpretatie ervan. Veelgebruikte analyses zijn onder andere:
– Parameter skew (mis. Bahr skew) untuk menilai tingkat 3D.
– Phase tensor untuk melihat arah dominan struktur tanpa terpengaruh static shift.
– Strike analysis untuk menentukan arah jurus (strike) dominan struktur 2D.
Als de data sterke 2D-kenmerken vertoont, wordt doorgaans een rotatie van de impedantietensor in de richting van de inslag uitgevoerd om de hoofdcomponenten (TE/TM) te verduidelijken. Als de 3D-kenmerken sterk aanwezig zijn, heeft een 3D-inversie de voorkeur.
-
Stap 6: MT-gegevensinversie (1D, 2D, 3D)
Interpretasi kuantitatif MT umumnya dilakukan melalui inversi , yaitu mencari model resistivitas yang menghasilkan respons data paling cocok. Inversi MT bersifat non-linear dan ill-posed, sehingga memerlukan regularisasi agar model tidak “liar”. Jenis inversi:
1. Inversi 1D : cocok untuk area berlapis (mis. sedimen mendatar). Cepat, tapi terbatas.
2. Inversi 2D : cocok untuk struktur memanjang seperti sesar, graben, atau sistem panas bumi berarah.
3. Inversi 3D : paling realistis untuk geologi kompleks, tetapi memerlukan data rapat, komputasi besar, dan QC ketat.
De inverse doelfunctie omvat doorgaans de afwijking van de data en de ruwheid van het model:
\[
\Phi = \Phi_d + \lambda \Phi_m
\]
waarbij λ de regularisatieparameter is. De keuze van λ is cruciaal: een te kleine waarde maakt het model te grof (overfitting), een te grote waarde maakt het model te glad en mist belangrijke geologische kenmerken.
-
Stap 7: Evaluatie van de inversieresultaten en de gevoeligheid
Na de inversie moet het model worden geëvalueerd, niet alleen visueel. Deze evaluatie omvat:
– Misfit (RMS) : apakah sesuai target (mis. RMS ~ 1–2 tergantung definisi error).
– Perbandingan kurva data vs respons model pada setiap stasiun.
– Uji resolusi : misalnya checkerboard test atau sensitivity map.
– Pengaruh prior model : mencoba beberapa model awal untuk melihat kestabilan hasil.
Jika model berubah drastis saat parameter kecil diubah, berarti interpretasi harus lebih hati-hati dan mungkin memerlukan data toevoeging.
-
Geologische interpretatie: Het verband tussen soortelijke weerstand en ondergrondse systemen
De soortelijke weerstand is geen directe indicator van het "gesteentetype", maar eerder een samengestelde fysieke respons. Er zijn echter wel een aantal algemene patronen die vaak als richtlijn kunnen dienen:
– Resistivitas rendah (konduktif) : lempung, alterasi hidrotermal (clay cap), zona fluida asin, grafit, atau mineral sulfida.
– Resistivitas tinggi (resistif) : batuan beku masif, batuan kering, zona silisifikasi, atau basement kristalin.
Bij geothermische exploratie laten klassieke modellen bijvoorbeeld vaak het volgende zien:
1. Clay cap konduktif di bagian atas
2. Reservoir lebih resistif di bawahnya
3. Zona upflow yang dikontrol struktur (patahan)
4. Sumber panas yang kadang tampak sebagai anomali resistif atau kompleks bergantung litologi dan temperatur
Interpretasi yang baik selalu dikaitkan dengan data lain: geologi permukaan, manifestasi termal, geokimia, gravitasi, seismik, atau data sumur.
-
Sluitend
De technieken voor de verwerking en interpretatie van MT-gegevens zijn met elkaar verbonden: van kwaliteitscontrole van tijdreeksen (QC), spectrale schatting, robuuste/remote referentieverwerking, statische verschuivingscorrectie en dimensionaliteitsanalyse tot 2D/3D-inversie en resolutie-evaluatie. Het succes van MT hangt niet alleen af van de inversiesoftware, maar ook van de kwaliteit van de metingen, het begrijpen van de ruis en de integratie ervan in de geologische context.
Met zijn gedisciplineerde workflow en op meerdere gegevens gebaseerde interpretatie is MT een krachtig instrument voor het in kaart brengen van weerstandsstructuren in de ondergrond en het ondersteunen van besluitvorming bij complexe exploratie- en geowetenschappelijke studies.
