MT数据处理和解释技术
彭达胡乱
大地电磁法（MT）是一种被动式地球物理勘探技术，它利用地球电场和磁场的自然变化来绘制地下电阻率分布图。电阻率与岩石类型、流体含量、温度以及断层或蚀变带等地质构造密切相关。因此，大地电磁法广泛应用于地热勘探、矿化作用、沉积盆地（油气）勘探以及构造和地壳研究。
然而，原始大地电磁数据很少能直接使用。它需要一系列处理步骤来提高信噪比，然后基于地质一致的电阻率模型进行解释。本文讨论了大地电磁数据的处理和解释技术，从数据采集和质量控制到阻抗处理，再到反演和地质解释。
-
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.
-
第一阶段：质量控制和预处理
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.
在这个阶段，通常会切断坏数据段，进行偏移校正，如果使用远程参考站，还会进行时间同步。
-
第二阶段：频域变换和频谱估计
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.
目标是在从高频（浅层）到低频（深层）的宽频率范围内获得准确的阻抗估计值。
-
第三阶段：稳健处理和远程参考
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.
此阶段的输出通常为：
– Zxy 和 Zyx 的 ρa 和 φ 曲线
– 误差线（标准差）
– 初级信号之间的相干性
-
第四阶段：降噪和失真处理（静态偏移）
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.
静态轮班处理可以通过以下方式完成：
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.
-
第五阶段：维度和打击分析
反演之前，确定地下结构是1D、2D还是3D至关重要。这会影响反演方法的选择和解释。常见的分析方法包括：
– 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.
如果数据呈现明显的二维特征，通常会进行沿走向的阻抗张量旋转，以明确主成分（TE/TM）。如果三维特征明显，则优选进行三维反演。
-
步骤 6：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.
逆目标函数通常考虑数据拟合误差和模型粗糙度：
\[
\Phi = \Phi_d + \lambda \Phi_m
\]
其中 λ 为正则化参数。λ 的选择至关重要：太小会导致模型过于粗糙（过拟合），太大则会导致模型过于平滑，从而遗漏重要的地质特征。
-
步骤 7：反演结果和灵敏度评估
反演之后，必须对模型进行评估，而不仅仅是目测。评估内容包括：
– 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 添加。
-
地质解释：电阻率与地下系统的关系
电阻率并非直接的“岩石类型”，而是一种综合的物理响应。然而，一些通用模式通常可作为指导原则：
– 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.
例如，在地热勘探中，经典模型通常显示：
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.
-
关闭
大地电磁（MT）数据处理和解释技术环环相扣：从时间序列质量控制（QC）、频谱估计、稳健/远程参考处理、静态偏移校正、维度分析，到二维/三维反演和分辨率评估。MT的成功不仅取决于反演软件，还取决于测量质量、噪声控制以及与地质背景的整合。
凭借其严谨的工作流程和基于多数据的解释，MT 是一种强大的工具，可用于绘制地下电阻率结构图，并辅助复杂的勘探和地球科学研究中的决策。
