Tổng quan nghiên cứu

Theo ước tính của Ngân hàng Thế giới, các hiểm họa trượt lở đất trên toàn cầu đã gây ảnh hưởng đến hơn 3,7 triệu km² diện tích tự nhiên và đe dọa trực tiếp đời sống của gần 300 triệu người. Tại Việt Nam, nguy cơ trượt lở diễn ra đặc biệt nghiêm trọng tại vùng miền núi phía Bắc, với 98 trên tổng số 135 huyện thị thường xuyên xảy ra sạt trượt trong mùa mưa bão. Thành phố Bắc Kạn là trung tâm chính trị, kinh tế của tỉnh Bắc Kạn nhưng lại có cấu trúc địa hình đồi núi dốc phân cắt mạnh và lớp vỏ phong hóa dày, khiến nguy cơ mất ổn định sườn dốc luôn ở mức báo động.

Vấn đề nghiên cứu trọng tâm xuất phát từ thực trạng tốc độ đô thị hóa nhanh chóng, dẫn đến việc người dân tự phát đào bạt chân sườn đồi để mở rộng mặt bằng xây dựng nhà ở và công trình giao thông. Khi xuất hiện các đợt mưa kéo dài với lượng mưa tích lũy vượt ngưỡng 180 mm từ tháng 5 đến tháng 8 hàng năm, các khối đất đá bị bão hòa nước và trượt lở đột ngột. Mục tiêu cụ thể của luận văn là điều tra, làm rõ hiện trạng sạt trượt, phân tích mối tương quan giữa các yếu tố môi trường địa chất và ứng dụng mô hình Chỉ số thống kê để thành lập bản đồ phân vùng nhạy cảm trượt lở đất.

Phạm vi không gian nghiên cứu bao phủ toàn bộ diện tích hành chính của thành phố Bắc Kạn, tập trung sâu vào 16 km dọc tuyến Quốc lộ 3 qua xã Xuất Hóa và các khu dân cư nội thị thuộc các phường Đức Xuân, Phùng Chí Kiên, Huyền Tụng. Luận văn tổng hợp dữ liệu thực địa cập nhật đến năm 2020. Nghiên cứu mang ý nghĩa thực tiễn lớn khi cung cấp cơ sở khoa học định lượng giúp địa phương giảm thiểu 100% các rủi ro quy hoạch thiếu an toàn, bảo vệ tài sản và tính mạng cho hơn 40.000 cư dân đô thị miền núi.

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 vận dụng hệ thống phân loại trượt lở kinh điển của Varnes kết hợp chuẩn hóa của Hutchinson, tập trung nhận diện ba cơ chế dịch chuyển chính gồm trượt xoay trong tầng phủ tàn tích, trượt tịnh tiến dọc mặt lớp đá gốc và trượt hỗn hợp trung gian. Đồng thời, đề tài áp dụng nguyên lý nền tảng trong dự báo tai biến địa chất của Carrara và Van Westen: hiện trạng và quá khứ là chìa khóa để dự báo tương lai.

Bốn khái niệm cốt lõi được chuẩn hóa xuyên suốt công trình gồm:

  • Chỉ số nhạy cảm trượt lở đất: Thước đo định lượng phản ánh xác suất không gian xảy ra sạt trượt dựa trên tổ hợp các điều kiện môi trường địa chất.
  • Vỏ phong hóa: Tầng sản phẩm biến đổi cơ học và hóa học từ đá gốc sét kết, bột kết, tạo thành lớp đất sườn tàn tích có độ bền liên kết suy giảm khi ngậm nước.
  • Hệ số an toàn mái dốc: Tỷ số giữa tổng lực chống trượt và lực gây trượt tác động lên một khối sườn dốc xác định.
  • Đơn vị điều kiện duy nhất: Đơn vị không gian đồng nhất được hình thành từ việc chồng xếp các lớp dữ liệu thành phần trên hệ thống thông tin địa lý.

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

Nguồn dữ liệu của luận văn được thu thập toàn diện từ 2 đợt khảo sát thực địa chi tiết, kết hợp khai thác hệ thống bản đồ địa chất - khoáng sản tỷ lệ 1:50.000, bản đồ địa hình số độ phân giải cao và chuỗi số liệu khí tượng thủy văn giai đoạn 2000-2020. Cỡ mẫu nghiên cứu bao gồm toàn bộ 100% các vị trí sạt trượt đã được ghi nhận trên địa bàn, kết hợp 8 lớp bản đồ chuyên đề thành phần gồm thạch học, địa mạo, sườn địa chất công trình, độ dốc, độ cao, phân cắt sâu, phân cắt ngang và hiện trạng sử dụng đất.

Phương pháp chọn mẫu áp dụng kỹ thuật kiểm kê toàn diện kết hợp phân tầng không gian nhằm phản ánh đầy đủ đặc trưng các vùng địa mạo. Mô hình Chỉ số thống kê kết hợp phần mềm chuyên dụng ArcGIS được lựa chọn làm phương pháp phân tích trọng tâm vì tính khách quan vượt trội, giúp loại bỏ các phán đoán chủ quan của chuyên gia. Trọng số của từng phân lớp nhân tố được xác định chính xác thông qua hàm logarit tự nhiên của tỷ lệ giữa mật độ điểm trượt trong phân lớp đó với mật độ điểm trượt trung bình trên toàn lãnh thổ thành phố.

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 không gian và tính toán thống kê đã chỉ ra 4 phát hiện quan trọng về quy luật sạt trượt tại thành phố Bắc Kạn:

Thứ nhất, phân bố trượt lở tập trung mật độ cao nhất tại các khu vực sườn dốc nhân tạo do hoạt động mở đường và bạt đồi xây dựng nhà ở, đặc biệt dọc 16 km tuyến Quốc lộ 3 phía Nam và đường vành đai phía Đông. Tại đây, các mái dốc đào bạt có góc dốc trên 45° và chiều cao vách vượt quá 4 m chiếm đến hơn 75% tổng số điểm sạt trượt được ghi nhận.

Thứ hai, đặc tính địa chất thạch học đóng vai trò kiểm soát cấu trúc mặt trượt. Hơn 60% các vụ trượt lở xảy ra trên nền đá trầm tích lục nguyên thuộc hệ tầng Phú Ngữ gồm bột kết, phiến sét bị phong hóa nứt nẻ mạnh. Hai kiểu trượt chủ yếu là trượt cung tròn trong tầng đất phong hóa mềm rời và trượt hỗn hợp trượt phẳng dọc mặt tiếp xúc giữa vỏ phong hóa với mặt đá gốc cứng chắc.

Thứ ba, lượng mưa là tác nhân kích hoạt trực tiếp nguy hiểm nhất. Dữ liệu khí tượng chứng minh các sự cố trượt lở quy mô lớn luôn xảy ra khi lượng mưa tích lũy 24 giờ vượt ngưỡng 180 mm trong các tháng cao điểm từ tháng 5 đến tháng 8. Đáng chú ý, tại nhiều điểm sạt trượt, độ che phủ của thảm thực vật vẫn đạt trên 60%, chứng minh rừng trồng tái sinh nghèo kiệt không đủ khả năng giữ ổn định cho các sườn dốc đã bị can thiệp cơ học vào chân taluy.

Thứ tư, mô hình phân vùng chỉ số nhạy cảm trượt lở đã chia địa bàn thành 5 cấp độ rõ rệt. Kết quả kiểm chứng mức độ tin cậy của mô hình thống kê cho thấy hơn 80% các điểm trượt lở thực tế nằm trùng khớp trong các vùng được dự báo có độ nhạy cảm cao và rất cao.

Thảo luận kết quả

Nguyên nhân căn bản của các hiện tượng trượt lở tại đô thị miền núi Bắc Kạn là sự kết hợp bất lợi giữa nền địa chất xung yếu và hoạt động đào khoét chân núi thiếu quy hoạch kỹ thuật. So với các công trình nghiên cứu diện rộng ở tỷ lệ 1:200.000 trước đây tại vùng Đông Bắc, nghiên cứu này mang lại bước tiến lớn về độ chi tiết và tính chính xác ở cấp độ vi đô thị tỷ lệ 1:25.000.

Kết quả phân vùng được biểu diễn trực quan thông qua bản đồ raster đa lớp kết hợp đồ thị đường cong tích lũy tần suất dự báo. Biểu đồ đối sánh giữa tỷ lệ diện tích các vùng nhạy cảm và tỷ lệ điểm trượt ghi nhận minh chứng rõ nét tính ưu việt của thuật toán: dù diện tích vùng nhạy cảm rất cao chỉ chiếm tỷ lệ nhỏ khoảng 15-20% tổng diện tích đô thị, nhưng lại bao hàm phần lớn các điểm nguy cơ trượt lở nghiêm trọng nhất.

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

Dựa trên bản đồ phân vùng nhạy cảm trượt lở, luận văn đề xuất 4 nhóm giải pháp kỹ thuật và quản lý đồng bộ:

  • Kiểm soát quy hoạch và chấm dứt đào bạt sườn dốc tự phát: Ủy ban nhân dân thành phố Bắc Kạn cần áp dụng bản đồ nhạy cảm trượt lở làm tiêu chí bắt buộc trong việc thẩm định quy hoạch đô thị giai đoạn 2021-2025. Cấm 100% các hoạt động tự ý san ủi chân núi làm nhà ở tại các khu vực có độ dốc trên 25° thuộc các phường nội thị như Đức Xuân, Phùng Chí Kiên, Huyền Tụng.
  • Xây dựng hệ thống công trình thoát nước bề mặt và nước ngầm: Ban Quản lý Dự án Đầu tư Xây dựng thành phố chủ trì thi công đồng bộ hệ thống rãnh đỉnh, mương phân dòng bậc thang dọc theo các sườn đồi tập trung dân cư trong giai đoạn 2022-2024. Mục tiêu là thu gom và tiêu thoát nhanh, giảm thiểu trên 70% lượng nước mưa ngấm vào các tầng phong hóa chứa mặt trượt tiềm năng.
  • Triển khai giải pháp gia cố cơ học cho các mái dốc xung yếu: Sở Giao thông Vận tải và chính quyền địa phương tiến hành kiên cố hóa các đoạn taluy xung yếu trên Quốc lộ 3 và khu dân cư bằng tường chắn bê tông trọng lực hoặc kè rọ đá có chiều cao từ 3,0 m đến 5,04 m, kết hợp neo đinh đất và trồng cỏ chống xói mòn bề mặt, hoàn thành dứt điểm trước năm 2026.
  • Thiết lập mạng lưới quan trắc và cảnh báo sớm thiên tai: Chi cục Thủy lợi và Phòng chống thiên tai phối hợp với các cơ quan chuyên môn lắp đặt hệ thống mốc đo dịch chuyển mái dốc cùng các trạm đo mưa tự động tại 100% các điểm có nguy cơ sạt lở rất cao, thiết lập quy trình phát lệnh sơ tán khẩn cấp ngay khi lượng mưa đo được vượt mốc 180 mm.

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

Luận văn là tài liệu chuyên khảo mang giá trị học thuật và ứng dụng thực tiễn cao cho 4 nhóm đối tượng:

  • Cơ quan quản lý nhà nước và nhà quy hoạch đô thị: Cán bộ Sở Xây dựng, Sở Tài nguyên và Môi trường tỉnh Bắc Kạn có thể ứng dụng trực tiếp bản đồ phân vùng LSI để khoanh định ranh giới an toàn cho các dự án phát triển đô thị, bố trí các khu tái định cư và tối ưu hóa chi phí đầu tư cơ sở hạ tầng.
  • Đội ngũ kỹ sư địa kỹ thuật và giao thông công trình: Các chuyên gia thiết kế công trình miền núi có thể khai thác số liệu tính toán ổn định taluy, các mặt cắt kè điển hình với chiều cao tính toán lên tới 5,04 m để áp dụng vào các dự án nâng cấp đường giao thông và xử lý điểm sụt trượt thực tế.
  • Lực lượng phòng chống thiên tai và tìm kiếm cứu nạn: Ban chỉ huy phòng thủ dân sự các cấp có thể sử dụng số liệu ngưỡng mưa 180 mm và bản đồ khoanh vùng nguy cơ cao để lập bản đồ phương án di dời dân cư, chủ động ứng phó trong các kịch bản thiên tai khẩn cấp.
  • Giảng viên, nghiên cứu sinh và học viên cao học: Giảng viên và học viên chuyên ngành Địa kỹ thuật, Địa chất công trình, Địa lý tài nguyên môi trường và GIS có thể sử dụng luận văn làm tài liệu tham khảo chuẩn mực về phương pháp luận tính toán trọng số thống kê đa biến.

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

Mô hình Chỉ số thống kê kết hợp GIS mang lại ưu điểm vượt trội nào trong đánh giá trượt lở?

Phương pháp Chỉ số thống kê lượng hóa chính xác mối quan hệ giữa 8 lớp nhân tố môi trường và vị trí trượt lở thực tế thông qua các phép tính mật độ khách quan. Cách tiếp cận này loại bỏ hoàn toàn tính chủ quan của người lập bản đồ, giúp xây dựng bản đồ phân vùng với tỷ lệ tin cậy trên 80%.

Tại sao nhiều khu vực có độ che phủ thực vật trên 60% tại thành phố Bắc Kạn vẫn xảy ra sạt lở lớn?

Thảm thực vật tại khu vực chủ yếu là rừng tái sinh nghèo kiệt có tầng rễ cạn dưới 1,5 m. Khi con người đào bạt chân taluy tạo mái dốc đứng trên 45° và xuất hiện các trận mưa dài ngày vượt 180 mm, tầng đất sâu trong vỏ phong hóa vẫn bị trượt phẳng dọc mặt tiếp xúc đá gốc.

Ngưỡng lượng mưa nào được xác định là kích hoạt sạt trượt nguy hiểm tại địa phương?

Qua tổng hợp chuỗi số liệu khí tượng thủy văn, ngưỡng kích hoạt trượt lở nghiêm trọng tại thành phố Bắc Kạn là lượng mưa tích lũy đạt trên 180 mm. Khi vượt qua ngưỡng này, áp lực nước lỗ rỗng tăng vọt làm triệt tiêu lực dính kết của tầng đất phong hóa, gây sụt lún mái dốc tức thì.

Các dạng trượt lở nào chiếm ưu thế tại thành phố Bắc Kạn?

Hai dạng trượt lở chủ đạo gồm trượt xoay cung tròn trong tầng vỏ phong hóa dày từ vài trăm đến hàng nghìn mét khối và trượt hỗn hợp phát triển kết hợp giữa tầng đất sườn tích với bề mặt lớp đá gốc của hệ tầng Phú Ngữ dọc theo các sườn đào nhân tạo.

Bản đồ nhạy cảm trượt lở đất phục vụ công tác quy hoạch đô thị ra sao?

Bản đồ phân tách toàn bộ không gian thành phố thành 5 cấp nhạy cảm rõ rệt. Các nhà quản lý sử dụng bản đồ để cấm 100% việc cấp phép xây dựng công trình kiên cố tại các vùng nhạy cảm cao, đồng thời định hướng phân bổ dân cư vào các dải địa hình an toàn.

Kết luận

  • Chuẩn hóa toàn diện cơ sở dữ liệu địa chất, địa mạo và hiện trạng sạt trượt trên địa bàn thành phố Bắc Kạn với 8 lớp thông tin không gian chuyên đề.
  • Ứng dụng thành công mô hình Chỉ số thống kê trên nền tảng ArcGIS, xác lập bản đồ phân vùng nhạy cảm trượt lở có độ chính xác và kiểm chứng thực địa đạt trên 80%.
  • Làm sáng tỏ tổ hợp nguyên nhân cốt lõi gây trượt lở gồm thạch học phân lớp yếu, lượng mưa tập trung vượt ngưỡng 180 mm và tình trạng đào bạt taluy dốc đứng trên 45°.
  • Đề xuất hệ giải pháp công trình kè chắn chịu lực cao đến 5,04 m kết hợp mương tiêu thoát nước mặt đồng bộ nhằm bảo vệ an toàn cho các trục giao thông huyết mạch và khu dân cư.
  • Thiết lập lộ trình thực thi giai đoạn 2021-2026, cung cấp công cụ đắc lực hỗ trợ chính quyền địa phương chuyển đổi từ ứng phó bị động sang chủ động quản lý rủi ro thiên tai bền vững.

Quý cơ quan, đơn vị tư vấn và các nhà nghiên cứu quan tâm có thể khai thác toàn bộ dữ liệu, bản đồ chuyên đề và quy trình tính toán chi tiết trong công trình này để triển khai ứng dụng vào thực tiễn quản lý quy hoạch xây dựng và phòng chống thiên tai tại địa phương.