Tổng quan luận án

Công trình khoa học "Tích hợp viễn thám và hệ thông tin địa lý nghiên cứu đánh giá tai biến trượt lở đất phục vụ giảm thiểu thiệt hại tại tỉnh Bắc Kạn" được thực hiện bởi nghiên cứu sinh Nguyễn Đình Tài dưới sự hướng dẫn khoa học của PGS. Nguyễn Ngọc Thạch. Luận án thuộc chuyên ngành Bản đồ, viễn thám và hệ thông tin địa lý (mã số: 62 44 02 14), bảo vệ học vị Tiến sĩ Địa lý tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội vào năm 2017. Đề tài được triển khai với sự hỗ trợ từ nhiệm vụ hợp tác quốc tế về khoa học và công nghệ theo Nghị định thư Việt Nam – Ấn Độ mang mã số 42/2009/HĐ-NĐT ("Tăng cường năng lực nghiên cứu, đào tạo về viễn thám và Hệ thông tin địa lý trong việc nghiên cứu, quản lý tai biến lũ lụt, lũ quét và trượt lở đất, nghiên cứu điển hình tại Vĩnh Phúc, Bắc Kạn").

Tính cấp thiết và khoảng trống nghiên cứu

Trượt lở đất là quá trình di chuyển đất đá trên sườn dốc dưới tác động trực tiếp của trọng lực. Đây là dạng tai biến diễn ra bất ngờ trong khoảng thời gian ngắn, có khả năng dịch chuyển hàng nghìn mét khối vật liệu sườn dốc, gây phá hủy cơ sở hạ tầng, vùi lấp khu dân cư và đe dọa trực tiếp tính mạng con người. Thống kê giai đoạn 2002–2012 trên phạm vi toàn cầu cho thấy trượt lở đất làm chết gần 9.000 người mỗi năm với thiệt hại kinh tế hàng tỷ USD tại các quốc gia như Ấn Độ, Trung Quốc, Pakistan, Philippines. Tại Việt Nam, giai đoạn 2000–2014 ghi nhận lũ quét và trượt lở đất khiến 646 người chết và mất tích, 351 người bị thương, làm hư hại và vùi lấp hơn 100.000 ngôi nhà cùng hơn 75.000 ha lúa, hoa màu.

Tại tỉnh Bắc Kạn, hiện tượng trượt lở đất xảy ra thường xuyên trong mùa mưa lũ. Các sự kiện điển hình được ghi nhận gồm vụ trượt lở ngày 04/7/2009 tại thôn Khên Lền (xã Công Bằng, huyện Pác Nặm) làm 13 người chết và mất tích; trận mưa lớn đầu tháng 8/2015 gây gần 400 điểm sạt lở với khối lượng đất đá phải thu dọn trên 30.000 m³; các sự cố vùi lấp đất canh tác và nhà dân do bãi thải khai thác khoáng sản tại huyện Chợ Đồn vào tháng 8/2016; sạt lở taluy dọc tuyến cao tốc Thái Nguyên – Bắc Kạn.

Khoảng trống nghiên cứu được tác giả xác định:

  • Các nghiên cứu trượt lở đất ở Việt Nam trước đây chủ yếu triển khai ở diện rộng, tỷ lệ nhỏ và mang tính dự báo định tính; thiếu các công trình điều tra chi tiết cấp tỉnh để hỗ trợ trực tiếp công tác phòng chống thiên tai.
  • Phương pháp thành lập bản đồ nguy cơ trượt lở chưa có sự thống nhất về thuật toán tích hợp giữa viễn thám (Remote Sensing - VT) và Hệ thông tin địa lý (GIS), đồng thời thường thiếu quy trình kiểm chứng độ chính xác định lượng của mô hình.
  • Ứng dụng ảnh vệ tinh đa thời gian, đa độ phân giải kết hợp phân tích thứ bậc và thống kê không gian chưa được khai thác đồng bộ trên địa bàn vùng núi Bắc Kạn.

Mục tiêu và nội dung nghiên cứu

Mục tiêu tổng quát của luận án là làm rõ vai trò, mức độ ảnh hưởng của các yếu tố tự nhiên và nhân sinh liên quan đến trượt lở đất; đồng thời xây dựng bản đồ nguy cơ trượt lở đất trên cơ sở tích hợp viễn thám và GIS nhằm phục vụ công tác giảm thiểu thiệt hại tại tỉnh Bắc Kạn.

Nội dung nghiên cứu gồm 5 nhiệm vụ chính:

  1. Tổng quan và xác lập cơ sở lý luận, phương pháp nghiên cứu về tích hợp viễn thám và GIS trong đánh giá nguy cơ trượt lở đất.
  2. Phân tích các yếu tố tự nhiên và nhân sinh ảnh hưởng đến quá trình phát sinh trượt lở đất trên địa bàn tỉnh Bắc Kạn.
  3. Ứng dụng kỹ thuật viễn thám và GIS để xác định, thành lập bản đồ hiện trạng trượt lở đất.
  4. Xử lý tích hợp không gian đa chỉ tiêu để thành lập bản đồ nguy cơ trượt lở đất tỉnh Bắc Kạn ở tỷ lệ 1:100.000.
  5. Đề xuất định hướng giải pháp giảm thiểu thiệt hại do tai biến trượt lở đất cho địa phương.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Các yếu tố tự nhiên (địa chất, địa hình, khí hậu thủy văn, lớp phủ thực vật) và các hoạt động khai thác lãnh thổ nhân sinh liên quan đến hiện trạng và nguy cơ phát sinh trượt lở đất.
  • Phạm vi khoa học: Đánh giá nguy cơ trượt lở đất bằng công nghệ viễn thám kết hợp hệ thông tin địa lý và các mô hình không gian định lượng.
  • Phạm vi không gian: Toàn bộ ranh giới hành chính tỉnh Bắc Kạn ở tỷ lệ nghiên cứu 1:100.000.
  • Giả thuyết nghiên cứu: Vận dụng nguyên tắc đồng dạng địa chất - địa mạo (quá khứ và hiện tại là chìa khóa cho tương lai), xem nguy cơ trượt lở đất là hàm số đa biến: $\text{TLĐ} = f(x_1, x_2, x_3, ..., x_n)$, trong đó $x_i$ là các biến số tự nhiên và nhân sinh.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trên thế giới

Tác giả tổng quan sự phát triển của các công trình quốc tế theo hai giai đoạn chính:

  1. Nghiên cứu hiện trạng bằng viễn thám:

    • Sử dụng ảnh hàng không và phân tích lập thể để xác định vết trượt (Mantovani và cộng sự, 1996), áp dụng phổ biến ở tỷ lệ lớn hơn 1:25.000.
    • Nhận diện vết trượt từ ảnh vệ tinh quang học: Petley và cộng sự (2002) dùng ảnh Landsat ETM+ ở vùng Himalaya; Marcelino và cộng sự (2009) sử dụng ảnh SPOT tăng cường độ phân giải lên 7 m để đạt tỷ lệ phát hiện 70%; Nichol và Wong (2005) dùng ảnh độ phân giải cao và siêu cao (SPOT-5 HRG, IKONOS, GeoEye).
    • Ứng dụng viễn thám chủ động: Singhroy và cộng sự (1998) tích hợp ảnh RADAR/SAR và InSAR; McKean và Roering (2004) dùng dữ liệu độ cao từ LiDAR (độ phân giải 1,5–10 m) để đo độ nhám bề mặt địa hình nhằm lập bản đồ khối trượt dưới thảm phủ thực vật dày.
    • Các kỹ thuật trắc địa hỗ trợ: quan trắc GNSS/CORS/RTK và giao thoa mặt đất IBIS (Tarchi và cộng sự).
  2. Đánh giá nguy cơ trượt lở đất bằng GIS:

    • Carrara (1983) khởi xướng việc phân tích dữ liệu raster trong GIS để đánh giá tính nhạy cảm trượt lở.
    • Soeters và Van Westen (1996) phân loại các phương pháp tiếp cận thành hai nhóm: trực tiếp (địa mạo, địa kỹ thuật - phù hợp tỷ lệ lớn) và gián tiếp (trọng số kinh nghiệm, thống kê - phù hợp tỷ lệ trung bình và nhỏ).
    • Nhóm trọng số kinh nghiệm nổi bật với phương pháp phân tích thứ bậc (Analytic Hierarchy Process - AHP) dựa trên ma trận so sánh cặp của Thomas Saaty (1980, 1986).
    • Nhóm thống kê định lượng: mô hình giá trị thông tin (Information Value Method) khởi xướng bởi Yin và Yan (1988), hoàn thiện bởi van Westen (1997); mô hình trọng số bằng chứng (Weight of Evidence); phân tích hồi quy Logistic; mạng thần kinh nhân tạo (Artificial Neural Network - ANN); cây quyết định và logic mờ (Fuzzy Logic).

Các hướng nghiên cứu tại Việt Nam và tại tỉnh Bắc Kạn

  • Nghiên cứu cấp công trình và tuyến giao thông: Điều tra nứt - trượt đồi Khau Cả và đồi Khí Tượng, thị xã Sơn La (Nguyễn Đức Mạnh, 1997); điều tra đoạn Đèo Gió - Nà Phặc dọc Quốc lộ 3 (Sở KH&CN Bắc Kạn, 2008; Uông Đình Khanh, 2007); điều tra sạt lở khu vực thị xã Bắc Kạn và Chợ Rã (Đỗ Minh Đức, 2005).
  • Nghiên cứu phân vùng quy mô vùng và cấp tỉnh: Các nghiên cứu tai biến địa chất dải ven biển miền Trung, dọc đường Hồ Chí Minh và khu vực Tây Bắc (Nguyễn Trọng Yêm, 2003; Lê Bá Thảo; Đỗ Tuyết); ứng dụng AHP đánh giá trượt lở tại Hòa Bình (Nguyễn Ngọc Thạch, 2004) và Thừa Thiên Huế (Trần Tuấn Anh, 2010); ứng dụng trọng số bằng chứng tại Bắc Giang (Nguyễn Viết Lương, 2012); ứng dụng SINMAP tại Quảng Bình (Nguyễn Kim Cương, 2010). Đề án điều tra phân vùng cảnh báo trượt lở đất đá miền núi của Viện Khoa học Địa chất và Khoáng sản (2013).

Vị trí và khoảng trống luận án lựa chọn giải quyết

Luận án kết hợp đồng thời phương pháp phân tích thống kê đơn biến (Giá trị thông tin - Information Value) để xác định thứ hạng cho từng cấp phụ trong mỗi lớp yếu tố, và phương pháp phân tích đa tiêu chí AHP của Saaty để xác định trọng số giữa 11 lớp yếu tố gây trượt. Hướng tiếp cận này khắc phục tính chủ quan của phương pháp chuyên gia thuần túy, phù hợp cho quy mô nghiên cứu cấp tỉnh ở tỷ lệ 1:100.000, đồng thời tích hợp quy trình kiểm chứng độ chính xác mô hình.


Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết và khái niệm sử dụng

  • Khái niệm trượt lở: Sử dụng định nghĩa của Varnes (1978), coi trượt lở đất là hiện tượng địa chất xảy ra trên sườn dốc khi đất, đá hoặc mảnh vụn di chuyển xuống phía dưới dưới tác động của trọng lực.
  • Phân loại trượt lở theo cơ chế chuyển động (Varnes, 1978): Rơi (fall), đổ/lật (topple), trượt xoay (rotational slide), trượt tịnh tiến (translational slide), sụt ngang (lateral spread), dòng chảy (flow) và trượt hỗn hợp (complex).
  • Phân loại theo quy mô kích thước khối trượt (van Schalkwyk, 1998):
    • Rất nhỏ: $< 10\text{ m}^2$
    • Nhỏ: $10 - 1.000\text{ m}^2$
    • Trung bình: $1.000 - 100.000\text{ m}^2$
    • Lớn: $100.000 - 1.000.000\text{ m}^2$
    • Rất lớn: $> 1.000.000\text{ m}^2$
  • Mô hình cân bằng vật lý: Sự trượt lở xảy ra khi tỷ số giữa lực kháng cắt chống trượt ($R$) và lực gây trượt ($T$) giảm xuống dưới ngưỡng an toàn (hệ số ổn định $F_s = R / T \le 1$).

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

  1. Khảo sát thực địa và điều tra xã hội học:
    • Thực hiện qua 3 đợt khảo sát (tháng 5/2010, tháng 6/2011, tháng 10/2012) cùng chuyên gia Ấn Độ.
    • Xây dựng 140 mẫu chìa khóa giải đoán; định vị điểm trượt bằng máy GPS cầm tay Garmin 78s; đo đạc góc dốc và hướng cắm của đá; phỏng vấn người dân bản địa về thời gian xuất hiện, quy mô và diễn biến khối trượt.
  2. Phương pháp viễn thám:
    • Thu thập và xử lý 16 cảnh ảnh vệ tinh quang học đa thời gian (2008–2015), đa độ phân giải gồm GeoEye-1 (độ phân giải 0,4–1,65 m), SPOT-4, SPOT-5 (độ phân giải 2,5–10 m).
    • Tiến hành giải đoán bằng mắt kết hợp mô hình số độ cao (DEM) và chìa khóa thực địa để thành lập bản đồ hiện trạng trượt lở.
    • Xử lý tự động tách chiết cấu trúc dạng tuyến (lineament) bằng các bộ lọc không gian (Sobel).
  3. Phương pháp GIS và phân tích không gian:
    • Xây dựng mô hình DEM từ bản đồ địa hình tỷ lệ 1:50.000.
    • Tính toán các chỉ tiêu trắc lượng hình thái: độ dốc, hướng dốc, độ chia cắt sâu (CCS), độ chia cắt ngang (CCN), vùng đệm (buffer) mạng lưới đường giao thông, nội suy phân bố lượng mưa trung bình năm.
  4. Phương pháp Giá trị thông tin ($W_i$):
    • Đánh giá mức độ ảnh hưởng của từng phân lớp thông tin $i$ dựa trên mật độ điểm ảnh trượt lở theo công thức của Yin và Yan (1988): $$W_i = \ln \left( \frac{\text{DenClass}}{\text{DensMap}} \right) = \ln \left( \frac{Npix(X_i) / Npix(N_i)}{\sum Npix(X_i) / \sum Npix(N_i)} \right)$$ Trong đó: $Npix(X_i)$ là số điểm ảnh trượt lở trong phân lớp $i$; $Npix(N_i)$ là tổng số điểm ảnh của phân lớp $i$; $\sum Npix(X_i)$ là tổng số điểm ảnh trượt lở trên toàn tỉnh; $\sum Npix(N_i)$ là tổng số điểm ảnh của toàn bộ vùng nghiên cứu.
  5. Phương pháp Phân tích thứ bậc (AHP):
    • Lập ma trận so sánh cặp giữa 11 lớp yếu tố gây trượt theo thang điểm 1–9 của Thomas Saaty.
    • Tính toán vector trọng số riêng ($\lambda_{max}$), chỉ số nhất quán $CI = (\lambda_{max} - n) / (n - 1)$ và tỷ số nhất quán $CR = CI / RI$. Giá trị $CR \le 0,10$ (10%) bảo đảm tính nhất quán của ma trận trọng số.
Kích thước ma trận ($n$) 1 2 3 4 5 6 7 8 9 10 11
Chỉ số ngẫu nhiên ($RI$) 0 0 0,58 0,90 1,12 1,24 1,32 1,41 1,45 1,49 1,51

Nguồn cơ sở dữ liệu

  • Bản đồ địa hình tỷ lệ 1:50.000 (nguồn: Cục Bản đồ).
  • Bản đồ địa chất và khoáng sản tỷ lệ 1:200.000.
  • Tư liệu viễn thám: 16 cảnh ảnh vệ tinh SPOT-4, SPOT-5, GeoEye-1 giai đoạn 2008–2015.
  • Chuỗi số liệu khí tượng thủy văn đo đạc tại các trạm khí tượng trên địa bàn tỉnh Bắc Kạn.
  • Cơ sở dữ liệu điều tra thực địa 140 điểm chìa khóa giải đoán và tọa độ các vết trượt.

Nội dung chính theo từng chương

Chương 1: Cơ sở lý luận và phương pháp nghiên cứu trượt lở đất tỉnh Bắc Kạn

Chương 1 hệ thống hóa cơ sở lý luận về cơ chế hình thành trượt lở đất và tổng quan tình hình nghiên cứu trong nước, quốc tế. Luận án tập trung phân tích mô hình vật lý của khối trượt sườn dốc và các yếu tố chi phối gồm địa chất, vỏ phong hóa, cấu trúc đứt gãy (lineament), độ dốc, lượng mưa và tác động nhân sinh.

Chương này xác lập quy trình nghiên cứu tích hợp viễn thám và GIS qua 5 giai đoạn:

  1. Xác định và số hóa 11 lớp biến số gây trượt.
  2. Thành lập bản đồ hiện trạng trượt lở từ ảnh vệ tinh đa thời gian.
  3. Tính toán thứ hạng các phân lớp bằng phương pháp giá trị thông tin ($W_i$) và trọng số lớp bằng phương pháp AHP.
  4. Chồng xếp không gian tính chỉ số nguy cơ trượt lở đất ($LSI$) và phân vùng nguy cơ theo phương pháp Natural Breaks, sau đó kiểm chứng độ chính xác.
  5. Đề xuất giải pháp quy hoạch và giảm thiểu thiệt hại.

Chương 2: Các yếu tố ảnh hưởng tới trượt lở đất tỉnh Bắc Kạn

Chương 2 phân tích đặc điểm tự nhiên và kinh tế - xã hội của tỉnh Bắc Kạn gắn với 11 biến số gây trượt lở:

  • Yếu tố địa chất: Phân tích đặc điểm thạch học các phân vị địa tầng, thành phần khoáng vật; đặc điểm cấu trúc vỏ phong hóa (chiều dày, thành phần hạt, độ dính sét); cấu trúc kiến tạo, mật độ đứt gãy và yếu tố dạng tuyến (lineament) được lọc bằng thuật toán Sobel; tính toán góc lệch giữa hướng cắm của mặt lớp đá với hướng dốc sườn địa hình.
  • Yếu tố địa hình, địa mạo: Phân tích trắc lượng hình thái từ DEM gồm độ cao tuyệt đối, độ dốc địa hình, độ chia cắt sâu (CCS) và độ chia cắt ngang (CCN), phân loại các bề mặt địa mạo (bóc mòn, xâm thực, tích tụ, karst).
  • Yếu tố khí hậu thủy văn và lớp phủ: Phân tích lượng mưa trung bình năm, tỷ lệ mưa theo mùa và mạng lưới địa chất thủy văn; hiện trạng lớp phủ thực vật năm 2014 phân loại theo cấu trúc tán rừng và hệ rễ.
  • Yếu tố nhân sinh: Đánh giá mật độ phân bố dân cư, mạng lưới hạ tầng giao thông (xây dựng các vùng đệm buffer quanh đường giao thông); hoạt động khai thác khoáng sản (vị trí các mỏ khai thác, bãi thải quặng chì kẽm, sắt, mangan); các hoạt động bạt núi, xẻ taluy mở đường làm mất chân sườn dốc tự nhiên.

Chương 3: Hiện trạng, nguy cơ trượt lở đất tỉnh Bắc Kạn và các giải pháp giảm thiểu

Chương 3 trình bày các kết quả thực nghiệm chính của luận án:

  • Hiện trạng trượt lở đất: Thành lập bản đồ hiện trạng trượt lở đất tỉnh Bắc Kạn từ giải đoán 16 ảnh vệ tinh đa thời gian, thống kê số lượng và quy mô diện tích các vết trượt phân bố theo từng đơn vị hành chính cấp huyện. Phân tích mối quan hệ giữa các điểm trượt lở đã xảy ra với 11 yếu tố thành phần.
  • Đánh giá nguy cơ trượt lở đất:
    • Tính toán giá trị $W_i$ cho từng phân lớp của 11 yếu tố gây trượt.
    • Xây dựng ma trận so sánh cặp AHP để xác định trọng số cho 11 lớp yếu tố, kiểm tra tỷ số nhất quán $CR < 0,10$.
    • Tính toán chỉ số nguy cơ trượt lở đất ($LSI$ - Landslide Susceptibility Index) bằng công thức nhân chập không gian: $$LSI = \sum_{j=1}^{11} \left( \omega_j \times W_{ij} \right)$$ Trong đó: $\omega_j$ là trọng số AHP của yếu tố $j$; $W_{ij}$ là giá trị thông tin của phân lớp $i$ trong yếu tố $j$.
    • Ứng dụng phương pháp Natural Breaks để chia giá trị $LSI$ thành 5 cấp nguy cơ: rất thấp, thấp, trung bình, cao và rất cao. Thành lập bản đồ nguy cơ trượt lở đất tỉnh Bắc Kạn tỷ lệ 1:100.000.
    • Kiểm chứng độ chính xác của mô hình dựa trên tần suất xuất hiện các vết trượt thực tế trong từng phân vùng nguy cơ.
  • Giải pháp giảm thiểu: Đề xuất định hướng quy hoạch sử dụng đất theo các vùng nguy cơ, lập danh mục các xã có nguy cơ cao và rất cao, đưa ra các biện pháp công trình (kè đá, rãnh thoát nước mặt, gia cố taluy) và phi công trình (cảnh báo sớm, quản lý bãi thải mỏ, bảo vệ lớp phủ rừng phòng hộ).

Kết quả và những đóng góp mới

Luận điểm bảo vệ và đóng góp mới của luận án

  1. Luận điểm 1: Tai biến trượt lở đất tại tỉnh Bắc Kạn chịu sự chi phối tổng hợp của các yếu tố tự nhiên và nhân sinh, trong đó các hoạt động nhân sinh (đặc biệt là việc bạt núi xẻ taluy mở đường giao thông và hoạt động đổ bãi thải khai thác khoáng sản) đóng vai trò ưu thế trong việc kích hoạt làm mất ổn định sườn dốc.
  2. Luận điểm 2: Việc tích hợp công nghệ viễn thám và GIS kết hợp phương pháp giá trị thông tin ($W_i$) và phương pháp phân tích thứ bậc (AHP) cho phép xây dựng bản đồ nguy cơ trượt lở đất tỷ lệ 1:100.000 đạt độ chính xác định lượng cao, có thể kiểm chứng bằng thực địa và số liệu thống kê.

Ý nghĩa khoa học và thực tiễn

  • Ý nghĩa khoa học: Xác lập quy trình công nghệ chuẩn hóa trong việc kết hợp giải đoán ảnh vệ tinh đa thời gian với các công cụ phân tích không gian GIS; bổ sung luận cứ khoa học về mối quan hệ giữa góc lệch hướng cắm của đá - hướng dốc sườn và mật độ lineament đối với cơ chế trượt lở đất ở vùng núi phía Bắc.
  • Ý nghĩa thực tiễn: Cung cấp bộ cơ sở dữ liệu số và bản đồ hiện trạng, bản đồ nguy cơ trượt lở đất tỷ lệ 1:100.000 phục vụ trực tiếp cho UBND tỉnh Bắc Kạn, Sở Tài nguyên và Môi trường, Sở Giao thông Vận tải trong công tác quy hoạch mạng lưới giao thông, cấp phép khai thác khoáng sản và bố trí các điểm dân cư an toàn.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của đề tài

  • Việc phát hiện các vết trượt lở có kích thước nhỏ gặp hạn chế do độ phân giải của một số cảnh ảnh quang học và hiện tượng thực vật che phủ nhanh sau mùa mưa, đòi hỏi phải kết hợp giải đoán bằng mắt và kiểm tra thực địa.
  • Mô hình đánh giá nguy cơ là mô hình tĩnh theo không gian, yếu tố thời gian và cường độ mưa theo từng trận mưa cực đoan chưa được mô phỏng theo thời gian thực.
  • Số lượng trạm đo mưa trên địa bàn tỉnh còn thưa, giá trị lượng mưa phân bố chủ yếu dựa vào phương pháp nội suy không gian.

Hướng nghiên cứu tiếp

  • Ứng dụng tư liệu viễn thám radar giao thoa (InSAR) và dữ liệu quét laser mặt đất/LiDAR độ phân giải cao để giám sát các biến dạng sườn dốc liên tục theo thời gian.
  • Tích hợp số liệu dự báo mưa định lượng theo thời gian thực vào GIS để nâng cấp hệ thống thành mô hình cảnh báo sớm tai biến trượt lở đất ở cấp độ thời gian thực.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên môn cho các nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và giảng viên: Chuyên ngành Bản đồ, Viễn thám, Hệ thông tin địa lý, Địa lý tự nhiên, Địa chất công trình và Quản lý tài nguyên môi trường có thể tham khảo quy trình tích hợp mô hình AHP và Giá trị thông tin ($W_i$) trong phân tích không gian đa chỉ tiêu.
  • Cơ quan quản lý nhà nước tại địa phương: Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn, Sở Xây dựng, Sở Giao thông Vận tải tỉnh Bắc Kạn có thể sử dụng trực tiếp bản đồ nguy cơ trượt lở đất và danh mục các xã trọng điểm nguy cơ cao để lập kế hoạch di dời dân cư và thiết kế công trình gia cố sườn dốc.
  • Kỹ sư quy hoạch hạ tầng miền núi: Tham khảo các phân tích về đới đệm an toàn quanh các trục đường giao thông (buffer) và sự phân bố bãi thải khoáng sản để giảm thiểu tác động làm mất chân sườn dốc.

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

1. Luận án sử dụng những loại tư liệu viễn thám và dữ liệu không gian nào?

Luận án sử dụng 16 cảnh ảnh vệ tinh đa thời gian (2008–2015) gồm GeoEye-1 (độ phân giải 0,4–1,65 m), SPOT-4 và SPOT-5 (độ phân giải 2,5–10 m); bản đồ địa hình tỷ lệ 1:50.000 của Cục Bản đồ để xây dựng mô hình số độ cao (DEM); bản đồ địa chất tỷ lệ 1:200.000; số liệu quan trắc mưa tại các trạm khí tượng và 140 mẫu chìa khóa giải đoán thực địa được định vị bằng máy GPS Garmin 78s.

2. Mô hình tích hợp không gian trong luận án bao gồm 11 yếu tố gây trượt nào?

11 lớp thông tin được lựa chọn làm biến số gồm:

  1. Thạch học
  2. Vỏ phong hóa
  3. Góc lệch giữa phương vị hướng cắm của đá với hướng đổ của địa hình
  4. Mật độ lineament (cấu trúc dạng tuyến)
  5. Đặc điểm địa mạo
  6. Độ dốc địa hình
  7. Độ chia cắt sâu (CCS)
  8. Độ chia cắt ngang (CCN)
  9. Lượng mưa trung bình năm
  10. Hiện trạng lớp phủ thực vật
  11. Vùng đệm khoảng cách đến đường giao thông (buffer giao thông)

3. Sự kết hợp giữa phương pháp Giá trị thông tin và AHP được thực hiện như thế nào?

Phương pháp Giá trị thông tin ($W_i$) được sử dụng để đánh giá phân dị theo "chiều ngang" bằng cách tính toán khách quan thứ hạng của từng phân lớp phụ dựa trên tỷ lệ mật độ trượt lở thực tế so với mật độ trung bình toàn tỉnh. Sau đó, phương pháp AHP của Thomas Saaty được sử dụng để đánh giá phân dị theo "chiều dọc" thông qua ma trận so sánh cặp nhằm gán trọng số ($\omega_j$) cho 11 lớp yếu tố lớn, bảo đảm tỷ số nhất quán $CR < 0,10$. Chỉ số nguy cơ trượt lở tổng hợp ($LSI$) là tổng tích phân giữa $\omega_j$ và $W_{ij}$.

4. Luận án phân chia nguy cơ trượt lở đất tại Bắc Kạn thành các cấp độ nào?

Dựa trên giá trị chỉ số $LSI$ thu được từ phép chồng xếp không gian, luận án sử dụng thuật toán phân nhóm Natural Breaks để chia toàn bộ diện tích tỉnh Bắc Kạn thành 5 cấp nguy cơ trượt lở đất:

  • Nguy cơ rất thấp
  • Nguy cơ thấp
  • Nguy cơ trung bình
  • Nguy cơ cao
  • Nguy cơ rất cao

Kết luận

Luận án Tiến sĩ Địa lý của tác giả Nguyễn Đình Tài đã xây dựng thành công quy trình công nghệ tích hợp viễn thám và GIS để đánh giá hiện trạng và nguy cơ trượt lở đất tại tỉnh Bắc Kạn ở tỷ lệ 1:100.000. Bằng việc kết hợp phương pháp thống kê Giá trị thông tin và phương pháp phân tích thứ bậc AHP trên 11 biến số tự nhiên - nhân sinh, công trình đã định lượng hóa mức độ nhạy cảm trượt lở đất và kiểm chứng mô hình với độ chính xác cao. Kết quả nghiên cứu khẳng định vai trò kích hoạt nổi bật của các hoạt động nhân sinh đối với tai biến sườn dốc trong giai đoạn hiện nay, đồng thời cung cấp cơ sở khoa học và hệ thống bản đồ phục vụ công tác quy hoạch, cảnh báo và giảm thiểu thiệt hại thiên tai tại địa phương.