Técnicas de procesamento e interpretación de datos de MT
Pendahuluan
O método magnetotelúrico (MT) é unha técnica xeofísica pasiva que utiliza as variacións naturais nos campos eléctricos e magnéticos da Terra para mapear a distribución da resistividade do subsolo. A resistividade está estreitamente relacionada co tipo de rocha, o contido de fluído, a temperatura e as estruturas xeolóxicas como as fallas ou as zonas de alteración. Polo tanto, o MT úsase amplamente para a exploración xeotérmica, a mineralización, as concas sedimentarias (hidrocarburos) e os estudos tectónicos e da codia terrestre.
Non obstante, os datos brutos de MT raramente se poden utilizar directamente. Requiren unha serie de pasos de procesamento para mellorar a relación sinal-ruído, seguidos dunha interpretación baseada nun modelo de resistividade xeoloxicamente consistente. Este artigo analiza as técnicas de procesamento e interpretación de datos de MT, desde a adquisición e o control de calidade ata o procesamento de impedancias, pasando pola inversión e a interpretación xeolóxica.
-
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} e Z_{xy} \\ Z_{yx} e 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: Control de calidade (CC) e preprocesamento
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.
Nesta fase, os segmentos de datos defectuosos adoitan cortarse, realízase unha corrección de desprazamento e realízase a sincronización horaria se se usa unha estación de referencia remota.
-
Fase 2: Transformación do dominio da frecuencia e estimación espectral
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.
O obxectivo é obter estimacións precisas da impedancia nun amplo rango de frecuencias, desde frecuencias altas (superficiais) ata frecuencias baixas (profundas).
-
Fase 3: Procesamento robusto e referencia remota
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.
A saída desta etapa adoita ser:
– Curvas ρa e φ para Zxy e Zyx
– Barra de erro (desviación estándar)
– Coherencia entre os sinais primarios
-
Fase 4: Eliminación de ruído e manexo da distorsión (desprazamento estático)
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.
A xestión de quendas estáticas pódese facer mediante:
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: Dimensionalidade e análise de golpes
Antes da inversión, é importante determinar se a estrutura do subsolo é unidimensional, bidimensional ou tridimensional. Isto inflúe na escolla do método de inversión e na súa interpretación. Algunhas análises comúns inclúen:
– 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.
Se os datos presentan características 2D fortes, normalmente realízase unha rotación tensorial de impedancia dirixida polo impacto para clarificar os compoñentes principais (TE/TM). Se as características 3D son fortes, prefírese unha inversión 3D.
-
Paso 6: Inversión de datos 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.
A función obxectivo inversa normalmente incorpora o desaxuste dos datos e a rugosidade do modelo:
\[
\Phi = \Phi_d + \lambda \Phi_m
\]
onde λ é o parámetro de regularización. A escolla de λ é crucial: demasiado pequeno fai que o modelo sexa demasiado groso (sobreaxuste), demasiado grande fai que o modelo sexa demasiado suave e pasa por alto características xeolóxicas importantes.
-
Paso 7: Avaliación dos resultados da inversión e da sensibilidade
Despois da inversión, o modelo debe ser avaliado, non só visualmente. Esta avaliación inclúe:
– 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 adición.
-
Interpretación xeolóxica: Relación da resistividade cos sistemas subterráneos
A resistividade non é un "tipo de rocha" directo, senón unha resposta física composta. Non obstante, algúns patróns xerais adoitan servir como guía:
– 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.
Na exploración xeotérmica, por exemplo, os modelos clásicos adoitan mostrar:
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.
-
Peche
As técnicas de procesamento e interpretación de datos de TA están interconectadas: desde o control de calidade (QC) de series temporais, a estimación espectral, o procesamento de referencias robusto/remoto, a corrección de desprazamento estático, a análise de dimensionalidade ata a inversión 2D/3D e a avaliación da resolución. O éxito da TA depende non só do software de inversión, senón tamén da calidade das medicións, da comprensión do ruído e da súa integración co contexto xeolóxico.
Co seu fluxo de traballo disciplinado e a súa interpretación baseada en múltiples datos, a MT é unha ferramenta poderosa para mapear estruturas de resistividade subterránea e axudar na toma de decisións en estudos complexos de exploración e xeociencia.
