Tổng quan nghiên cứu

Công tác đo bay địa vật lý đóng vai trò chiến lược trong việc điều tra cơ bản địa chất và thăm dò khoáng sản tại Việt Nam. Khối lượng dữ liệu thu thập từ các đề án bay đo từ - phổ gamma tỷ lệ 1/50.000 rất phong phú nhưng cũng đặt ra thách thức lớn về xử lý và chiết xuất thông tin. Khi số lượng tham số đầu vào tăng cao, việc sử dụng đồng thời toàn bộ dữ liệu thường gây nhòe thông tin do nhiễu địa chất, làm giảm độ tin cậy của kết quả dự báo.

Luận văn thạc sĩ khoa học chuyên ngành Vật lý địa cầu (Mã số: 60440111) của học viên Phạm Bảo Chi thực hiện năm 2015 tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội đã giải quyết trực tiếp bài toán này. Đề tài tập trung nghiên cứu, hoàn thiện và ứng dụng phương pháp Tần suất – Nhận dạng mới nhằm tối ưu hóa việc phân loại và dự báo khoáng sản. Phạm vi nghiên cứu thực nghiệm được triển khai trên diện tích khoảng 600 km2 tại khu vực các huyện Ninh Phước và Thuận Nam, tỉnh Ninh Thuận (giới hạn từ 11°21'34'' đến 11°33'00'' vĩ độ Bắc và 108°44'49'' đến 109°00'25'' kinh độ Đông).

Nghiên cứu mang ý nghĩa thực tiễn sâu sắc khi cung cấp một công cụ toán học tự động hóa cao, giúp giảm thiểu 25% đến 30% chi phí khảo sát mặt đất không cần thiết. Phương pháp cho phép khoanh định chính xác các đới khoáng hóa thiếc, vonfram và sulfur tiềm năng, nâng cao hiệu quả kinh tế - kỹ thuật cho các đề án địa vật lý hàng không tại Việt Nam.

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 nền tảng lý thuyết nhận dạng mẫu trong địa vật lý và lý thuyết xử lý tổ hợp số liệu thống kê đa chiều. Bản chất của nhận dạng địa vật lý là bài toán phân loại các điểm quan sát vào các lớp đối tượng địa chất tương ứng dựa trên tập hợp dấu hiệu trường đo được.

Khung lý thuyết vận dụng 2 mô hình toán học cốt lõi:

  1. Lý thuyết phân phối xác suất: Sử dụng hàm mật độ xác suất và luật phân phối chuẩn Gaussian để mô tả các đại lượng ngẫu nhiên địa vật lý như hàm lượng phóng xạ, độ từ cảm và mật độ đất đá.
  2. Thuật toán phân tích tần suất Griffiths – Vinni: Đánh giá chất lượng của từng chủng loại thông tin thông qua tần suất xuất hiện đồng thời của các dấu hiệu trên đối tượng mẫu, từ đó xác định tỷ trọng thông tin tương đối.

Các khái niệm chính được làm rõ bao gồm: Ma trận thông tin nhị phân (mã hóa 0 và 1 qua cửa sổ quét giá trị đặc trưng), Tỷ trọng thông tin tương đối, Ngưỡng tỷ trọng tích lũy thông tin và Chỉ số đồng dạng phục vụ việc ra quyết định.

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

Nguồn dữ liệu đầu vào của nghiên cứu kế thừa từ Đề án Bay đo từ - phổ gamma tỷ lệ 1/50.000 vùng Phan Rang - Nha Trang thực hiện năm 2003. Bộ dữ liệu khảo sát bao gồm 11 tham số đặc trưng: hàm lượng Uran, Thori, Kali, Kênh tổng, hệ số F và 6 tỷ số hàm lượng phóng xạ liên quan.

Cỡ mẫu và phương pháp chọn mẫu:

  • Mẫu chuẩn đại diện: Vùng Nhị Hà (diện tích ký hiệu số 6 trên bản đồ phân vùng) được chọn làm đối tượng mẫu chuẩn có tiền đề địa chất rõ ràng về quặng thiếc - vonfram - sulfur loại B, nằm tại tọa độ 11°27'33'' Bắc và 108°53'55'' Đông.
  • Mẫu đối sánh: 5 đối tượng triển vọng trong khu vực (từ số 1 đến số 5) và toàn bộ các điểm nút của mạng lưới nội suy đều trên diện tích 600 km2.

Quy trình phân tích được thực hiện qua các bước:

  1. Xây dựng đường cong biến phân thực nghiệm để xác định khoảng giá trị đặc trưng cho 11 tham số của mẫu chuẩn.
  2. Chuyển đổi dữ liệu thực tế sang ma trận thông tin logic nhị phân.
  3. Áp dụng thuật toán Griffiths – Vinni để sàng lọc tập hợp tham số tối ưu với ngưỡng tích lũy thông tin đạt từ 80% trở lên.
  4. Quét cửa sổ di động đa giác trên mạng lưới nội suy để tính toán chỉ số đồng dạng trên toàn bộ diện tích thông qua phần mềm QTS, tự động khoanh định ranh giới quặng mà không cần biết trước đối tượng đối sánh.

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 địa vật lý máy bay vùng Ninh Thuận đã mang lại 3 phát hiện khoa học quan trọng:

  • Sàng lọc tổ hợp thông tin tối ưu: Từ 11 tham số trường xạ ban đầu, phương pháp Tần suất – Nhận dạng mới đã xác định được 6 tính chất đặc trưng có tỷ trọng thông tin cao nhất bao gồm: tỷ số K/Tg (tỷ trọng đạt 0.49927), hàm lượng Kali, hệ số F, tỷ số U/Th, Th/Tg và U/Tg. Tập hợp 6 tính chất này tích lũy trên 80% lượng thông tin hữu ích của đối tượng mẫu, cho phép loại bỏ 5 tham số kém đặc trưng, giúp giảm thiểu 45% nhiễu thông tin nền.

  • Đánh giá định lượng mức độ đồng dạng: Khi đối sánh 5 vùng triển vọng với mẫu chuẩn Nhị Hà, Đối tượng số 1 đạt chỉ số đồng dạng cao vượt trội lên tới 90.4%. Các đối tượng còn lại có mức độ đồng dạng dao động từ 48.6% đến 72.3%, chứng minh Đối tượng số 1 có cấu trúc trường địa vật lý tương đồng nhất với mỏ quặng Nhị Hà.

  • Tự động xác định và hiệu chỉnh đường biên quặng: Thuật toán quét cửa sổ di động đã khoanh định chính xác ranh giới các vùng đồng dạng trên diện tích 600 km2. Kết quả giúp tái lập và thu hẹp đường biên của các vùng triển vọng ban đầu khoảng 15% đến 20%, loại bỏ các diện tích dị thường giả do địa hình hoặc phong hóa bề mặt gây ra.

Thảo luận kết quả

Các kết quả nghiên cứu được minh chứng và trực quan hóa rõ nét thông qua hệ thống bảng biểu và bản đồ chuyên đề:

  • Bảng ma trận đối sánh thông tin thể hiện chi tiết sự phân bố giá trị nhị phân giữa các đối tượng.
  • Sơ đồ đường đồng mức chỉ số đồng dạng dạng lưới 2D/3D cho phép các nhà địa chất quan sát trực quan sự biến thiên không gian của đới khoáng hóa.

Nguyên nhân cốt lõi giúp phương pháp Tần suất – Nhận dạng mới đạt độ chính xác cao là việc kết hợp chặt chẽ giữa đặc trưng biến thiên hàm lượng xạ và bản chất địa chất nhiệt dịch của mỏ quặng thiếc - vonfram. Các tỷ số xạ Kali/Kênh tổng và Uran/Thori phản ánh trực tiếp quá trình biến đổi đá vây quanh và kali hóa đi kèm khoáng hóa quặng.

So với các phương pháp nhận dạng truyền thống trong bộ công cụ COSCAD hoặc hồi quy đa biến vốn đòi hỏi phải có sẵn ranh giới đối tượng so sánh, phương pháp mới mang tính đột phá khi tự động dò tìm và khoanh định đối tượng chưa biết trên toàn vùng khảo sát.

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

Dựa trên kết quả nghiên cứu, 4 giải pháp chiến lược được đề xuất nhằm nâng cao hiệu quả công tác thăm dò khoáng sản:

  • Ứng dụng chuẩn hóa phần mềm: Liên đoàn Địa vật lý Địa chất và các đơn vị thăm dò thuộc Tổng cục Địa chất và Khoáng sản Việt Nam cần đưa quy trình phân tích Tần suất – Nhận dạng mới (chương trình QTS, QTSM) vào quy chuẩn xử lý tài liệu bay đo giai đoạn 2026-2028, hướng tới mục tiêu tự động hóa 100% công tác phân tích sơ bộ dữ liệu từ - phổ gamma.

  • Triển khai kiểm tra thực địa trọng điểm: Sở Tài nguyên và Môi trường tỉnh Ninh Thuận phối hợp với các doanh nghiệp khai khoáng ưu tiên cấp phép và thực hiện công tác kiểm tra mặt đất, khoan thăm dò chi tiết tại Đối tượng số 1 (đạt chỉ số đồng dạng 90.4%) theo mạng lưới 100m x 50m trong thời gian 6 đến 12 tháng tới.

  • Tích hợp đa phương pháp địa vật lý: Nhóm nghiên cứu khuyến nghị tích hợp thêm dữ liệu trọng lực tỷ lệ 1/100.000 và ảnh viễn thám vệ tinh vào ma trận thông tin nhận dạng để nâng ngưỡng tin cậy tích lũy thông tin từ 80% lên mức 90%, nâng cao độ sâu nghiên cứu.

  • Xây dựng ngân hàng mẫu chuẩn địa vật lý: Các viện nghiên cứu và trường đại học cần phối hợp thu thập, chuẩn hóa cơ sở dữ liệu mẫu chuẩn cho các loại hình khoáng sản chiến lược (như đất hiếm, bauxit, titan sa khoáng) trên toàn quốc trong vòng 24 tháng nhằm phục vụ nhận dạng tự động trên quy mô vùng.

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

Luận văn là tài liệu tham khảo giá trị cho 4 nhóm đối tượng chính:

  • Kỹ sư và chuyên gia địa vật lý - địa chất: Nắm vững quy trình xử lý dữ liệu bay đo, áp dụng thuật toán Griffiths – Vinni để lọc bỏ thông tin nhiễu, nâng cao độ chính xác khi lập bản đồ dự báo khoáng sản tại các vùng mỏ mới.

  • Cơ quan quản lý nhà nước về khoáng sản: Sử dụng các sơ đồ khoanh vùng triển vọng để quy hoạch phân vùng thăm dò, tối ưu hóa phân bổ ngân sách nhà nước và rút ngắn 20% thời gian đánh giá tiềm năng tài nguyên.

  • Giảng viên, nghiên cứu sinh và học viên ngành Khoa học Trái đất: Sử dụng công trình làm tài liệu giảng dạy chuyên sâu về môn học Xử lý số liệu địa vật lý, Lý thuyết nhận dạng mẫu và Thống kê địa chất ứng dụng.

  • Chuyên gia phân tích dữ liệu không gian và GIS: Tham khảo thuật toán quét cửa sổ di động và ma trận nhị phân để phát triển các module mở rộng trên nền tảng phần mềm mã nguồn mở phục vụ viễn thám và trắc địa.

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

  • Điểm khác biệt lớn nhất giữa phương pháp Tần suất – Nhận dạng mới và phương pháp cũ là gì? Phương pháp mới khắc phục hoàn toàn nhược điểm của phương pháp cũ nhờ khả năng tự động khoanh vùng và xác định đường biên đối tượng đồng dạng trên toàn bộ diện tích mạng lưới mà không cần biết trước ranh giới của các đối tượng đối sánh.

  • Tại sao lại chọn 6 tính chất đặc trưng thay vì sử dụng toàn bộ 11 tham số đo được? Việc sử dụng toàn bộ 11 tham số sẽ đưa cả các biến số có lượng tin thấp vào mô hình, gây nhiễu dữ liệu. Chọn 6 tính chất có tỷ trọng thông tin tích lũy đạt 80% giúp tối ưu hóa thuật toán và tăng độ chính xác phân loại.

  • Đối tượng số 1 tại Ninh Thuận có tiềm năng quặng hóa cụ thể ra sao? Đối tượng số 1 có chỉ số đồng dạng 90.4% so với mỏ mẫu Nhị Hà, đồng thời có các chỉ số phóng xạ Kali và hệ số F dị thường rất cao, là dấu hiệu địa vật lý điển hình của khoáng hóa thiếc, vonfram và sulfur nhiệt dịch.

  • Phương pháp này có thể áp dụng cho các tài liệu đo đạc khác ngoài phổ gamma không? Hoàn toàn có thể. Thuật toán có tính mở cao, xử lý tốt ma trận dữ liệu từ trường, trọng lực, điện trở suất hoặc tài liệu địa hóa mặt đất với số lượng tham số đầu vào không hạn chế.

  • Phần mềm QTS xử lý dữ liệu trên mạng lưới nội suy theo cơ chế nào? Phần mềm sử dụng cửa sổ quét di động (hình tròn, vuông hoặc elip) trượt đều qua từng điểm nút trên mạng lưới, tính toán chỉ số đồng dạng tại tâm cửa sổ và xuất ra tệp dữ liệu không gian để vẽ bản đồ đẳng trị.

Kết luận

  • Luận văn đã hoàn thiện thành công cơ sở toán học và quy trình công nghệ của phương pháp Tần suất – Nhận dạng mới phục vụ xử lý tài liệu địa vật lý hàng không.
  • Nghiên cứu đã tối ưu hóa việc lựa chọn 6 tham số trường xạ then chốt, loại bỏ hiệu quả các nhiễu địa chất cục bộ trên diện tích nghiên cứu 600 km2.
  • Kết quả thực nghiệm tại Ninh Thuận đã xác định chuẩn xác Đối tượng số 1 có độ đồng dạng 90.4% với mẫu chuẩn quặng Nhị Hà và tái lập thành công ranh giới mỏ.
  • Thuật toán được tự động hóa hoàn toàn qua phần mềm QTS, mở ra hướng ứng dụng rộng rãi cho công tác bay đo thăm dò địa vật lý trên toàn quốc giai đoạn 2026-2030.
  • Các đơn vị điều tra địa chất nên sớm áp dụng phương pháp này để đẩy nhanh tiến độ khoanh vùng và nâng cao hiệu quả thăm dò khoáng sản ẩn sâu.