Teknik Pemprosesan dan Interpretasi Data MT
pengenalan
Kaedah Magnetotelurik (MT) ialah teknik geofizik pasif yang menggunakan variasi semula jadi dalam medan elektrik dan magnet Bumi untuk memetakan taburan kerintangan bawah permukaan. Kerintangan berkait rapat dengan jenis batuan, kandungan bendalir, suhu dan struktur geologi seperti sesar atau zon perubahan. Oleh itu, MT digunakan secara meluas untuk penerokaan geoterma, mineralisasi, lembangan sedimen (hidrokarbon) dan kajian tektonik dan kerak bumi.
Walau bagaimanapun, data MT mentah jarang boleh digunakan secara langsung. Ia memerlukan satu siri langkah pemprosesan untuk meningkatkan nisbah isyarat-ke-hingar, diikuti dengan interpretasi berdasarkan model kerintangan yang konsisten secara geologi. Artikel ini membincangkan teknik pemprosesan dan interpretasi data MT, daripada pemerolehan dan kawalan kualiti kepada pemprosesan impedans, melalui penyongsangan dan interpretasi geologi.
-
Pangkalan Data MT: Medan E, Medan H dan Tensor Impedans
Dalam tinjauan MT, sensor merekodkan komponen medan elektrik (Ex, Ey) dan komponen medan magnet (Hx, Hy, kadangkala Hz) sebagai fungsi masa. Daripada hubungan antara medan elektrik dan magnet dalam domain frekuensi, tensor impedans diperoleh:
\[
\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}
\]
Unsur-unsur tensor ini kemudiannya diterbitkan menjadi parameter interpretasi utama seperti kerintangan ketara (ρa) dan fasa (φ). Secara amnya, data yang baik mempamerkan trend yang lancar dan konsisten dalam ρa dan φ merentasi frekuensi, dan mempunyai bar ralat yang realistik.
-
Peringkat 1: Kawalan Kualiti (QC) dan Pra-pemprosesan
Pemprosesan MT yang baik bermula dengan QC, bermula dengan data lapangan. Peringkat ini merangkumi:
1. Pemeriksaan siri masa: mencari lonjakan, hanyutan, ketepuan sensor atau gangguan berkala.
2. Bunyi budaya: sumber bunyi yang biasa ialah talian kuasa 50/60 Hz, kereta api, aktiviti perindustrian, pagar elektrik dan telekomunikasi.
3. Keadaan elektrod dan sentuhan pembumian: rintangan sentuhan yang tinggi akan memburukkan kualiti Ex/Ey.
4. Orientasi dan kedudukan sensor: ralat azimut sensor boleh menyebabkan herotan tafsiran, terutamanya dalam kajian struktur arah.
Pada peringkat ini, segmen data yang rosak biasanya dipotong, pembetulan ofset dilakukan, dan penyegerakan masa dilakukan jika menggunakan stesen rujukan jauh.
-
Peringkat 2: Transformasi Domain Frekuensi dan Penganggaran Spektrum
Oleh kerana MT menganalisis hubungan antara E dan H dalam domain frekuensi, siri masa diubah menggunakan teknik spektrum seperti Transformasi Fourier. Data kemudiannya dibahagikan kepada tetingkap (segmen) untuk memastikan kestabilan statistik. Beberapa teknik yang biasa digunakan termasuk:
– Pemasangan tingkap dan tirus (cth. Hanning) untuk mengurangkan kebocoran spektrum.
– Membuat purata antara tetingkap untuk mendapatkan anggaran spektrum yang lebih kukuh.
– Spektrum kuasa silang untuk menentukan hubungan E dengan H.
Matlamatnya adalah untuk mendapatkan anggaran impedans yang tepat pada julat frekuensi yang luas, daripada frekuensi tinggi (cetek) hingga frekuensi rendah (dalam).
-
Peringkat 3: Pemprosesan yang Mantap dan Rujukan Jauh
Salah satu cabaran MT ialah hingar yang sering berkorelasi dalam saluran tertentu. Oleh itu, teknik pemprosesan yang mantap telah dibangunkan untuk mengurangkan pengaruh outlier dan segmen data yang buruk. Kaedah yang mantap menggunakan pemberat lelaran untuk menyekat sumbangan data yang tidak konsisten.
Satu lagi teknik yang sangat penting ialah rujukan jauh (RR). Konsepnya ialah medan magnet diukur di stesen yang jauh dari sumber hingar tempatan. Dengan menghubungkan E dan H di stesen utama dengan H di stesen RR, kesan hingar tempatan terhadap H dapat dikurangkan. RR amat berkesan di kawasan yang bising, contohnya berhampiran kawasan perumahan atau infrastruktur.
Output dari peringkat ini biasanya:
– lengkung ρa dan φ untuk Zxy dan Zyx
– Bar ralat (sisihan piawai)
– Kekoherenan antara isyarat utama
-
Peringkat 4: Penyah-bising dan Pengendalian Herotan (Anjakan Statik)
Selain hingar, MT sering dipengaruhi oleh anjakan statik, iaitu anjakan menegak dalam lengkung kerintangan ketara disebabkan oleh heterogeniti cetek (contohnya, lapisan tanah liat nipis, kerikil, atau keadaan elektrod). Anjakan statik tidak mengubah fasa dengan ketara, tetapi mengalihkan ρa ke atas atau ke bawah secara berganda.
Pengendalian anjakan statik boleh dilakukan melalui:
1. Pembetulan berasaskan data TDEM/CSAMT sebagai kawalan kerintangan cetek.
2. Songsangan bersama MT-TDEM untuk menjadikan model cetek lebih terikat.
3. Pendekatan songsangan yang mantap untuk anjakan, contohnya membenarkan parameter anjakan setiap stesen.
Di samping itu, pengesanan penyimpangan frekuensi juga dijalankan: titik-titik tertentu pada lengkung yang menyimpang secara mendadak biasanya dibuang atau diberikan ralat yang lebih besar.
-
Peringkat 5: Dimensi dan Analisis Hentaman
Sebelum penyongsangan, adalah penting untuk menentukan sama ada struktur bawah permukaan adalah 1D, 2D atau 3D. Ini mempengaruhi pilihan kaedah dan tafsiran penyongsangan. Analisis biasa termasuk:
– Parameter condong (cth. Bahr condong) untuk menilai tahap 3D.
– Tensor fasa untuk melihat arah dominan struktur tanpa terjejas oleh anjakan statik.
– Analisis hentakan untuk menentukan arah hentakan dominan bagi struktur 2D.
Jika data menunjukkan ciri-ciri 2D yang kuat, putaran tensor impedans arah strike biasanya dilakukan untuk menjelaskan komponen utama (TE/TM). Jika ciri-ciri 3D adalah kuat, penyongsangan 3D adalah lebih baik.
-
Langkah 6: Penyongsangan Data MT (1D, 2D, 3D)
Tafsiran kuantitatif MT biasanya dilakukan melalui penyongsangan, yang melibatkan pencarian model kerintangan yang paling sesuai dengan tindak balas data. Penyongsangan MT adalah tidak linear dan tidak tepat, memerlukan regularisasi untuk mengelakkan model daripada menjadi "liar". Jenis penyongsangan:
1. Songsangan 1D: sesuai untuk kawasan berlapis (cth., sedimen mendatar). Cepat, tetapi terhad.
2. Songsangan 2D: sesuai untuk struktur memanjang seperti sesar, graben atau sistem geoterma berarah.
3. Songsangan 3D: paling realistik untuk geologi kompleks, tetapi memerlukan data yang padat, pengiraan yang besar dan QC yang ketat.
Fungsi objektif songsang biasanya menggabungkan ketidakpadanan data dan kekasaran model:
\[
\Phi = \Phi_d + \lambda \Phi_m
\]
di mana λ ialah parameter regularisasi. Pemilihan λ adalah penting: terlalu kecil menjadikan model terlalu kasar (overfitting), terlalu besar menjadikan model terlalu licin dan terlepas ciri-ciri geologi yang penting.
-
Langkah 7: Penilaian Keputusan dan Kepekaan Penyongsangan
Selepas penyongsangan, model mesti dinilai, bukan hanya secara visual. Penilaian ini merangkumi:
– Tidak Padan (RMS): sama ada ia sepadan dengan sasaran (cth. RMS ~ 1–2 bergantung pada definisi ralat).
– Perbandingan lengkung data vs tindak balas model di setiap stesen.
– Ujian resolusi: contohnya ujian papan dam atau peta kepekaan.
– Pengaruh model terdahulu: cuba beberapa model awal untuk melihat kestabilan keputusan.
Jika model berubah secara drastik apabila parameter kecil diubah, ini bermakna interpretasi harus lebih berhati-hati dan mungkin memerlukan data tambahan.
-
Tafsiran Geologi: Menghubungkan Kerintangan dengan Sistem Bawah Permukaan
Kerintangan bukanlah "jenis batuan" secara langsung, tetapi sebaliknya tindak balas fizikal komposit. Walau bagaimanapun, beberapa corak umum sering menjadi garis panduan:
– Kerintangan rendah (konduktif): tanah liat, pengubahan hidroterma (tudung tanah liat), zon bendalir garam, grafit atau mineral sulfida.
– Kerintangan tinggi (rintangan): batuan igneus masif, batuan kering, zon silisifikasi atau dasar kristal.
Dalam penerokaan geoterma, sebagai contoh, model klasik sering menunjukkan:
1. Penutup tanah liat konduktif di atas
2. Takungan lebih berrintangan di bawahnya
3. Zon aliran naik terkawal struktur (kerosakan)
4. Sumber haba yang kadangkala kelihatan sebagai anomali resistif atau kompleks bergantung pada litologi dan suhu.
Tafsiran yang baik sentiasa dikaitkan dengan data lain: geologi permukaan, manifestasi terma, geokimia, graviti, seismik atau data telaga.
-
penutup
Teknik pemprosesan dan interpretasi data MT saling berkaitan: daripada kawalan kualiti siri masa (QC), anggaran spektrum, pemprosesan rujukan teguh/jauh, pembetulan anjakan statik, analisis dimensi, kepada penilaian penyongsangan dan resolusi 2D/3D. Kejayaan MT bukan sahaja bergantung pada perisian penyongsangan, tetapi juga pada kualiti pengukuran, pemahaman hingar dan penyepaduannya dengan konteks geologi.
Dengan aliran kerja yang berdisiplin dan interpretasi berasaskan berbilang data, MT merupakan alat yang ampuh untuk memetakan struktur kerintangan bawah permukaan dan membantu membuat keputusan dalam penerokaan kompleks dan kajian geosains.
