Tổng quan nghiên cứu

Khu vực miền núi phía Bắc Việt Nam luôn là điểm nóng về thiên tai địa chất, trong đó tỉnh Yên Bái chịu ảnh hưởng nặng nề với địa hình có độ cao trung bình khoảng 600m so với mực nước biển và chiếm đến 67,56% diện tích vùng núi cao toàn tỉnh. Nghiên cứu tập trung vào phạm vi không gian trọng điểm gồm 3 đơn vị hành chính là xã Minh Bảo, xã Nam Cường và phường Yên Ninh thuộc thành phố Yên Bái, với tổng diện tích tự nhiên là 24,22 km2. Khu vực này nằm trong vùng khí hậu nhiệt đới gió mùa có lượng mưa trung bình hàng năm đạt khoảng 1.200 mm, độ ẩm từ 83% đến 87%, kết hợp với địa hình phân cắt mạnh và lớp vỏ phong hóa dày từ 5m đến hơn 50m.

Vấn đề nghiên cứu cốt lõi xuất phát từ thực trạng tốc độ đô thị hóa nhanh chóng, hoạt động bạt sườn đào taluy làm nhà và mở đường giao thông đã làm mất chân sườn dốc tự nhiên, kết hợp với các đợt mưa lớn kéo dài kích hoạt hàng chục vụ trượt lở đất mỗi năm, đe dọa trực tiếp đến tính mạng cùng tài sản của người dân. Mục tiêu cụ thể của luận văn là kết hợp phân tích ảnh viễn thám độ phân giải cao, khảo sát thực địa và ứng dụng công nghệ hệ thống thông tin địa lý (GIS) cùng mô hình cơ học tiền nghiệm (deterministic model) để xây dựng bản đồ phân vùng nguy cơ trượt lở đất tỷ lệ lớn 1/10.000.

Ý nghĩa của nghiên cứu thể hiện qua việc lượng hóa chỉ số an toàn mái dốc, chuyển đổi từ các phương pháp đánh giá định tính diện rộng sang mô hình dự báo định lượng theo pixel không gian. Kết quả nghiên cứu cung cấp cơ sở khoa học tin cậy giúp chính quyền địa phương nâng cao năng lực cảnh báo sớm, tối ưu hóa quy hoạch sử dụng đất đô thị và giảm thiểu từ 60% đến 80% thiệt hại do tai biến địa chất gây ra trong tương lai.

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 sự tích hợp của hai nền tảng lý thuyết cơ bản trong địa kỹ thuật và thủy văn học:

  • Lý thuyết cân bằng giới hạn sườn dốc vô tận (Infinite Slope Model): Dựa trên hệ thống phân loại trượt lở của Varnes (1984) và nguyên lý cơ học đất của Skempton & Hutchinson (1969). Mô hình xác định hệ số ổn định sườn dốc ($F_s$) thông qua tỷ số giữa tổng lực kháng cắt giữ trượt và lực gây trượt dọc theo mặt trượt phẳng song song với bề mặt địa hình.
  • Mô hình thủy văn địa hình TOPMODEL (Beven & Kirkby): Sử dụng để mô phỏng sự tích tụ dòng chảy ngầm và xác định chỉ số bão hòa nước cục bộ ($m$) trong tầng đất dưới tác động của các trận mưa cực đại.

Các khái niệm chuyên ngành then chốt bao gồm:

  • Hệ số an toàn ($F_s$): Tỷ số định lượng đánh giá mức độ ổn định mái dốc; khi $F_s \le 1,0$ sườn dốc rơi vào trạng thái mất ổn định và trượt lở xảy ra.
  • Chỉ số bão hòa ($m = h/D$): Tỷ lệ giữa chiều cao cột nước ngầm bão hòa ($h$) và tổng bề dày lớp đất phong hóa ($D$).
  • Lực dính kết tổng hợp ($C_s + C_r$): Sự kết hợp giữa lực dính tự nhiên của đất đá ($C_s$) và lực dính bổ sung do mạng lưới rễ cây tạo ra ($C_r$).
  • Vỏ phong hóa (VPH): Tầng sản phẩm bở rời tại chỗ phát triển trên các hệ tầng đá biến chất Sông Hồng, Ngòi Chi và Núi Voi.

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

Nguồn dữ liệu đầu vào của luận văn được thu thập đồng bộ từ nhiều nguồn:

  • Dữ liệu viễn thám và mô hình số độ cao (DEM) nội suy từ bản đồ địa hình tỷ lệ 1/10.000 để trích xuất bản đồ độ dốc ($\theta$), hướng dốc và diện tích lưu vực đơn vị ($a$).
  • Chuỗi số liệu lượng mưa ngày cực đại quan trắc qua hơn 30 năm tại trạm khí tượng Yên Bái, được tính toán ngoại suy theo phân bố Gumbel và Pearson Type III với chu kỳ lặp lại 20 năm (lượng mưa dự báo đạt mức trên 150 - 200 mm/ngày).
  • Bộ dữ liệu phân tích cơ lý đất từ hơn 50 mẫu thí nghiệm hiện trường và trong phòng, xác định các chỉ tiêu: góc ma sát trong ($\varphi$), dung trọng khô ($\gamma_d$), dung trọng bão hòa ($\gamma_{sat}$) và lực dính kết ($C_s$).

Phương pháp chọn mẫu thí nghiệm được thực hiện theo nguyên tắc ngẫu nhiên phân tầng dựa trên 5 kiểu vỏ phong hóa chính. Lý do lựa chọn mô hình tiền nghiệm tích hợp GIS là khả năng tính toán cơ học chi tiết từng ô lưới không gian (grid cell), phản ánh trực tiếp mối quan hệ vật lý giữa đặc tính đất đá, chế độ thủy văn và tải trọng thực vật, khắc phục triệt để tính chủ quan của các phương pháp thống kê hoặc chỉ số trọng số kinh nghiệm truyền thống.

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

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

Nghiên cứu đã làm sáng tỏ cấu trúc không gian và quy luật phát sinh tai biến trượt lở đất trên địa bàn 24,22 km2 của thành phố Yên Bái qua các phát hiện cụ thể:

  • Đặc điểm thạch học và cấu trúc vỏ phong hóa: Chiếm tỷ lệ áp đảo trong vùng là kiểu vỏ Sialferit trên đá trầm tích lục nguyên và metapelit (SAF) với 44,67% tổng diện tích, tiếp theo là vỏ hỗn hợp Sialit - Sialferit (SA-SAF) chiếm 23,86%, vỏ Sialferit - Ferosialit (SAF-FSA) chiếm 20,69%, vỏ Ferosialit (FSA) chiếm 5,10%, và vỏ Sialit (SA) chiếm 0,18%. Tầng phong hóa bở rời phát triển trên hệ tầng Ngòi Chi có chiều dày lớn từ 5m đến 15m, nhiều nơi đạt tới 40m - 60m, có thành phần giàu sét, bột, xen kẹp thạch anh vỡ vụn, khiến độ bền liên kết nội tại rất thấp.
  • Mối tương quan giữa độ dốc và thảm phủ thực vật: Toàn bộ các điểm trượt lở thực tế ghi nhận qua điều tra đều phân bố ở sườn có độ dốc trên 20 độ. Tại xã Nam Cường và phường Yên Ninh, rừng trồng keo tai tượng đơn loài và đồi chè chiếm phần lớn diện tích; hệ thống rễ cây keo chỉ cung cấp lực gia cố $C_r$ dao động từ 1,5 đến 3,0 kN/m2, thấp hơn nhiều so với rừng tự nhiên đa tầng tán (đạt trên 5,0 kN/m2).
  • Hiện trạng và quy mô các khối trượt: Đã điều tra và lập danh mục hàng chục điểm trượt lở điển hình. Quy mô khối trượt tại xã Nam Cường chủ yếu từ 30 m2 đến 60 m2 thuộc dạng trượt chảy bùn đất và đổ lở taluy. Tại xã Minh Bảo, quy mô trượt mở rộng từ 30 m2 đến hơn 500 m2. Bản đồ phân vùng nguy cơ bằng GIS đã nhóm lãnh thổ thành 4 cấp độ: Rất cao, Cao, Trung bình và Thấp, trong đó vùng nguy cơ Cao và Rất cao chiếm khoảng 25% đến 35% diện tích các sườn đồi dốc.

Thảo luận kết quả

Cơ chế gây trượt chính tại khu vực nghiên cứu được xác định là sự suy giảm đột ngột cường độ kháng cắt khi đất đá bão hòa nước. Vào mùa mưa lũ, lượng mưa ngày cực đại vượt ngưỡng thấm làm gia tăng nhanh chóng áp suất nước lỗ rỗng và đẩy chỉ số bão hòa $m$ tiến gần tới giá trị 1,0. Kết hợp với việc nhân dân tự ý bạt chân taluy để mở rộng diện tích xây dựng nhà ở tại các khu vực sườn dốc thuộc phường Yên Ninh và xã Nam Cường đã kích hoạt thế mất cân bằng trọng lực nghiêm trọng.

Kết quả dự báo từ mô hình tiền nghiệm có độ tin cậy vượt trội, với tỷ lệ trùng khớp giữa các điểm trượt lở quan sát thực tế và các vùng được dự báo có nguy cơ cao đạt trên 85%. Dữ liệu nghiên cứu được biểu diễn trực quan thông qua hệ thống bản đồ số chuyên đề đa lớp trong GIS, cho phép chồng lớp giữa bản đồ độ dốc, thảm phủ, hệ số dẫn nước $K$ và mô hình số độ cao DEM. So sánh với các nghiên cứu tại các vùng núi lân cận, phương pháp phân tích độ ổn định sườn dốc vô tận kết hợp mô hình thủy văn mang lại độ chi tiết cao, phục vụ hiệu quả cho cấp quản lý đô thị cơ sở thay vì chỉ mang tính định hướng vùng như các bản đồ tỷ lệ 1/100.000 trước đây.

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

Nhằm chủ động phòng tránh và giảm thiểu thiệt hại do trượt lở đất tại thành phố Yên Bái, luận văn đề xuất 4 nhóm giải pháp chiến lược và công trình cụ thể:

  • Hoàn thiện quy hoạch không gian và quản lý xây dựng sườn dốc: Ban hành quy chế nghiêm cấm các hoạt động đào bạt chân taluy tự phát có góc dốc lớn hơn 25 độ; thiết lập hành lang an toàn tối thiểu 30m đến 50m tính từ chân các sườn đồi có nguy cơ trượt lở cao. (Chủ thể: UBND thành phố Yên Bái, Phòng Quản lý Đô thị; Timeline: Hoàn thành trong giai đoạn 2015 - 2017; Target metric: Kiểm soát 100% hoạt động xây dựng trên đất dốc).
  • Xây dựng hệ thống công trình thoát nước và gia cố cơ học: Thi công rãnh đỉnh đón nước mặt sườn đồi, rãnh phân lưu và tường chắn rọ đá có tầng lọc ngược tại các điểm xung yếu dọc các trục đường dân sinh ở phường Yên Ninh và xã Minh Bảo. (Chủ thể: Ban Quản lý Dự án Đầu tư Xây dựng thành phố; Timeline: Định kỳ triển khai trước tháng 5 hàng năm; Target metric: Nâng hệ số an toàn mái dốc $F_s \ge 1,30$).
  • Tái cấu trúc thảm phủ sinh học chống xói mòn: Giảm thiểu việc khai thác trắng rừng trồng keo tai tượng theo chu kỳ ngắn; chuyển đổi sang mô hình nông lâm kết hợp, trồng xen cây bản địa rễ sâu (như lát hoa, lim xẹt, cỏ vetiver) trên các sườn dốc trên 20 độ để nâng cao lực dính rễ cây $C_r$ thêm 2,0 - 4,0 kN/m2. (Chủ thể: Hạt Kiểm lâm thành phố và các hộ dân nhận khoán đất rừng; Timeline: Thực hiện liên tục trong 3 - 5 năm; Target metric: Nâng độ che phủ rừng phòng hộ ổn định lên trên 60%).
  • Thiết lập mạng lưới quan trắc và cảnh báo thiên tai sớm: Lắp đặt các trạm đo mưa tự động tại các xã Minh Bảo, Nam Cường và mốc đo chuyển vị sườn dốc; xây dựng ngưỡng cảnh báo khi lượng mưa tích lũy 24 giờ vượt quá 100 mm hoặc lượng mưa 3 ngày liên tiếp vượt quá 200 mm để kịp thời sơ tán dân cư. (Chủ thể: Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn các cấp; Timeline: Vận hành từ năm 2015; Target metric: Đảm bảo 100% người dân vùng nguy cơ nhận tin báo trước 2 - 4 giờ).

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

Tài liệu luận văn là nguồn tư liệu khoa học giá trị cho 4 nhóm đối tượng chính:

  • Cơ quan quản lý nhà nước và chính quyền địa phương: Cán bộ Sở Tài nguyên và Môi trường, Sở Xây dựng, UBND thành phố Yên Bái sử dụng bản đồ phân vùng nguy cơ tỷ lệ 1/10.000 để rà soát quy hoạch đô thị, cấp phép xây dựng và lập phương án di dời dân cư vùng nguy hiểm.
  • Kỹ sư địa kỹ thuật và tư vấn thiết kế hạ tầng: Các đơn vị khảo sát, thiết kế công trình giao thông, thủy lợi và dân dụng ứng dụng các thông số cơ lý đất ($c$, $\varphi$, $\gamma$) và công thức tính toán $F_s$ để thiết kế mái dốc, tường chắn đạt chuẩn an toàn.
  • Giảng viên, học viên cao học và nhà nghiên cứu: Các chuyên gia trong ngành Khoa học Môi trường, Quản lý Môi trường, Địa chất ứng dụng khai thác quy trình tích hợp giữa mô hình thủy văn TOPMODEL, phương pháp cơ học sườn dốc và công nghệ GIS làm tài liệu giảng dạy, nghiên cứu chuyên sâu.
  • Các tổ chức phi chính phủ và quản lý dự án cộng đồng: Các điều phối viên dự án giảm nhẹ rủi ro thiên tai và thích ứng biến đổi khí hậu tham khảo khung đánh giá để xây dựng các mô hình sinh kế bền vững và truyền thông nâng cao nhận thức cộng đồng miền núi.

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

Mô hình tiền nghiệm (deterministic model) có ưu điểm gì so với các phương pháp thống kê trong đánh giá trượt lở?
Mô hình tiền nghiệm tính toán hệ số an toàn ($F_s$) dựa trên các định luật vật lý và cơ học đất thực tế thay vì dựa vào tần suất thống kê quá khứ. Phương pháp này cho phép đánh giá chính xác trạng thái mái dốc theo từng điều kiện mưa giả định với độ tin cậy vượt trên 85% ở tỷ lệ bản đồ lớn 1/10.000.

Tại sao thành phố Yên Bái có nguy cơ trượt lở đất cao dù độ cao trung bình chỉ khoảng 70m?
Nguy cơ cao bắt nguồn từ sự kết hợp của lớp vỏ phong hóa dày từ 5m đến hơn 40m trên nền đá biến chất bở rời của hệ tầng Ngòi Chi, độ dốc địa hình cục bộ thường vượt 20 độ, lượng mưa năm lớn trên 1.200 mm và áp lực đào đắp sườn đồi do đô thị hóa nhanh.

Rừng trồng keo tai tượng có vai trò như thế nào đối với việc ổn định sườn dốc?
Keo tai tượng giúp hạn chế xói mòn mặt nhưng bộ rễ ăn nông chỉ tạo lực dính kết $C_r$ khoảng 1,5 đến 3,0 kN/m2. Khi khai thác trắng theo chu kỳ ngắn 5 - 7 năm, đất bị phơi ải và mất hoàn toàn lực giữ cơ học, làm nguy cơ trượt lở tăng cao gấp nhiều lần so với rừng tự nhiên.

Các thông số cơ lý đất chính nào quyết định trực tiếp đến hệ số an toàn mái dốc trong nghiên cứu?
Ba thông số cốt lõi gồm lực dính kết của đất ($C_s$), góc ma sát trong ($\varphi$) và dung trọng bão hòa ($\gamma_{sat}$). Khi nước mưa thấm sâu làm đất chuyển sang trạng thái bão hòa ($m \to 1$), áp lực nước lỗ rỗng tăng lên làm triệt tiêu lực kháng ma sát, khiến hệ số $F_s$ tụt xuống dưới 1,0.

Chính quyền cơ sở có thể khai thác trực tiếp bản đồ nguy cơ của luận văn vào công việc nào?
Bản đồ phân vùng tỷ lệ 1/10.000 là công cụ khoanh định chính xác các ô đất có nguy cơ rất cao và cao tại 3 xã, phường nghiên cứu. Chính quyền sử dụng trực tiếp để cắm biển cảnh báo, dừng cấp phép xây dựng tại các sườn dốc mất an toàn và lập phương án sơ tán dân cư khi mưa bão.

Kết luận

  • Xây dựng thành công cơ sở dữ liệu số địa không gian và bản đồ hiện trạng trượt lở đất chi tiết trên diện tích 24,22 km2 thuộc địa bàn thành phố Yên Bái.
  • Ứng dụng hiệu quả mô hình cơ học sườn dốc vô tận kết hợp GIS để thành lập bản đồ phân vùng nguy cơ trượt lở đất tỷ lệ 1/10.000 với độ chính xác đạt trên 85%.
  • Lượng hóa vai trò của các yếu tố tự nhiên, xác định nhóm vỏ phong hóa SAF chiếm 44,67% và lượng mưa cực đại chu kỳ 20 năm là hai tác nhân tự nhiên kích hoạt trượt lở hàng đầu.
  • Chỉ rõ tác động tiêu cực của các hoạt động nhân sinh như đào bạt chân sườn đồi và chuyển đổi thảm thực vật kém bền vững đến độ ổn định mái dốc.
  • Đề xuất hệ thống 4 nhóm giải pháp kỹ thuật, quy hoạch, sinh học và quản lý có tính khả thi cao nhằm bảo vệ an toàn tính mạng và tài sản cho cộng đồng.

Đóng góp chính của luận văn là đã hoàn thiện quy trình công nghệ tích hợp viễn thám, GIS và mô hình tiền nghiệm trong dự báo tai biến trượt lở đất ở tỷ lệ lớn tại Việt Nam. Trong giai đoạn 2015 - 2020, nghiên cứu cần được tiếp tục mở rộng áp dụng cho toàn bộ các xã, phường còn lại của thành phố Yên Bái và tích hợp hệ thống cảm biến đo mưa tự động theo thời gian thực. Hãy chủ động ứng dụng kết quả phân vùng nguy cơ trượt lở đất vào quy hoạch phát triển kinh tế - xã hội ngay hôm nay để xây dựng một đô thị Yên Bái an toàn, bền vững trước biến đổi khí hậu.