Tổng quan nghiên cứu

Phương pháp phổ gamma hàng không là một trong những công nghệ địa vật lý máy bay nòng cốt được triển khai mạnh mẽ tại Việt Nam từ thập niên 1980 đến nay. Với khả năng khảo sát nhanh chóng trên diện tích rộng hàng ngàn kilomet vuông ở tỷ lệ 1:25.000 và 1:50.000, phương pháp này đóng vai trò then chốt trong công tác thành lập bản đồ địa chất, tìm kiếm trực tiếp quặng phóng xạ Uran cũng như phát hiện các diện tích triển vọng khoáng sản kim loại quý như vàng, thiếc, molipden, bauxit và sa khoáng ven biển.

Tuy nhiên, thực tế xử lý số liệu địa vật lý máy bay bộc lộ một thách thức kỹ thuật lớn: các bản đồ phân bố dị thường phổ gamma truyền thống chỉ thể hiện các điểm dị thường đơn lẻ nằm rải rác. Trong khi đó, các đối tượng quặng hóa và biến đổi nhiệt dịch trong vỏ Trái Đất luôn phân bố tập trung theo các đới, diện tích có cấu trúc địa chất đồng nhất tạo thành các cụm dị thường. Việc chỉ dựa vào các điểm dị thường đơn lẻ khiến công tác kiểm tra thực địa mặt đất gặp nhiều khó khăn, thiếu tính định lượng và phụ thuộc lớn vào đánh giá cảm quan.

Luận văn thạc sĩ khoa học chuyên ngành Địa Vật lý của tác giả Cao Văn Chỉnh, dưới sự hướng dẫn khoa học của Phó Giáo sư 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, đã tập trung giải quyết bài toán cấp thiết này. Mục tiêu trọng tâm của nghiên cứu là áp dụng đồng bộ phương pháp mã hóa - phân loại dị thường đơn và phương pháp đánh giá - phân loại cụm dị thường vào xử lý, minh giải tài liệu phổ gamma hàng không thực tế tại khu vực Tuy Hòa và vùng Đông Bắc tỉnh Đắk Lắk. Phạm vi nghiên cứu thuộc tọa độ 12°39'N - 13°05'N và 108°39'E - 108°54'E, bao phủ diện tích hơn 800 km² với cấu trúc địa chất phức tạp. Đóng góp của luận văn giúp nâng cao độ tin cậy trong khoanh định mục tiêu khoáng sản lên trên 85%, đồng thời tiết kiệm hơn 60% kinh phí và thời gian cho công tác kiểm tra thực địa chi tiết trên mặt đấ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 vững chắc của vật lý hạt nhân và địa hóa phóng xạ tự nhiên. Cơ sở vật lý của phương pháp dựa trên hiện tượng phân rã hạt nhân tự nhiên của ba nguyên tố chính trong vỏ Trái Đất: Uran (U-238), Thôri (Th-232) và đồng vị phóng xạ của Kali (K-40). Khi các hạt nhân này phân rã, bức xạ gamma năng lượng cao được phát ra với các mức năng lượng đặc trưng. Kỹ thuật phổ gamma hàng không ghi nhận tín hiệu thông qua phổ kế nhấp nháy đa kênh, chia thành một kênh tích phân tổng và ba cửa sổ năng lượng vi phân riêng biệt:

  • Cửa sổ Kali: năng lượng từ 1,35 MeV đến 1,55 MeV (đỉnh hấp thụ 1,46 MeV).
  • Cửa sổ Uran: năng lượng từ 1,55 MeV đến 1,85 MeV (đỉnh hấp thụ 1,76 MeV của đồng vị con Bismuth-214).
  • Cửa sổ Thôri: năng lượng từ 2,40 MeV đến 2,80 MeV (đỉnh hấp thụ 2,62 MeV của đồng vị con Thallium-208).

Về cơ sở địa chất - thạch học, sự phân bố hàm lượng các nguyên tố phóng xạ tuân theo quy luật tiến hóa magma của A. Vinogradov. Trong các đá magma nguyên sinh, hàm lượng các nguyên tố phóng xạ tỷ lệ thuận với hàm lượng SiO2:

  • Đá axit (granitoid): Hàm lượng trung bình của Uran là 3,50 × 10⁻⁴%, Thôri là 18,00 × 10⁻⁴%, Kali là 3,34%, tỷ số Th/U đạt 5,10.
  • Đá trung tính (diorit): Hàm lượng Uran là 1,80 × 10⁻⁴%, Thôri là 7,00 × 10⁻⁴%, Kali là 2,31%, tỷ số Th/U đạt 3,90.
  • Đá bazơ (gabro, bazan): Hàm lượng Uran giảm xuống 0,50 × 10⁻⁴%, Thôri là 5,00 × 10⁻⁴%, Kali là 0,88%, tỷ số Th/U đạt 3,70.
  • Đá siêu bazơ (peridotit): Hàm lượng Uran chỉ đạt 0,003 × 10⁻⁴%, Thôri đạt 0,005 × 10⁻⁴%, Kali đạt 0,03%, tỷ số Th/U đạt 3,70.

Trong các quá trình biến chất nhiệt dịch tạo quặng, mối quan hệ tỷ lệ cân bằng tự nhiên này bị phá vỡ, dẫn tới hiện tượng tích tụ hoặc suy giảm cục bộ của từng nguyên tố, tạo ra các dị thường xạ địa hóa có ý nghĩa tìm kiếm khoáng sản.

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

Nghiên cứu sử dụng toàn bộ tập dữ liệu bay đo từ - phổ gamma hàng không tỷ lệ 1:50.000 tại vùng Tuy Hòa, trải dài qua các huyện Ea Kar, M'Đrăk, Krông Năng (Đắk Lắk) và phần phía Tây tỉnh Phú Yên. Cỡ mẫu nghiên cứu bao gồm hàng chục tuyến bay đo với hơn 100 điểm dị thường đơn lẻ được phát hiện và nhóm thành các cụm dị thường theo phân vùng địa chất. Phương pháp chọn mẫu là chọn mẫu toàn bộ các điểm dị thường vượt ngưỡng phông xạ địa hóa khu vực để đảm bảo tính đại diện và loại bỏ hoàn toàn các sai số ngẫu nhiên.

Quy trình xử lý dữ liệu được chuẩn hóa qua 3 bước nghiêm ngặt:

  1. Xử lý số liệu sơ cấp: San phẳng đường đo (làm trơn) bằng phương pháp trung bình trượt 3 điểm trên cả 4 kênh năng lượng và kênh độ cao nhằm loại bỏ thăng giáng thống kê. Tiếp đó, hiệu chỉnh phông dư tự nhiên (bức xạ vũ trụ và radon khí quyển), hiệu chỉnh độ cao bay thực tế về độ cao chuẩn theo quy luật hàm mũ giảm dần exp(-μh) bằng phương pháp bình phương cực tiểu, và hiệu chỉnh tán xạ Compton giữa các kênh năng lượng vi phân.
  2. Tính chuyển tham số địa vật lý: Chuyển đổi tốc độ đếm xung (cps) sang cường độ phóng xạ gamma tổng (µR/h) và hàm lượng nồng độ định lượng của Uran (ppm), Thôri (ppm), Kali (%). Toàn bộ số liệu được liên kết tuyến kiểm tra đầu và cuối ngày bay theo mô hình tuyến tính để triệt tiêu sai số trôi dạt thiết bị.
  3. Thuật toán phân loại cụm dị thường: Đánh giá tập hợp các dị thường đơn trong cùng một cụm thông qua đường cong mật độ phân bố xác suất để tìm giá trị cực đại có tần suất xuất hiện lớn nhất (Mode). Đồng thời, tính toán bổ sung 3 hệ số tương quan tuyến tính Pearson giữa các cặp nguyên tố (RU-Th, RU-K, RK-Th) nhằm xác định mức độ đồng nhất hoặc phá vỡ cân bằng địa hóa của môi trường đất đá.

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 bay đo khu vực Tuy Hòa – Đông Bắc Đắk Lắk đã mang lại nhiều kết quả định lượng rõ ràng:

Thứ nhất, kết quả tính toán sai số đo đạc thực tế bằng chương trình tự động hóa QSS trên các tuyến bay lặp khẳng định độ chính xác cao của tập số liệu. Cụ thể, sai số bình phương trung bình thực tế của kênh tổng TC chỉ đạt 0,308 µR/h (thấp hơn giới hạn cho phép 0,35 µR/h); sai số kênh Kali đạt 0,192% (thấp hơn giới hạn 0,25%); sai số kênh Uran đạt 0,392 ppm (thấp hơn giới hạn 0,40 ppm); và sai số kênh Thôri đạt 0,910 ppm (thấp hơn giới hạn 1,00 ppm). Các sai số thực tế đều thấp hơn mức trần quy định từ 2,0% đến 23,2%, đảm bảo độ tin cậy tuyệt đối cho việc phân tích chuyên sâu.

Thứ hai, nghiên cứu đã chuẩn hóa mã hóa dị thường đơn thành chuỗi ký tự gồm 6 vị trí (2 chữ cái và 4 chữ số) dựa trên 6 chỉ tiêu: độ rộng nửa biên độ T(1/2), gia số cường độ ΔJ, cường độ bức xạ tương đối Ji, tỷ số ΔTh/ΔU, tỷ số ΔU/ΔK và chỉ số nhiều thành phần F = (U × K) / Th². Trên cơ sở này, hơn 100 dị thường đơn được phân loại chính xác vào 7 nhóm bản chất phóng xạ: nhóm Uran (U), nhóm Thôri (Th), nhóm Kali (K), nhóm Uran-Kali (U-K), nhóm Thôri-Kali (Th-K), nhóm Thôri-Uran (Th-U) và nhóm hỗn hợp (H-H).

Thứ ba, việc áp dụng phương pháp đánh giá cụm dị thường đã tổng hợp thành công các cụm dị thường tiêu biểu trong khu vực. Điển hình như cụm dị thường số 18 được mã hóa định lượng là CB.3213 kết hợp với hệ số tương quan RU/Th đạt 0,4197. Giá trị hệ số tương quan thấp này phản ánh rõ nét sự phá vỡ tương quan xạ địa hóa thông thường, chỉ ra đới biến đổi nhiệt dịch giàu Kali cục bộ trên nền đá magma xâm nhập.

Thảo luận kết quả

Các kết quả nghiên cứu sau khi được xử lý bằng chương trình QLK đã được thể hiện trực quan trên hệ thống bản đồ địa vật lý tỷ lệ 1:50.000, bao gồm bản đồ các tuyến bay thực tế, bản đồ cường độ tổng gamma với tiết diện đường đẳng trị 1 µR/h, bản đồ hàm lượng Kali tiết diện 0,5%, bản đồ hàm lượng Uran tiết diện 1 ppm, bản đồ hàm lượng Thôri tiết diện 3 ppm, bản đồ tỷ số Th/U và bản đồ chỉ số nhiều thành phần F.

Việc phân tích không gian các cụm dị thường kết hợp với bản đồ địa chất khu vực cho thấy sự tương quan chặt chẽ giữa các cụm dị thường phóng xạ với các thành tạo địa chất đặc trưng. Tại các diện tích phân bố đá biến chất cổ thuộc hệ tầng Khâm Đức (PR2-3 kđ) và các khối xâm nhập granitoid phức hệ Bến Giằng (P1 bg), phức hệ Đèo Cả (K đc), các cụm dị thường mang bản chất Uran-Kali và Thôri-Uran xuất hiện với mật độ cao, chỉ số F vượt ngưỡng 1,5. Đây là dấu hiệu xạ địa hóa điển hình của các đới biến đổi nhiệt dịch liên quan trực tiếp đến quặng hóa vàng, bạc, chì sa khoáng và các khoáng sản kim loại hiếm với quy mô dự báo sa khoáng đạt trên 300.000 tấn.

Ngược lại, tại các vùng phủ bazan trẻ thuộc hệ tầng Đại Nga (N2 đn) và trầm tích Đệ tứ, trường phóng xạ có tính chất phân bố đồng nhất, hệ số tương quan giữa các nguyên tố duy trì ở mức cao (R > 0,75), cho thấy môi trường địa chất ổn định không có biểu hiện khoáng hóa. So sánh với các phương pháp tách trường truyền thống trước đây, phương pháp đánh giá cụm dị thường đã loại bỏ các tín hiệu nhiễu cục bộ, làm nổi bật bản chất địa chất thực sự của các nguồn gây dị thường.

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

Dựa trên các kết quả xử lý và phân tích tài liệu thực tế, nghiên cứu đề xuất 4 nhóm giải pháp mang tính khả thi cao:

Thứ nhất, chuẩn hóa và tự động hóa quy trình xử lý dữ liệu: Trung tâm Địa vật lý máy bay và Tổng cục Địa chất và Khoáng sản Việt Nam cần tích hợp trọn bộ thuật toán mã hóa 6 chỉ tiêu và phân loại cụm dị thường vào các phần mềm xử lý chuyên dụng như QLK, COSCAD và ER Mapper. Việc tự động hóa hoàn toàn chu trình tính toán sẽ giúp giảm thiểu 40% thời gian xử lý dữ liệu thô trong các đề án bay đo giai đoạn 2024-2026.

Thứ hai, định hướng kiểm tra thực địa mặt đất trọng tâm: Các đơn vị điều tra địa chất cần ưu tiên bố trí nguồn lực khoan thăm dò, đo phổ gamma mặt đất và lấy mẫu hóa địa tại các cụm dị thường có bản chất Uran-Kali (U-K) và Thôri-Uran (Th-U) sở hữu hệ số tương quan R < 0,5. Thời gian thực hiện kiểm tra chi tiết nên tập trung trong 12 đến 24 tháng tới tại khu vực M'Đrăk (Đắk Lắk) và Sông Hinh (Phú Yên) để nhanh chóng xác định ranh giới thân quặng.

Thứ ba, thiết lập hệ thống cơ sở dữ liệu GIS xạ địa hóa quốc gia: Bộ Tài nguyên và Môi trường cần chỉ đạo xây dựng cơ sở dữ liệu địa vật lý số tích hợp đa tỷ lệ, số hóa 100% tài liệu bay đo phổ gamma từ năm 1980 đến nay trước năm 2028. Cơ sở dữ liệu này không chỉ phục vụ định hướng khai thác khoáng sản bền vững mà còn hỗ trợ đắc lực công tác quan trắc môi trường phóng xạ sinh thái trên toàn quốc.

Thứ tư, ứng dụng công nghệ trí tuệ nhân tạo và học máy: Các viện nghiên cứu và trường đại học chuyên ngành cần đẩy mạnh nghiên cứu kết hợp phương pháp phân tích cụm dị thường với các thuật toán nhận dạng mẫu tự điều chỉnh và mạng nơ-ron nhân tạo. Mục tiêu nâng độ chính xác dự báo diện tích triển vọng khoáng sản lên 90-95% trong giai đoạn 2025-2030.

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

Nội dung và phương pháp luận của luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cao cho nhiều nhóm đối tượng:

  • Các nhà nghiên cứu và nghiên cứu sinh chuyên ngành Địa vật lý, Địa chất khoáng sản: Tài liệu cung cấp cơ sở lý thuyết sâu sắc, các công thức toán học hiệu chỉnh bức xạ và thuật toán thống kê phân tích trường phóng xạ đa biến để phát triển các công trình nghiên cứu chuyên sâu.
  • Kỹ sư và chuyên viên địa vật lý tại các liên đoàn địa chất, công ty khảo sát tài nguyên: Cung cấp tài liệu hướng dẫn kỹ thuật chi tiết về quy trình đúc kết số liệu bay đo, phương pháp tính sai số bằng phần mềm QSS và quy chuẩn thành lập các bản đồ phổ gamma tỷ lệ 1:50.000.
  • Giảng viên và sinh viên các trường đại học khối Khoa học Tự nhiên, Mỏ - Địa chất: Là tài liệu tham khảo và học liệu giảng dạy giá trị cho các học phần Địa vật lý máy bay, Thăm dò phóng xạ, Xử lý số liệu địa vật lý và Địa hóa ứng dụng.
  • Doanh nghiệp khai thác khoáng sản và tổ chức quản lý môi trường: Giúp nắm bắt phương pháp đánh giá định lượng các đới biến đổi quặng hóa triển vọng, tối ưu hóa chi phí đầu tư thăm dò và đánh giá hiện trạng phông bức xạ môi trường tại khu vực miền Trung - Tây Nguyên.

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

1. Phương pháp phổ gamma hàng không ghi nhận những nguyên tố nào?

Phương pháp tập trung đo đạc bức xạ gamma phát ra từ ba nguyên tố phóng xạ tự nhiên gồm Uran (qua đỉnh phổ Bismuth-214 ở 1,76 MeV), Thôri (qua đỉnh phổ Thallium-208 ở 2,62 MeV) và Kali (qua đồng vị K-40 ở 1,46 MeV). Thông qua 3 cửa sổ năng lượng này, các nhà địa vật lý xác định chính xác hàm lượng nồng độ của từng nguyên tố trong đất đá.

2. Vì sao cần chuyển từ đánh giá dị thường đơn sang đánh giá cụm dị thường?

Bản đồ dị thường đơn chỉ thể hiện hàng ngàn điểm rời rạc khiến việc lựa chọn vị trí kiểm tra mặt đất gặp nhiều khó khăn và mang tính chủ quan. Đánh giá theo cụm giúp gom các dị thường có cùng nguồn gốc địa chất, xác định các tham số đặc trưng chung (Mode) và làm nổi bật bản chất của cả một đới biến đổi quặng hóa trên quy mô diện tích.

3. Chuỗi mã hóa dị thường phổ gamma gồm những thông số nào?

Chuỗi mã hóa gồm 6 ký tự (2 chữ cái và 4 chữ số) biểu thị cho 6 chỉ tiêu: độ rộng nửa biên độ T(1/2) (mức A, B, C), gia số cường độ ΔJ (mức A, B, C), cường độ tương đối Ji (mức 1 đến 6), tỷ số ΔTh/ΔU (mức 1, 2, 3), tỷ số ΔU/ΔK (mức 1, 2) và chỉ số nhiều thành phần F (mức 1, 2, 3).

4. Hệ số tương quan giữa các nguyên tố phóng xạ mang ý nghĩa gì trong thăm dò khoáng sản?

Trong điều kiện địa chất ổn định, hàm lượng U, Th, K có mối tương quan tuyến tính chặt chẽ (hệ số tương quan R tiến gần 1). Khi xảy ra quá trình nhiệt dịch tạo quặng làm biến đổi đất đá, tương quan này bị phá vỡ khiến R giảm mạnh về gần 0. Do đó, hệ số tương quan thấp là dấu hiệu then chốt để phát hiện các đới quặng hóa.

5. Sai số đo đạc trong xử lý phổ gamma hàng không được kiểm soát ra sao?

Sai số được tính toán tự động qua kết quả đo lặp trên các tuyến sản xuất và tuyến kiểm tra bằng công thức sai số bình phương trung bình. Trong nghiên cứu này, sai số thực tế kênh tổng TC đạt 0,308 µR/h (chuẩn 0,35 µR/h), kênh K đạt 0,192% (chuẩn 0,25%), kênh U đạt 0,392 ppm (chuẩn 0,40 ppm) và kênh Th đạt 0,910 ppm (chuẩn 1,00 ppm), hoàn toàn đáp ứng quy phạm kỹ thuật hiện hành.

Kết luận

  • Luận văn đã hệ thống hóa và làm sáng tỏ cơ sở khoa học - địa chất của phương pháp phổ gamma hàng không trong điều tra địa chất và tìm kiếm khoáng sản.
  • Xây dựng hoàn chỉnh quy trình xử lý, hiệu chỉnh số liệu bay đo và tính toán sai số đo đạc đạt độ chính xác cao hơn từ 2% đến 23% so với tiêu chuẩn ngành.
  • Chuẩn hóa phương pháp mã hóa 6 chỉ tiêu cho dị thường đơn và phát triển phương pháp đánh giá cụm dị thường kết hợp hệ số tương quan Pearson giữa các cặp nguyên tố phóng xạ.
  • Ứng dụng thành công vào thực tế vùng Tuy Hòa – Đông Bắc Đắk Lắk, khoanh định chính xác các cụm dị thường triển vọng khoáng sản vàng, sa khoáng và kim loại hiếm liên quan đến các phức hệ magma xâm nhập.
  • Để khai thác tối đa giá trị khoa học của đề tài, các cơ quan chuyên môn và doanh nghiệp địa chất cần sớm áp dụng quy trình đánh giá cụm dị thường vào các đề án bay đo thực tế, kết hợp kiểm tra chi tiết mặt đất nhằm phát hiện các mỏ khoáng sản mới phục vụ phát triển kinh tế đất nước.