Tổng quan nghiên cứu

Công tác đo bay địa vật lý tỷ lệ lớn (từ - phổ gamma hàng không) tại Việt Nam đã trải qua hơn 30 năm hình thành và phát triển, đóng vai trò nền tảng trong công tác điều tra cơ bản địa chất và thăm dò tài nguyên. Tuy nhiên, thách thức lớn nhất trong thực tiễn là khâu xử lý và minh giải khối lượng dữ liệu khổng lồ. Trong mỗi đề án bay đo, số lượng dị thường điểm đơn lẻ thường lên tới hàng nghìn điểm, trong khi công tác kiểm tra thực địa mặt đất chỉ có thể đáp ứng khoảng 5% đến 10% các vị trí trọng điểm do hạn chế về kinh phí và địa hình hiểm trở.

Vấn đề nghiên cứu trọng tâm của luận văn là khắc phục tính chủ quan trong phân tích trực quan truyền thống bằng việc ứng dụng các thuật toán xử lý thống kê đa biến định lượng. Mục tiêu cụ thể là nghiên cứu cơ sở khoa học và áp dụng thử nghiệm hai phương pháp tiên tiến: Phương pháp Hệ số tương quan và Phương pháp Tần suất - Nhận dạng nhằm nâng cao độ chính xác trong phân loại cụm dị thường và dự báo khoáng sản.

Phạm vi nghiên cứu được triển khai thực nghiệm trên diện tích khoảng 250 km2, bao gồm 3 xã thuộc huyện Diên Khánh và 5 xã thuộc huyện Cam Lâm, tỉnh Khánh Hòa (khu vực Đông Nam thành phố Nha Trang), dựa trên nguồn số liệu bay đo tỷ lệ 1/50.000 hoàn thành năm 2015. Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc khi giúp tăng độ tin cậy khoanh định ranh giới các đới biến đổi nhiệt dịch chứa quặng (vàng, thiếc, đa kim) lên trên 85%, đồng thời giảm thiểu khoảng 30% chi phí khảo sát kiểm tra mặt đất không cần thiết.

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 xử lý tổ hợp số liệu địa vật lý và lý thuyết nhận dạng mẫu (Pattern Recognition). Trong địa vật lý thăm dò, các trường quan sát luôn chứa đựng thành phần ngẫu nhiên do nhiễu địa hình, thảm phủ và sai số thiết bị. Do đó, mô hình thống kê đóng vai trò quyết định trong việc phân loại trường xạ thành hai nhóm đối tượng chính: nhóm nhận dạng có mẫu chuẩn (áp dụng khi đã biết trước đặc trưng thống kê của mỏ hoặc đới khoáng hóa) và nhóm nhận dạng không có mẫu chuẩn (phân chia trường tự động dựa trên mức độ đồng nhất).

Lý thuyết phân bố trường xạ địa hóa tự nhiên của ba nguyên tố phóng xạ chủ đạo gồm Uran (U), Thori (Th) và Kali (K) là tiền đề vật lý cốt lõi. Trong điều kiện địa chất ổn định, tỷ lệ hàm lượng giữa các nguyên tố này duy trì mối tương quan tuyến tính chặt chẽ. Khi xuất hiện các quá trình biến đổi nhiệt dịch hoặc khoáng hóa quặng, sự cân bằng này bị phá vỡ cục bộ. Luận văn vận dụng 5 khái niệm chuyên ngành then chốt: cụm dị thường phổ gamma, hệ số tương quan tuyến tính Pearson ($R$), chỉ số nhiều thành phần ($F = U \cdot Th / K$), tỷ trọng thông tin ($P_m$) theo thuật toán Griffiths - Vinni, và chỉ số đồng dạng ($P^*_m$) giữa đối tượng nghiên cứu với các cấu trúc quặng chuẩn.

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

Nguồn dữ liệu đầu vào sử dụng bộ bản đồ trường phóng xạ tỷ lệ 1/50.000 gồm 4 kênh số liệu gốc: bản đồ hàm lượng Uran, Thori, Kali và bản đồ cường độ phóng xạ kênh tổng (Tc), kết hợp các tài liệu kiểm chứng mặt đất từ đề án bay đo vùng Phan Rang - Nha Trang. Các tham số trung gian được tính toán bổ sung gồm $U/Th$, $Th/K$, $U/K$ và tỷ lệ liều lượng từng thành phần so với kênh tổng.

Về quy mô mẫu và phương pháp chọn mẫu, nghiên cứu thiết lập lưới điểm đều trên toàn bộ diện tích 250 km2 với hơn 1.200 nút mạng quan sát. Phương pháp lấy mẫu theo cơ chế quét cửa sổ trượt (sliding window) di động liên tục với tâm cửa sổ là từng điểm nút mạng. Lý do lựa chọn phương pháp này là nhằm bảo đảm tính liên tục của trường dữ liệu, đồng thời triệt tiêu hiệu quả các dị thường giả cục bộ.

Quy trình phân tích kết hợp hai kỹ thuật:

  1. Thuật toán hệ số tương quan ($R_{U/Th}, R_{U/K}, R_{Th/K}$) để xác lập 3 lớp dữ liệu không gian mới phản ánh sự phá hủy tương quan trường xạ.
  2. Phương pháp Tần suất - Nhận dạng đối sánh với 2 đối tượng mẫu chuẩn đã biết rõ bản chất địa chất là mỏ Trà Năng và Ma Ty Du Long. Toàn bộ quy trình tính toán được lập trình tự động hóa và xử lý hoàn tất trong giai đoạn nghiên cứu năm 2015.

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

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

Quá trình xử lý thử nghiệm số liệu thực tế tại khu vực Diên Khánh và Cam Lâm đã mang lại các kết quả định lượng cụ thể:

Thứ nhất, phân tích hệ số tương quan trên toàn diện tích đã xác lập được 3 hệ lớp bản đồ đẳng trị $R_{U/Th}, R_{U/K}$ và $R_{Th/K}$. Tại các diện tích địa chất đồng nhất, hệ số tương quan đạt giá trị cao từ 0,75 đến 0,92. Ngược lại, nghiên cứu đã phát hiện và khoanh định chính xác 4 đới biến đổi trường xạ địa hóa cục bộ có hệ số tương quan sụt giảm mạnh xuống dưới 0,35, phản ánh sự tái phân bố các nguyên tố phóng xạ do hoạt động magma - nhiệt dịch.

Thứ hai, phương pháp Tần suất - Nhận dạng với ngưỡng tỷ trọng thông tin $P_m \ge 75%$ đã sàng lọc từ 8 thông số đầu vào ban đầu để chọn ra 4 tham số có lượng tin cao nhất, chiếm 82,4% tổng lượng thông tin hữu ích của đối tượng mẫu.

Thứ ba, kết quả nhận dạng đối sánh theo mẫu Trà Năng và mẫu Ma Ty Du Long đã phân định rõ rệt diện tích có triển vọng khoáng sản cao với chỉ số đồng dạng $P^*_m > 70%$. Đặc biệt, tại đới Núi Go.B, ranh giới triển vọng quặng hóa vàng (Au), thiếc (Sn) và sunfur đa kim được khoanh vùng với mức độ trùng khớp lên đến 88% so với tài liệu kiểm tra địa chất mặt đất thực tế.

+-----------------------------------------------------------------------------------+
|               KẾT QUẢ PHÂN TÍCH TỔ HỢP ĐỊA VẬT LÝ MÁY BAY KHU VỰC NGHIÊN CỨU      |
+----------------------+--------------------+---------------------+-----------------+
| Đối tượng phân tích  | Phương pháp áp dụng| Giá trị tham số     | Đánh giá        |
+----------------------+--------------------+---------------------+-----------------+
| Vùng nền ổn định     | Hệ số tương quan R | R = 0,75 - 0,92     | Đồng nhất       |
| Đới biến đổi nhiệt   | Hệ số tương quan R | R < 0,35            | Dị thường cao   |
| Mẫu chuẩn Trà Năng   | Tần suất - Nhận dạng| P*m > 75%          | Đồng dạng quặng |
| Vùng quặng Núi Go.B  | Phân tích tổ hợp   | Khớp 88% thực địa   | Triển vọng Au-Sn|
+----------------------+--------------------+---------------------+-----------------+

Thảo luận kết quả

Các phát hiện nghiên cứu được trực quan hóa thông qua hệ thống bản đồ số và đồ thị phân tán hồi quy đa biến trên nền tảng hệ thông tin địa lý GIS và MapInfo. Biểu đồ phân tán cho thấy sự phân tách rõ rệt giữa các điểm phông môi trường và các điểm dị thường quặng.

Nguyên nhân căn bản của sự suy giảm hệ số tương quan tại các đới triển vọng là do tác động của dung dịch nhiệt dịch chứa quặng xâm tán làm giàu bất thường nguyên tố Uran hoặc Kali trong khi Thori tương đối trơ và cố định. So sánh với phương pháp truyền thống xử lý trên bộ phần mềm CODCAD trước đây, phương pháp Hệ số tương quan kết hợp Tần suất - Nhận dạng cho phép vạch rõ ranh giới biên chi tiết hơn 35%, không bị giới hạn bởi việc gộp cụm chủ quan. Kết quả này chứng minh tính ưu việt vượt trội của việc số hóa tự động trong xử lý tài liệu địa vật lý máy bay tại Việt Nam.

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

Dựa trên kết quả thực nghiệm thành công, luận văn đề xuất 4 nhóm giải pháp cụ thể:

  1. Chuẩn hóa quy trình xử lý dữ liệu: Đề nghị Tổng cục Địa chất và Khoáng sản Việt Nam ban hành hướng dẫn kỹ thuật áp dụng quy trình xử lý 2 bước (quét hệ số tương quan sơ bộ kết hợp nhận dạng tần suất chi tiết) cho 100% các đề án bay đo từ - phổ gamma tỷ lệ 1/50.000 và 1/25.000 trên toàn quốc, thực hiện trong giai đoạn 2016-2020.
  2. Xây dựng ngân hàng mẫu chuẩn số hóa: Bộ Tài nguyên và Môi trường cần chủ trì xây dựng cơ sở dữ liệu số về đặc trưng thống kê trường xạ cho ít nhất 20 mỏ khoáng sản điển hình (vàng, uran, đất hiếm, thiếc, sa khoáng) tại Việt Nam trước năm 2022 nhằm cung cấp thư viện mẫu chuẩn chính xác cho các thuật toán nhận dạng.
  3. Nâng cấp phần mềm chuyên dụng: Các nhóm nghiên cứu tại Trường Đại học Khoa học Tự nhiên tiếp tục module hóa thuật toán Tần suất - Nhận dạng thành phần mềm thương mại tích hợp trực tiếp với các nền tảng GIS hiện đại, đặt mục tiêu rút ngắn thời gian xử lý dữ liệu từ 30 ngày xuống dưới 7 ngày làm việc.
  4. Tối ưu hóa kinh phí khảo sát thực địa: Các liên đoàn địa chất cần sử dụng bản đồ chỉ số đồng dạng $P^*_m > 75%$ làm tiêu chí tiên quyết để định vị mạng lưới hào, khoan thăm dò và đo trọng sa, giúp nâng tỷ lệ phát hiện quặng thực tế lên trên 40% và tiết kiệm 25% ngân sách nhà nước.

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

Nội dung luận văn mang giá trị ứng dụng cao và là tài liệu tham khảo thiết thực cho 4 nhóm đối tượng:

  • Kỹ sư và chuyên gia địa vật lý: Khai thác quy trình xử lý số liệu tự động, sử dụng bộ công cụ hệ số tương quan để phân loại cụm dị thường và xác định ranh giới đứt gãy kiến tạo ẩn.
  • Học viên cao học và nghiên cứu sinh ngành Vật lý địa cầu, Địa chất học: Tiếp cận phương pháp luận nghiên cứu khoa học chuẩn mực, khung lý thuyết nhận dạng mẫu và các thuật toán thống kê đa biến ứng dụng trong khoa học Trái Đất.
  • Cơ quan quản lý nhà nước về địa chất và khoáng sản: Tham khảo cơ sở định lượng để thẩm định, phê duyệt và đánh giá hiệu quả kinh tế - kỹ thuật của các dự án bay đo địa vật lý hàng không.
  • Chuyên gia phát triển phần mềm khoa học: Ứng dụng các thuật toán giải mã ma trận thông tin logic Griffiths - Vinni và thuật toán cửa sổ trượt để xây dựng các phần mềm khai phá dữ liệu không gian.

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

Hỏi: Hệ số tương quan trong phân tích phổ gamma hàng không có ý nghĩa vật lý như thế nào? Đáp: Hệ số tương quan phản ánh mức độ gắn kết địa hóa giữa các nguyên tố U, Th, K. Khi đất đá ổn định, hệ số này tiến gần đến 1. Khi có đới biến đổi quặng hóa nhiệt dịch làm dịch chuyển cục bộ một nguyên tố, hệ số tương quan sẽ sụt giảm mạnh xuống gần 0, giúp phát hiện đới chứa quặng.

Hỏi: Điểm khác biệt cốt lõi giữa nhận dạng có mẫu chuẩn và không mẫu chuẩn là gì? Đáp: Nhận dạng có mẫu chuẩn so sánh dữ liệu khảo sát với một mỏ quặng đã biết rõ để tìm đối tượng tương đồng. Nhận dạng không mẫu chuẩn tự động gom cụm các điểm quan sát có cùng đặc trưng thống kê thành các lớp trường riêng biệt khi chưa có thông tin kiểm chứng trước.

Hỏi: Chỉ số đồng dạng $P^*_m$ trong phương pháp Tần suất - Nhận dạng được xác định ra sao? Đáp: $P^_m$ được tính toán bằng cách đối sánh ma trận thông tin nhị phân của đối tượng nghiên cứu với tập hợp các tham số đặc trưng nhất của đối tượng mẫu. Giá trị $P^_m \ge 75%$ khẳng định đối tượng nghiên cứu có cấu trúc trường xạ tương đồng với mỏ chuẩn.

Hỏi: Tại sao khu vực Diên Khánh và Cam Lâm (Khánh Hòa) được chọn để thử nghiệm? Đáp: Khu vực này có cấu trúc địa chất phức tạp, phát triển nhiều khối xâm nhập magma và đã có tài liệu bay đo tỷ lệ 1/50.000 cùng các số liệu kiểm chứng thực tế tại mỏ Núi Go.B, tạo điều kiện lý tưởng để đối sánh độ chính xác của thuật toán.

Hỏi: Phương pháp này có thể áp dụng cho các khoáng sản không chứa phóng xạ không? Đáp: Có thể áp dụng hiệu quả. Nhiều loại khoáng sản như vàng, đồng, chì, kẽm, thiếc thường đi kèm với quá trình biến đổi nhiệt dịch làm giàu Kali hoặc biến đổi tỷ lệ Thori/Uran, tạo thành dấu hiệu trường xạ chỉ thị gián tiếp rất rõ ràng.

Kết luận

  • Luận văn đã hệ thống hóa toàn diện cơ sở lý thuyết xử lý tổ hợp số liệu địa vật lý và nhận dạng mẫu trong thăm dò khoáng sản.
  • Xây dựng thành công quy trình ứng dụng phương pháp Hệ số tương quan đa biến, phát hiện 4 đới biến đổi trường xạ địa hóa đặc trưng tại khu vực thử nghiệm.
  • Ứng dụng xuất sắc phương pháp Tần suất - Nhận dạng với 2 mẫu chuẩn Trà Năng và Ma Ty Du Long, sàng lọc 4 tham số có lượng tin chiếm 82,4%.
  • Khoanh định chuẩn xác ranh giới đới triển vọng quặng hóa đa kim Núi Go.B với độ trùng khớp 88% so với tài liệu thực tế.
  • Khẳng định giá trị thực tiễn trong việc giảm 25% đến 30% chi phí khảo sát mặt đất và nâng cao năng lực chuyển đổi số ngành địa chất.

Để tiếp tục hoàn thiện, các nghiên cứu tiếp theo cần mở rộng tích hợp thêm dữ liệu từ hàng không và trọng lực mô hình 3D. Độc giả và các đơn vị chuyên môn quan tâm có thể khai thác trực tiếp giải thuật và bộ số liệu mẫu được trình bày chi tiết trong luận văn để ứng dụng vào thực tế sản xuất.