Tổng quan nghiên cứu

Địa vật lý máy bay là một trong những phương pháp thăm dò hiện đại, đóng vai trò then chốt trong công tác điều tra cơ bản địa chất và đánh giá tiềm năng khoáng sản tại Việt Nam. Với khả năng thu thập dữ liệu nhanh chóng trên diện tích rộng lớn, các đề án bay đo từ - phổ gamma đã tích lũy khối lượng thông tin khổng lồ. Tuy nhiên, việc khai thác nguồn tài nguyên số liệu này đang đối mặt với rào cản lớn do tính chất phân bố ngẫu nhiên của các trường vật lý, ảnh hưởng của nhiễu địa hình và sự giới hạn về số lượng tham số đầu vào trong các thuật toán phân tích truyền thống.

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 tác giả Nguyễn Thị Thu Hiền, thực hiện dưới sự hướng dẫn khoa học của PGS. Võ Thanh Quỳnh tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội năm 2016, tập trung giải quyết bài toán cấp thiết trên. Mục tiêu nghiên cứu trọng tâm là xây dựng, hoàn thiện và áp dụng thử nghiệm phương pháp Khoảng cách – Tần suất – Nhận dạng mới nhằm nâng cao độ chính xác trong phân loại cấu trúc địa chất và khoanh vùng triển vọng quặng hóa.

Nghiên cứu được triển khai thực nghiệm trên nguồn số liệu bay đo phổ gamma tỷ lệ 1/50.000 thuộc khu vực xã Khánh Thượng, huyện Khánh Vĩnh, tỉnh Khánh Hòa. Địa bàn có diện tích tự nhiên 209,42 km2 trong tổng thể diện tích huyện Khánh Vĩnh là 1.165 km2 với mật độ dân cư khoảng 29 người/km2. Bằng việc xử lý tổ hợp 15 thông số phóng xạ sơ cấp và biến đổi trung gian, đề tài đã xác lập quy trình tối ưu hóa việc chọn lọc 8 đặc trưng tiêu biểu nhất, mang lại hệ số đồng dạng nhận dạng đạt từ 69,58% đến 82,30% cho các đối tượng kiểm chứng. Kết quả nghiên cứu khẳng định tiềm năng to lớn của phương pháp trong việc gia tăng 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

Khung lý thuyết của đề tài được xây dựng dựa trên sự kết hợp chặt chẽ giữa lý thuyết nhận dạng thống kê có đối tượng mẫu và mô hình xử lý tổ hợp số liệu địa vật lý hiện đại. Luận văn kế thừa và phát triển 2 nền tảng toán học - địa vật lý cốt lõi:

Thứ nhất là lý thuyết phân tích khoảng cách khái quát của Paguonov. Phương pháp này đánh giá mức độ chứa thông tin của từng dấu hiệu địa vật lý thông qua bình phương khoảng cách khái quát giữa trọng tâm các đám mây dữ liệu của hai lớp đối tượng mẫu đối nghịch trong không gian tham số đa chiều. Công thức tính toán dựa trên độ chênh lệch giá trị trung bình chuẩn hóa theo tổng phương sai của từng tính chất giữa tập mẫu quặng và tập mẫu không quặng.

Thứ hai là phương pháp phân tích tần suất nhận dạng theo thuật toán Griffiths - Vinni, được PGS.TS Võ Thanh Quỳnh phát triển từ năm 2007. Phương pháp sử dụng tần suất xuất hiện đồng thời của các cặp dấu hiệu để xác định lượng thông tin tương đối, từ đó thiết lập tỷ trọng thông tin tích lũy nhằm chọn lọc tổ hợp tham số có chất lượng cao nhất.

Các khái niệm cơ bản được chuẩn hóa trong nghiên cứu bao gồm: ma trận thông tin logic nhị phân, cửa sổ quét giá trị đặc trưng, chỉ số đồng dạng và tỷ trọng thông tin tương đối. Nhờ việc chuyển đổi các thang đo vật lý liên tục sang không gian logic có kiểm soát, thuật toán mới triệt tiêu được sự khác biệt về đơn vị đo lường giữa các trường thành phần.

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

Nguồn dữ liệu sử dụng trong luận văn được trích xuất từ Đề án Bay đo từ - phổ gamma tỷ lệ 1/50.000 và đo vẽ trọng lực vùng Phan Rang - Nha Trang do Liên đoàn Vật lý Địa chất hoàn thành năm 2003. Tập dữ liệu sơ cấp bao gồm 4 kênh trường xạ: hàm lượng Uran, hàm lượng Thori, hàm lượng Kali và cường độ bức xạ kênh tổng. Từ dữ liệu gốc, 11 tham số biến đổi trung gian không thứ nguyên được tính toán bổ sung, bao gồm các tỷ số hàm lượng, chỉ số F và các tích số chuẩn hóa trên kênh tổng, hình thành tập hợp 15 tính chất nghiên cứu đồng bộ.

Cỡ mẫu phân tích bao gồm 9 đối tượng dị thường phóng xạ đã được khoanh định trước đó. Trong đó, đối tượng số 6 mang tiềm năng quặng vàng loại B được lựa chọn làm đối tượng mẫu chuẩn, đối sánh trực tiếp với 1 diện tích mẫu không quặng và 8 đối tượng khảo sát còn lại. Lý do lựa chọn phương pháp Khoảng cách – Tần suất – Nhận dạng là nhằm khắc phục triệt để hiện tượng thừa tham số nhưng thiếu thông tin hữu ích trong các phương pháp phân loại trường cổ điển.

Quy trình nghiên cứu 4 bước được triển khai nghiêm ngặt:

  • Bước 1: Nội suy lưới điểm đồng đều và xác định khoảng giá trị đặc trưng của 15 tính chất trên đối tượng mẫu.
  • Bước 2: Thiết lập ma trận thông tin nhị phân chuẩn hóa cho các đối tượng mẫu và đối tượng đối sánh.
  • Bước 3: Đánh giá khoảng cách khái quát Paguonov, sắp xếp lượng tin giảm dần và xác định số lượng tính chất tối ưu dựa trên điểm uốn của đường cong tỷ trọng thông tin.
  • Bước 4: Tính toán hệ số đồng dạng tại từng nút lưới hoặc từng đối tượng để phân loại và khoanh định ranh giới địa chất.

Toàn bộ quá trình tính toán, lập trình thử nghiệm và kiểm chứng mô hình được thực hiện liên tục trong khoảng thời gian từ năm 2015 đến đầu năm 2016.

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

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

Phân tích định lượng trên tập số liệu thực tế vùng Khánh Thượng đã đem lại các phát hiện khoa học quan trọng:

Thứ nhất, việc đánh giá khoảng cách khái quát giữa đối tượng mẫu quặng và đối tượng không quặng trên 15 tính chất phóng xạ đã cho thấy sự phân hóa thông tin rất rõ rệt. Tính chất Thori thể hiện khoảng cách khái quát cao nhất đạt 526,88; tiếp theo là tỷ số Uran trên Thori đạt 4,46; trong khi nhiều tính chất khác có giá trị dưới 1,0. Dựa trên đạo hàm bậc hai của đường cong tỷ trọng thông tin tích lũy, nghiên cứu đã chọn lọc chính xác 8 tính chất có khả năng phân biệt mạnh nhất gồm: U/Th, K/Tổng, K, U/Tổng, Tổng, UK/Tổng, ThK/Tổng và Th, giúp giảm 46,67% số chiều không gian dữ liệu đầu vào mà vẫn bảo toàn lượng thông tin cốt lõi.

Thứ hai, kết quả tính toán hệ số đồng dạng cho 8 đối tượng đối sánh với mẫu chuẩn số 6 đạt mức độ tương đồng rất cao. Cụ thể: đối tượng 2 đạt 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 đồng trung bình trên toàn bộ 8 đối tượng đạt 76,08%, khẳng định tính đồng nhất cao về bản chất quặng hóa của cụm dị thường này.

Thứ ba, trong bài toán thử nghiệm phân định ranh giới địa chất giữa hệ tầng Đơn Dương và pha 2 phức hệ Định Quán, phương pháp quét cửa sổ dịch chuyển đã xây dựng thành công sơ đồ đường đồng mức hệ số đồng dạng. Ranh giới tiếp xúc địa chất được tái lập có độ trùng khớp trên 85% so với tài liệu thực địa, đồng thời phát hiện thêm các điểm uốn nếp kiến tạo cục bộ chưa từng được thể hiện trên bản đồ truyền thống.

Thảo luận kết quả

Nguyên nhân cốt lõi giúp phương pháp Khoảng cách – Tần suất – Nhận dạng đạt độ chính xác vượt trội nằm ở cơ chế kết hợp hai cấp lọc thông tin. Phân tích khoảng cách khái quát đóng vai trò loại bỏ các biến số gây nhiễu, trong khi phân tích tần suất logic hóa dữ liệu giúp bảo toàn cấu trúc phân bố xác suất của đối tượng mẫu. Khi so sánh với các kỹ thuật phân tích thành phần chính hoặc lọc trường đơn biến, phương pháp mới cung cấp chỉ số đồng dạng định lượng cụ thể, hạn chế tối đa yếu tố chủ quan của chuyên gia phân tích.

Về mặt trực quan hóa, toàn bộ hệ số đồng dạng tính toán được thể hiện hiệu quả dưới dạng bản đồ đường đồng mức đẳng trị kết hợp bảng ma trận thông tin nhị phân. Khi xếp chồng lớp bản đồ đẳng trị hệ số đồng dạng lên bản đồ địa chất tỷ lệ 1/50.000, các đới biến đổi nhiệt dịch liên quan đến khoáng hóa vàng và ranh giới thạch học xuất hiện vô cùng sắc nét, mở ra cách tiếp cận trực quan và chuẩn xác cho công tác biên tập bản đồ địa chất hàng không.

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

Dựa trên kết quả nghiên cứu thực nghiệm thành công, luận văn đưa ra 4 giải pháp và khuyến nghị hành động thiết thực:

  1. Chuẩn hóa quy trình tiền xử lý số liệu phổ gamma máy bay: Bộ Tài nguyên và Môi trường cùng Tổng cục Địa chất và Khoáng sản cần ban hành hướng dẫn kỹ thuật tích hợp phân tích 15 tham số phóng xạ và lọc không gian logic nhị phân cho 100% các đề án bay đo tỷ lệ 1/50.000 và 1/25.000 trong giai đoạn 2026 - 2030.
  2. Xây dựng phần mềm xử lý chuyên dụng tích hợp thuật toán mới: Các đơn vị nghiên cứu như Viện Địa chất, Viện Vật lý Địa cầu cần phối hợp phát triển module phần mềm tự động hóa tính toán chỉ số đồng dạng trên nền tảng GIS, đặt mục tiêu rút ngắn khoảng 40% thời gian xử lý và tổng hợp dữ liệu địa vật lý hàng không.
  3. Mở rộng địa bàn thử nghiệm tại các cấu trúc địa chất phức tạp: Giao Liên đoàn Vật lý Địa chất chủ trì áp dụng phương pháp tại các đới khoáng hóa ẩn sâu ở Tây Bắc và Trung Trung Bộ, hướng tới mục tiêu nâng tỷ lệ xác thực qua khoan kiểm chứng đạt trên 80% đối với các dị thường dự báo đến năm 2028.
  4. Nâng cao năng lực chuyên môn cho đội ngũ kỹ sư: Tổ chức các chương trình tập huấn chuyên sâu về thuật toán nhận dạng thống kê và xử lý số liệu tổ hợp cho khoảng 150 cán bộ kỹ thuật địa vật lý tại các viện nghiên cứu và trường đại học trong vòng 24 tháng tới.

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

Công trình nghiên cứu mang giá trị học thuật và ứng dụng thực tiễn cao, đặc biệt hữu ích cho 4 nhóm đối tượng sau:

  • Các nhà khoa học và giảng viên ngành Vật lý địa cầu, Địa chất thăm dò: Cung cấp tài liệu tham khảo giá trị về phương pháp luận xử lý dữ liệu ngẫu nhiên, mô hình hóa toán - địa vật lý và thuật toán phân loại trường hiện đại.
  • Kỹ sư địa vật lý và chuyên viên xử lý dữ liệu tại các Liên đoàn địa chất: Áp dụng trực tiếp quy trình 4 bước để xử lý tài liệu bay đo phổ gamma, nâng cao độ chính xác khi khoanh vùng quặng và phân chia các thành tạo địa chất.
  • Học viên cao học và nghiên cứu sinh chuyên ngành Khoa học Trái đất: Sử dụng làm tài liệu nghiên cứu chuyên đề về xử lý tổ hợp thông tin, kỹ thuật trích xuất đặc trưng và lập trình thuật toán nhận dạng mẫu.
  • Cán bộ quản lý quy hoạch khoáng sản tại cơ quan nhà nước: Nắm bắt cơ sở khoa học định lượng để thẩm định báo cáo kết quả bay đo, phục vụ công tác quy hoạch thăm dò và khai thác khoáng sản bền vững.

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

Phương pháp Khoảng cách – Tần suất – Nhận dạng có điểm gì ưu việt hơn các phương pháp thống kê truyền thống? Phương pháp kết hợp được ưu điểm của thuật toán khoảng cách khái quát để loại bỏ biến số nhiễu và thuật toán tần suất để định lượng hóa mức độ tương đồng thông qua chỉ số đơn nhất, giúp phân loại chính xác các tập dữ liệu đa biến có độ phân tán lớn.

Tại sao nghiên cứu lại chọn lọc 8 tính chất trong tổng số 15 tính chất ban đầu? Dựa trên nguyên lý cực trị đạo hàm bậc hai của đường cong tỷ trọng thông tin tích lũy, 8 tính chất này tập trung phần lớn lượng thông tin phân biệt giữa mẫu quặng và không quặng, giúp tối ưu hóa khối lượng tính toán mà không làm giảm độ chính xác nhận dạng.

Chỉ số đồng dạng đạt từ 69,58% đến 82,30% tại vùng Khánh Thượng phản ánh điều gì? Kết quả này khẳng định 8 đối tượng đối sánh có cấu trúc trường xạ tương đồng rất cao với đối tượng mẫu quặng số 6, chứng minh khu vực nghiên cứu có tiềm năng lớn về quặng vàng nhiệt dịch phù hợp với các kết luận địa chất thực địa.

Phương pháp này có thể áp dụng cho các dạng tài liệu địa vật lý khác không? Về mặt toán học, thuật toán hoàn toàn có thể mở rộng xử lý cho dữ liệu từ trường hàng không, trọng lực, địa hóa hoặc tích hợp đa phương pháp nhằm phân loại đối tượng địa chất trong không gian nhiều chiều.

Yêu cầu dữ liệu đầu vào để áp dụng thành công phương pháp này là gì? Cần có mạng lưới điểm đo nội suy đồng đều, tối thiểu 2 tập mẫu đối nghịch đã xác định rõ bản chất địa chất qua khoan thăm dò, cùng tập hợp các tham số biến đổi trung gian được chuẩn hóa đầy đủ.

Kết luận

  • Luận văn đã hệ thống hóa và phát triển thành công phương pháp Khoảng cách – Tần suất – Nhận dạng mới phục vụ phân tích tổ hợp số liệu địa vật lý máy bay.
  • Xác lập quy trình định lượng chọn lọc 8 tính chất đặc trưng nhất từ 15 thông số phóng xạ, giúp giảm thiểu đáng kể số chiều dữ liệu và loại bỏ nhiễu hiệu quả.
  • Đánh giá thử nghiệm thành công 9 đối tượng dị thường tại xã Khánh Thượng với hệ số đồng dạng đạt từ 69,58% đến 82,30%, khẳng định tiềm năng triển vọng quặng vàng loại B.
  • Chứng minh tính khả thi vượt trội của phương pháp trong việc phân định chính xác ranh giới tiếp xúc giữa hệ tầng Đơn Dương và phức hệ Định Quán.
  • Đóng góp luận cứ khoa học quan trọng cho công tác điều tra cơ bản địa chất và mở ra hướng ứng dụng tự động hóa phân tích dữ liệu địa vật lý hàng không tại Việt Nam giai đoạn 2026 - 2030. Các đơn vị chuyên môn và nhà nghiên cứu quan tâm có thể tiếp cận trích xuất dữ liệu, khai thác quy trình và nhân rộng mô hình thuật toán vào các đề án thăm dò thực địa.