Tổng quan nghiên cứu

Tai biến trượt lở đất tại các vùng núi cao Tây Bắc luôn là mối hiểm họa địa chất nghiêm trọng, đe dọa trực tiếp đến tính mạng con người và phá hủy các công trình hạ tầng kinh tế - xã hội. Thống kê thiên tai năm 2013 tại huyện Bắc Yên, tỉnh Sơn La cho thấy các đợt mưa lũ cực đoan đã kích hoạt hàng loạt vụ sạt lở, làm hư hỏng 918 ngôi nhà, sập đổ hoàn toàn nền nhà của 48 hộ dân, phá hủy 348,267 ha hoa màu và gây ách tắc 17 điểm giao thông huyết mạch với khối lượng đất đá tràn lấp hơn 4.040 m³, ước tính tổng thiệt hại vượt mức 4,6 tỷ đồng.

Vấn đề nghiên cứu trọng tâm xuất phát từ đặc thù tự nhiên của huyện Bắc Yên với địa hình chia cắt sâu sắc, biên độ độ cao dao động lớn từ 120 m tại lòng hồ Sông Đà lên tới 2.982 m ở đỉnh Phù Sa Phin (độ cao trung bình 1.400 m), sườn dốc hiểm trở và lượng mưa tập trung cao điểm từ tháng 6 đến tháng 9 chiếm tới 85% tổng lượng mưa cả năm. Mục tiêu nghiên cứu cụ thể của đề tài là phân tích định lượng vai trò của các nhân tố tự nhiên - nhân sinh, xây dựng cơ sở dữ liệu không gian và thành lập bản đồ phân vùng nguy cơ tai biến trượt lở đất tỷ lệ 1:50.000 nhằm phục vụ công tác cảnh báo sớm và quy hoạch an toàn lãnh thổ.

Phạm vi không gian nghiên cứu bao trùm toàn bộ diện tích tự nhiên 110.371 ha của huyện Bắc Yên (chiếm 7,78% diện tích tỉnh Sơn La) gồm 16 đơn vị hành chính xã, thị trấn với 58.701 người dân sinh sống. Khung thời gian dữ liệu tập trung thu thập giai đoạn 2010–2014 kết hợp ảnh viễn thám năm 2012 và khảo sát thực địa năm 2014. Ý nghĩa của nghiên cứu thể hiện qua việc cung cấp giải pháp chuyển đổi từ ứng phó bị động sang phòng ngừa chủ động, giúp giảm thiểu tối đa rủi ro thiên tai và bảo vệ các tuyến giao thông huyết mạch như Quốc lộ 37 và Tỉnh lộ 112.

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 dựa trên ba trụ cột lý thuyết chuyên ngành vững chắc:

  1. Thuyết cân bằng giới hạn cơ học đất đá (Mohr-Coulomb): Giải thích cơ chế phá hủy mái dốc khi lực kéo trượt T thắng lực kháng cắt giữ trượt G. Trạng thái mất an toàn xảy ra khi hệ số ổn định F = G / T nhỏ hơn 1 dưới tác động của trọng lực khối trượt, áp lực nước thủy động và áp lực thủy tĩnh bão hòa trong đất đá.
  2. Phương pháp phân tích thứ bậc đa tiêu chuẩn Saaty (AHP - Analytic Hierarchy Process): Mô hình bán định lượng cho phép thiết lập ma trận so sánh cặp đôi giữa các nhân tố kích hoạt và nhân tố tiền đề để chuẩn hóa trọng số một cách khách quan với chỉ số nhất quán CR đạt độ tin cậy khoa học.
  3. Khung phân tích quản lý tổng hợp lưu vực sườn dốc (theo mô hình nghiên cứu của trường ITC - Hà Lan): Đánh giá trượt lở như một hàm đa biến của các hợp phần tự nhiên và nhân tác trong không gian địa lý lưu vực sông suối.

Hệ thống khái niệm học thuật cốt lõi trong luận văn bao gồm: Tai biến trượt lở đất (Landslide hazard), Hệ số ổn định sườn dốc (Factor of safety), Vùng đệm không gian (Buffer zone), và Chỉ số nhạy cảm trượt lở (Landslide Susceptibility Index).

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

Nguồn dữ liệu đầu vào phục vụ phân tích bao gồm: ảnh vệ tinh đa phổ Landsat độ phân giải không gian 30 m thu nhận năm 2012, mô hình độ cao số DEM nội suy từ bản đồ địa hình đường bình độ khoảng cao đều 10 m, bản đồ thổ nhưỡng tỷ lệ 1:100.000, bản đồ địa chất - kiến tạo, dữ liệu lượng mưa trung bình năm (1.300–2.000 mm) và hồ sơ lưu trữ thiệt hại thiên tai thực tế từ Ban chỉ huy Phòng chống lụt bão huyện Bắc Yên.

Cỡ mẫu nghiên cứu bao phủ toàn diện 110.371 ha lãnh thổ thông qua mạng lưới ô lưới pixel 30x30 m đồng nhất trên Hệ thống thông tin địa lý (GIS). Phương pháp chọn mẫu kiểm chứng thực địa được thực hiện bằng cách khoanh vùng và định vị tọa độ GPS tại các vết trượt lở lịch sử trên 4 xã trọng điểm (Mường Khoa, Song Pe, Xín Vàng, Chim Vàn) cùng toàn bộ hành lang 22 tuyến giao thông nông thôn và các tuyến Quốc lộ, Tỉnh lộ.

Lý do lựa chọn phương pháp phân tích ma trận chéo (Crossing matrix) kết hợp AHP đa tiêu chí trên nền tảng GIS là khả năng tích hợp đồng thời 9 lớp dữ liệu chuyên đề đa nguồn. Cách tiếp cận này loại bỏ tính chủ quan cục bộ của chuyên gia, xử lý định lượng mối tương quan phi tuyến giữa các điểm trượt lở thực tế với từng cấp độ dốc, thành phần thạch học và khoảng cách đứt gãy trong suốt giai đoạn thực hiện đề tài từ năm 2012 đến 2014.

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

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

Quá trình xử lý và chồng xếp 9 lớp dữ liệu chuyên đề đã mang lại 4 phát hiện quan trọng có giá trị khoa học và thực tiễn cao:

  1. Độ dốc và cấu trúc thổ nhưỡng là hai nhân tố chi phối mạnh mẽ nhất: Sườn dốc trên 25° tập trung phần lớn diện tích phân bố đất đỏ vàng trên đá sét biến chất Fs (chiếm 30% diện tích, khoảng 32.980 ha) và đất vàng nhạt trên đá cát Fq (chiếm 29% diện tích, khoảng 31.880 ha). Đây là những loại đất có tầng mỏng đến trung bình (30–100 cm), giàu khoáng vật sét dẻo và độ rỗng lớn, chiếm tới hơn 65% tổng số điểm trượt lở ghi nhận trên toàn huyện.
  2. Hành lang đứt gãy kiến tạo và mạng lưới giao thông hình thành các đới nguy cơ tập trung: Vùng đệm bán kính 100–300 m xung quanh các đứt gãy hoạt động và các tuyến đường Quốc lộ 37, Tỉnh lộ 112 có mật độ tai biến sạt lở cao hơn 3,2 lần so với các vùng sườn dốc tự nhiên chưa bị chia cắt. Việc bạt chân taluy cơ giới hóa kết hợp vùng dập vỡ kiến tạo làm giảm lực dính kết C xuống dưới ngưỡng giới hạn.
  3. Lượng mưa cực đoan đóng vai trò kích hoạt tức thời: Vùng có lượng mưa trung bình năm trên 1.900 mm ghi nhận nguy cơ trượt lở cao gấp 2,8 lần so với khu vực lượng mưa dưới 1.300 mm/năm. Đáng chú ý, 85% các vụ trượt lở diễn ra dồn dập vào mùa mưa lũ tập trung từ tháng 6 đến tháng 9 hàng năm.
  4. Phân vùng rủi ro thực tế: Toàn bộ 110.371 ha diện tích tự nhiên huyện Bắc Yên được phân định thành 5 cấp độ nguy cơ, trong đó vùng có nguy cơ cao và rất cao chiếm hơn 38% tổng diện tích, tập trung dày đặc tại các xã Mường Khoa, Song Pe, Xín Vàng và Chim Vàn.

Thảo luận kết quả

Cơ chế trượt lở tại huyện Bắc Yên được giải thích rõ nét qua sự kết hợp giữa trọng lực tự nhiên và sự suy giảm sức kháng cắt của tầng đất mặt. Khi những trận mưa lớn kéo dài trên 100 mm ngấm sâu vào lớp vỏ phong hóa Fs và Fq, khối lượng thể tích đất đá tăng mạnh từ 2,0 g/cm³ lên trên 2,3 g/cm³, đồng thời sinh ra áp lực nước thủy động đẩy trôi các liên kết thạch học kém bền vững dọc theo sườn dốc. So sánh với các nghiên cứu tương đồng tại vùng hồ thủy điện Sơn La và Hòa Bình của các nhà khoa học trước đây, kết quả nghiên cứu này hoàn toàn đồng nhất về quy luật phân bố nhưng có độ chi tiết không gian vượt trội nhờ tích hợp ảnh vệ tinh Landsat 30 m và mô hình DEM địa phương.

Toàn bộ kết quả định lượng được trình bày trực quan qua bản đồ phân vùng nguy cơ trượt lở tỷ lệ 1:50.000, kết hợp biểu đồ hình tròn phân bố tỷ lệ diện tích theo 5 cấp nguy cơ và bảng ma trận trọng số AHP 9x9 với chỉ số nhất quán CR đạt dưới 0,08 (thỏa mãn tiêu chuẩn khoa học quốc tế).

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

Dựa trên kết quả phân vùng nguy cơ tai biến trượt lở đất, 4 nhóm giải pháp chiến lược khả thi được đề xuất triển khai cụ thể:

  • Thiết lập mạng lưới quan trắc và cảnh báo sớm tự động: Lắp đặt 05 trạm đo mưa tự động và mốc quan trắc chuyển vị mái dốc tại các điểm xung yếu dọc tuyến Quốc lộ 37 và Tỉnh lộ 112. Target metric: Cung cấp bản tin cảnh báo trước từ 2 đến 6 giờ khi lượng mưa đạt ngưỡng kích hoạt 100 mm/24h. Timeline: Giai đoạn 2015–2018. Chủ thể thực hiện: Ban chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn huyện Bắc Yên phối hợp với Sở Tài nguyên và Môi trường tỉnh Sơn La.
  • Quy hoạch di dời và bố trí tái định cư chủ động: Tổ chức di dời khẩn cấp 48 hộ dân nằm trong vùng sạt lở nguy hiểm đã sụt lún nền nhà và lập phương án tái định cư an toàn cho 200 hộ dân nằm trong đới nguy cơ rất cao tại các xã Mường Khoa, Song Pe, Chim Vàn. Target metric: Giảm 100% rủi ro thương vong về người do sạt lở đất. Timeline: Lộ trình 2015–2020. Chủ thể thực hiện: Phòng Nông nghiệp và Phát triển nông thôn cùng UBND các xã liên quan.
  • Gia cố công trình hạ tầng giao thông và thoát nước sườn dốc: Áp dụng các giải pháp công trình như tường chắn rọ đá có lưới thép mạ kẽm, neo đất tầng sâu và hệ thống rãnh thoát nước mặt dạng bậc thang tại 17 điểm đen thường xuyên bị ách tắc. Target metric: Giảm 80% khối lượng đất đá sụt trượt tràn ra lòng đường (xuống dưới mức 800 m³/năm). Timeline: Giai đoạn 2015–2017. Chủ thể thực hiện: Sở Giao thông Vận tải tỉnh Sơn La.
  • Phục hồi thảm phủ thực vật và quản lý nương rẫy bền vững: Khoanh nuôi và trồng bổ sung 3.000 ha rừng phòng hộ đầu nguồn sông Đà, nâng tỷ lệ che phủ rừng toàn huyện từ 42% lên 55%. Khuyến khích chuyển đổi 1.500 ha nương rẫy trên sườn dốc trên 25° sang mô hình nông lâm kết hợp cây ăn quả có rễ cọc bám sâu. Target metric: Tăng độ liên kết cơ học của đất mặt thêm 35%. Timeline: Lộ trình 2015–2025. Chủ thể thực hiện: Hạt Kiểm lâm huyện Bắc Yên cùng cộng đồng dân tộc bản địa.

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

Nội dung nghiên cứu và hệ thống dữ liệu số hóa của luận văn mang lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng sau:

  • Cơ quan quản lý nhà nước và chính quyền địa phương cấp huyện/tỉnh: Sử dụng bản đồ nguy cơ tỷ lệ 1:50.000 để lồng ghép vào quy hoạch sử dụng đất, lập kế hoạch ngân sách phòng chống thiên tai và xây dựng bản đồ sơ tán dân cư chi tiết cho 16 xã, thị trấn.
  • Đơn vị tư vấn, khảo sát và xây dựng công trình giao thông - thủy lợi: Ứng dụng tập dữ liệu phân tích vùng đệm đứt gãy và hệ số dốc để thiết kế độ dốc taluy an toàn, tính toán tải trọng móng và bố trí công trình tiêu thoát nước khi cải tạo các tuyến đường miền núi như Quốc lộ 37.
  • Nhà nghiên cứu, giảng viên và học viên cao học chuyên ngành Quản lý tài nguyên môi trường, Địa lý, Địa chất, Viễn thám - GIS: Tham khảo quy trình tích hợp mô hình AHP - GIS, bộ khóa giải đoán ảnh vệ tinh Landsat 30 m và thuật toán ma trận chéo trong đánh giá các dạng tai biến thiên nhiên tương tự.
  • Các tổ chức phi chính phủ và cơ quan phát triển quốc tế: Khai thác kết quả đánh giá tính tổn thương cộng đồng tại các xã vùng cao như Xín Vàng, Chim Vàn để xây dựng các dự án sinh kế thích ứng biến đổi khí hậu và quản lý rủi ro thiên tai dựa vào cộng đồng.

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

Tại sao công nghệ Viễn thám và GIS lại có ưu thế vượt trội trong đánh giá trượt lở đất? Viễn thám cho phép cập nhật định kỳ biến động thảm phủ thực vật và phát hiện các vết trượt lở trên diện tích rộng lớn 110.371 ha mà khảo sát mặt đất khó tiếp cận toàn diện. Kết hợp với GIS, các mô hình không gian có thể chuẩn hóa và tích hợp đồng thời 9 lớp dữ liệu đa nguồn để tính toán chỉ số nguy cơ chính xác tới từng ô pixel 30 m.

Những yếu tố tự nhiên nào đóng vai trò quyết định nhất đến tai biến trượt lở tại huyện Bắc Yên? Yếu tố quyết định hàng đầu là độ dốc địa hình trên 25° kết hợp với lớp thổ nhưỡng đất đỏ vàng biến chất Fs (chiếm 30% diện tích) có tầng phong hóa dày giàu sét. Khi chịu tác động của lượng mưa tập trung trên 1.500 mm/năm (chiếm 85% trong mùa mưa), độ bền kháng cắt của đất giảm mạnh, trực tiếp gây mất ổn định sườn dốc.

Làm thế nào để xác định trọng số giữa các yếu tố gây trượt lở một cách khách quan? Nghiên cứu áp dụng phương pháp phân tích thứ bậc AHP của Saaty, thiết lập ma trận so sánh cặp đôi dựa trên tần suất xuất hiện thực tế của các vết trượt lở lịch sử trên từng lớp chuyên đề. Trọng số chỉ được chấp nhận khi tỷ số nhất quán CR nhỏ hơn 0,1, đảm bảo tính logic và độ tin cậy toán học.

Bản đồ nguy cơ trượt lở đất có thể ứng dụng trực tiếp vào quy hoạch xây dựng như thế nào? Bản đồ phân định rõ 5 cấp độ rủi ro từ rất thấp đến rất cao. Các cơ quan quy hoạch có thể sử dụng trực tiếp ranh giới này để cấm xây dựng công trình kiên cố hoặc dân cư tại vùng nguy cơ cao (chiếm trên 38% diện tích), đồng thời lựa chọn các khu vực an toàn để bố trí tái định cư bền vững.

Hoạt động canh tác nương rẫy của con người ảnh hưởng như thế nào đến độ ổn định sườn dốc? Việc đốt rừng làm nương rẫy trên 45.463 ha đất đồi dốc làm suy giảm độ che phủ rừng xuống còn 42%, phá hủy hệ thống rễ cây liên kết đất mặt. Khi mất đi tầng tán cản nước mưa, nước mặt chảy tràn xói mòn dữ dội và ngấm nhanh vào tầng ngầm, đẩy nhanh quá trình trượt lở đất đá.

Kết luận

  • Luận văn đã xác lập cơ sở khoa học vững chắc về cơ chế hình thành và quy luật phân bố không gian của tai biến trượt lở đất tại huyện Bắc Yên, tỉnh Sơn La.
  • Tích hợp thành công 9 lớp thông tin chuyên đề (độ dốc, địa chất, địa mạo, thổ nhưỡng, lượng mưa, thảm phủ, thủy văn, đứt gãy, giao thông) trên nền tảng công nghệ Viễn thám và GIS.
  • Thành lập hoàn chỉnh bản đồ phân vùng nguy cơ tai biến trượt lở đất tỷ lệ 1:50.000, nhận diện chính xác 38% diện tích tự nhiên thuộc vùng nguy cơ cao và rất cao.
  • Đóng góp giải pháp khoa học thực tiễn trong việc chuyển đổi phương thức ứng phó thiên tai từ cứu trợ thụ động sang cảnh báo sớm và quy hoạch phòng ngừa chủ động.
  • Lộ trình triển khai các giải pháp công trình và phi công trình được phân kỳ chi tiết theo các giai đoạn 2015–2018, 2015–2020 và tầm nhìn 2025.

Chính quyền địa phương và các cơ quan quản lý chuyên ngành cần sớm tiếp nhận và số hóa bộ cơ sở dữ liệu này vào hệ thống quản lý đô thị và quy hoạch nông thôn nhằm bảo vệ an toàn tính mạng, tài sản của nhân dân và thúc đẩy phát triển kinh tế - xã hội bền vững vùng cao Tây Bắc.