Tổng quan nghiên cứu

Ứng dụng công nghệ viễn thám và hệ thống thông tin địa lý (GIS) trong điều tra địa chất hiện đại mang lại hiệu quả vượt bậc, giúp tiết kiệm từ 40% đến 60% chi phí khảo sát và rút ngắn 50% thời gian thực địa so với các phương pháp truyền thống. Tuy nhiên, tại các khu vực nhiệt đới gió mùa như Việt Nam, mức độ che phủ của thảm thực vật thường vượt quá 70% diện tích tự nhiên kết hợp với tầng vỏ phong hóa dày từ 5 đến 15 mét tạo ra rào cản kỹ thuật rất lớn trong việc giải đoán trực tiếp các thành tạo thạch học và cấu trúc kiến tạo.

Đề tài tập trung giải quyết bài toán xác định các thông số kỹ thuật cốt lõi trên tư liệu ảnh vệ tinh quang học đa phổ nhằm nhận dạng chính xác các nhóm đá chính và hệ thống đứt gãy kiến tạo. Mục tiêu cụ thể là xây dựng quy trình tích hợp giữa phân tích ảnh số tự động và giải đoán chuyên gia, phục vụ công tác thành lập sơ đồ địa chất cấu trúc ảnh tỷ lệ 1:50.000 tại khu vực huyện Cao Lộc, tỉnh Lạng Sơn.

Phạm vi không gian nghiên cứu bao trùm toàn bộ diện tích tự nhiên của huyện Cao Lộc trên khung dữ liệu ảnh vệ tinh Landsat 8 OLI/TIRS và ASTER giai đoạn 2014 – 2015. Kết quả nghiên cứu nâng độ phân giải không gian từ 30 mét lên 15 mét nhờ kỹ thuật trộn kênh toàn sắc, hỗ trợ khoanh định thành công 5 nhóm đá thạch học và xác lập mạng lưới đứt gãy đới phá hủy với độ trùng khớp thực tế trên 80%, đóng góp cơ sở dữ liệu quan trọng cho công tác cảnh báo tai biến trượt lở đất và quy hoạch khai thác 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 bức xạ điện từ và đặc trưng phổ phản xạ của vật chất địa chất. Theo các nguyên lý quang phổ kinh điển, mỗi loại khoáng vật có dải hấp thụ đặc trưng riêng biệt: ion Fe3+ trong hematit tạo vùng hấp thụ mạnh ở dải sóng khả kiến từ 0,4 đến 0,7 micromet; ion cacbonat (CO32-) trong đá vôi calcite tạo dải hấp thụ sâu ở khoảng 1,8 đến 2,4 micromet; trong khi các khoáng vật sét như kaolinite và montmorillonite thể hiện mức độ hấp thụ rõ nét ở bước sóng 1,4 và 1,9 micromet do dao động của nhóm hydroxyl (OH-).

Để giải quyết vấn đề thảm thực vật che phủ, mô hình phân giải pixel hỗn hợp tuyến tính (Linear Spectral Mixture Model - LSMM) được ứng dụng thông qua chỉ số thực vật chuẩn hóa NDVI (Normalized Difference Vegetation Index). Mô hình tách biệt giá trị phản xạ của 3 hợp phần cơ bản: đất, nước và thực vật trên từng pixel ảnh. Bên cạnh đó, khung nghiên cứu vận dụng 4 khái niệm kỹ thuật cốt lõi bao gồm: Tone ảnh (chia thành 10 đến 12 cấp độ xám), Cấu trúc ảnh (Texture biểu thị độ nhám, mịn hay dạng khảm), Cấu trúc dạng tuyến (Photolineament phản ánh khe nứt kiến tạo) và Phương pháp phân tích thành phần chính (PCA) nhằm giảm chiều dữ liệu quang phổ.

+-------------------------------------------------------------------------+
|                  QUY TRÌNH XỬ LÝ VIỄN THÁM ĐỊA CHẤT                     |
+-------------------------------------------------------------------------+
                                     |
                                     v
     +---------------------------------------------------------------+
     |              TIỀN XỬ LÝ DỮ LIỆU ĐA PHỔ (11 BANDS)             |
     | - Hiệu chỉnh bức xạ TOA & Áp suất khí quyển                  |
     | - Chuyển đổi tọa độ UTM WGS84 -> VN2000 (Bộ tham số chuẩn)   |
     | - Nâng độ phân giải lên 15m với Pan Band 8                    |
     +---------------------------------------------------------------+
                                     |
                                     v
     +---------------------------------------------------------------+
     |             PHÂN TÍCH QUANG PHỔ & GIẢI ĐOÁN ẢNH               |
     | - Phân tầng thảm phủ NDVI (Ngưỡng 0.2 - 0.6)                  |
     | - Tổ hợp màu tối ưu: Landsat 8 (PCA 624, RGB 572)             |
     | - Tách chiết Photolineament từ tỷ số kênh (6, 4/5, 3/4, 2)    |
     | - Phân loại có kiểm định SAM (Spectral Angle Mapper)          |
     +---------------------------------------------------------------+
                                     |
                                     v
     +---------------------------------------------------------------+
     |           KIỂM CHỨNG & THÀNH LẬP BẢN ĐỒ CẤU TRÚC ẢNH          |
     | - Đối sánh 15 điểm kiểm tra thực địa                          |
     | - Đối chiếu bản đồ địa chất tỷ lệ 1:50.000                    |
     | -> Hoàn thiện Sơ đồ địa chất cấu trúc ảnh tỷ lệ 1:50.000      |
     +---------------------------------------------------------------+

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

Nguồn dữ liệu đầu vào sử dụng cảnh ảnh Landsat 8 mã hiệu LC81260452014364LGN01 thu nhận ngày 30/12/2014 gồm 11 kênh phổ với dải quét 185 x 180 km và ảnh ASTER mã hiệu ASTB151022033507 chụp ngày 22/10/2015 phủ diện tích 60 x 60 km. Mẫu dữ liệu thực địa gồm 15 điểm lộ kiểm chứng tiêu biểu đại diện cho các phân vị địa tầng và đới dập vỡ kiến tạo trên địa bàn huyện Cao Lộc. Phương pháp chọn mẫu dựa trên mạng lưới phân bố thạch học từ các mảnh bản đồ địa chất tỷ lệ 1:50.000 hiện có kết hợp khảo sát tuyến địa mạo thực tế.

Lý do lựa chọn phương pháp phân tích số là khả năng tối ưu hóa giá trị phản xạ quang học đa băng tần. Quy trình phân tích gồm các bước chặt chẽ:

  • Tiền xử lý hiệu chỉnh bức xạ tầng trên khí quyển (TOA) và áp suất khí quyển;
  • Nắn chỉnh tọa độ từ UTM WGS84 sang VN2000 bằng bộ 7 tham số chuẩn (dx = 191.45032835);
  • Nâng cao độ phân giải không gian từ 30 mét lên 15 mét bằng thuật toán tích hợp kênh Panchromatic (Band 8);
  • Phân tích chỉ số thực vật NDVI;
  • Tổ hợp màu và phân tích thành phần chính PCA;
  • Phân loại có kiểm định bằng thuật toán góc phổ (Spectral Angle Mapper - SAM) và phân tích mật độ Photolineament. Toàn bộ quy trình được thực hiện trên các nền tảng phần mềm ENVI 5.3, MapInfo 11.0 và Global Mapper 19 trong giai đoạn nghiên cứu từ năm 2016 đến 2018.

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

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

Nghiên cứu đã xác lập các chỉ số định lượng cụ thể để nhận dạng chính xác các nhóm đá và cấu trúc kiến tạo tại Cao Lộc, Lạng Sơn:

  • Phân vùng che phủ thực vật (NDVI): Kết quả tính toán chỉ số NDVI trên ảnh Landsat 8 cho thấy khoảng 65% diện tích có chỉ số NDVI > 0,6 (thảm phủ rừng nhiệt đới dày đặc); 25% diện tích có chỉ số từ 0,2 đến 0,5 (cây bụi, nương rẫy); và khoảng 10% diện tích có chỉ số NDVI < 0,2 (đá lộ, đất trống, đô thị). Dải giá trị này định hướng phân lập phương pháp giải đoán trực tiếp cho 10% diện tích lộ và giải đoán gián tiếp qua địa mạo cho 90% diện tích còn lại.
  • Tối ưu hóa tổ hợp màu và phân tích thành phần chính (PCA): Tổ hợp màu PCA 624 và RGB 572 trên ảnh Landsat 8 cùng tổ hợp PCA 213 trên ảnh ASTER cho phép bóc tách rõ nét 5 nhóm đá chính. Độ tương phản màu sắc giữa đá vôi carbonat (màu lam sáng, cấu trúc dạng khảm) và đá phun trào mafic (màu xanh đậm, cấu trúc đốm đặc) đạt mức phân biệt thị giác trên 85%.
  • Chiết xuất mật độ đứt gãy photolineament: Phân tích tỷ số kênh 6, 4/5, 3/4, 2 trên Landsat 8 đã khoanh định chính xác hệ thống đứt gãy phương tây bắc - đông nam và á kinh tuyến. Các khu vực có mật độ photolineament cao trên 2,5 km/km2 trùng khớp hoàn toàn với các đới milonit và dăm kết kiến tạo trên bản đồ địa chất thực địa.
  • Nâng cao chất lượng ảnh giải đoán: Quá trình trộn kênh toàn sắc 15m giúp tăng độ sắc nét biên giới địa chất lên 100% so với ảnh đa phổ 30m nguyên bản, cho phép nhận dạng các thể trầm tích bở rời thung lũng sông suối với sai số hình học dưới 12 mét.
Nhóm thạch học Tổ hợp màu tối ưu (Landsat 8 / ASTER) Cấu trúc ảnh (Texture) Giá trị NDVI trung bình Dấu hiệu địa mạo đặc trưng
Đá carbonat (Đá vôi) PCA 624 (Lam sáng) / RGB 312 (Vàng nhạt) Dạng khảm, gồ ghề, dạng móc 0,45 - 0,65 Địa hình karst, vách đứng hiểm trở
Phun trào mafic & Tuf RGB 572 (Lục thẫm) / PCA 213 (Tím sẫm) Dạng đốm mịn, khối đặc 0,35 - 0,50 Sườn dốc thoải, đỉnh tròn
Phun trào axit - trung tính PCA 624 (Hồng nhạt) / PCA 213 (Cam sáng) Dạng khối thô, phân nhánh 0,50 - 0,70 Đỉnh nhọn, phân cắt sâu
Trầm tích lục nguyên RGB 572 (Vàng nâu) / RGB 312 (Xanh rêu) Phân lớp rõ nét, dạng dải 0,40 - 0,60 Sườn phân bậc, khe rãnh song song
Trầm tích bở rời (Aluvia) RGB 432 (Xám trắng) / RGB 543 (Hồng tươi) Mịn, đồng nhất 0,15 - 0,30 Bãi bồi, lòng sông, thung lũng bằng

Thảo luận kết quả

Khả năng phân biệt vượt trội của các kênh sóng ngắn hồng ngoại (SWIR - Band 6 và Band 7 trên Landsat 8) bắt nguồn từ sự hấp thụ năng lượng đặc thù của các liên kết phân tử Al-OH, Mg-OH và gốc CO32- trong khoáng vật tạo đá. Khi áp dụng thuật toán PCA 624, các biến sai quang phổ được cô lập vào các thành phần chính, loại bỏ hiện tượng dư thừa thông tin giữa các dải phổ tương quan cao.

So sánh với các nghiên cứu viễn thám tại các vùng hoang mạc khô cằn ở Ai Cập hay Iran (nơi tỷ lệ lộ đá gốc đạt hơn 90% và có thể phân loại khoáng vật trực tiếp bằng phương pháp SAM với độ chính xác trên 92%), nghiên cứu tại Cao Lộc đã chứng minh tính hiệu quả của phương pháp tiếp cận đa chỉ số thích ứng khí hậu nhiệt đới. Việc kết hợp chỉ số NDVI để loại trừ nền thực vật và phân tích cấu trúc địa mạo dòng chảy đã nâng độ tin cậy của sơ đồ địa chất cấu trúc ảnh lên tương đương với các phương pháp đo vẽ truyền thống.

Các dữ liệu giải đoán được trực quan hóa tối ưu thông qua sơ đồ vector hóa đới nứt nẻ kiến tạo, ma trận nhầm lẫn (Confusion Matrix) kiểm chứng thực địa và biểu đồ phân bố giá trị phổ phản xạ 11 kênh, tạo lập bộ công cụ kiểm tra tiếp biên bản đồ địa chất tỷ lệ 1:50.000 đạt độ tin cậy 90%.

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

  • Chuẩn hóa quy trình xử lý viễn thám địa chất: Cục Địa chất và Khoáng sản Việt Nam cần chủ trì ban hành bộ quy chuẩn kỹ thuật quốc gia về tiền xử lý dữ liệu viễn thám quang học (bắt buộc hiệu chỉnh bức xạ TOA và nắn chuyển hệ tọa độ VN2000) áp dụng cho 100% các đề án đo vẽ địa chất trong vòng 12 tháng tới.
  • Tích hợp công nghệ đa cảm biến (Multi-sensor Fusion): Giai đoạn 2024 – 2026, các cơ quan nghiên cứu địa chất cần mở rộng kết hợp dữ liệu quang học Landsat 8, Sentinel-2 với dữ liệu radar khẩu độ tổng hợp (SAR) như Sentinel-1 hoặc ALOS PALSAR. Công nghệ giao thoa radar sẽ xuyên qua tầng phủ thực vật, nâng tỷ lệ phát hiện đứt gãy ẩn sâu lên trên 90%.
  • Xây dựng thư viện phổ phản xạ mặt đất quốc gia: Bộ Tài nguyên và Môi trường cần phân bổ ngân sách trong 2 năm tới để đo đạc và số hóa phổ phản xạ thực địa của ít nhất 50 nhóm đất đá đặc thù tại các tỉnh miền núi phía Bắc, tạo cơ sở dữ liệu mẫu chuẩn cho thuật toán phân loại SAM.
  • Ứng dụng cảnh báo sớm tai biến địa chất: Ủy ban nhân dân và Sở Tài nguyên và Môi trường tỉnh Lạng Sơn cần ứng dụng ngay sơ đồ mật độ photolineament và các đới dập vỡ kiến tạo từ luận văn để khoanh vùng các điểm có nguy cơ trượt lở đất đá cao dọc tuyến Quốc lộ 1A và Quốc lộ 4B, duy trì chu kỳ cập nhật dữ liệu 6 tháng/lần trước mùa mưa bão.

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

  • Học viên cao học và nhà nghiên cứu chuyên ngành Địa chất, Viễn thám và GIS: Cung cấp khung phương pháp luận chi tiết từ tiền xử lý dữ liệu, trích xuất chỉ số NDVI đến kỹ thuật nâng cao độ phân giải không gian 15m bằng kênh Panchromatic phục vụ nghiên cứu học thuật.
  • Kỹ sư địa chất tại các Liên đoàn Bản đồ và Doanh nghiệp thăm dò khoáng sản: Ứng dụng trực tiếp các bộ thông số tổ hợp màu (PCA 624, RGB 572) để giải đoán nhanh ranh giới thạch học và kiểm tra, ghép nối tiếp biên các mảnh bản đồ địa chất tỷ lệ 1:50.000.
  • Cán bộ quản lý tại Sở Tài nguyên & Môi trường, Chi cục Phòng chống Thiên tai: Sử dụng sơ đồ cấu trúc ảnh và bản đồ mật độ đứt gãy để đánh giá mức độ bất ổn định sườn dốc, phục vụ quy hoạch không gian an toàn cho các khu dân cư vùng núi.
  • Chuyên gia tư vấn thiết kế công trình hạ tầng và giao thông: Tham khảo dữ liệu khoanh vùng các đới phá hủy dăm kết và milonit để khảo sát địa chất công trình, tối ưu hóa hướng tuyến và tính toán ổn định nền móng cho các dự án giao thông, thủy lợi.

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

Tại sao cần thực hiện kỹ thuật Pan-sharpening kết hợp kênh 8 trên ảnh Landsat 8?
Kỹ thuật này sử dụng kênh toàn sắc Band 8 có độ phân giải không gian 15 mét để nâng cao độ nét của 7 kênh đa phổ (vốn có độ phân giải 30 mét). Quá trình này tăng gấp đôi độ chi tiết hình học trên ảnh tổ hợp, giúp nhận diện rõ các thể địa chất quy mô nhỏ và ranh giới phân lớp trầm tích mà không làm sai lệch đặc tính phổ phản xạ gốc.

Chỉ số thực vật NDVI hỗ trợ thế nào trong giải đoán địa chất ở vùng nhiệt đới?
Chỉ số NDVI phân loại chính xác mức độ che phủ của thực vật theo các ngưỡng định lượng (từ dưới 0,2 đối với đất đá trơ trọi đến trên 0,6 đối với rừng rậm). Việc xác định tỷ lệ che phủ giúp chuyên gia linh hoạt chuyển đổi giữa phương pháp giải đoán phổ trực tiếp ở vùng lộ và giải đoán gián tiếp qua mạng lưới sông suối, địa mạo ở vùng rừng dày.

Ưu thế nổi bật của phương pháp phân loại góc phổ (Spectral Angle Mapper - SAM) là gì?
Thuật toán SAM so sánh độ tương đồng giữa phổ của từng pixel ảnh với phổ mẫu dựa trên góc vector n chiều thay vì cường độ bức xạ. Điều này giúp loại bỏ triệt để ảnh hưởng của góc chiếu sáng Mặt Trời và hiện tượng bóng đổ địa hình vùng núi hiểm trở, mang lại độ chính xác phân loại thạch học đạt trên 80%.

Tổ hợp màu nào hiệu quả nhất để nhận diện đá vôi carbonat và đứt gãy kiến tạo?
Đối với ảnh Landsat 8, tổ hợp màu thành phần chính PCA 624 và ảnh tỷ số kênh 6, 4/5, 3/4, 2 mang lại hiệu quả cao nhất. Đá vôi hiển thị màu lam sáng dạng khảm nổi bật, trong khi các đới phá hủy dập vỡ kiến tạo biểu hiện qua các dải dị thường màu sắc tuyến tính kéo dài với độ tương phản sắc nét.

Quy trình xác định thông số ảnh trong luận văn có thể chuyển giao cho vùng khác không?
Hoàn toàn có thể chuyển giao cho các địa phương có điều kiện tự nhiên tương đồng tại khu vực Đông Bắc và Tây Bắc Việt Nam. Quy trình tiền xử lý, tính toán chỉ số NDVI và phân tích mật độ photolineament có tính chuẩn hóa cao, chỉ cần hiệu chỉnh lại giá trị góc thiên đỉnh Mặt Trời và bổ sung mẫu kiểm chứng thực địa tại địa phương mới.

Kết luận

  • Luận văn đã chuẩn hóa thành công quy trình xử lý ảnh vệ tinh quang học Landsat 8 và ASTER, nâng độ phân giải không gian từ 30 mét lên 15 mét phục vụ giải đoán địa chất.
  • Thiết lập bộ thông số tổ hợp màu tối ưu (PCA 624, RGB 572, PCA 213) giúp phân lập chính xác 5 nhóm đá thạch học chính với độ tin cậy trên 85%.
  • Xây dựng thành công phương pháp trích xuất mật độ photolineament để xác định mạng lưới đứt gãy và đới dập vỡ kiến tạo, trùng khớp với bản đồ địa chất thực tế tỷ lệ 1:50.000.
  • Giải quyết triệt để rào cản che phủ thực vật nhiệt đới (>70%) thông qua việc tích hợp mô hình phân giải pixel hỗn hợp NDVI và giải đoán gián tiếp qua địa mạo.
  • Thành lập hoàn chỉnh Sơ đồ địa chất cấu trúc ảnh tỷ lệ 1:50.000 huyện Cao Lộc, tỉnh Lạng Sơn, cung cấp công cụ đắc lực cho công tác cảnh báo tai biến trượt lở.

Đóng góp lớn nhất của công trình là mở ra hướng tiếp cận viễn thám hiện đại, chi phí thấp nhưng đạt độ chính xác cao cho ngành địa chất Việt Nam. Kế hoạch tiếp theo trong giai đoạn 2024 – 2025 là tích hợp ảnh radar siêu cao tần và triển khai áp dụng quy trình này cho toàn bộ các tỉnh vùng núi phía Bắc. Các cơ quan quản lý và đơn vị nghiên cứu hãy chủ động ứng dụng ngay giải pháp công nghệ này vào các đề án điều tra tai biến địa chất và quy hoạch khoáng sản nhằm tối ưu hóa nguồn lực và bảo đảm an toàn cộng đồng.