Tổng quan nghiên cứu
Lưu vực sông Sê San là một trong những mạng lưới thủy văn trọng điểm tại khu vực Tây Nguyên với tổng diện tích tự nhiên đạt 17.000 km², trong đó phần diện tích thuộc địa bàn tỉnh Kon Tum chiếm tới 11.659 km² cùng dòng chính dài 237 km. Với tổng lượng dòng chảy hàng năm vượt trên 10 tỷ m³, phân bố dòng chảy có sự phân hóa mùa cực đoan khi mùa lũ chỉ kéo dài 5 tháng nhưng chiếm từ 70% đến 75% tổng lượng nước năm, gây ra trung bình từ 4 đến 6 trận lũ lớn hàng năm. Dưới tác động kép của biến đổi khí hậu toàn cầu và các hoạt động xây dựng công trình bậc thang thủy điện quy mô lớn như Yaly và Pleikrông, hình thái lòng dẫn và bờ sông tại vùng thượng lưu Sê San liên tục chịu biến động mạnh mẽ, đe dọa trực tiếp đến an toàn dân sinh và cơ sở hạ tầng.
Mục tiêu cốt lõi của nghiên cứu là ứng dụng công nghệ viễn thám đa thời kỳ kết hợp hệ thống thông tin địa lý GIS nhằm xác lập sơ đồ biến động lòng dẫn, định lượng tốc độ xói lở và bồi tụ bờ sông, phân tích sự thay đổi lớp phủ bề mặt và dự báo xu thế biến đổi lòng dẫn trong tương lai. Phạm vi không gian của đề tài tập trung vào hai chi lưu huyết mạch là sông Đăk Bla (diện tích lưu vực 3.530 km², độ dốc 1,3% đến 1,7%) và sông Pô Kô (diện tích lưu vực 3.530 km², chiều dài 121 km) trên địa bàn tỉnh Kon Tum trong khung thời gian kéo dài 23 năm từ năm 1990 đến năm 2013. Kết quả nghiên cứu mang lại ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp cơ sở dữ liệu định lượng chính xác phục vụ công tác quy hoạch phòng chống thiên tai, bảo vệ nguồn tài nguyên đất cho hơn 453.206 người dân và đóng góp trực tiếp cho chương trình phát triển kinh tế - xã hội vùng Tây Nguyên thông qua các đề tài cấp nhà nước TN3/T16 và TN3/T20.
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 nền tảng lý thuyết tương tác bức xạ điện từ trong viễn thám và lý thuyết động lực học hình thái sông ngòi miền núi. Về mặt vật lý viễn thám, năng lượng bức xạ mặt trời khi chiếu tới mặt đất được phân tách thành ba thành phần cơ bản gồm năng lượng phản xạ, năng lượng hấp thụ và năng lượng truyền qua. Khả năng phản xạ phổ của các đối tượng tự nhiên biến thiên phụ thuộc vào bước sóng: chất diệp lục trong thực vật phản xạ mạnh ở dải sóng cận hồng ngoại từ 0,7 đến 1,1 µm, trong khi nước tự nhiên hấp thụ hầu như toàn bộ bức xạ ở dải hồng ngoại sóng ngắn. Nghiên cứu khai thác đặc tính nồng độ tạp chất lơ lửng trên 10 mg/l và độ sâu giới hạn 30 m để phân tách các ranh giới thủy văn.
Các khái niệm then chốt được chuẩn hóa bao gồm: Đặc trưng phản xạ phổ, Chỉ số khác biệt nước chuẩn hóa NDWI, Hệ thống phân tích đường bờ kỹ thuật số DSAS và Ma trận chuyển dịch lớp phủ sau phân loại. Kết hợp với mô hình thủy lực học dòng dẫn uốn khúc và quy luật xói mòn lòng sông do dòng chảy xiết vùng đồi núi có độ dốc từ 1,8% đến 3,3%, khung lý thuyết giúp định lượng chính xác sự thay đổi vị trí không gian của các dải bờ sông và bãi bồi qua các thời kỳ.
Phương pháp nghiên cứu
Nguồn dữ liệu đầu vào sử dụng chuỗi ảnh vệ tinh quang học Landsat gồm Landsat TM và ETM+ độ phân giải không gian 30 m thu nhận vào ba mốc thời gian tiêu biểu: tháng 1 năm 1990, năm 2002 và tháng 1 năm 2013. Toàn bộ ảnh được chuẩn hóa về hệ tọa độ UTM, phép chiếu WGS84 múi 48 Bắc. Phương pháp chọn mẫu ảnh được thực hiện nghiêm ngặt vào các tháng mùa khô (tháng 12 đến tháng 3 năm sau) với tỷ lệ mây che phủ dưới 10% nhằm loại trừ sai số do biến động mực nước theo mùa và đồng nhất điều kiện chiếu sáng mặt trời với độ ẩm trung bình không khí dao động từ 66% đến 78%.
Quy trình phân tích tích hợp xử lý ảnh số và phân tích không gian trong môi trường GIS:
- Xử lý ảnh và sửa lỗi sọc ảnh Landsat 7 ETM+ bằng các thuật toán nội suy chuyên dụng.
- Tách chiết ranh giới mặt nước tự động bằng tỷ số phổ kết hợp băng sóng Band 5, Band 7 và Band 2 kết hợp chỉ số NDWI nhằm phân tách chính xác ranh giới nước với các bãi cát sỏi có độ ẩm cao.
- Sử dụng thuật toán phân loại không giám sát K-Means phân chia diện tích lòng sông cực đại thành 4 nhóm lớp phủ chính: mặt nước, bãi cát sỏi, đất nông nghiệp và thảm thực vật cây bụi.
- Đo lường định lượng dịch chuyển bờ sông qua công cụ DSAS bằng cách tạo hệ thống các tuyến mặt cắt transect vuông góc với đường bờ cơ sở với khoảng cách đều nhau 50 m. Lý do lựa chọn phương pháp này là khả năng tự động hóa cao, hạn chế tối đa sai số chủ quan của người giải đoán và cho phép tính toán tốc độ biến động hàng năm với sai số vị trí dưới ± 15 m.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Phân tích chuỗi dữ liệu 23 năm đã chỉ ra các quy luật biến động hình thái lòng dẫn rất rõ nét trên hai hệ thống sông nghiên cứu:
Thứ nhất, trên sông Đăk Bla, hệ số uốn khúc cao đạt 2,03 cùng mật độ lưới sông 0,49 km/km² tạo nên các đoạn cong lõm liên tục bị xâm thực. Trong giai đoạn 1990 - 2002, bờ sông chịu xói lở chủ yếu ở các đoạn uốn khúc trung lưu với tốc độ trung bình 2,5 đến 4,8 m/năm. Bước sang giai đoạn 2002 - 2013, dưới tác động gia tăng của hoạt động khai thác cát sỏi và áp lực xây dựng hạ tầng đô thị thành phố Kon Tum, diện tích bãi bồi cát sỏi ven bờ suy giảm 18,5%, nhiều vị trí bờ lõm ghi nhận tốc độ sạt lở tăng vọt lên mức 6,2 m/năm, gây sạt trượt hàng nghìn mét khối đất bờ.
Thứ hai, trên sông Pô Kô, đoạn sông chảy qua xã Kroong thuộc huyện Sa Thầy nằm xen kẹp giữa hai hồ chứa thủy điện lớn Yaly và Pleikrông là điểm nóng xói lở nghiêm trọng nhất. Tốc độ sạt lở bờ sông tại khu vực này đã tăng từ mức 1,8 m/năm trong giai đoạn 1990 - 2002 lên mức 5,4 m/năm trong giai đoạn 2002 - 2013, tương đương mức tăng hơn 200% về cường độ xâm thực bờ.
Thứ ba, về biến động lớp phủ lòng sông cực đại, ma trận chuyển dịch giai đoạn 1990 - 2013 cho thấy có hơn 320 ha diện tích đất bãi bồi nông nghiệp và cây bụi ven bờ bị chuyển hóa thành mặt nước mở rộng hoặc bãi sỏi trơ trụi do các đợt lũ lịch sử với lưu lượng đỉnh vượt 2.500 m³/s gây xói trôi lớp đất mặt.
Thảo luận kết quả
Nguyên nhân cốt lõi dẫn đến sự gia tăng cường độ sạt lở bắt nguồn từ sự tương tác phức tạp giữa điều kiện tự nhiên khắc nghiệt và nhân tai. Địa hình lưu vực thượng nguồn dốc lớn trên 3,3%, lượng mưa trung bình năm tập trung cao từ 2.200 đến 2.520 mm kết hợp với nền địa chất phong hóa bở rời tạo tiền đề cho quá trình xói mòn. Đặc biệt, việc vận hành cắt đỉnh lũ và tích nước xả dòng đột ngột từ các công trình thủy điện đã làm thay đổi hoàn toàn chế độ thủy lực tự nhiên, tạo dòng chảy kiệt và dòng chảy xiết bất thường gây xói hàm ếch chân bờ.
Khi đối chiếu với các công trình nghiên cứu biến động sông Brahmaputra tại Ấn Độ hay sông Hồng tại đồng bằng Bắc Bộ, kết quả này hoàn toàn tương đồng về quy luật gia tốc xói lở khi lòng dẫn chịu tác động của các công trình ngăn dòng nhân tạo. Toàn bộ các phát hiện định lượng được thể hiện trực quan thông qua hệ thống 43 bản đồ chuyên đề biến động, sơ đồ mặt cắt transect và 9 bảng ma trận chuyển dịch lớp phủ, bảo đảm độ tin cậy của thuật toán phân loại đạt trên 85%.
Đề xuất và khuyến nghị
Dựa trên các bằng chứng khoa học và số liệu định lượng thu được, nghiên cứu đưa ra 4 nhóm giải pháp cấp bách nhằm giảm thiểu nguy cơ tai biến môi trường:
- Thiết lập hệ thống giám sát viễn thám và cảnh báo sớm: Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn tỉnh Kon Tum cần phối hợp với các cơ quan chuyên môn định kỳ 6 tháng một lần cập nhật dữ liệu ảnh vệ tinh quang học và radar đa thời kỳ, đặt mục tiêu giảm thiểu 40% thiệt hại kinh tế do sạt lở đất ven sông bất ngờ trong giai đoạn 2024 - 2028.
- Hoàn thiện quy chế vận hành liên hồ chứa thủy điện: Bộ Tài nguyên và Môi trường cùng Tập đoàn Điện lực Việt Nam cần tái cấu trúc quy trình xả nước điều tiết của các hồ thủy điện Yaly và Pleikrông, duy trì lưu lượng xả đáy ổn định nhằm khống chế tốc độ sạt lở bờ sông Pô Kô xuống dưới ngưỡng an toàn 2,0 m/năm trước năm 2026.
- Quy hoạch và siết chặt cấp phép khai thác khoáng sản lòng sông: Sở Tài nguyên và Môi trường tỉnh Kon Tum cần lập tức khoanh vùng cấm khai thác cát sỏi tại 15 điểm nóng xung yếu trên sông Đăk Bla, hoàn thành cắm mốc chỉ giới hành lang bảo vệ nguồn nước trong vòng 12 tháng.
- Khôi phục thảm thực vật bảo vệ bờ và gia cố kè sinh thái: Ủy ban nhân dân các huyện Sa Thầy và Đăk Hà cần triển khai chương trình trồng mới 500 ha rừng phòng hộ ven sông bằng các loài cây bản địa có hệ rễ sâu, hoàn thành giai đoạn 1 vào năm 2027 nhằm gia tăng độ ổn định mái dốc bờ sông thêm 35%.
Đối tượng nên tham khảo luận văn
Luận văn cung cấp nguồn tài liệu học thuật và thực tiễn giá trị cao cho 4 nhóm đối tượng chính:
- Cơ quan quản lý nhà nước về tài nguyên và quy hoạch xây dựng: Các Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển nông thôn có thể sử dụng bản đồ khoanh vùng nguy cơ sạt lở để lập quy hoạch sử dụng đất, cấp phép khai thác khoáng sản và bố trí các công trình hạ tầng kỹ thuật an toàn.
- Đơn vị quản lý, tư vấn thiết kế và vận hành công trình thủy điện: Các ban quản lý dự án thủy điện khai thác bộ dữ liệu biến động dòng dẫn để đánh giá tác động môi trường hạ lưu, tối ưu hóa quy trình điều tiết nước và thiết kế các công trình bảo vệ bờ đập.
- Các nhà khoa học, giảng viên và học viên sau đại học: Chuyên ngành Biến đổi khí hậu, Địa lý tài nguyên, Bản đồ - Viễn thám và GIS có thể kế thừa quy trình công nghệ kết hợp chỉ số NDWI, giải thuật K-Means và công cụ DSAS làm khung phương pháp chuẩn cho các lưu vực sông tương tự.
- Các tổ chức phát triển cộng đồng và ứng phó biến đổi khí hậu: Vận dụng số liệu thực địa để xây dựng các dự án sinh kế thích ứng bền vững, bảo vệ đất canh tác cho hơn 53% dân cư là đồng bào các dân tộc thiểu số tại khu vực Tây Nguyên.
Câu hỏi thường gặp
Vì sao ảnh vệ tinh Landsat đa thời kỳ lại được ưu tiên lựa chọn trong nghiên cứu biến động lòng sông miền núi? Chuỗi vệ tinh Landsat cung cấp tư liệu ảnh liên tục từ năm 1990 với độ phân giải không gian 30 m và hoàn toàn miễn phí. Khả năng bao quát toàn bộ lưu vực rộng 17.000 km² giúp khắc phục triệt để khó khăn trong việc tiếp cận địa hình đồi núi hiểm trở và tiết kiệm hàng tỷ đồng chi phí khảo sát thực địa.
Công cụ DSAS đo lường mức độ xói lở và bồi tụ bờ sông theo nguyên lý nào? Công cụ DSAS thiết lập các đường transect vuông góc với đường bờ cơ sở tại các khoảng cách đều nhau 50 m. Bằng cách tính toán giao điểm của transect với các đường bờ số hóa qua các năm 1990, 2002 và 2013, hệ thống tự động xuất ra các chỉ số dịch chuyển thực tế và tốc độ biến động tính theo đơn vị mét trên năm.
Sự xuất hiện của các hồ chứa thủy điện Yaly và Pleikrông tác động như thế nào đến lòng dẫn sông Pô Kô? Việc tích nước và xả lũ định kỳ làm đảo lộn quy luật dòng chảy tự nhiên. Lưu lượng nước thay đổi đột ngột làm gia tăng áp lực thủy động lực học lên mái dốc, thúc đẩy tốc độ sạt lở bờ sông tại khu vực xã Kroong tăng vọt lên mức 5,4 m/năm trong giai đoạn 2002 - 2013.
Làm thế nào để phân biệt chính xác ranh giới mặt nước với các bãi cát ẩm ven sông trên ảnh viễn thám? Nghiên cứu áp dụng đồng thời chỉ số khác biệt nước chuẩn hóa NDWI và tỷ số phổ kết hợp giữa các băng sóng hồng ngoại trung Band 5, Band 7 với băng xanh lục Band 2. Thuật toán này giúp triệt tiêu phản xạ phổ của cát ướt, bảo đảm độ chính xác phân tách mặt nước đạt trên 85%.
Kết quả nghiên cứu có thể ứng dụng trực tiếp cho công tác phòng chống thiên tai đô thị tại thành phố Kon Tum không? Có, hệ thống bản đồ nguy cơ và phương trình hồi quy dự báo cung cấp ranh giới sạt lở chính xác, giúp chính quyền thành phố Kon Tum chủ động cắm biển cảnh báo và lập phương án di dời kịp thời các hộ dân ven sông Đăk Bla trước mùa lũ hàng năm.
Kết luận
- Ứng dụng thành công công nghệ viễn thám Landsat đa thời kỳ và GIS để xây dựng quy trình giám sát biến động lòng sông miền núi với độ chính xác cao và chi phí tối ưu.
- Lượng hóa chi tiết tốc độ xói lở và bồi tụ trên hai nhánh sông Đăk Bla và Pô Kô qua chuỗi số liệu 23 năm (1990 - 2013).
- Chỉ rõ tác động gia tăng của các bậc thang thủy điện làm tốc độ sạt lở bờ sông Pô Kô tăng hơn 200% trong giai đoạn 2002 - 2013.
- Xây dựng hệ thống 43 bản đồ chuyên đề và 9 bảng ma trận chuyển dịch lớp phủ vùng lòng sông cực đại đạt độ tin cậy trên 85%.
- Đề xuất các nhóm giải pháp công trình và phi công trình khả thi nhằm bảo vệ an toàn sinh kế cho cộng đồng dân cư vùng thượng lưu Sê San.
Luận văn thạc sĩ là một công trình nghiên cứu hoàn chỉnh, đóng góp luận cứ khoa học quan trọng cho công tác quản trị lưu vực sông bền vững. Trong giai đoạn tiếp theo từ năm 2025 đến năm 2030, các đơn vị quản lý và nghiên cứu cần tiếp tục ứng dụng các thế hệ ảnh vệ tinh độ phân giải cao như Sentinel-2 kết hợp dữ liệu viễn thám radar để chủ động ứng phó với diễn biến phức tạp của biến đổi khí hậu trên toàn vùng Tây Nguyên.