Teknik Pengolahan dan Interpretasi Data MT
Pendahuluan
Metode Magnetotellurik (MT) adalah salah satu teknik geofisika pasif yang memanfaatkan variasi alami medan listrik dan medan magnet bumi untuk memetakan distribusi resistivitas bawah permukaan. Resistivitas berkaitan erat dengan jenis batuan, kandungan fluida, temperatur, dan struktur geologi seperti patahan atau zona alterasi. Karena itu, MT banyak dipakai untuk eksplorasi panas bumi, mineralisasi, cekungan sedimen (hidrokarbon), hingga studi tektonik dan kerak bumi.
Namun, data MT mentah jarang bisa langsung dipakai. Ia memerlukan rangkaian pengolahan (processing) untuk meningkatkan rasio sinyal terhadap derau (noise), lalu interpretasi berbasis model resistivitas yang konsisten dengan geologi. Artikel ini membahas teknik pengolahan dan interpretasi data MT, mulai dari tahap akuisisi, quality control, pemrosesan impedansi, hingga inversi dan penafsiran geologinya.
—
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.
—
Tahap 1: Quality Control (QC) dan Pra-pemrosesan
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.
Pada tahap ini biasanya dilakukan pemotongan segmen data yang buruk, koreksi offset, serta sinkronisasi waktu bila menggunakan stasiun remote reference.
—
Tahap 2: Transformasi Domain Frekuensi dan Estimasi Spektral
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.
Tujuannya adalah mendapatkan estimasi impedansi yang akurat pada rentang frekuensi luas, dari frekuensi tinggi (dangkal) hingga rendah (dalam).
—
Tahap 3: Robust Processing dan Remote Reference
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.
Output dari tahap ini biasanya berupa:
– Kurva ρa dan φ untuk Zxy dan Zyx
– Error bar (deviasi standar)
– Koherensi (coherence) antara sinyal-sinyal utama
—
Tahap 4: De-noising dan Penanganan Distorsi (Static Shift)
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.
Penanganan static shift dapat dilakukan melalui:
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.
—
Tahap 5: Analisis Dimensionalitas dan Strike
Sebelum inversi, penting menentukan apakah struktur bawah permukaan bersifat 1D, 2D, atau 3D. Hal ini memengaruhi pilihan metode inversi dan cara interpretasi. Analisis yang umum:
– 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.
Jika data menunjukkan karakter 2D kuat, biasanya dilakukan rotasi tensor impedansi ke arah strike agar komponen utama (TE/TM) lebih jelas. Jika karakter 3D kuat, inversi 3D lebih disarankan.
—
Tahap 6: Inversi Data MT (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.
Fungsi objektif inversi biasanya memadukan misfit data dan kekasaran model:
\[
\Phi = \Phi_d + \lambda \Phi_m
\]
di mana λ adalah parameter regularisasi. Pemilihan λ krusial: terlalu kecil membuat model terlalu kasar (overfitting), terlalu besar membuat model terlalu halus dan kehilangan fitur geologi penting.
—
Tahap 7: Evaluasi Hasil Inversi dan Sensitivitas
Setelah inversi, model harus dievaluasi, bukan hanya dilihat secara visual. Evaluasinya meliputi:
– 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 tambahan.
—
Interpretasi Geologi: Menghubungkan Resistivitas dengan Sistem Bawah Permukaan
Resistivitas bukan “jenis batuan” secara langsung, melainkan respon fisik gabungan. Namun, beberapa pola umum sering dijadikan pedoman:
– 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.
Dalam eksplorasi panas bumi, misalnya, model klasik sering menunjukkan:
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.
—
Penutup
Teknik pengolahan dan interpretasi data MT merupakan rangkaian yang saling terkait: mulai dari QC time series, estimasi spektral, robust/remote reference processing, koreksi static shift, analisis dimensionalitas, hingga inversi 2D/3D dan evaluasi resolusi. Keberhasilan MT tidak hanya bergantung pada software inversi, tetapi pada kualitas pengukuran, pemahaman noise, dan integrasi dengan konteks geologi.
Dengan workflow yang disiplin dan interpretasi yang berbasis multi-data, MT menjadi alat yang sangat kuat untuk memetakan struktur resistivitas bawah permukaan serta membantu pengambilan keputusan pada eksplorasi dan studi geosains yang kompleks.
