Tổng quan nghiên cứu

Trong công tác điều tra địa chất và thăm dò khoáng sản hiện đại, phương pháp phổ Gamma hàng không đóng vai trò trụ cột khi đảm nhận tới hơn 70% khối lượng khảo sát diện rộng ở tỷ lệ 1:25.000 đến 1:50.000. Tuy nhiên, việc xử lý nguồn dữ liệu bay đo khổng lồ thường gặp rào cản lớn do các phương pháp truyền thống chủ yếu phân tích định tính trên từng dị thường đơn lẻ, dẫn đến nguy cơ bỏ sót hoặc đánh giá thiếu chính xác các đới khoáng hóa ẩn sâu. Nhằm giải quyết triệt để vấn đề này, đề tài tập trung nghiên cứu, hoàn thiện và thử nghiệm hai thuật toán tiên tiến gồm phương pháp Hệ số tương quan và phương pháp Tần suất - Nhận dạng trên cơ sở thuật toán Griffiths-Vinni.

Mục tiêu cụ thể của công trình là xây dựng quy trình tự động hóa giải đoán các trường xạ địa hóa, phân loại chính xác các cụm dị thường và nâng cao độ tin cậy trong công tác khoanh vùng khoáng sản có ích. Phạm vi nghiên cứu được triển khai trên tập số liệu bay đo địa vật lý thực tế tại khu vực Nam Trung Bộ, bao phủ diện tích khoảng 300 km2 thuộc địa bàn 6 xã thuộc hai huyện Khánh Vĩnh và Diên Khánh, tỉnh Khánh Hòa (tọa độ 12°11'48" đến 12°20'02" vĩ độ Bắc và 108°53' đến 108°59' kinh độ Đông). Kết quả nghiên cứu có ý nghĩa khoa học và thực tiễn sâu sắc khi giúp tối ưu hóa 37,5% số lượng tham số đầu vào, nâng cao độ chính xác nhận dạng đối tượng quặng hóa vàng lên trên 80%, đồng thời tiết kiệm ước tính khoảng 25% đến 30% chi phí khảo sát kiểm tra mặt đất so với các quy trình phân tích thủ công trước đây.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên hai nền tảng lý thuyết cốt lõi của địa vật lý hạt nhân và địa hóa học phóng xạ:

Thứ nhất là lý thuyết phân rã phóng xạ tự nhiên của các đồng vị không bền vững, phát xạ lượng tử Gamma mang năng lượng đặc trưng cho từng hạt nhân nguyên tử. Cường độ bức xạ được đo đạc tập trung tại 3 cửa sổ năng lượng tiêu chuẩn: cửa sổ Kali (K) trong dải 1,35–1,55 MeV, cửa sổ Uranium (U) trong dải 1,55–1,85 MeV và cửa sổ Thorium (Th) trong dải 2,40–2,80 MeV. Các hệ số chuyển đổi chuẩn quy định 1% Kali trong đất đá tương ứng 313 Bq/kg của đồng vị 40K, 1 ppm Uranium tương ứng 12,35 Bq/kg của 238U và 1 ppm Thorium tương ứng 4,06 Bq/kg của 232Th.

Thứ hai là lý thuyết trường xạ địa hóa vỏ Trái Đất theo mô hình của A. Vinogradov. Trong các thành tạo đá magma nguyên sinh, hàm lượng các nguyên tố phóng xạ có mối tương quan chặt chẽ với hàm lượng SiO2, thể hiện qua tỷ số Th/U ổn định trong khoảng 3,7 đến 5,1 (đá axit đạt tỷ số Th/U khoảng 5,1; đá kiềm và bazơ dao động quanh 3,7 đến 3,9). Mô hình nghiên cứu vận dụng các khái niệm then chốt gồm: dị thường xạ địa hóa cục bộ, 8 tham số đặc trưng cụm dị thường, hàm mật độ phân bố xác suất và chỉ số đồng dạng P*m.

Phương pháp nghiên cứu

Nguồn dữ liệu của đề tài trích xuất từ đề án bay đo từ - phổ Gamma tỷ lệ 1:50.000 và trọng lực tỷ lệ 1:100.000 vùng Phan Rang – Nha Trang do Liên đoàn Vật lý Địa chất thực hiện, kết hợp cùng các công trình nghiên cứu chuyên sâu giai đoạn 2017–2020.

Về cỡ mẫu và phương pháp chọn mẫu, tác giả áp dụng phương pháp chọn mẫu điển hình có mục đích với cỡ mẫu gồm 8 cụm đối tượng địa chất - xạ địa hóa trên tổng diện tích 300 km2. Trong đó, 6 đối tượng là các diện tích có biểu hiện quặng vàng đã được kiểm chứng thực địa để làm mẫu đối sánh chuẩn và 2 đối tượng đóng vai trò mẫu kiểm định khách quan.

Lý do lựa chọn phương pháp phân tích dựa trên tính ưu việt của hai công cụ:

  • Phương pháp Hệ số tương quan đa biến (RU/Th, RU/K, RTh/K) trên cửa sổ quét di động: Cho phép phát hiện sự phá vỡ cân bằng xạ địa hóa tự nhiên tại các đới biến đổi nhiệt dịch kiến tạo, nơi hệ số tương quan R suy giảm mạnh so với trạng thái bình thường.
  • Phương pháp Tần suất - Nhận dạng Griffiths-Vinni kết hợp phần mềm QTS và chương trình xử lý tự động viết trên môi trường MATLAB: Xử lý hiệu quả ma trận thông tin nhị phân (0 và 1), tự động sàng lọc tổ hợp dấu hiệu có tỷ trọng thông tin tích lũy Pm vượt ngưỡng 80%.

Toàn bộ quy trình thu thập, lập trình thuật toán và thử nghiệm được hoàn thành trong lộ trình 12 tháng.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình phân tích thực nghiệm trên tài liệu phổ Gamma hàng không vùng nghiên cứu đã mang lại 4 phát hiện quan trọng:

Thứ nhất, phân tích hệ số tương quan di động trên 3 cặp thông số RU/Th, RU/K và RTh/K cho thấy sự phân hóa địa chất rất rõ nét. Tại các vùng nền địa chất ổn định, hệ số tương quan R đạt mức cao từ 0,72 đến 0,88. Ngược lại, tại các đới đứt gãy và biến đổi nhiệt dịch, giá trị R giảm mạnh xuống dưới 0,35 (giảm hơn 50% độ liên kết), tạo thành các dị thường ranh giới rõ rệt chỉ thị quặng hóa.

Thứ hai, thuật toán Griffiths-Vinni đã tối ưu hóa thành công ma trận dữ liệu từ 8 tham số ban đầu xuống còn 5 tham số đặc trưng nhất (gồm gia số cường độ ΔJ, độ rộng nửa biên độ T(1/2), tỷ số ΔTh/ΔU, tỷ số ΔU/ΔK và chỉ số đa thành phần F). Tổ hợp 5 tham số này đạt tỷ trọng thông tin tích lũy Pm = 85,4%, vượt 15,4% so với ngưỡng yêu cầu thông thường là 70%, giúp tinh giản 37,5% khối lượng tính toán phức tạp.

Thứ ba, kết quả xác định chỉ số đồng dạng Pm trên 7 đối tượng thử nghiệm cho thấy 5/6 đối tượng có triển vọng quặng vàng đạt giá trị Pm rất cao từ 0,78 đến 0,92, khẳng định mức độ tương đồng xạ địa hóa trên 80% so với mỏ chuẩn đã biết.

Thứ tư, phương pháp đã khoanh định chính xác 3 cụm dị thường xạ địa hóa mới có quy mô diện tích từ 1,8 km2 đến 4,5 km2, trùng khớp trên 90% với các cấu trúc kiến tạo sâu đã được định vị bằng tài liệu từ và trọng lực.

Thảo luận kết quả

Sự sụt giảm đột ngột của hệ số tương quan tuyến tính tại các cụm dị thường được giải thích bởi cơ chế nhiệt dịch hậu magma. Quá trình biến đổi magma đã dẫn đến sự tái phân bố cục bộ, làm giàu đột biến nguyên tố Kali hoặc Uranium và phá vỡ tỷ số nguyên thủy Th/U = 3,9–5,1 vốn có trong đá vây quanh.

So sánh với các phương pháp tách trường truyền thống trong bộ chương trình COSCAD hay các bộ lọc tuyến tính thông thường, phương pháp Tần suất - Nhận dạng cải tiến có ưu thế vượt trội khi tự động mở rộng không gian dự báo cho cả những khu vực chưa biết trước đối tượng đối sánh. Độ nhạy và khả năng phân lớp đối tượng tăng khoảng 22% so với cách phân loại định tính thủ công.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua đồ thị phân tán quan hệ đa biến X-Y, sơ đồ khối quy trình xử lý phần mềm QTS, bản đồ đẳng trị hệ số tương quan và bản đồ đồng mức chỉ số đồng dạng P*m(x,y). Việc trực quan hóa này giúp các nhà địa chất nhanh chóng phân định ranh giới tiếp xúc địa chất và phát hiện các thân quặng ẩn sâu dưới tầng phủ phong hóa dày.

Đề xuất và khuyến nghị

Dựa trên kết quả nghiên cứu lý thuyết và kiểm nghiệm thực tế, tác giả đưa ra 4 khuyến nghị then chốt:

  1. Nâng cấp và chuẩn hóa công cụ phần mềm: Hoàn thiện mã nguồn MATLAB và module nhận dạng QTS thành gói phần mềm thương mại tích hợp trực tiếp vào hệ thống GIS, hướng tới nâng cao tốc độ xử lý dữ liệu mạng lưới bay đo tỷ lệ 1:25.000 thêm 40% trong thời gian 6 tháng tới. Đơn vị thực hiện: Nhóm nghiên cứu thuộc Trường Đại học Khoa học Tự nhiên phối hợp cùng các viện chuyên ngành.
  2. Ban hành quy chuẩn phân loại cụm dị thường: Cục Địa chất và Khoáng sản Việt Nam cần sớm xây dựng thông tư hướng dẫn kỹ thuật đánh giá cụm dị thường phổ Gamma dựa trên bộ 8 chỉ tiêu định lượng (ΔJ, T(1/2), ΔTh/ΔU, ΔU/ΔK, JU, JTh, JK, F) cho 100% các đề án đo vẽ địa vật lý hàng không giai đoạn 2024–2026.
  3. Triển khai kiểm tra thực địa chi tiết: Đề nghị các đoàn địa chất khu vực Nam Trung Bộ tiến hành đo vẽ chi tiết mặt đất tỷ lệ 1:10.000 tại 3 cụm dị thường có chỉ số P*m > 0,80 ở huyện Khánh Vĩnh và Diên Khánh trong vòng 12 đến 18 tháng, nhằm xác định tọa độ các lỗ khoan thăm dò quặng vàng và kim loại hiếm.
  4. Đào tạo và chuyển giao công nghệ: Tổ chức 4 khóa tập huấn chuyên sâu về thuật toán Griffiths-Vinni và giải đoán phổ Gamma nâng cao cho ít nhất 60 kỹ sư địa vật lý trên toàn quốc trong năm 2025. Đơn vị chủ trì: Bộ môn Vật lý Địa cầu, Trường Đại học Khoa học Tự nhiên - ĐHQGHN.

Đối tượng nên tham khảo luận văn

Công trình luận văn là nguồn tài liệu giá trị cho 4 nhóm đối tượng chính:

  1. Kỹ sư và chuyên gia phân tích địa vật lý: Ứng dụng trực tiếp thuật toán hệ số tương quan di động và chương trình tính toán MATLAB để tự động hóa khâu xử lý tài liệu bay đo, giúp rút ngắn 35% thời gian xử lý dữ liệu thô tại thực địa.
  2. Nhà nghiên cứu và giảng viên chuyên ngành Địa chất - Địa vật lý: Khai thác cơ sở lý thuyết xạ địa hóa hạt nhân, bảng số liệu chuẩn Vinogradov và thuật toán nhận dạng mẫu để phục vụ công tác giảng dạy các học phần Địa vật lý máy bay, Xử lý số liệu và Thăm dò phóng xạ.
  3. Cán bộ quản lý tại các cơ quan tài nguyên và khoáng sản: Sử dụng bản đồ dự báo và chỉ số đồng dạng P*m để xây dựng quy hoạch cấp phép thăm dò, định hướng đầu tư điều tra khoáng sản chiến lược trên diện tích hơn 300 km2 tại khu vực Trung Bộ.
  4. Học viên cao học và sinh viên ngành Vật lý Địa cầu: Tham khảo phương pháp luận khoa học, cách tiếp cận xây dựng ma trận thông tin nhị phân và mô hình hóa toán học phục vụ thực hiện luận văn thạc sĩ và đồ án tốt nghiệp đạt chuẩn học thuật.

Câu hỏi thường gặp

  1. Phương pháp phổ Gamma hàng không hoạt động trên nguyên lý vật lý nào? Phương pháp đo cường độ bức xạ Gamma tự nhiên phát ra từ quá trình phân rã của Potassium (40K), Uranium (238U) và Thorium (232Th) trong đất đá. Máy thu đặt trên máy bay ghi nhận tín hiệu ở 3 cửa sổ năng lượng từ 1,35 đến 2,80 MeV, giúp xác định hàm lượng khoáng vật trên diện tích hàng trăm kilomet vuông với độ chính xác cao.

  2. Tại sao cần phân tích theo cụm dị thường thay vì các điểm dị thường đơn lẻ? Trong mỗi đề án bay đo tỷ lệ 1:50.000 có hàng nghìn điểm dị thường đơn lẻ khiến việc kiểm tra mặt đất tốn kém và dễ sai lệch. Việc gộp thành cụm dị thường giúp xác lập 8 tham số xạ địa hóa chung của cấu trúc địa chất, loại bỏ hơn 60% nhiễu ngẫu nhiên và tập trung kiểm chứng đúng các đới quặng tiềm năng.

  3. Thuật toán Griffiths-Vinni đóng vai trò gì trong việc xử lý dữ liệu? Thuật toán chuyển đổi toàn bộ số liệu đo phức tạp thành ma trận nhị phân 0 và 1 qua các cửa sổ quét đặc trưng, sau đó tính tỷ trọng thông tin tương đối Pm. Phương pháp này giúp tự động chọn lọc tổ hợp dấu hiệu mang hơn 80% dung lượng thông tin, từ đó tính chỉ số đồng dạng P*m với độ tin cậy vượt trội.

  4. Hệ số tương quan R giữa các nguyên tố phóng xạ chỉ thị điều gì? Ở các thành tạo đá nguyên thủy, hệ số tương quan giữa U, Th và K thường đạt mức cao từ 0,70 đến 0,90. Khi có đứt gãy kiến tạo hoặc biến đổi nhiệt dịch làm biến đổi khoáng vật, liên kết này bị phá vỡ khiến R giảm sâu xuống dưới 0,35, đóng vai trò là dấu hiệu nhận diện trực tiếp các đới quặng.

  5. Hiệu quả thực tế của phương pháp tại khu vực Khánh Hòa được chứng minh như thế nào? Thử nghiệm trên 300 km2 tại Khánh Vĩnh và Diên Khánh đã phát hiện 5/6 cụm dị thường triển vọng vàng đạt chỉ số đồng dạng P*m từ 0,78 đến 0,92. Kết quả này giúp khoanh vùng chuẩn xác các thân quặng, tiết kiệm khoảng 25% đến 30% chi phí khảo sát thực địa cho các đơn vị thăm dò.

Kết luận

  • Hệ thống hóa hoàn chỉnh cơ sở lý thuyết vật lý hạt nhân và địa hóa phóng xạ, thiết lập chuẩn xác 3 cửa sổ năng lượng đặc trưng của Potassium, Uranium và Thorium.
  • Xây dựng và kiểm chứng thành công phương pháp Hệ số tương quan đa biến di động (RU/Th, RU/K, RTh/K) nhằm định vị chính xác các đới đứt gãy và biến đổi nhiệt dịch.
  • Hoàn thiện thuật toán Tần suất - Nhận dạng Griffiths-Vinni trên phần mềm QTS, tối ưu hóa ma trận tham số đạt tỷ trọng thông tin tích lũy 85,4%.
  • Đánh giá thực nghiệm xuất sắc trên diện tích 300 km2 tại Khánh Hòa, nhận dạng chính xác 6 cụm dị thường quặng vàng triển vọng với chỉ số đồng dạng P*m đạt từ 0,78 đến 0,92.
  • Phát triển bộ chương trình tính toán tự động bằng ngôn ngữ MATLAB, tạo bước đột phá trong việc chuyển dịch từ giải đoán định tính sang phân tích định lượng hiện đại.

Công trình luận văn của tác giả Tạ Quang Hân dưới sự hướng dẫn của PGS. Võ Thanh Quỳnh là đóng góp học thuật và công nghệ xuất sắc cho ngành Vật lý Địa cầu Việt Nam. Kế hoạch hoàn thiện và chuyển giao bộ công cụ phân tích tự động trên toàn quốc dự kiến hoàn thành trong giai đoạn 2024–2026. Các tổ chức khoa học và doanh nghiệp địa chất quan tâm có thể khai thác tài liệu luận văn để ứng dụng ngay vào các đề án điều tra, đánh giá triển vọng khoáng sản thực tế.