Tổng quan nghiên cứu

Công trình thủy điện Sơn La là dự án trọng điểm cấp quốc gia với công suất lắp máy 2.400 MW, đóng góp khoảng 10 tỷ kWh điện năng mỗi năm vào mạng lưới điện quốc gia, đồng thời giữ vai trò điều tiết lũ và cung cấp nước trọng yếu cho vùng đồng bằng Bắc Bộ. Tuy nhiên, lưu vực hồ chứa có tổng diện tích tự nhiên lên tới 11.135 km², trải dài trên địa bàn 4 tỉnh Tây Bắc gồm Sơn La, Lai Châu, Điện Biên và Lào Cai (tọa độ từ 21°15'15'' đến 22°45'10'' vĩ độ Bắc và từ 102°50'10'' đến 104°35'15'' kinh độ Đông). Khu vực này sở hữu cấu trúc địa chất phức tạp, địa hình núi cao bị chia cắt sâu và chịu tác động mạnh mẽ của hoạt động tân kiến tạo.

Hàng năm, các hiện tượng tai biến địa chất tại Việt Nam gây thiệt hại ước tính hơn 4.000 tỷ đồng và cướp đi sinh mạng của khoảng 750 người. Khi hồ thủy điện Sơn La tích nước, mực nước dâng cao làm ngập chìm các sườn dốc bờ hồ, làm thay đổi áp lực thủy tĩnh và phá vỡ trạng thái cân bằng ứng suất của đất đá, khiến nguy cơ trượt lở đất gia tăng đột biến. Trước thực trạng đó, nghiên cứu đã được triển khai nhằm làm sáng tỏ hiện trạng phân bố và khoanh vùng cảnh báo nguy cơ trượt lở đất thông qua việc điều tra 392 vị trí trượt lở thực tế. Kết quả nghiên cứu mang lại giá trị khoa học to lớn trong việc hiện đại hóa phương pháp đánh giá tai biến địa chất, đồng thời cung cấp cơ sở dữ liệu thực tiễn tin cậy phục vụ công tác quy hoạch không gian lãnh thổ, bảo vệ an toàn công trình hồ đập và ổn định sinh kế cho hơn 500.000 người dân sinh sống tại lưu vực Tây Bắc.

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 địa chất động lực công trình hiện đại kết hợp lý thuyết phân tích không gian địa lý. Khung lý thuyết xem xét trượt lở đất như một hệ thống mở phức hợp, chịu sự chi phối đồng thời của các nhân tố nội sinh (hoạt động đứt gãy, biến dạng kiến tạo trẻ) và ngoại sinh (lực trọng trường, lượng mưa, mức độ phong hóa, xói mòn dòng chảy). Ba mô hình cốt lõi được tích hợp bao gồm: Mô hình định lượng chỉ số nhạy cảm trượt lở đất (Landslide Susceptibility Index - LSI), Mô hình phân tích thứ bậc đa tiêu chí (Analytic Hierarchy Process - AHP) do Thomas L. Saaty đề xuất năm 1994, và Mô hình cấu trúc vỏ phong hóa nhiệt đới.

Trong khuôn khổ nghiên cứu, 4 khái niệm nền tảng được chuẩn hóa:

  • Tai biến địa chất (Geological Hazard): Các quá trình tự nhiên hoặc nhân sinh có tiềm năng gây tổn hại nghiêm trọng đến tính mạng, tài sản và môi trường sinh thái.
  • Khối trượt lở đất (Landslide mass): Sự dịch chuyển trọng lực nhanh chóng của tầng đất đá bở rời xuống chân sườn dốc theo mặt trượt xác định, tính toán thể tích theo hình học bán elipxoid với công thức $V = \frac{1}{6} \pi \times D \times R \times S$ (trong đó $D$ là chiều dài, $R$ là chiều rộng, $S$ là chiều sâu lớp trượt).
  • Đới động lực đứt gãy: Vùng dập vỡ kiến tạo có mật độ nứt nẻ cao bao quanh các đứt gãy hoạt động, làm suy giảm mạnh cường độ kháng cắt của đá gốc.
  • Chỉ số nhạy cảm trượt lở (LSI): Tổng điểm trọng số tích hợp từ các lớp thông tin không gian biểu thị khả năng phát sinh tai biến tại từng ô lưới phân tích.

Hệ thống phân loại thạch học trong khu vực nghiên cứu được chuẩn hóa thành 5 nhóm đá chính: trầm tích bở rời Đệ Tứ, trầm tích gắn kết, đá biến chất, đá magma axit và đá cacbonat (đá vôi), kết hợp thang điểm đánh giá nhạy cảm từ 1 đến 9 bậc.

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

Nghiên cứu kết hợp đồng bộ nguồn dữ liệu đa nguồn gồm: ảnh vệ tinh quang học độ phân giải cao SPOT-5 (độ phân giải không gian 2,5m - 5m) và VNREDSat-1 (độ phân giải 2,5m kênh toàn sắc và 10m kênh đa phổ) thuộc đề tài mã số VT/UD-03/13-15; hệ thống bản đồ địa hình, địa chất, thủy văn, thảm phủ thực vật tỷ lệ 1/50.000; và chuỗi số liệu khí tượng thủy văn giai đoạn 2013-2015.

Cỡ mẫu nghiên cứu bao gồm 392 điểm trượt lở đất được định vị tọa độ GPS chính xác trên toàn bộ phạm vi 11.135 km². Phương pháp chọn mẫu áp dụng là phương pháp điều tra toàn diện kết hợp giải đoán ảnh viễn thám đa thời gian và kiểm chứng thực địa tại các khu vực trọng điểm như Mường Lay, Sìn Hồ, Thuận Châu, Mường La.

Lý do lựa chọn phương pháp phân tích: Việc tích hợp công nghệ Viễn thám và GIS cùng thuật toán AHP là giải pháp tối ưu cho vùng nghiên cứu diện rộng có địa hình hiểm trở khó tiếp cận trực tiếp. Viễn thám cho phép chiết xuất nhanh hiện trạng trượt lở, mạng lưới đứt gãy lineamen và biến động thảm phủ. Hệ thống thông tin địa lý (GIS) đóng vai trò chuẩn hóa hệ tọa độ không gian, số hóa 11 lớp bản đồ chuyên đề (độ dốc, chia cắt sâu, chia cắt ngang, lượng mưa, thạch học, vỏ phong hóa, nước ngầm, mật độ đứt gãy, đới động lực đứt gãy, thảm phủ thực vật, mật độ giao thông) và thực hiện phép chồng lớp trọng số đại số bản đồ chính xác đến từng pixel kích thước 30m x 30m.

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

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

Tổng hợp phân tích không gian trên 11.135 km² lưu vực hồ thủy điện Sơn La đã xác lập 392 khối trượt lở đất với quy mô dao động từ dưới 200 m³ đến trên 100.000 m³. Các phát hiện then chốt được lượng hóa cụ thể như sau:

  • Mối quan hệ phi tuyến giữa trượt lở đất và độ dốc địa hình: Tai biến trượt lở tập trung cao nhất ở cấp độ dốc 25°-35° với 125 khối trượt (chiếm 31,89% tổng số khối trượt trên diện tích 3.703 km²), tiếp theo là cấp độ dốc 15°-25° với 106 khối trượt (chiếm 27,04% trên diện tích 3.255 km²). Tại cấp độ dốc 35°-45°, ghi nhận 81 khối trượt trên diện tích 1.737 km². Ngược lại, cấp độ dốc >45° chỉ ghi nhận 9 khối trượt (chiếm 2,30% trên diện tích 293 km²) do cơ chế tại đây chuyển dịch chủ yếu sang dạng đổ lở đá gốc thay vì trượt lở tầng vỏ phong hóa.
  • Sự khống chế của mật độ chia cắt sâu địa hình: Mật độ chia cắt sâu trung bình từ 268 đến 536 m/km² là khu vực phát sinh trượt lở mãnh liệt nhất với 256 khối trượt (chiếm tỷ lệ áp đảo 65,31% tổng số khối trượt) phân bố trên diện tích 6.866 km². Cấp chia cắt sâu <268 m/km² ghi nhận 92 khối trượt trên diện tích 3.413 km², trong khi cấp chia cắt sâu >1072 m/km² hầu như không ghi nhận khối trượt nào (0 khối trượt trên 3 km²).
  • Mức độ tập trung theo mật độ chia cắt ngang: Cấp chia cắt ngang 360-700 m/km² có số lượng khối trượt lớn nhất với 147 vị trí (chiếm 37,50% trên 3.448 km²). Cấp chia cắt ngang <360 m/km² có 78 vị trí trên 3.523 km². Đáng chú ý, cấp chia cắt ngang rất cao >1400 m/km² dù chỉ chiếm diện tích 582 km² nhưng xuất hiện tới 24 khối trượt quy mô lớn, cho thấy mật độ tập trung tai biến cục bộ rất nghiêm trọng tại các vùng sườn thung lũng sông Đà.
  • Quy luật phân bố dạng tuyến theo cấu trúc kiến tạo: 392 điểm trượt lở không phân bố ngẫu nhiên mà kết cụm thành các dải tuyến tính rõ nét theo hai phương cấu trúc chính: phương á kinh tuyến (dải Phong Thổ - Mường Chà dọc thung lũng Nậm Nay và Quốc lộ 12) và phương Tây Bắc - Đông Nam (dải Quỳnh Nhai - Mường La, dải Lai Châu - Tân Uyên dọc tuyến Quốc lộ 32 và đới đứt gãy Hoàng Liên Sơn).

Thảo luận kết quả

Cơ chế phát sinh trượt lở đất tại lưu vực hồ thủy điện Sơn La là kết quả của sự tương tác bất lợi giữa cấu trúc địa chất nội sinh và tác động ngoại cảnh. Khi hồ chứa tích nước đến mực nước dâng bình thường, áp lực nước lỗ rỗng trong tầng đất đá ven bờ tăng cao, làm bão hòa tầng sét dẻo của vỏ phong hóa và hạ thấp hệ số ổn định trượt $K_s$ xuống dưới 1,0. Sự phân dị địa hình với độ chênh cao địa hình từ 500m đến hơn 3.000m tạo nên thế năng trọng lực cực lớn dọc theo các sườn dốc bị cắt xẻ bởi đứt gãy trẻ.

So sánh với các nghiên cứu trước đây tại khu vực Tây Bắc (như các khảo sát của Viện Địa chất và Viện Hàn lâm Khoa học và Công nghệ Việt Nam), nghiên cứu này đã khắc phục được tính chất định tính rời rạc nhờ ứng dụng mô hình AHP và dữ liệu ảnh vệ tinh VNREDSat-1. Việc lượng hóa trọng số cho thấy yếu tố độ dốc sườn giữ vai trò tối quan trọng (trọng số đạt 9 điểm trong ma trận so sánh cặp), theo sau là yếu tố mật độ đới động lực đứt gãy và lượng mưa năm.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua bảng ma trận tương quan đa nhân tố kết hợp biểu đồ phân bố tần suất tích lũy (Cumulative Frequency Curve). Kết quả phân vùng cảnh báo nguy cơ trên nền bản đồ GIS đã tách biệt rõ ràng 5 cấp độ nguy cơ: Rất thấp, Thấp, Trung bình, Cao và Rất cao. Các vùng có nguy cơ Cao và Rất cao chiếm khoảng 22,4% tổng diện tích lưu vực nhưng chứa tới 78,5% các điểm trượt lở đã và đang hoạt động, chứng minh độ chính xác và độ tin cậy vượt trội của mô hình không gian.

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

Dựa trên kết quả phân vùng nguy cơ trượt lở đất tại lưu vực hồ thủy điện Sơn La, 4 nhóm giải pháp kỹ thuật và quản lý trọng tâm được đề xuất:

  • Xây dựng hệ thống quan trắc và cảnh báo sớm tự động: Thiết lập mạng lưới trạm quan trắc địa kỹ thuật tự động trang bị cảm biến đo dịch chuyển sâu (inclinometer), đo áp lực nước lỗ rỗng (piezometer) và trạm đo mưa tự động tại 100% các điểm có nguy cơ trượt lở Rất cao dọc tuyến Quốc lộ 12, Quốc lộ 6 và khu vực sườn đập thủy điện Sơn La. Thời gian triển khai hoàn thành trước tháng 12/2026, do Bộ Giao thông Vận tải phối hợp cùng Tập đoàn Điện lực Việt Nam (EVN) chủ trì thực hiện.
  • Khoanh định ranh giới và số hóa bản đồ cảnh báo chi tiết: Chuyển giao bộ bản đồ phân vùng nguy cơ trượt lở đất tỷ lệ 1/25.000 cho Ủy ban nhân dân 4 tỉnh (Sơn La, Lai Châu, Điện Biên, Lào Cai) nhằm tích hợp vào đồ án quy hoạch sử dụng đất và quy hoạch dân cư nông thôn. Mục tiêu hoàn thành số hóa 100% dữ liệu không gian trong giai đoạn 2026-2027 do Sở Tài nguyên và Môi trường các địa phương điều phối.
  • Phục hồi thảm phủ thực vật và công trình sinh học phòng thái: Triển khai chương trình trồng rừng phòng hộ đầu nguồn và phủ xanh đất trống đồi núi trọc, nâng tỷ lệ che phủ rừng tại các lưu vực sườn dốc từ 25° đến 45° lên trên 70% vào năm 2030. Sở Nông nghiệp và Phát triển Nông thôn các tỉnh Tây Bắc phối hợp cùng các Ban Quản lý rừng phòng hộ chịu trách nhiệm giám sát và thực hiện.
  • Áp dụng giải pháp công trình gia cố sườn dốc và di dời tái định cư: Thi công kết cấu tường chắn rọ đá, neo đất đá tầng sâu kết hợp khoan thoát nước ngầm ngang tại 45 điểm có nguy cơ sạt trượt trực tiếp đe dọa các cụm dân cư tại huyện Mường Chà, Mường Lay và Quỳnh Nhai. Đồng thời, tổ chức di dời khẩn cấp 100% các hộ dân nằm trong vùng trượt lở nguy hiểm trước mùa mưa bão năm 2027 do Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn cấp tỉnh phụ trách.

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

Nội dung và sản phẩm khoa học của luận văn mang giá trị ứng dụng cao cho 4 nhóm đối tượng chuyên môn:

  • Cơ quan quản lý nhà nước và chính quyền địa phương: Ủy ban nhân dân và các Sở Tài nguyên và Môi trường, Giao thông Vận tải, Xây dựng tại các tỉnh miền núi phía Bắc sử dụng bản đồ phân vùng nguy cơ để lập quy hoạch không gian, bố trí điểm tái định cư an toàn và cấp phép xây dựng công trình hạ tầng.
  • Chủ đầu tư và đơn vị vận hành công trình thủy điện: Tập đoàn Điện lực Việt Nam (EVN) cùng Ban Quản lý Thủy điện Sơn La - Lai Châu ứng dụng dữ liệu để quan trắc ổn định sườn bờ hồ, dự báo bồi lắng lòng hồ do đất đá trượt lở và điều tiết quy trình xả lũ an toàn.
  • Các viện nghiên cứu, trường đại học và chuyên gia học thuật: Giảng viên, học viên cao học và sinh viên các ngành Bản đồ, Viễn thám, GIS, Địa mạo, Địa chất công trình khai thác phương pháp luận kết hợp mô hình AHP và viễn thám độ phân giải cao để phục vụ nghiên cứu giảng dạy.
  • Cán bộ quản lý rủi ro thiên tai và cứu hộ cứu nạn: Lực lượng phòng chống thiên tai cấp huyện, xã sử dụng cơ sở dữ liệu để xây dựng kịch bản ứng phó tình huống khẩn cấp, thiết lập phương án sơ tán dân cư cục bộ khi xảy ra mưa lũ cực đoan kéo dài.

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

Công nghệ viễn thám và GIS mang lại ưu thế gì vượt trội trong nghiên cứu trượt lở đất?

Viễn thám và GIS cho phép quét bao quát toàn bộ 11.135 km² địa hình hiểm trở mà phương pháp khảo sát truyền thống không thể thực hiện toàn diện. Ảnh SPOT-5 và VNREDSat-1 giúp định vị chính xác 392 điểm trượt lở, trong khi GIS tích hợp đồng bộ 11 lớp thông tin không gian để tự động hóa tính toán chỉ số nhạy cảm LSI với độ chính xác cao.

Tại sao cấp độ dốc từ 25° đến 35° lại ghi nhận số lượng điểm trượt lở cao nhất?

Cấp độ dốc 25°-35° chiếm tới 125 khối trượt (31,89%) vì đây là ngưỡng địa hình tích tụ tầng vỏ phong hóa bở rời đủ dày nhưng đồng thời có góc dốc tạo ra lực trượt lớn dọc sườn. Tại độ dốc này, chỉ cần lượng mưa đạt trên 1.600 mm/năm kết hợp tác động dòng chảy là đủ kích hoạt hiện tượng mất cân bằng trọng lực.

Vì sao khu vực có độ dốc trên 45° lại có rất ít điểm trượt lở được ghi nhận?

Tại cấp độ dốc >45°, nghiên cứu chỉ ghi nhận 9 khối trượt (2,30% trên diện tích 293 km²). Nguyên nhân là do góc sườn quá dốc khiến lớp vỏ phong hóa và trầm tích bở rời không thể tích tụ dày mà bị xói mòn ngay, quá trình phá hủy sườn tại đây diễn ra dưới dạng đổ lở đá gốc cục bộ thay vì trượt lở khối lớn.

Việc tích nước hồ thủy điện Sơn La tác động như thế nào đến sự ổn định của sườn dốc?

Khi lòng hồ chứa 11.135 km² lưu vực tích nước, mực nước dâng làm ngập chân các sườn dốc ven sông Đà. Hiện tượng ngâm nước làm bão hòa tầng đất đá chứa khoáng vật sét, làm giảm góc ma sát trong và lực dính kết, đồng thời tạo ra áp lực thủy động khi hạ mực nước hồ nhanh, kích hoạt các pha trượt lở bờ hồ liên tục.

Phương pháp phân tích thứ bậc AHP của Saaty được ứng dụng như thế nào trong đề tài?

Phương pháp AHP sử dụng ma trận so sánh cặp thông minh theo thang điểm từ 1 đến 9 để xác định trọng số khách quan cho từng nhân tố phát sinh trượt lở. Trọng số tính toán dựa trên vector riêng nguyên lý Eigen giúp loại bỏ sự chủ quan, lượng hóa chính xác vai trò áp đảo của độ dốc (9 điểm) và mật độ đứt gãy trong mô hình tổng thể.

Kết luận

  • Luận văn đã thiết lập thành công cơ sở dữ liệu số hóa hoàn chỉnh với 392 khối trượt lở đất phân bố trên diện tích 11.135 km² thuộc lưu vực hồ thủy điện Sơn La.
  • Khẳng định tính ưu việt của công nghệ tích hợp ảnh viễn thám độ phân giải cao VNREDSat-1, SPOT-5 và hệ thống thông tin địa lý (GIS) trong giải đoán tai biến địa chất.
  • Xác lập ma trận trọng số thứ bậc AHP khoa học, lượng hóa vai trò chi phối của độ dốc địa hình (đặc biệt ở dải 25°-35°) và đới động lực đứt gãy kiến tạo Tây Bắc - Đông Nam.
  • Xây dựng thành công bản đồ phân vùng cảnh báo nguy cơ trượt lở đất tỷ lệ 1/50.000, phân tách chính xác các khu vực nguy cơ Cao và Rất cao chiếm 22,4% diện tích.
  • Đề xuất hệ thống giải pháp công trình và phi công trình khả thi nhằm bảo vệ an toàn tuyệt đối cho công trình thủy điện trọng điểm quốc gia và cộng đồng dân sinh.

Trong giai đoạn 2026-2030, các cơ quan ban ngành cần tiếp tục ứng dụng và nâng cấp mô hình cảnh báo thời gian thực phục vụ quy hoạch bền vững vùng Tây Bắc. Hãy liên hệ và khai thác ngay nguồn dữ liệu từ công trình luận văn thạc sĩ này để nâng cao hiệu quả quản trị rủi ro thiên tai tại địa phương!