Tổng quan nghiên cứu

Trong hơn 25 năm qua, công tác bay đo địa vật lý tỷ lệ lớn (từ và phổ gamma hàng không) tại Việt Nam đã phát triển mạnh mẽ, thu thập khối lượng dữ liệu khổng lồ phục vụ công tác điều tra cơ bản địa chất và tìm kiếm khoáng sản. Tuy nhiên, việc xử lý và khai thác nguồn tài liệu phong phú này đang đối mặt với bài toán nghẽn thông tin: khối lượng dữ liệu đa chiều rất lớn nhưng số lượng tham số đầu vào của các phần mềm phân tích nhận dạng truyền thống thường bị giới hạn. Đồng thời, việc đưa toàn bộ các kênh dữ liệu vào nhận dạng không qua chọn lọc làm suy giảm độ chính xác khoảng 20% đến 30% do tác động của các tín hiệu nhiễu.

Trước thách thức đó, nghiên cứu được thực hiện nhằm xây dựng và áp dụng quy trình đánh giá, lựa chọn thông tin tối ưu trong xử lý số liệu địa vật lý máy bay. Mục tiêu cụ thể là xác định các tổ hợp thông tin đặc trưng có lượng tin cao nhất, giảm thiểu sai số nhận dạng và nâng cao độ tin cậy trong khoanh định diện tích triển vọng quặng.

Phạm vi không gian của đề tài bao phủ toàn bộ vùng Tây Nam Tuy Hòa với tổng diện tích 4.200 km2, thuộc địa bàn ba tỉnh Phú Yên, Đắk Lắk và Gia Lai - Kon Tum, dựa trên tập dữ liệu bay đo tỷ lệ 1:25.000 hoàn thành năm 1988 cùng các kiểm chứng thực địa năm 2012. Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc nâng tỷ lệ dự báo chính xác các điểm quặng tiềm năng lên trên 92%, giúp tiết kiệm ước tính 35% chi phí và thời gian khảo sát thực địa cho các đề án điều tra khoáng sản.

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 thông tin thống kê. Mô hình phân tích dựa trên giả định các trường vật lý do đối tượng địa chất tạo ra mang tính ngẫu nhiên do ảnh hưởng của môi trường và nhiễu địa hình. Khung lý thuyết tích hợp 2 thuật toán cốt lõi:

Thứ nhất, lý thuyết phân tích tần suất theo thuật toán Griffiths - Vinni xác định tỉ trọng thông tin tương đối của từng tính chất dựa trên ma trận tần suất xuất hiện đồng thời. Thuật toán thiết lập chỉ số thông tin tích lũy Pm, cho phép chọn lọc tập hợp m tính chất chủ đạo khi Pm đạt ngưỡng từ 80% trở lên.

Thứ hai, lý thuyết phân tích khoảng cách khái quát theo thuật toán Paguônôv đo lường mức độ phân tách giữa 2 lớp đối tượng đối nghịch (quặng và không quặng) 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. Kết hợp với việc tìm cực trị của đạo hàm bậc 2 đường cong thông tin tích lũy Qh theo h, thuật toán tách biệt chính xác tập thông tin tin cậy cao và loại bỏ các biến có sai số lớn hơn 15%.

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

Nguồn dữ liệu đầu vào bao gồm hệ thống bản đồ đo bay từ và phổ gamma tỷ lệ 1:25.000 với 6 kênh trường chính: Cường độ trường từ Delta T, Tổng cường độ Gamma, hàm lượng Kali (K), Uran (U), Thori (Th) và các tỷ số tương quan khoáng hóa.

Về cỡ mẫu và phương pháp chọn mẫu, nghiên cứu đã tiến hành thu thập và xử lý toàn bộ 120 cụm dị thường phổ gamma trên diện tích 4.200 km2. Phương pháp chọn mẫu chuẩn được tiến hành dựa trên 28 cụm dị thường đã có kết quả kiểm tra mặt đất và công tác khoan thăm dò chi tiết. Trong đó, cụm 68 (Xuân Sơn - Suối Cái với trữ lượng dự báo khoảng 5.000 kg vàng) và cụm 38 (triển vọng thiếc, wolfram) được chọn làm đối tượng mẫu quặng chuẩn; cụm 88 và cụm 55 được chọn làm mẫu đối nghịch không quặng.

Lý do lựa chọn phương pháp phân tích: Việc ứng dụng 2 phần mềm chuyên dụng tự lập trình là QTS (Tần suất - Nhận dạng) và QKC (Khoảng cách - Tần suất - Nhận dạng) giúp khắc phục triệt để nhược điểm xử lý thủ công, cho phép nạp số lượng tham số đầu vào tùy ý và tự động hóa quy trình phân lớp trong mốc thời gian hoàn thiện luận văn năm 2012.

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

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

Quá trình áp dụng hệ phương pháp đánh giá và lựa chọn thông tin vào khu vực Tây Nam Tuy Hòa đã mang lại 4 kết quả then chốt:

  1. Khẳng định tính ưu việt của phân tích theo cụm: Việc gom hàng nghìn điểm dị thường đơn lẻ thành 120 cụm dị thường quy chuẩn đã giúp giảm 75% độ phân tán dữ liệu và phản ánh chính xác bản chất trường phóng xạ của các cấu trúc địa chất.
  2. Lựa chọn tối ưu tập thông tin đặc trưng: Qua phân tích thuật toán Paguônôv và Griffiths - Vinni, tổ hợp 5 tham số có lượng tin cao nhất gồm hàm lượng K, tỷ số U/Th, Th/K, Delta T và hệ số tương quan R(U,Th) đã giữ lại hơn 85% tổng lượng thông tin hữu ích của toàn vùng.
  3. Hiệu quả vượt trội của chương trình QTS: Từ 2 cụm mẫu chuẩn, chương trình Tần suất - Nhận dạng QTS đã khoanh định được 23 cụm đồng dạng. Đối chiếu với 9 cụm đã kiểm tra mặt đất, có tới 8 cụm (đạt tỷ lệ chính xác 88,89%) thuộc nhóm triển vọng loại 1 và loại 2, chỉ có 1 cụm ít triển vọng.
  4. Độ tin cậy cao của chương trình QKC: Phương pháp Khoảng cách - Tần suất - Nhận dạng xác định được 21 cụm đồng dạng với quặng. Khi đối chiếu với 13 cụm đã kiểm tra thực địa, có 12 cụm (6 cụm loại 1 và 6 cụm loại 2, đạt độ tin cậy 92,31%) khẳng định có biểu hiện khoáng hóa vàng và kim loại hiếm.

Thảo luận kết quả

Kết quả phân tích cho thấy việc lọc thông tin trước khi nhận dạng đóng vai trò quyết định đến độ chính xác của mô hình địa vật lý. Trong các nghiên cứu truyền thống áp dụng hệ thống phân tích phổ COSCAD hay phần mềm ERMAPPER, việc nạp đồng thời toàn bộ biến trường thường dẫn đến hiện tượng trôi trọng tâm nhận dạng, khiến sai số phân loại lên tới 25% đến 30%. Ngược lại, phương pháp kết hợp QKC đã tối ưu hóa ranh giới phân tách lớp nhờ điểm uốn của đạo hàm bậc 2 trên đường cong tích lũy thông tin.

Về mặt trực quan hóa, toàn bộ dữ liệu được tổng hợp qua bảng ma trận đối sánh chỉ số tương đồng P*m và biểu đồ phân bố tần suất 6 kênh phổ. Bản đồ giải đoán cấu trúc và sơ đồ dự báo triển vọng khoáng sản tỷ lệ 1:50.000 thể hiện rõ các cụm dị thường triển vọng cao tập trung dọc theo các đới đứt gãy kiến tạo phương Tây Bắc - Đông Nam và ranh giới tiếp xúc của granitoid phức hệ Định Quán, Đèo Cả. Điều này hoàn toàn phù hợp với quy luật sinh khoáng nội sinh khu vực Nam Trung Bộ.

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

Dựa trên kết quả nghiên cứu thực nghiệm, luận văn đề xuất 4 giải pháp cụ thể nhằm hoàn thiện công tác xử lý số liệu địa vật lý:

  1. Tích hợp module đánh giá thông tin vào phần mềm chuyên dụng: Cục Địa chất và Khoáng sản Việt Nam cần chủ trì tích hợp thuật toán QTS và QKC vào các bộ công cụ phân tích hiện hành trước năm 2027, đặt mục tiêu tự động hóa 100% khâu tiền xử lý dữ liệu bay đo.
  2. Triển khai kiểm tra thực địa chi tiết cho 12 cụm dị thường ưu tiên: Liên đoàn Địa chất Trung Trung Bộ cần lập kế hoạch khảo sát thực địa giai đoạn 2026 - 2028 đối với 12 cụm dị thường loại 1 và loại 2 đã được khoanh định tại khu vực Tây Nam Tuy Hòa, hướng tới mục tiêu gia tăng 20% trữ lượng tài nguyên dự báo cấp 333 cho vàng và thiếc.
  3. Mở rộng thuật toán cho các dạng trường địa vật lý khác: Các viện nghiên cứu địa chất cần tiến hành thử nghiệm mô hình đối với số liệu từ trọng lực, điện trở và phân cực kích thích trong vòng 18 tháng tới, nhằm tăng độ chính xác phân giải mặt cắt sâu lên trên 90%.
  4. Chuẩn hóa cơ sở dữ liệu địa vật lý - địa hóa quốc gia: Bộ Tài nguyên và Môi trường cần ban hành quy chuẩn kỹ thuật số hóa dữ liệu bay đo tỷ lệ 1:25.000 và 1:50.000, hoàn thành chỉ tiêu 100% diện tích các vùng trọng điểm kinh tế khoáng sản trước năm 2030.

Đố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 chuyên môn sau:

  1. Kỹ sư địa vật lý và chuyên viên đo vẽ bản đồ: Nắm vững quy trình xử lý dữ liệu phổ gamma hàng không, ứng dụng trực tiếp thuật toán QKC để nâng cao 30% hiệu suất khoanh vùng dị thường.
  2. Cán bộ quản lý tại các cơ quan quản lý khoáng sản: Sử dụng kết quả dự báo 21 cụm dị thường triển vọng để định hướng cấp phép thăm dò, giảm 40% rủi ro đầu tư vốn ban đầu.
  3. Giảng viên và học viên cao học chuyên ngành Địa chất - Địa vật lý: Khai thác làm học liệu giảng dạy về lý thuyết nhận dạng mẫu, phương pháp thống kê đa biến với 2 thuật toán tiêu chuẩn Griffiths - Vinni và Paguônôv.
  4. Các doanh nghiệp khai thác và thăm dò khoáng sản: Áp dụng sơ đồ dự báo triển vọng tại vùng Tây Nam Tuy Hòa để tối ưu hóa vị trí đặt các lỗ khoan thăm dò sâu từ 100 đến 300 mét.

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

Phương pháp phân tích tần suất theo Griffiths - Vinni hoạt động như thế nào?

Thuật toán chuyển đổi ma trận dữ liệu thực địa thành các giá trị logic nhị phân 0 và 1 dựa trên khoảng giá trị đặc trưng. Sau đó, công thức tính toán tần suất xuất hiện đồng thời giữa các tính chất để xếp hạng tỉ trọng thông tin tương đối, giúp chọn lọc chính xác các biến đại diện chiếm trên 80% lượng tin tổng thể.

Tại sao cần chuyển từ xử lý dị thường đơn lẻ sang đánh giá theo cụm dị thường?

Trong thực tế mỗi đề án bay đo có tới hàng nghìn dị thường đơn lẻ phân bố dày đặc, gây khó khăn lớn cho công tác kiểm tra thực địa. Việc gom thành 120 cụm dị thường đồng nhất về bản chất phóng xạ giúp khái quát hóa trường vật lý, nâng độ chính xác nhận dạng đối tượng quặng lên hơn 90%.

Sự khác biệt cơ bản giữa 2 chương trình QTS và QKC là gì?

Chương trình QTS sử dụng thuần túy thuật toán tần suất trên mẫu chuẩn để tính chỉ số tương đồng P*m. Trong khi đó, chương trình QKC kết hợp cả khoảng cách khái quát Paguônôv giữa 2 tập mẫu đối nghịch (quặng và không quặng) cùng phân tích tần suất, cho phép phân tách ranh giới nhận dạng chặt chẽ hơn.

Luận văn đã kiểm chứng kết quả dự báo tại khu vực nào?

Nghiên cứu thử nghiệm trên diện tích 4.200 km2 thuộc vùng Tây Nam Tuy Hòa. Kết quả đã dự báo chính xác 21 cụm dị thường đồng dạng quặng, trong đó đối chứng thực địa 13 cụm đạt độ tin cậy thực tế 92,31% với các mỏ vàng Xuân Sơn và khoáng hóa thiếc.

Mô hình nhận dạng này có thể áp dụng cho các dữ liệu địa vật lý khác không?

Hoàn toàn có thể. Về mặt nguyên lý toán học, thuật toán không bị giới hạn số chiều tham số đầu vào nên có thể áp dụng hiệu quả cho dữ liệu trọng lực, từ trường, địa hóa và viễn thám trong các bài toán phân lớp địa chất phức tạp.

Kết luận

  • Luận văn đã giải quyết thành công bài toán nghẽn thông tin trong xử lý số liệu địa vật lý máy bay bằng việc ứng dụng 2 thuật toán chọn lọc thông tin tiên tiến Griffiths - Vinni và Paguônôv.
  • Xây dựng hoàn chỉnh 2 chương trình tính toán QTS và QKC, cho phép tự động hóa quy trình phân tích nhận dạng mẫu với số lượng biến đầu vào không giới hạn.
  • Đã phân loại thành công 120 cụm dị thường phổ gamma vùng Tây Nam Tuy Hòa (4.200 km2), khoanh định chính xác 21 cụm triển vọng khoáng sản với độ tin cậy thực tế đạt trên 92%.
  • Đóng góp luận cứ khoa học quan trọng cho việc định hướng công tác điều tra cơ bản địa chất và thăm dò khoáng sản vàng, thiếc giai đoạn 2026 - 2030.
  • Hãy tham khảo và ứng dụng ngay quy trình phân tích nhận dạng này để tối ưu hóa độ chính xác và tiết kiệm chi phí cho các đề án địa vật lý của đơn vị bạn!