Các kỹ thuật xử lý và diễn giải dữ liệu MT
Giới thiệu
Phương pháp đo điện từ (MT) là một kỹ thuật địa vật lý thụ động sử dụng các biến thiên tự nhiên trong điện trường và từ trường của Trái đất để lập bản đồ phân bố điện trở suất dưới lòng đất. Điện trở suất có liên quan chặt chẽ đến loại đá, hàm lượng chất lỏng, nhiệt độ và các cấu trúc địa chất như đứt gãy hoặc vùng biến đổi địa chất. Do đó, MT được sử dụng rộng rãi trong thăm dò địa nhiệt, khoáng hóa, các lưu vực trầm tích (hydrocarbon) và các nghiên cứu về kiến tạo và vỏ Trái đất.
Tuy nhiên, dữ liệu MT thô hiếm khi có thể sử dụng trực tiếp. Nó đòi hỏi một loạt các bước xử lý để cải thiện tỷ lệ tín hiệu trên nhiễu, tiếp theo là diễn giải dựa trên mô hình điện trở suất phù hợp về mặt địa chất. Bài viết này thảo luận về các kỹ thuật xử lý và diễn giải dữ liệu MT, từ thu thập và kiểm soát chất lượng đến xử lý trở kháng, thông qua phép nghịch đảo và diễn giải địa chất.
-
Cơ sở dữ liệu MT: Trường E, Trường H và Tensor trở kháng
Trong các khảo sát MT, các cảm biến ghi lại các thành phần điện trường (Ex, Ey) và các thành phần từ trường (Hx, Hy, đôi khi là Hz) theo thời gian. Từ mối quan hệ giữa điện trường và từ trường trong miền tần số, người ta thu được tenxơ trở kháng:
\[
\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}
\]
Các phần tử tenxơ này sau đó được chuyển đổi thành các tham số diễn giải quan trọng như điện trở suất biểu kiến (ρa) và pha (φ). Nhìn chung, dữ liệu tốt thể hiện xu hướng mượt mà và nhất quán của ρa và φ trên các tần số, và có sai số thực tế.
-
Giai đoạn 1: Kiểm soát chất lượng (QC) và tiền xử lý
Quá trình xử lý MT tốt bắt đầu từ khâu kiểm soát chất lượng, bắt đầu từ dữ liệu thực địa. Giai đoạn này bao gồm:
1. Kiểm tra chuỗi thời gian: tìm kiếm các đỉnh nhọn, sự trôi lệch, hiện tượng bão hòa cảm biến hoặc các nhiễu loạn định kỳ.
2. Tiếng ồn văn hóa: các nguồn tiếng ồn phổ biến là đường dây điện 50/60 Hz, tàu hỏa, hoạt động công nghiệp, hàng rào điện và viễn thông.
3. Tình trạng của điện cực và các điểm tiếp xúc nối đất: điện trở tiếp xúc cao sẽ làm giảm chất lượng của Ex/Ey.
4. Hướng và vị trí cảm biến: sai số góc phương vị của cảm biến có thể gây ra sự sai lệch trong việc diễn giải kết quả, đặc biệt là trong các nghiên cứu cấu trúc định hướng.
Ở giai đoạn này, các đoạn dữ liệu lỗi thường bị cắt bỏ, việc hiệu chỉnh độ lệch được thực hiện và đồng bộ hóa thời gian được tiến hành nếu sử dụng trạm tham chiếu từ xa.
-
Giai đoạn 2: Biến đổi miền tần số và ước lượng phổ
Vì MT phân tích mối quan hệ giữa E và H trong miền tần số, chuỗi thời gian được biến đổi bằng các kỹ thuật phổ như biến đổi Fourier. Dữ liệu sau đó được chia thành các cửa sổ (phân đoạn) để đảm bảo tính ổn định thống kê. Một số kỹ thuật thường được sử dụng bao gồm:
– Sử dụng cửa sổ và bộ lọc làm mờ (ví dụ: Hanning) để giảm hiện tượng rò rỉ phổ.
– Lấy trung bình giữa các cửa sổ để thu được ước lượng phổ chính xác hơn.
– Phân tích phổ công suất chéo để thiết lập mối quan hệ giữa E và H.
Mục tiêu là thu được các ước tính trở kháng chính xác trên một dải tần số rộng, từ tần số cao (nông) đến tần số thấp (sâu).
-
Giai đoạn 3: Xử lý mạnh mẽ và tham chiếu từ xa
Một trong những thách thức của MT là nhiễu thường tương quan trong một kênh cụ thể. Do đó, các kỹ thuật xử lý mạnh mẽ đã được phát triển để giảm ảnh hưởng của các giá trị ngoại lai và các đoạn dữ liệu xấu. Các phương pháp mạnh mẽ sử dụng trọng số lặp để triệt tiêu sự đóng góp của dữ liệu không nhất quán.
Một kỹ thuật rất quan trọng khác là tham chiếu từ xa (RR). Khái niệm là từ trường được đo tại một trạm cách xa nguồn nhiễu cục bộ. Bằng cách tương quan E và H tại trạm chính với H tại trạm RR, ảnh hưởng của nhiễu cục bộ lên H có thể được giảm thiểu. RR đặc biệt hiệu quả ở những khu vực nhiều nhiễu, ví dụ như gần khu dân cư hoặc cơ sở hạ tầng.
Thông thường, kết quả thu được từ giai đoạn này là:
– Đường cong ρa và φ cho Zxy và Zyx
– Sai số chuẩn (độ lệch chuẩn)
– Sự nhất quán giữa các tín hiệu chính
-
Giai đoạn 4: Khử nhiễu và xử lý méo tiếng (Dịch chuyển tĩnh)
Ngoài nhiễu, MT thường bị ảnh hưởng bởi sự dịch chuyển tĩnh, là sự dịch chuyển theo chiều dọc của đường cong điện trở suất biểu kiến do các dị chất nông (ví dụ: lớp đất sét mỏng, sỏi hoặc điều kiện điện cực). Sự dịch chuyển tĩnh không làm thay đổi pha một cách đáng kể, nhưng làm dịch chuyển ρa lên hoặc xuống theo cấp số nhân.
Việc xử lý sự thay đổi tĩnh có thể được thực hiện thông qua:
1. Hiệu chỉnh dựa trên dữ liệu TDEM/CSAMT để kiểm soát điện trở suất nông.
2. Kết hợp đảo ngược MT-TDEM để làm cho mô hình nông hơn và có độ chính xác cao hơn.
3. Phương pháp đảo ngược mạnh mẽ để dịch chuyển, ví dụ như cho phép các tham số dịch chuyển cho từng trạm.
Ngoài ra, việc phát hiện các giá trị ngoại lệ về tần số cũng được thực hiện: một số điểm trên đường cong có độ lệch lớn thường bị loại bỏ hoặc được gán giá trị sai số lớn hơn.
-
Giai đoạn 5: Phân tích chiều không gian và cú đánh
Trước khi thực hiện phép nghịch đảo, điều quan trọng là phải xác định xem cấu trúc dưới lòng đất là 1D, 2D hay 3D. Điều này ảnh hưởng đến việc lựa chọn phương pháp nghịch đảo và cách diễn giải kết quả. Các phân tích phổ biến bao gồm:
– Các tham số độ lệch (ví dụ: độ lệch Bahr) để đánh giá mức độ 3D.
– Sử dụng tenxơ pha để xác định hướng chủ đạo của cấu trúc mà không bị ảnh hưởng bởi sự dịch chuyển tĩnh.
– Phân tích phương vị để xác định hướng phương vị chủ đạo của cấu trúc 2D.
Nếu dữ liệu thể hiện rõ các đặc điểm 2D, phương pháp xoay tensor trở kháng theo hướng phương vị thường được thực hiện để làm rõ các thành phần chính (TE/TM). Nếu các đặc điểm 3D mạnh, phương pháp đảo ngược 3D được ưu tiên hơn.
-
Bước 6: Đảo ngược dữ liệu MT (1D, 2D, 3D)
Việc diễn giải định lượng MT thường được thực hiện thông qua phép nghịch đảo, bao gồm việc tìm ra mô hình điện trở suất phù hợp nhất với dữ liệu phản hồi. Phép nghịch đảo MT là phi tuyến tính và không ổn định, đòi hỏi phải có sự điều chỉnh để ngăn mô hình trở nên "không kiểm soát". Các loại phép nghịch đảo:
1. Phương pháp đảo ngược 1D: thích hợp cho các khu vực có cấu trúc phân lớp (ví dụ: trầm tích nằm ngang). Nhanh nhưng có hạn chế.
2. Phương pháp nghịch đảo 2D: thích hợp cho các cấu trúc kéo dài như đứt gãy, rãnh sụt hoặc hệ thống địa nhiệt định hướng.
3. Phương pháp đảo ngược 3D: phương pháp thực tế nhất cho địa chất phức tạp, nhưng đòi hỏi dữ liệu dày đặc, khối lượng tính toán lớn và kiểm soát chất lượng nghiêm ngặt.
Hàm mục tiêu nghịch đảo thường kết hợp sai lệch dữ liệu và độ không chính xác của mô hình:
\[
\Phi = \Phi_d + \lambda \Phi_m
\]
trong đó λ là tham số điều chỉnh. Việc lựa chọn λ rất quan trọng: quá nhỏ sẽ làm cho mô hình quá thô (quá khớp), quá lớn sẽ làm cho mô hình quá trơn tru và bỏ sót các đặc điểm địa chất quan trọng.
-
Bước 7: Đánh giá kết quả đảo ngược và độ nhạy
Sau khi đảo ngược mô hình, cần phải đánh giá nó, không chỉ bằng mắt thường. Việc đánh giá này bao gồm:
– Sai số trung bình bình phương (RMS): mức độ phù hợp với mục tiêu (ví dụ: RMS ~ 1–2 tùy thuộc vào định nghĩa về sai số).
– So sánh đường cong dữ liệu với phản hồi của mô hình tại mỗi trạm.
– Kiểm tra độ phân giải: ví dụ như kiểm tra bàn cờ hoặc bản đồ độ nhạy.
– Ảnh hưởng của các mô hình trước đó: hãy thử một vài mô hình ban đầu để xem tính ổn định của kết quả.
Nếu mô hình thay đổi đáng kể khi một tham số nhỏ được thay đổi, điều đó có nghĩa là việc diễn giải cần thận trọng hơn và có thể cần thêm dữ liệu.
-
Giải thích địa chất: Mối liên hệ giữa điện trở suất và các hệ thống dưới lòng đất
Điện trở suất không phải là một "loại đá" trực tiếp, mà là một phản ứng vật lý tổng hợp. Tuy nhiên, một số quy luật chung thường được dùng làm hướng dẫn:
– Điện trở suất thấp (dẫn điện): đất sét, biến đổi thủy nhiệt (lớp đất sét phủ), vùng dung dịch muối, than chì hoặc khoáng chất sunfua.
– Điện trở suất cao (điện trở): đá magma khối lớn, đá khô, vùng silic hóa hoặc nền móng tinh thể.
Ví dụ, trong thăm dò địa nhiệt, các mô hình cổ điển thường cho thấy:
1. Nắp đất sét dẫn điện ở phía trên
2. Phần dưới của bể chứa có điện trở lớn hơn.
3. Vùng dòng chảy ngược do cấu trúc kiểm soát (đứt gãy)
4. Các nguồn nhiệt đôi khi xuất hiện dưới dạng các dị thường điện trở hoặc phức tạp tùy thuộc vào thành phần đá và nhiệt độ.
Việc giải thích tốt luôn gắn liền với các dữ liệu khác: địa chất bề mặt, biểu hiện nhiệt, địa hóa học, trọng lực, địa chấn hoặc dữ liệu giếng khoan.
-
Đóng cửa
Các kỹ thuật xử lý và diễn giải dữ liệu MT có mối liên hệ mật thiết với nhau: từ kiểm soát chất lượng chuỗi thời gian (QC), ước lượng phổ, xử lý tham chiếu từ xa/mạnh mẽ, hiệu chỉnh dịch chuyển tĩnh, phân tích chiều, đến đảo ngược 2D/3D và đánh giá độ phân giải. Sự thành công của MT không chỉ phụ thuộc vào phần mềm đảo ngược mà còn phụ thuộc vào chất lượng của các phép đo, hiểu biết về nhiễu và tích hợp nó với bối cảnh địa chất.
Với quy trình làm việc bài bản và khả năng diễn giải dựa trên nhiều nguồn dữ liệu, MT là một công cụ mạnh mẽ để lập bản đồ cấu trúc điện trở suất dưới lòng đất và hỗ trợ ra quyết định trong các nghiên cứu thăm dò và địa chất phức tạp.
