Tổng quan nghiên cứu

Tài nguyên nước mặt đóng vai trò huyết mạch trong việc duy trì an ninh lương thực, an ninh năng lượng và phát triển kinh tế xã hội tại Việt Nam. Dù sở hữu mạng lưới sông ngòi dày đặc với tổng lượng dòng chảy mặt đạt khoảng 830 đến 847 tỷ m³/năm, Việt Nam vẫn đối mặt với nghịch lý khan hiếm nguồn nước khi lượng nước mặt bình quân đầu người chỉ đạt khoảng 3.840 m³/người/năm, thấp hơn tiêu chuẩn an toàn 4.000 m³/người/năm theo khuyến nghị quốc tế. Đặc biệt, khoảng 63% tổng lượng dòng chảy bắt nguồn từ ngoài lãnh thổ, khiến nguồn nước nội địa chịu nhiều rủi ro và phụ thuộc lớn vào các quốc gia thượng nguồn.

Lưu vực sông Vu Gia – Thu Bồn có diện tích tự nhiên khoảng 10.350 km², trải rộng trên địa bàn thành phố Đà Nẵng và tỉnh Quảng Nam, là một trong chín hệ thống sông lớn nhất cả nước. Khu vực này sở hữu tiềm năng thủy văn dồi dào nhưng phân bố cực kỳ không đồng đều theo không gian và thời gian. Tình trạng biến đổi khí hậu, nước biển dâng kết hợp với tốc độ đô thị hóa nhanh chóng đã làm gia tăng tần suất các hiện tượng thiên tai cực đoan như lũ quét, lũ lụt cục bộ trong mùa mưa và hạn hán, xâm nhập mặn nghiêm trọng trong mùa khô.

Luận văn thạc sĩ kỹ thuật của tác giả Nguyễn Mạnh Cường tập trung giải quyết bài toán định lượng dòng chảy mặt thông qua việc nghiên cứu ứng dụng công nghệ viễn thám đa phổ kết hợp với hệ thống thông tin địa lý GIS và mô hình thủy văn SWAT. Phạm vi nghiên cứu bao quát toàn bộ lưu vực Vu Gia – Thu Bồn với 33 tiểu lưu vực được phân định chi tiết. Kết quả của công trình mang lại giá trị thực tiễn to lớn, cung cấp cơ sở dữ liệu số hóa và công cụ mô phỏng chuẩn xác giúp nâng cao hiệu quả quản lý, điều tiết và bảo vệ bền vững nguồn tài nguyên nước mặt cho vùng kinh tế trọng điểm miền Trung.

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 trên nền tảng kết hợp giữa lý thuyết chu trình thủy văn tự nhiên và các mô hình không gian số địa lý hiện đại. Khung lý thuyết cốt lõi của đề tài bao gồm phương trình cân bằng nước tổng thể, mô hình thủy văn phân phối bán cấu trúc SWAT và lý thuyết phân loại ảnh viễn thám đa phổ có kiểm định. Trong đó, mô hình SWAT thực hiện tính toán sự vận động của nước dựa trên phương trình cân bằng ẩm với các thành phần chính gồm lượng mưa ngày, lượng dòng chảy mặt, lượng bốc thoát hơi nước thực tế, lượng nước thấm vào tầng đất chưa bão hòa và lượng dòng ngầm bổ cập.

Đề tài vận dụng đồng bộ năm khái niệm chuyên ngành cốt lõi:

  • Mô hình số độ cao DEM nhằm biểu diễn hình thái địa hình và trích xuất mạng lưới dòng chảy tự nhiên.
  • Đơn vị phản hồi thủy văn HRU đại diện cho các vùng đồng nhất về thổ nhưỡng, lớp phủ thực vật và độ dốc.
  • Hệ số phân loại Kappa dùng để kiểm soát và đánh giá độ chính xác của quá trình giải đoán ảnh viễn thám.
  • Lưu lượng dòng chảy mặt biểu thị qua thể tích nước di chuyển qua mặt cắt ngang lòng dẫn trong một đơn vị thời gian.
  • Chỉ số phản xạ phổ của các đối tượng tự nhiên làm căn cứ phân tách ranh giới lớp phủ bề mặt.

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

Nghiên cứu sử dụng nguồn dữ liệu không gian và thuộc tính đa nguồn được chuẩn hóa đồng bộ. Nguồn tư liệu viễn thám sử dụng ảnh vệ tinh Landsat với độ phân giải không gian 30m x 30m chụp vào tháng 2 năm 2010. Dữ liệu mô hình số độ cao DEM được trích xuất từ dữ liệu vệ tinh SRTM độ phân giải 30m của NASA. Bản đồ thổ nhưỡng được kế thừa từ cơ sở dữ liệu đất toàn cầu của FAO kết hợp bản đồ đất tỷ lệ 1:100.000 do Bộ Nông nghiệp và Phát triển Nông thôn xây dựng, phân chia rõ đặc tính lý hóa của 2 tầng đất chính là 0 đến 30 cm và 30 đến 100 cm. Dữ liệu khí tượng thủy văn được tích hợp từ hệ thống trạm quan trắc mặt đất và nguồn dữ liệu tái phân tích toàn cầu CFSR.

Quy trình phân tích dữ liệu được triển khai theo ba giai đoạn:

  • Giai đoạn 1: Tiền xử lý ảnh Landsat gồm hiệu chỉnh bức xạ, hiệu chỉnh khí quyển và hiệu chỉnh địa hình theo mô hình Minnaert.
  • Giai đoạn 2: Phân loại lớp phủ bề mặt có kiểm định bằng thuật toán xác suất cực đại Maximum Likelihood, giúp phân tách chính xác các kiểu sử dụng đất tự nhiên và nhân tạo.
  • Giai đoạn 3: Tích hợp các lớp thông tin DEM, thổ nhưỡng, lớp phủ và dữ liệu khí tượng vào mô hình ArcSWAT trên môi trường ArcGIS để tiến hành phân chia 33 tiểu lưu vực và hàng trăm đơn vị HRU, từ đó chạy mô phỏng lưu lượng dòng chảy theo từng tháng.

Phương pháp mô hình hóa ArcSWAT được lựa chọn vì khả năng mô phỏng liên tục theo thời gian, phù hợp với các lưu vực có địa hình phức tạp, diện tích trên 10.000 km² và điều kiện trạm đo thủy văn thực địa còn thưa thớt.

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

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

Nghiên cứu đã mang lại các kết quả định lượng cụ thể về lớp phủ và chế độ thủy văn lưu vực Vu Gia – Thu Bồn:

  • Thứ nhất, quá trình giải đoán ảnh vệ tinh Landsat bằng thuật toán Maximum Likelihood đã phân loại thành công các trạng thái lớp phủ chính với độ chính xác tổng thể cao, đạt hệ số Kappa trên 0,85. Kết quả này đảm bảo độ tin cậy tuyệt đối khi đưa thông tin lớp phủ vào mô hình thủy văn.

  • Thứ hai, mô hình đã số hóa và phân định toàn bộ lưu vực 10.350 km² thành 33 tiểu lưu vực thủy văn độc lập, phản ánh chính xác mạng lưới sông suối và hướng dốc dòng chảy tự nhiên từ vùng núi cao phía Tây sang vùng đồng bằng duyên hải phía Đông.

  • Thứ ba, kết quả mô phỏng thủy văn đã định lượng rõ rệt tính quy luật của dòng chảy mùa khô. Các tháng 2, 3, 4 và 5 là giai đoạn kiệt nước nghiêm trọng nhất trên toàn lưu vực, trong đó lưu lượng dòng chảy bình quân chạm đáy vào khoảng tháng 1 và tháng 2 với giá trị dòng chảy nhỏ nhất ghi nhận ở mức xấp xỉ 4,8 m³/s tại các tiểu lưu vực vùng núi cao.

  • Thứ tư, chế độ dòng chảy mùa lũ tập trung cao độ vào các tháng 8, 9, 10, 11 và 12. Dòng chảy đạt đỉnh lũ lịch sử vào các tháng 9, 10 và 11 với giá trị lưu lượng mô phỏng cực đại lên tới 226,5 m³/s. Tổng lượng dòng chảy trong giai đoạn mùa mưa này chiếm tới khoảng 70% đến 85% tổng lượng nước mặt của cả năm.

Thảo luận kết quả

Sự chênh lệch sâu sắc giữa lưu lượng đỉnh lũ 226,5 m³/s và lưu lượng đáy kiệt 4,8 m³/s giải thích nguyên nhân căn bản dẫn đến tình trạng lũ quét dồn dập vào cuối năm và hạn hán, nhiễm mặn gay gắt vào đầu mùa hè tại hạ lưu Đà Nẵng – Hội An. Nguyên nhân chính bắt nguồn từ sự kết hợp giữa địa hình đón gió của sườn Đông dãy Trường Sơn tạo nên các trung tâm mưa đặc biệt lớn như Ba Na, Trà Mi với lượng mưa hằng năm vượt mức 2.600 mm đến trên 5.000 mm, cùng với độ dốc sườn núi lớn khiến thời gian tập trung dòng chảy diễn ra rất nhanh.

Khi so sánh với các phương pháp thủy văn truyền thống vốn chỉ tính toán được lưu lượng tại các điểm trạm quan trắc cố định, việc ứng dụng mô hình SWAT tích hợp GIS đã chứng minh tính ưu việt vượt trội trong việc phân bố lưu lượng chi tiết trên từng tiểu lưu vực không gian. Trong thực tế, toàn bộ kết quả mô phỏng được trực quan hóa qua hệ thống bản đồ số chuyên đề và biểu đồ phân phối dòng chảy 12 tháng của 33 tiểu lưu vực. Cách trình bày dạng ma trận bảng biểu kết hợp bản đồ GIS phân vùng dòng chảy giúp các nhà quản lý dễ dàng nhận diện những tiểu lưu vực có nguy cơ lũ quét cao hoặc các vùng dễ bị cạn kiệt nguồn nước.

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

Dựa trên kết quả định lượng dòng chảy, nghiên cứu đưa ra 4 nhóm giải pháp chiến lược nhằm tối ưu hóa công tác quản trị tài nguyên nước:

  • Thiết lập hệ thống giám sát và dự báo dòng chảy thời gian thực: Ứng dụng công nghệ ảnh viễn thám chu kỳ chụp lặp 16 ngày kết hợp các trạm quan trắc tự động để cập nhật biến động lớp phủ và dự báo lưu lượng dòng chảy. Mục tiêu giảm sai số dự báo lũ và hạn hán xuống dưới 10%. Thời gian thực hiện từ năm 2024 đến năm 2026, do Sở Tài nguyên và Môi trường thành phố Đà Nẵng và tỉnh Quảng Nam chủ trì.

  • Hoàn thiện quy trình vận hành liên hồ chứa thủy điện và thủy lợi: Xây dựng biểu đồ điều tiết nước thông minh cho hệ thống hồ chứa bậc thang ở thượng nguồn dựa trên dữ liệu mô phỏng của 33 tiểu lưu vực. Đảm bảo duy trì lưu lượng xả tối thiểu trong mùa cạn đạt từ 15 đến 20 m³/s nhằm đẩy mặn cho các trạm bơm hạ du. Thời gian thực hiện giai đoạn 2025–2027, phối hợp giữa Bộ Tài nguyên và Môi trường cùng Tập đoàn Điện lực Việt Nam.

  • Phục hồi và bảo vệ rừng phòng hộ đầu nguồn: Triển khai các dự án trồng rừng bản địa và khoanh nuôi tái sinh rừng tại các tiểu lưu vực vùng núi có độ dốc cao, phấn đấu nâng tỷ lệ che phủ rừng toàn lưu vực lên trên 55% vào năm 2030. Giải pháp này giúp gia tăng khả năng giữ nước ngầm và nâng mô đun dòng chảy mùa kiệt thêm 30% đến 50%, do Chi cục Kiểm lâm địa phương và các ban quản lý rừng phòng hộ thực hiện.

  • Xây dựng nền tảng cơ sở dữ liệu GIS mở phục vụ quản lý tổng hợp: Số hóa toàn diện dữ liệu thổ nhưỡng, địa hình DEM 30m, lớp phủ và khí tượng thành một nền tảng WebGIS dùng chung cho 100% các cơ quan quản lý đô thị, nông nghiệp và phòng chống thiên tai trước năm 2028.

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

Công trình nghiên cứu mang tính ứng dụng liên ngành cao, là tài liệu tham khảo giá trị cho 4 nhóm đối tượng:

  • Cơ quan quản lý nhà nước về tài nguyên và môi trường: Các đơn vị như Sở Tài nguyên và Môi trường Đà Nẵng, Quảng Nam, Cục Quản lý Tài nguyên nước có thể khai thác quy trình mô phỏng để hoạch định quy hoạch lưu vực sông, cấp phép khai thác nước mặt và xây dựng kịch bản ứng phó biến đổi khí hậu.

  • Đơn vị quản lý vận hành công trình thủy lợi, thủy điện: Các kỹ sư vận hành hồ chứa thủy điện sông Tranh, A Vương, Đak Mi có thể sử dụng số liệu dòng chảy mô phỏng để tính toán dung tích đón lũ an toàn và lập kế hoạch xả nước phục vụ sản xuất nông nghiệp trong mùa khô.

  • Giảng viên, nhà khoa học chuyên ngành Đo đạc Bản đồ, Viễn thám và GIS: Đề tài cung cấp tài liệu tham khảo chi tiết về quy trình tiền xử lý ảnh vệ tinh quang học, kỹ thuật phân loại Maximum Likelihood và phương pháp tích hợp dữ liệu raster đa nguồn vào mô hình ArcSWAT.

  • Học viên cao học và sinh viên ngành Quản lý Môi trường, Kỹ thuật Tài nguyên nước: Luận văn là bản mẫu học thuật chuẩn xác về phương pháp xây dựng phương trình cân bằng nước, phân tích đơn vị thủy văn HRU và kỹ năng trình bày nghiên cứu định lượng.

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

  • Tại sao việc tích hợp viễn thám và GIS lại vượt trội hơn phương pháp quan trắc thủy văn truyền thống? Phương pháp truyền thống chỉ đo đạc cục bộ tại một số trạm thủy văn ven sông, tốn kém chi phí và khó bao quát các vùng núi non hiểm trở. Viễn thám và GIS cung cấp ảnh chụp diện rộng, mô hình hóa số độ cao DEM 30m và tính toán liên tục dòng chảy cho toàn bộ 10.350 km² lưu vực mà không bị giới hạn bởi địa hình.

  • Dữ liệu ảnh Landsat đóng vai trò gì trong mô hình tính toán lượng nước SWAT? Ảnh vệ tinh Landsat sau khi hiệu chỉnh khí quyển được sử dụng để phân loại các trạng thái lớp phủ bề mặt như rừng, đất nông nghiệp, mặt nước và đô thị. Mỗi loại lớp phủ có hệ số thấm và khả năng giữ nước khác nhau, là thông số đầu vào quyết định tính toán dòng chảy mặt trong mô hình.

  • Lưu lượng nước mặt tại lưu vực Vu Gia – Thu Bồn phân bố chênh lệch như thế nào giữa hai mùa? Nghiên cứu chứng minh sự tương phản sâu sắc: mùa khô từ tháng 2 đến tháng 5 có dòng chảy chạm đáy ở mức khoảng 4,8 m³/s gây thiếu nước hạ du; trong khi mùa mưa từ tháng 8 đến tháng 12 có dòng chảy đỉnh lũ lên tới 226,5 m³/s, chiếm tới trên 70% tổng lượng nước năm.

  • Mô hình số độ cao DEM 30m được ứng dụng như thế nào trong đề tài? Dữ liệu DEM SRTM 30m được xử lý trên môi trường GIS để tự động hóa trích xuất mạng lưới sông suối, xác định hướng dốc dòng chảy, ranh giới lưu vực tổng thể và phân chia chi tiết 33 tiểu lưu vực thủy văn độc lập.

  • Hệ số Kappa có ý nghĩa gì đối với kết quả nghiên cứu của luận văn? Hệ số Kappa là thước đo thống kê độ chính xác của bản đồ phân loại lớp phủ từ ảnh viễn thám. Với giá trị đạt trên 0,85, kết quả phân loại đảm bảo độ tin cậy tuyệt đối, loại bỏ các sai số ngẫu nhiên trước khi đưa dữ liệu vào mô phỏng dòng chảy SWAT.

Kết luận

  • Luận văn đã xây dựng thành công quy trình công nghệ tích hợp ảnh viễn thám Landsat, mô hình số độ cao DEM 30m và hệ thống thông tin địa lý GIS vào mô hình thủy văn SWAT.
  • Nghiên cứu đã định lượng hóa bức tranh thủy văn chi tiết trên 33 tiểu lưu vực của toàn bộ diện tích 10.350 km² lưu vực Vu Gia – Thu Bồn.
  • Làm sáng tỏ quy luật biến động dòng chảy cực đoan với sự chênh lệch lớn giữa mùa kiệt đạt mức xấp xỉ 4,8 m³/s và mùa lũ đạt đỉnh tới 226,5 m³/s.
  • Cung cấp cơ sở khoa học tin cậy phục vụ công tác quy hoạch liên hồ chứa, bảo vệ rừng phòng hộ và cảnh báo thiên tai tại miền Trung.
  • Đề xuất lộ trình mở rộng ứng dụng mô hình kết hợp dữ liệu vệ tinh đa thời gian để dự báo dòng chảy thích ứng với biến đổi khí hậu trong giai đoạn 2025–2030.

Hãy tiếp tục ứng dụng và mở rộng phương pháp tiếp cận không gian số này vào các lưu vực sông trọng điểm khác trên toàn quốc nhằm bảo đảm vững chắc an ninh nguồn nước quốc gia.