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 xuyên biên giới trọng yếu của lưu vực sông Mê Kông với tổng diện tích tự nhiên đạt 17.000 km2, trong đó phần thượng lưu chảy qua địa phận tỉnh Kon Tum chiếm khoảng 11.659 km2. Dưới tác động phức hợp của biến đổi khí hậu toàn cầu, chế độ mưa lũ cực đoan và tốc độ gia tăng nhanh chóng của các dự án thủy điện bậc thang như Plei Krông hay Yaly, hiện tượng biến động hình thái lòng dẫn và sạt lở bờ sông tại hai chi lưu chính là sông Đăk Bla và sông Pô Kô đã diễn ra hết sức phức tạp. Vấn đề đặt ra là các phương pháp khảo sát trắc địa truyền thống trên địa hình đồi núi dốc đứng của Tây Nguyên thường gặp trở ngại lớn về chi phí và khả năng tiếp cận thực địa liên tục, dẫn đến sự thiếu hụt dữ liệu không gian dài hạn phục vụ công tác quy hoạch.

Luận văn thạc sĩ chuyên ngành Biến đổi khí hậu với đề tài nghiên cứu biến động lòng sông thượng lưu sông Sê San đã được thực hiện nhằm giải quyết trọn vẹn thách thức trên. Mục tiêu cốt lõi của nghiên cứu là xác lập sơ đồ biến động lòng sông, lượng hóa mức độ sạt lở, bồi tụ bờ sông và phân tích sự chuyển dịch lớp phủ vùng ngập lũ cực đại trong chuỗi thời gian kéo dài 23 năm từ năm 1990 đến năm 2013. Phạm vi không gian nghiên cứu bao phủ 23 xã, phường trọng điểm dọc theo hai chi lưu Đăk Bla và Pô Kô thuộc tỉnh Kon Tum. Về mặt ý nghĩa khoa học và thực tiễn, nghiên cứu đã xây dựng thành công bộ dữ liệu giám sát số hóa đạt độ tin cậy phân loại tổng thể 92,74%, giúp cơ quan quản lý khoanh vùng chính xác 100% các đoạn sông xung yếu, giảm thiểu rủi ro thiên tai và bảo vệ an toàn sinh kế cho hơn 453.000 người dân trên địa bàn tỉnh.

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 vững chắc của lý thuyết tương tác bức xạ điện từ trong viễn thám quang học và lý thuyết động lực học hình thái sông ngòi (Fluvial Geomorphology). Về mặt vật lý viễn thám, năng lượng mặt trời khi chiếu xuống bề mặt trái đất sẽ 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 thấu quang. Đối tượng mặt nước tự nhiên có đặc tính hấp thụ gần như hoàn toàn sóng điện từ ở dải cận hồng ngoại (NIR) và hồng ngoại sóng ngắn (SWIR), trong khi thảm thực vật và thổ nhưỡng lại có độ phản xạ rất mạnh ở các dải sóng này. Đây là cơ sở khoa học để phân tách ranh giới thủy văn với nền đất xung quanh.

Bên cạnh đó, nghiên cứu vận dụng các khái niệm chuyên ngành then chốt bao gồm:

  • Chỉ số tỷ số phổ cải biên $(b5+b7)/b2$: Tận dụng kênh hồng ngoại sóng ngắn Band 5, Band 7 và kênh màu lục Band 2 của ảnh vệ tinh Landsat nhằm tối ưu hóa độ tương phản phổ giữa nước, đất trống và thực vật.
  • Hệ thống phân tích đường bờ số (Digital Shoreline Analysis System - DSAS): Phương pháp toán học xây dựng các tuyến cắt ngang vuông góc để tính toán khoảng cách dịch chuyển và tốc độ xói lở, bồi tụ đường bờ qua các thời kỳ.
  • Ma trận biến động lớp phủ (Transition Matrix): Mô hình so sánh sau phân loại nhằm xác định cụ thể diện tích chuyển dịch giữa các trạng thái lớp phủ bề mặt trong lòng sông.

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

Cơ sở dữ liệu nghiên cứu sử dụng cỡ mẫu gồm 9 cảnh ảnh vệ tinh Landsat 5 TM và Landsat 7 ETM+ có độ phân giải không gian 30m, bao gồm các mốc thời gian tiêu biểu: 1990, 1994, 1997, 2002, 2004, 2008, 2010 và 2013. Tiêu chuẩn chọn mẫu dữ liệu được thực hiện nghiêm ngặt theo các tiêu chí: toàn bộ ảnh được chụp vào mùa khô (từ tháng 12 đến tháng 3 năm sau) để loại trừ sai số ngập úng mùa lũ, độ che phủ mây dưới 10%, và được quy chuẩn đồng nhất về hệ tọa độ UTM WGS84 Zone 49N. Đối với các ảnh Landsat 7 sau năm 2003 bị lỗi sọc quét (SLC-off), kỹ thuật bù khoảng trống dữ liệu trong phần mềm ENVI đã được áp dụng để khôi phục 100% diện mạo không gian ảnh.

Lý do lựa chọn chuỗi phương pháp phân tích tích hợp bắt nguồn từ tính ưu việt của kỹ thuật chiết tách phổ kết hợp không gian. Ảnh tỷ số $(b5+b7)/b2$ cho thấy 95% giá trị phổ tập trung trong khoảng từ 0,3 đến 6 với độ lệch chuẩn 0,8727, giúp thuật toán phân ngưỡng lát cắt mật độ (Density Slicing) tách mặt nước ở miền giá trị từ 0,17 đến 0,96 một cách chuẩn xác. Sau đó, thuật toán phân nhóm tự động K-means của MacQueen kết hợp khoảng cách Euclidean được ứng dụng để phân loại 4 nhóm lớp phủ chính mà không bị chi phối bởi định kiến chủ quan. Tốc độ biến động đường bờ được lượng hóa thông qua công thức $R = D_n / T$ (với $D_n$ là khoảng cách dịch chuyển và $T$ là khoảng thời gian giữa hai thời kỳ). Quá trình kiểm định chất lượng sử dụng ảnh vệ tinh độ phân giải cao GeoEye 0,4m cùng số liệu điều tra thực địa năm 2013, xác nhận độ chính xác tổng thể đạt 92,74% và hệ số Kappa vượt mức 0,85.

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à phân tích chuỗi tư liệu viễn thám đa thời kỳ giai đoạn 1990 - 2013 đã mang lại 4 phát hiện mang tính đột phá về quy luật biến động lòng dẫn thượng lưu sông Sê San:

  • Sự phân hóa mạnh mẽ về cường độ sạt lở và bồi tụ giữa hai chi lưu: Sông Đăk Bla có chiều dài thượng lưu khoảng 74 km với hệ số uốn khúc cao đạt 2,03 và độ dốc lòng dẫn 1,3% - 1,7%, thể hiện xu thế xói lở ngang rất mạnh tại các đoạn khúc uốn qua thành phố Kon Tum, tốc độ sạt lở cục bộ vượt ngưỡng 12 m/năm. Trong khi đó, sông Pô Kô với chiều dài 121 km và độ dốc đoạn thượng nguồn lên tới 3,3% lại có hình thái sạt lở dọc theo các thung lũng hẹp dạng chữ V với mức độ xói lở sâu.
  • Tốc độ dịch chuyển đường bờ gia tăng đột biến trong giai đoạn 2002 - 2013: So với giai đoạn 1990 - 2002, tốc độ xói lở bờ sông trong giai đoạn 2002 - 2013 đã tăng trung bình từ 1,5 đến 2,2 lần trên nhiều phân đoạn. Điểm nóng sạt lở nghiêm trọng nhất được ghi nhận tại khu vực xã Krong (nằm kẹp giữa hai lòng hồ thủy điện lớn Plei Krông và Yaly), nơi đường bờ bị dịch chuyển ngang hơn 45 mét trong suốt chu kỳ nghiên cứu.
  • Chuyển dịch cơ cấu lớp phủ lòng sông cực đại diễn ra phức tạp: Kết quả ma trận biến động chỉ ra diện tích mặt nước và bãi bồi cát sỏi liên tục hoán đổi vị trí cho nhau. Hơn 25% diện tích thảm thực vật cây bụi ven bờ đã bị xói trôi hoặc chuyển đổi sang đất canh tác nông nghiệp ngắn ngày, làm suy giảm khả năng giữ đất tự nhiên của thảm phủ sinh thái.
  • Quy luật hình thành bãi bồi ngầm và cồn cát: Quá trình bồi tụ diễn ra không đồng đều, tập trung chủ yếu ở phần bờ lõm thoải và hạ lưu các công trình nhân tạo, làm thay đổi dòng chảy chủ lưu và gián tiếp thúc đẩy sạt lở mạnh ở bờ đối diện.

Thảo luận kết quả

Nguyên nhân dẫn đến bức tranh biến động trên là sự kết hợp giữa các yếu tố tự nhiên nội tại và tác động nhân sinh mạnh mẽ. Về mặt tự nhiên, địa hình Kon Tum có tính phân dị cao, lượng mưa trung bình hàng năm lên tới 2.121 mm (có năm đạt đỉnh 2.520 mm), tập trung tới 70% - 75% lượng dòng chảy năm vào mùa lũ (tổng lượng dòng chảy toàn hệ thống đạt 8,61 tỷ m3/năm), trong khi mùa kiệt chỉ chiếm 3% - 5%. Sự chênh lệch lưu lượng cực đoan này tạo ra động lực dòng chảy cắt xẻ đáy và khoét chân bờ vô cùng dữ dội.

Về mặt nhân sinh, hoạt động khai thác cát, sỏi rầm rộ trên sông Đăk Bla đã làm hạ thấp cao trình đáy sông, phá vỡ thế cân bằng động của lòng dẫn. Đồng thời, việc tích nước và vận hành xả lũ của các hồ chứa thủy điện trên lưu vực sông Pô Kô đã làm biến đổi đột ngột trường thủy lực hạ lưu đập. Khi so sánh với các nghiên cứu của Zhengyi Yao trên sông Hoàng Hà hay Archana Sarkar trên sông Brahmaputra, kết quả tại sông Sê San thể hiện rõ tính đặc thù của các lưu vực sông miền núi nhiệt đới có độ dốc lớn kết hợp bậc thang thủy điện. Toàn bộ các phát hiện này được minh chứng trực quan thông qua bảng ma trận chuyển đổi không gian lớp phủ và hệ thống biểu đồ phân bố tốc độ dịch chuyển bờ dọc theo hàng trăm tuyến transect từ thượng nguồn về hạ lưu.

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

Nhằm kiểm soát tình trạng xói lở bờ sông, bảo vệ cơ sở hạ tầng và phát triển bền vững lưu vực thượng lưu sông Sê San, nghiên cứu đề xuất 4 nhóm giải pháp hành động cụ thể:

  • Hoàn thiện và ban hành hành lang thoát lũ, nghiêm cấm khai thác cát sỏi trái phép: Ủy ban Nhân dân tỉnh Kon Tum phối hợp với Sở Tài nguyên và Môi trường tiến hành cắm mốc chỉ giới hành lang bảo vệ nguồn nước cho 100% các đoạn bờ sông có nguy cơ sạt lở cao trước quý 4 năm 2025; đình chỉ toàn bộ các bãi khai thác khoáng sản không phép tại sông Đăk Bla.
  • Tối ưu hóa quy trình vận hành liên hồ chứa thủy điện: Tập đoàn Điện lực Việt Nam và các chủ hồ chứa (Plei Krông, Yaly) cần phối hợp chặt chẽ với cơ quan khí tượng thủy văn để xây dựng kịch bản điều tiết dòng chảy xả lũ khoa học, duy trì dòng chảy sinh thái tối thiểu, phấn đấu giảm 35% áp lực thủy động lực dòng chảy tác động trực tiếp vào chân bờ sông trong giai đoạn 2024 - 2030.
  • Xây dựng hệ thống công trình kè mềm kết hợp phục hồi sinh thái rừng ven sông: Sở Nông nghiệp và Phát triển Nông thôn tỉnh Kon Tum chủ trì triển khai chương trình trồng mới 500 ha rừng phòng hộ bán ngập và cây giữ đất bản địa có bộ rễ sâu tại các bãi bồi ven sông Pô Kô và Đăk Bla, hoàn thành phủ xanh 85% diện tích bờ xung yếu trước năm 2027.
  • Thiết lập hệ thống quan trắc và cảnh báo biến động lòng sông tự động bằng công nghệ viễn thám - GIS: Trung tâm Viễn thám Quốc gia phối hợp với các viện nghiên cứu xây dựng trạm giám sát không gian định kỳ 6 tháng/lần dựa trên ảnh vệ tinh đa thời kỳ, cung cấp thông tin dự báo xói lở sớm cho toàn bộ 23 xã, phường ven sông từ năm 2025.

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

Công trình nghiên cứu mang giá trị thực tiễn và học thuật sâu sắc, đặc biệt hữu ích cho 4 nhóm đối tượng sau:

  • Cơ quan quản lý nhà nước và chính quyền địa phương (Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển Nông thôn tỉnh Kon Tum): Cung cấp bộ bản đồ số hóa chi tiết và dữ liệu định lượng 23 năm để lập quy hoạch sử dụng đất ven sông, quy hoạch dân cư và cấp phép khai thác tài nguyên nước an toàn.
  • Ban quản lý các dự án thủy điện và thủy lợi trên lưu vực sông Sê San: Cung cấp cơ sở khoa học về động lực dòng chảy và diễn biến lòng dẫn để tính toán bồi lắng lòng hồ, bảo vệ an toàn chân đập và các công trình phụ trợ hạ du.
  • Chuyên gia nghiên cứu, giảng viên, học viên cao học và sinh viên chuyên ngành Biến đổi khí hậu, Viễn thám, GIS, Địa lý - Địa mạo: Sử dụng làm tài liệu tham khảo chuẩn mực về quy trình công nghệ chiết tách đường bờ bằng ảnh tỷ số $(b5+b7)/b2$, thuật toán K-means và mô hình DSAS.
  • Các doanh nghiệp tư vấn xây dựng công trình thủy công và hạ tầng giao thông ven sông: Khai thác số liệu về tốc độ xói lở tối đa theo từng mét/năm để thiết kế kết cấu kè bảo vệ, mố cầu và tuyến đường ven sông đạt độ bền vững tối ưu.

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

Tại sao nghiên cứu lại ưu tiên sử dụng ảnh vệ tinh Landsat độ phân giải 30m thay vì ảnh hàng không độ phân giải cao?
Ảnh Landsat sở hữu nguồn dữ liệu lịch sử liên tục, miễn phí và có tính lặp lại định kỳ cao từ năm 1990 đến năm 2013 với diện tích bao phủ rộng lớn 185 km x 170 km mỗi cảnh. Điều này giúp tối ưu hóa chi phí nghiên cứu và cho phép đánh giá tổng thể biến động lòng sông trên toàn lưu vực 11.659 km2 mà ảnh máy bay khó đáp ứng được do kinh phí đắt đỏ và tần suất bay chụp thưa thớt.

Công thức tỷ số phổ $(b5+b7)/b2$ có ưu điểm gì vượt trội trong việc chiết tách đường bờ sông?
Công thức này kết hợp hai kênh hồng ngoại sóng ngắn (Band 5 và Band 7) vốn bị nước hấp thụ cực mạnh cùng với kênh màu lục (Band 2) có độ phản xạ nước nhất định. Nhờ đó, tỷ số phổ tạo ra độ tương phản sắc nét giữa mặt nước và cả hai đối tượng bao quanh là đất trống và thảm thực vật dày, giúp loại bỏ nhiễu phổ và nâng độ chính xác phân ngưỡng lên 92,74%.

Đâu là yếu tố chính khiến đoạn sông Pô Kô chảy qua xã Krong bị sạt lở nghiêm trọng nhất?
Đoạn sông này chịu tác động kép từ điều kiện tự nhiên và nhân tạo: sông nằm trong thung lũng hẹp có nền địa chất trầm tích bở rời, đồng thời vị trí địa lý nằm kẹp giữa hai hồ thủy điện Plei Krông và Yaly khiến chế độ thủy lực dòng chảy bị dao động mạnh mỗi khi các nhà máy tích nước hoặc xả lũ đột ngột.

Hệ thống phân tích đường bờ DSAS xử lý sai số tại các khúc uốn cong phức tạp như thế nào?
Khi các tuyến cắt ngang (transect) giao cắt đường bờ tại nhiều điểm trên khúc uốn cong hoặc giao với đường cơ sở (baseline) ở góc nhọn, phần mềm sẽ áp dụng thuật toán tạo đường đệm baseline mịn và tự động lựa chọn điểm giao gần nhất. Kỹ thuật này giúp loại bỏ hoàn toàn các lỗi hình học và đảm bảo tính nhất quán của dữ liệu thống kê khoảng cách.

Quy trình nghiên cứu trong luận văn có thể chuyển giao cho các lưu vực sông khác tại Việt Nam không?
Quy trình tích hợp ảnh viễn thám đa thời kỳ, chỉ số tỷ số phổ và mô hình phân loại K-means hoàn toàn có thể nhân rộng hiệu quả cho các lưu vực sông lớn khác tại Việt Nam như sông Hồng, sông Thu Bồn hay sông Đồng Nai nhằm phục vụ giám sát tai biến xói lở định kỳ với chi phí thấp.

Kết luận

  • Luận văn đã thiết lập thành công chuỗi dữ liệu không gian 23 năm (1990 - 2013) phản ánh toàn diện bức tranh biến động lòng dẫn thượng lưu sông Sê San với độ chính xác kiểm định đạt 92,74%.
  • Xác định rõ sự khác biệt hình thái biến động giữa sông Đăk Bla (thiên về xói lở ngang, mở rộng khúc uốn) và sông Pô Kô (thiên về cắt xẻ dọc thung lũng dốc và chịu tác động thủy điện).
  • Chứng minh tính hiệu quả vượt trội của việc kết hợp ảnh tỷ số $(b5+b7)/b2$, thuật toán phân cụm K-means và công cụ DSAS trong nghiên cứu địa mạo thủy văn tại các địa bàn đồi núi hiểm trở.
  • Cung cấp luận cứ khoa học thực chứng để chính quyền tỉnh Kon Tum quy hoạch hành lang thoát lũ, kiểm soát khai thác cát sỏi và bảo vệ an toàn cho 23 xã, phường ven sông.
  • Mở ra định hướng ứng dụng công nghệ viễn thám thời gian thực kết hợp trí tuệ nhân tạo để tự động hóa công tác cảnh báo thiên tai sạt lở bờ sông trong giai đoạn 2025 - 2030.