Tổng quan nghiên cứu

Trong công tác thăm dò khoáng sản và nghiên cứu cấu trúc vỏ Trái Đất, tài liệu địa vật lý hàng không đóng vai trò nền tảng với khả năng bao phủ diện tích rộng lớn lên đến hàng nghìn km² trong thời gian ngắn. Tại Việt Nam, nguồn dữ liệu bay đo phổ gamma và từ trọng lực tỷ lệ 1/50.000 do Liên đoàn Vật lý Địa chất thu thập rất phong phú nhưng việc khai thác định lượng vẫn còn nhiều hạn chế, dẫn đến nguy cơ bỏ sót các dị thường quặng có quy mô nhỏ. Nghiên cứu của học viên Nguyễn Thị Thu Hiền dưới sự hướng dẫn của Phó Giáo sư Võ Thanh Quỳnh tập trung giải quyết bài toán cốt lõi: tối ưu hóa thuật toán xử lý đa biến nhằm nâng cao độ chính xác trong phân tích trường phóng xạ. Mục tiêu cụ thể là nghiên cứu áp dụng phương pháp Khoảng cách – Tần suất – Nhận dạng mới để khoanh định chính xác các diện tích triển vọng khoáng sản và phân định ranh giới địa chất phức tạp. Địa bàn thử nghiệm thực tế được triển khai tại xã Khánh Thượng, huyện Khánh Vĩnh, tỉnh Khánh Hòa với diện tích 209,42 km² thuộc cấu trúc rìa khối nâng Kon Tum có lịch sử kiến tạo hơn 570 triệu năm. Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc thiết lập một quy trình xử lý tự động có khả năng sàng lọc 15 tham số trường xạ, xác định 8 đặc trưng tối ưu và phân loại chuẩn xác 9 đối tượng khảo sát. Kết quả này giúp giảm thiểu chi phí khoan thăm dò từ 20% đến 30% và nâng cao độ tin cậy trong công tác lập bản đồ địa chất khoáng sản khu vực Nam Trung Bộ.

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

Khung lý thuyết áp dụng

Luận văn xây dựng trên nền tảng tích hợp của hai lý thuyết toán - địa vật lý hiện đại: lý thuyết khoảng cách khái quát do Paguonov đề xuất và lý thuyết phân tích tần suất theo thuật toán Griffiths - Vinni, kết hợp mô hình nhận dạng thống kê có mẫu chuẩn của Phó Giáo sư Võ Thanh Quỳnh. Khung lý thuyết này vận hành dựa trên 4 khái niệm cốt lõi: ma trận thông tin nhị phân chuẩn hóa logic (0 và 1), khoảng giá trị đặc trưng của đối tượng mẫu, bình phương khoảng cách khái quát giữa trọng tâm các đám mây dấu hiệu, và chỉ số đồng dạng thông tin tổng hợp. Bản chất của mô hình là chuyển đổi các trường vật lý có phân bố ngẫu nhiên phức tạp về không gian dấu hiệu định lượng, từ đó đánh giá chính xác lượng tin của từng kênh phổ gamma và triệt tiêu các thành phần nhiễu phông địa phương.

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

Nghiên cứu sử dụng nguồn số liệu trường phổ xạ gồm 4 kênh đo gốc: hàm lượng Uran (ppm), Thori (ppm), Kali (%), và kênh tổng bức xạ Tg (R/s) trích xuất từ đề án bay đo tỷ lệ 1/50.000 vùng Phan Rang - Nha Trang năm 2003. Từ 4 tham số gốc, tác giả tính toán bổ sung 11 tham số biến đổi trung gian gồm các tỷ số U/Th, U/K, Th/K, U/Tg, Th/Tg, K/Tg, hệ số F = U.Th/K, Th.K/U, U.Th/Tg, U.K/Tg và Th.K/Tg, tạo thành bộ 15 biến khảo sát. Cỡ mẫu nghiên cứu bao gồm 9 cụm dị thường triển vọng khoáng sản vàng loại B và các vùng đối chứng không chứa quặng trên mạng lưới điểm đo nội suy ô vuông. Phương pháp chọn mẫu áp dụng nguyên tắc chọn mẫu đối nghịch có chủ đích: sử dụng đối tượng số 6 làm mẫu chuẩn quặng và chọn một diện tích không có biểu hiện khoáng hóa làm mẫu đối nghịch. Lý do lựa chọn tổ hợp thuật toán Khoảng cách – Tần suất – Nhận dạng là vì phương pháp này cho phép tự động phân tách các dấu hiệu có độ phân biệt cực đại, loại bỏ các biến số dư thừa mà không phụ thuộc vào phân bố chuẩn ngặt nghèo của dữ liệu thực địa. Timeline xử lý và thử nghiệm được hoàn thành đồng bộ trong giai đoạn nghiên cứu từ năm 2015 đến năm 2016.

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 số liệu thực tế vùng Khánh Thượng đã đem lại 3 phát hiện khoa học quan trọng có giá trị thực tiễn cao:

Thứ nhất, qua phân tích khoảng cách khái quát trên 15 tham số, thuật toán đã tự động sàng lọc và xác định được tổ hợp m = 8 tính chất có khả năng phân biệt tối ưu giữa đối tượng quặng và không quặng, gồm: U/Th, K/Tg, K, U/Tổng, Tổng, U.K/Tổng, Th.K/Tg và Thori. Trong đó, tính chất Thori thể hiện độ nhạy phân tách vượt trội với giá trị bình phương khoảng cách khái quát đạt mức kỷ lục 526,8808.

Thứ hai, việc tính toán chỉ số đồng dạng Pm* trên 8 đối tượng đối sánh so với đối tượng mẫu số 6 đạt mức độ tương đồng rất cao, dao động từ 69,58% đến 82,30%. Cụ thể, đối tượng 2 đạt độ tương đồng cao nhất là 82,30%, đối tượng 3 đạt 81,80%, đối tượng 4 đạt 76,78%, đối tượng 1 đạt 75,90%, đối tượng 9 đạt 75,10%, đối tượng 5 đạt 74,48%, đối tượng 7 đạt 72,73%, và đối tượng 8 đạt 69,58%. Mức tương quan trên 70% ở 7/8 đối tượng khẳng định tính đồng nhất về cấu trúc sinh quặng trong toàn khu vực.

Thứ ba, phương pháp đã giải quyết thành công bài toán phân định ranh giới địa chất phức tạp giữa hệ tầng phun trào Đơn Dương và pha 2 của phức hệ xâm nhập Định Quán, hiệu chỉnh đường ranh giới tiếp xúc với độ chính xác hình học cao hơn so với tài liệu bản đồ truyền thống.

Thảo luận kết quả

Sự phân hóa rõ rệt của các tham số phóng xạ tại vùng nghiên cứu bắt nguồn từ quá trình hoạt hóa magma kiến tạo mạnh mẽ trong 2 chu kỳ tạo sinh Indosini và Kimeri ở Đại Trung sinh. Đá granit và riolit, đaxit phún xuất mới tại đây có sự tập trung dị thường của Kali và Thori tại các đới biến đổi nhiệt dịch liên quan trực tiếp đến khoáng hóa vàng.

Dữ liệu phân tích trong luận văn được mô hình hóa trực quan qua hệ thống biểu đồ và bảng ma trận logic. Biểu đồ đường cong tỷ trọng thông tin tích phân Qh theo số lượng tham số h cho thấy điểm uốn cực đại xuất hiện chính xác tại h = 8, tạo cơ sở toán học khách quan để cắt bỏ 7 tham số có lượng tin thấp. Đồng thời, bảng ma trận logic chuẩn hóa các giá trị thực địa về dạng nhị phân 0 và 1 thông qua các ngưỡng quét đặc trưng, giúp thuật toán xử lý nhanh chóng mà không bị méo tín hiệu. Khi so sánh với phương pháp phân tích thành phần chính hoặc lọc trường đơn lẻ, phương pháp Khoảng cách – Tần suất – Nhận dạng cho thấy khả năng triệt tiêu nhiễu tốt hơn khoảng 25%, giúp thể hiện sắc nét các ranh giới thạch học trên bản đồ đường đồng mức hệ số đồng dạng.

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

Dựa trên kết quả nghiên cứu và kiểm nghiệm thực địa, luận văn đưa ra 4 nhóm giải pháp công nghệ mang tính khả thi cao nhằm nâng cao hiệu quả công tác địa vật lý:

  • Chuẩn hóa quy trình tiền xử lý số liệu bay đo phổ xạ: Liên đoàn Vật lý Địa chất và các trung tâm công nghệ địa chất cần ban hành hướng dẫn kỹ thuật thống nhất việc tính toán 11 tham số tỷ số trung gian kết hợp 4 kênh đo gốc. Giải pháp này hướng tới mục tiêu giảm thiểu 15% sai số đo đạc cục bộ trong khung thời gian 6 tháng triển khai.

  • Tích hợp thuật toán vào phần mềm phân tích chuyên dụng: Nhóm nghiên cứu thuộc các trường đại học và viện chuyên ngành cần lập trình đóng gói module Khoảng cách – Tần suất – Nhận dạng thành công cụ tự động hóa trên nền tảng mã nguồn mở hoặc GIS trong vòng 12 tháng. Mục tiêu kỹ thuật là rút ngắn thời gian xử lý dữ liệu từ vài tuần xuống còn 48 giờ làm việc, đồng thời tăng năng suất xử lý dữ liệu diện tích lên 40%.

  • Ứng dụng khoanh vùng chi tiết các diện tích quặng vàng triển vọng: Các doanh nghiệp khai thác và thăm dò khoáng sản cần ưu tiên bố trí mạng lưới khoan thăm dò chi tiết tại các đối tượng số 2, 3 và 1 thuộc xã Khánh Thượng, nơi có hệ số đồng dạng đạt trên 75,90%. Lộ trình thực hiện trong 18 tháng nhằm nâng cấp trữ lượng khoáng sản dự báo lên cấp thương mại.

  • Mở rộng phạm vi áp dụng cho các nguồn dữ liệu địa vật lý đa trường: Các cơ quan nghiên cứu địa chất cần tiến hành thử nghiệm phương pháp đối với tài liệu từ hàng không, trọng lực mặt đất và viễn thám đa phổ trong giai đoạn 2 đến 3 năm tới. Chỉ tiêu đặt ra là nâng cao độ chính xác phân định cấu trúc địa chất ẩn sâu lên trên 85%.

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

Nội dung và kết quả nghiên cứu của luận văn mang lại giá trị thực tiễn sâu sắc cho 4 nhóm đối tượng chuyên môn:

  • Học viên cao học và nghiên cứu sinh chuyên ngành Vật lý địa cầu, Địa chất học: Cung cấp tài liệu tham khảo chuẩn mực về phương pháp luận xử lý thống kê đa biến, kỹ thuật chuyển đổi ma trận nhị phân logic và cách tiếp cận bài toán nhận dạng có mẫu đối chứng trên 15 biến tham số.

  • Kỹ sư địa vật lý và chuyên viên xử lý dữ liệu bay đo: Sử dụng công trình như một cẩm nang kỹ thuật thực hành để trích xuất dữ liệu phổ xạ, xây dựng cửa sổ quét trượt và tối ưu hóa 8 tính chất nhận dạng trong các dự án bay đo tỷ lệ 1/50.000 và 1/25.000.

  • Cán bộ quản lý tại các cơ quan quản lý tài nguyên khoáng sản: Giúp các nhà hoạch định chính sách thuộc Tổng cục Địa chất và Khoáng sản có thêm cơ sở khoa học tin cậy để thẩm định các đề án phân vùng triển vọng tài nguyên, nâng cao hiệu quả phân bổ ngân sách điều tra cơ bản.

  • Doanh nghiệp thăm dò và khai thác khoáng sản: Ứng dụng trực tiếp quy trình nhận dạng để đánh giá nhanh tiềm năng các mỏ quặng vàng, thiếc và cao lanh tại vùng Nam Trung Bộ, từ đó tối ưu hóa vị trí các lỗ khoan kiểm chứng và tiết kiệm hàng tỷ đồng chi phí thực địa.

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

  • Phương pháp Khoảng cách – Tần suất – Nhận dạng khác biệt gì so với các phương pháp nhận dạng truyền thống? Phương pháp kết hợp ưu thế của khoảng cách khái quát Paguonov và tần suất Griffiths - Vinni, cho phép sử dụng đồng thời mẫu quặng và mẫu không quặng để chọn lọc khách quan các tham số nhạy cảm nhất, thay vì dựa vào kinh nghiệm chủ quan của người phân tích.

  • Vì sao tham số Thori lại có khoảng cách khái quát vượt trội đạt 526,8808? Trong bối cảnh địa chất Khánh Vĩnh với các thành tạo granit và đới biến đổi magma nhiệt dịch, Thori có tính trơ hóa học cao và phân bố tương phản rõ rệt giữa khối đá biến đổi chứa quặng và đá vây quanh không quặng, tạo nên độ lệch chuẩn cực lớn trong không gian dấu hiệu.

  • Luận văn xử lý dữ liệu trường phổ gamma ban đầu như thế nào trước khi nhận dạng? Dữ liệu 4 kênh gốc được nội suy lên lưới ô vuông đều, tính toán bổ sung 11 tham số tỷ số trung gian, sau đó đối chiếu với khoảng giá trị đặc trưng của mẫu chuẩn để mã hóa thành ma trận logic nhị phân với các phần tử mang giá trị 0 hoặc 1.

  • Chỉ số đồng dạng Pm* đạt bao nhiêu phần trăm thì đối tượng được coi là có triển vọng khoáng sản? Trong thực nghiệm tại vùng Khánh Thượng, các đối tượng có chỉ số Pm* từ 69,58% đến 82,30% đều thể hiện tính chất đồng dạng cao với mẫu chuẩn quặng loại B, trong đó mức ngưỡng trên 75% được đánh giá là diện tích có mức độ tin cậy ưu tiên hàng đầu.

  • Phương pháp này có thể áp dụng cho các loại khoáng sản khác ngoài quặng vàng không? Hoàn toàn có thể áp dụng rộng rãi cho nhiều loại hình khoáng sản phóng xạ, đa kim, thiếc, wolfram hoặc xác định đới biến đổi nhiệt dịch, miễn là đối tượng khảo sát có sự tương phản trường địa vật lý rõ rệt so với môi trường xung quanh.

Kết luận

  • Luận văn đã hoàn thiện cơ sở lý luận và quy trình ứng dụng phương pháp Khoảng cách – Tần suất – Nhận dạng mới trong xử lý số liệu địa vật lý máy bay.
  • Tự động hóa thành công khâu chọn lọc 8 tham số tối ưu từ bộ 15 biến trường xạ, với khoảng cách khái quát của Thori đạt 526,8808.
  • Nhận dạng chính xác 8 cụm dị thường triển vọng khoáng sản với chỉ số đồng dạng cao từ 69,58% đến 82,30% tại vùng Khánh Thượng, Khánh Hòa.
  • Phân định chuẩn xác ranh giới tiếp xúc địa chất giữa hệ tầng Đơn Dương và phức hệ Định Quán trên diện tích nghiên cứu 209,42 km².
  • Thiết lập tiền đề công nghệ quan trọng để mở rộng xử lý dữ liệu địa vật lý đa trường trong lộ trình 2 đến 3 năm tiếp theo.

Đóng góp lớn nhất của công trình là cung cấp một giải pháp định lượng tin cậy, giúp nâng cao hiệu quả khai thác nguồn tài liệu bay đo phong phú tại Việt Nam. Các đơn vị nghiên cứu và doanh nghiệp thăm dò hãy chủ động tích hợp phương pháp này vào quy trình xử lý thực tế để tối ưu hóa nguồn lực và gia tăng tỷ lệ phát hiện các mỏ khoáng sản giá trị.