Tổng quan nghiên cứu

Xử lý và giải đoán tài liệu địa vật lý máy bay đóng vai trò then chốt trong công tác điều tra địa chất cơ bản và tìm kiếm khoáng sản quy mô vùng. Nghiên cứu tập trung giải quyết bài toán tối ưu hóa nguồn dữ liệu đo bay từ - phổ gamma tỷ lệ 1/50.000 trên diện tích khoảng 600 km2 tại khu vực Ninh Thuận (bao gồm huyện Ninh Phước và huyện Thuận Nam, giới hạn trong tọa độ 11°21'34'' đến 11°33'00'' vĩ độ Bắc và 108°44'49'' đến 109°00'25'' kinh độ Đông). Vấn đề cốt lõi đặt ra là khối lượng thông tin thu thập từ các tuyến bay đo rất phong phú với 11 tham số trường xạ, nhưng việc đưa đồng thời tất cả các tham số vào phân tích thường gây ra hiện tượng nhòe thông tin, tăng sai số và làm giảm độ tin cậy của các vùng dự báo triển vọng.

Mục tiêu cụ thể của luận văn là nghiên cứu, hoàn thiện và áp dụng phương pháp Tần suất – Nhận dạng mới để lựa chọn tổ hợp thông tin đặc trưng có lượng tin cao, từ đó tự động khoanh định và chính xác hóa ranh giới các diện tích có triển vọng khoáng hóa thiếc, vonfram và sulfur loại B. Dữ liệu thực nghiệm được kế thừa từ Đề án Bay đo từ - phổ gamma tỷ lệ 1/50.000 và đo vẽ trọng lực tỷ lệ 1/100.000 vùng Phan Rang - Nha Trang, đối sánh trực tiếp với đối tượng mẫu mỏ Nhị Hà. Ý nghĩa thực tiễn của công trình thể hiện qua việc nâng cao độ chính xác nhận dạng dị thường lên trên 85%, giúp giảm thiểu từ 30% đến 40% chi phí khảo sát thực địa và định hướng chính xác các vị trí khoan kiểm chứng 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 lý thuyết nhận dạng mẫu kết hợp lý thuyết xác suất thống kê ứng dụng trong địa vật lý. Hai hệ thống lý thuyết chủ đạo được tích hợp gồm:

  • 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 tần suất xuất hiện đồng thời các đặc trưng của đối tượng nghiên cứu.
  • Lý thuyết khoảng cách khái quát Paguonov: Đánh giá độ phân tách giữa hai tập mẫu đối nghịch (quặng và không quặng) trong không gian dấu hiệu đa chiều.

Bên cạnh đó, nghiên cứu vận dụng các mô hình phân phối xác suất lý thuyết, đặc biệt là luật phân phối chuẩn Gauss đối với các trường phóng xạ và phân phối log-normal đối với hàm lượng nguyên tố. Các khái niệm toán học và địa vật lý then chốt bao gồm:

  • Đường cong biến phân: Đường cong mật độ phân bố thực nghiệm dùng để xác định khoảng giá trị đặc trưng của từng trường đo đạc.
  • Ma trận thông tin nhị phân: Bảng dữ liệu chuyển đổi các giá trị trường thực tế thành các phần tử logic 0 hoặc 1 thông qua cửa sổ quét giá trị đặc trưng.
  • Tỷ trọng thông tin tương đối: Đại lượng định lượng phản ánh mức độ đóng góp của từng thông số vào việc nhận diện đối tượng.
  • Chỉ số đồng dạng: Tiêu chí định lượng xác định mức độ tương đồng giữa diện tích khảo sát với mẫu chuẩn mỏ khoáng sản.

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

Nguồn dữ liệu sử dụng trong luận văn gồm tập hợp các bản đồ hàm lượng phóng xạ Uran, Thori, Kali, kênh tổng và trường tỉ số F được thu thập từ các tuyến bay đo đồng bộ. Mẫu chuẩn được thiết lập tại diện tích mỏ Nhị Hà (vùng triển vọng số 6) có tọa độ 11°27′33″ vĩ độ Bắc và 108°53′55″ kinh độ Đông, kết hợp tập mẫu kiểm chứng gồm 5 cụm dị thường đối sánh trên toàn vùng nghiên cứu 600 km2.

Phương pháp chọn mẫu dựa trên nguyên tắc địa chất - địa vật lý thực nghiệm: lựa chọn các diện tích đã có kết quả kiểm tra khoan mặt đất xác nhận sự tồn tại của khoáng hóa thiếc - vonfram để làm tập mẫu dương tính, và các diện tích đá vây quanh không biến đổi làm tập mẫu âm tính. Toàn bộ dữ liệu trường khảo sát được nội suy lên mạng lưới nút đều ô vuông bằng hệ thống phần mềm chuyên dụng COSCAD.

Lý do lựa chọn phương pháp Tần suất – Nhận dạng mới là nhờ khả năng giải quyết triệt để bài toán đa trị của trường địa vật lý và khắc phục giới hạn số lượng tham số đầu vào của các thuật toán truyền thống. Điểm đột phá của phương pháp là việc ứng dụng kỹ thuật cửa sổ quét di động đa hình dạng (tròn, chữ nhật, elip) trượt trên toàn bộ mạng lưới nút đều, cho phép nhận dạng và khoanh định ranh giới dị thường tự động ngay cả khi chưa biết trước vị trí các đối tượng đối sánh. Toàn bộ quy trình tính toán được thực thi tự động qua phần mềm QTS và QTSM trong giai đoạn nghiên cứu từ năm 2011 đến 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 phân tích thực nghiệm trên tài liệu địa vật lý máy bay vùng Ninh Thuận đã đem lại các kết quả định lượng cụ thể:

  • Sàng lọc và tối ưu hóa tổ hợp thông tin: Từ 11 tính chất xạ ban đầu (U, Th, K, Tg, F, U/Th, U/K, U/Tg, Th/Tg, K/Tg, Th/K), thuật toán Griffiths – Vinni đã xác định được 6 tính chất có lượng tin cao nhất đạt ngưỡng tỷ trọng tích lũy 80%, bao gồm K/Tg (tỷ trọng 0.49927), K, F, U/Th, Th/Tg và U/Tg. Việc loại bỏ 5 tính chất có lượng tin thấp giúp tinh gọn 45.4% dữ liệu đầu vào mà không làm mất đi các thông tin bản chất của trường quặng.
  • Đánh giá chỉ số đồng dạng định lượng: Phân tích đối sánh 5 đối tượng triển vọng trong vùng với mẫu chuẩn Nhị Hà đã xác định được sự phân hóa rõ rệt về mức độ tương đồng khoáng hóa. Chỉ số đồng dạng đạt giá trị phân hóa từ khoảng 65% đến trên 90%, trong đó đối tượng 1 và đối tượng 2 đạt mức tương đồng cao nhất với trên 88%, khẳng định tiềm năng chứa quặng thiếc - vonfram tương đương mỏ chuẩn.
  • Tự động khoanh định và hiệu chỉnh ranh giới mỏ: Ứng dụng kỹ thuật cửa sổ quét di động trên mạng lưới nút nội suy đã xác lập bản đồ đường đồng mức chỉ số đồng dạng, giúp tái khoanh định chính xác ranh giới vùng triển vọng Nhị Hà và các diện tích lân cận. Kết quả giúp thu hẹp 25% diện tích dị thường mờ nhạt so với phương pháp khoanh vẽ thủ công, tập trung tối đa vào đới biến đổi nhiệt dịch giàu khoáng hóa.

Thảo luận kết quả

Các kết quả nhận dạng phản ánh chính xác bản chất trường xạ địa chất của khu vực nghiên cứu. Sự biến thiên mạnh của các tỷ số K/Tg và U/Th tương ứng với các đới biến đổi nhiệt dịch greisen hóa và felspat hóa – môi trường sinh khoáng đặc trưng của thiếc và vonfram. Tương quan nghịch giữa các nguyên tố phóng xạ tại các cụm dị thường khẳng định sự xáo trộn địa hóa mạnh mẽ do hoạt động magma xâm nhập.

So sánh với phương pháp phân tích hồi quy đa thức truyền thống hoặc phương pháp nhận dạng không có mẫu chuẩn K-means, phương pháp Tần suất – Nhận dạng mới thể hiện tính ưu việt vượt trội:

  • Không bị phụ thuộc vào giả định tuyến tính của các trường số liệu.
  • Loại bỏ hoàn toàn sự can thiệp chủ quan của người giải đoán khi lựa chọn dấu hiệu nhận dạng.
  • Tránh được hiện tượng bão hòa thông tin khi xử lý các tập dữ liệu lớn.

Dữ liệu nghiên cứu được trực quan hóa tối ưu thông qua hệ thống bảng ma trận thông tin nhị phân, đồ thị biến phân xác định khoảng đặc trưng và các bản đồ đường đồng mức chỉ số đồng dạng. Cách trình bày này giúp các nhà địa chất dễ dàng nhận diện cấu trúc dị thường không gian 2D và đánh giá trực tiếp mức độ tin cậy của từng đới khoáng hóa.

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

Nhằm phát huy tối đa giá trị thực tiễn của phương pháp Tần suất – Nhận dạng mới trong công tác điều tra tài nguyên khoáng sản, các giải pháp và khuyến nghị cụ thể được đề xuất như sau:

  • Chuẩn hóa quy trình tiền xử lý số liệu địa vật lý máy bay: Ban hành quy chuẩn kỹ thuật về nội suy lưới đều và lọc nhiễu địa chất cho 100% dữ liệu đo bay từ - phổ gamma hiện có trên phạm vi toàn quốc. Hoạt động này cần được thực hiện bởi Liên đoàn Vật lý Địa chất trong giai đoạn 2026-2027 nhằm tạo nguồn dữ liệu đầu vào đồng nhất.
  • Tích hợp phần mềm QTS và QTSM vào hệ thống quản lý dữ liệu địa chất số quốc gia: Nâng cấp và chuyển giao gói phần mềm nhận dạng tự động cho Cục Địa chất và Khoáng sản Việt Nam, hướng tới mục tiêu rút ngắn 50% thời gian xử lý và tổng hợp tài liệu dự báo khoáng sản trong vòng 12 tháng tới.
  • Triển khai công tác kiểm tra thực địa và khoan thăm dò chi tiết: Ưu tiên bố trí các tuyến đo địa vật lý mặt đất và công trình khoan kiểm chứng tại các diện tích đạt chỉ số đồng dạng trên 85% tại khu vực Ninh Thuận. Các đoàn địa chất và doanh nghiệp khai khoáng cần đặt mục tiêu xác minh trữ lượng từ 2 đến 3 thân quặng tiềm năng trong khung thời gian 2027-2028.
  • Mở rộng ứng dụng tích hợp đa phương pháp địa vật lý: Kết hợp đồng thời tài liệu phổ gamma với tài liệu từ hàng không, trọng lực và ảnh viễn thám trong mô hình nhận dạng nhằm nâng cao độ chính xác dự báo không gian 3D lên trên 90%. Nhiệm vụ này do các viện nghiên cứu chuyên ngành và các trường đại học khối khoa học trái đất chủ trì thực hiện giai đoạn 2027-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ị thực tiễn và học thuật cho 4 nhóm đối tượng chính:

  • Cán bộ nghiên cứu và chuyên gia địa vật lý: Khai thác quy trình toán học trong việc tối ưu hóa lượng tin và thuật toán nhận dạng mẫu để phát triển các thuật toán xử lý tín hiệu địa vật lý tiên tiến.
  • Kỹ sư địa chất và cán bộ kỹ thuật tại các liên đoàn điều tra khoáng sản: Áp dụng trực tiếp quy trình xử lý của phần mềm QTS vào việc phân tích các tập bản đồ dị thường phóng xạ, rút ngắn thời gian khoanh vùng triển vọng trên các diện tích bay đo tỷ lệ 1/50.000 và 1/25.000.
  • Doanh nghiệp thăm dò và khai thác khoáng sản: Sử dụng các kết quả khoanh vùng định lượng tại Ninh Thuận để tối ưu hóa vị trí các lỗ khoan thăm dò, giảm thiểu rủi ro đầu tư và tiết kiệm từ 30% đến 40% chi phí khảo sát ban đầu.
  • Giảng viên, học viên cao học và sinh viên ngành Địa chất - Địa vật lý: Sử dụng công trình làm tài liệu tham khảo chuẩn mực về phương pháp kết hợp toán thống kê hiện đại trong giải quyết các bài toán khoa học trái đất thực tế.

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

Điểm khác biệt cốt lõi giữa phương pháp Tần suất – Nhận dạng mới và phương pháp truyền thống là gì?

Phương pháp mới tích hợp thuật toán cửa sổ quét di động trên mạng lưới nút đều nội suy, cho phép phát hiện và tự động khoanh định ranh giới các đối tượng đồng dạng trên toàn vùng ngay cả khi chưa biết trước vị trí của chúng, thay vì chỉ so sánh từng cặp đối tượng đơn lẻ như trước đây.

Tại sao lại chỉ lựa chọn 6 tham số trong tổng số 11 tham số xạ để tiến hành nhận dạng?

Theo định luật phân tích tần suất Griffiths – Vinni, 6 tham số được chọn (K/Tg, K, F, U/Th, Th/Tg, U/Tg) đã chiếm trên 80% tổng lượng tin đặc trưng của mỏ chuẩn. Việc loại bỏ 5 tham số dư thừa giúp triệt tiêu nhiễu địa chất và ngăn chặn hiện tượng làm mờ các dấu hiệu khoáng hóa quan trọng.

Ma trận thông tin nhị phân được xây dựng dựa trên nguyên lý nào?

Mỗi giá trị trường thực tế được đối sánh với khoảng đặc trưng xác định từ đường cong biến phân của mẫu chuẩn. Nếu giá trị khảo sát nằm trong khoảng đặc trưng của đối tượng quặng, phần tử ma trận sẽ nhận giá trị 1; nếu nằm ngoài khoảng đặc trưng, phần tử sẽ nhận giá trị 0.

Vùng quặng Nhị Hà đóng vai trò như thế nào trong nghiên cứu thực nghiệm?

Vùng mỏ Nhị Hà (triển vọng loại B về thiếc và vonfram) được sử dụng làm đối tượng mẫu chuẩn có độ tin cậy cao. Toàn bộ các khoảng giá trị đặc trưng của 11 thông số xạ và tỷ trọng thông tin được trích xuất từ vùng này để làm chuẩn đối sánh cho toàn bộ diện tích 600 km2.

Phương pháp này có thể áp dụng cho các loại hình khoáng sản khác không?

Phương pháp hoàn toàn có thể mở rộng cho mọi loại hình khoáng sản khác (vàng, đồng, sắt, đất hiếm) hoặc bài toán đo vẽ địa chất. Người dùng chỉ cần thay đổi tập mẫu chuẩn tương ứng với các trường địa vật lý đặc trưng (từ trường, điện trở suất hoặc trọng lực).

Kết luận

  • Luận văn đã hoàn thiện thành công cơ sở lý thuyết và thuật toán của phương pháp Tần suất – Nhận dạng mới, giải quyết triệt để bài toán khoanh định tự động dị thường khi chưa biết trước đối tượng đối sánh.
  • Phân tích và tinh chọn thành công tổ hợp 6 thông số xạ có lượng tin chiếm trên 80%, giúp giảm thiểu 45.4% chiều dữ liệu dư thừa và nâng cao chất lượng nhận dạng.
  • Ứng dụng thực nghiệm trên diện tích 600 km2 tại Ninh Thuận đã xác định chính xác các diện tích triển vọng khoáng sản thiếc - vonfram với chỉ số đồng dạng đạt trên 88% so với mỏ chuẩn Nhị Hà.
  • Tự động hóa toàn bộ quy trình tính toán thông qua phần mềm chuyên dụng QTS, hỗ trợ tái lập ranh giới mỏ chính xác hơn 25% so với phương pháp thủ công truyền thống.
  • Đóng góp quan trọng vào công tác điều tra cơ bản tài nguyên khoáng sản, mở ra hướng ứng dụng rộng rãi cho các đề án bay đo từ - phổ gamma trên toàn quốc trong giai đoạn 2026-2030.

Hãy áp dụng ngay phương pháp Tần suất – Nhận dạng mới vào công tác xử lý số liệu địa vật lý để tối ưu hóa chi phí thăm dò và nâng cao độ chính xác trong dự báo tài nguyên khoáng sản.