Tổng quan nghiên cứu
Lưu vực sông Nậm Mức đóng vai trò huyết mạch trong mạng lưới thủy văn Tây Bắc, với chiều dài dòng chính đạt 165 km và cung cấp tổng lượng nước hàng năm ước tính khoảng 1,62 tỷ m³ cho hệ thống sông Đà. Công trình thủy điện Nậm Mức có công suất lắp máy 44 MW chính thức phát điện từ tháng 5 năm 2015, giữ nhiệm vụ trọng yếu trong cung cấp năng lượng và điều tiết nguồn nước cho tỉnh Điện Biên. Tuy nhiên, dưới tác động của biến đổi khí hậu toàn cầu, các hiện tượng thời tiết cực đoan như mưa lớn tập trung và hạn hán kéo dài xuất hiện với tần suất ngày càng dày đặc, làm gia tăng nghiêm trọng hiện tượng xói mòn sườn dốc và bồi lắng lòng hồ.
Nghiên cứu tập trung giải quyết bài toán định lượng tác động của biến đổi khí hậu đến dòng chảy mặt và lượng bùn cát gia nhập hồ thủy điện Nậm Mức. Mục tiêu cụ thể là thiết lập công cụ mô phỏng thủy văn phân bố độ chính xác cao, nhận diện quy luật biến động dòng chảy theo mùa, từ đó dự báo rủi ro bồi lắng dung tích hữu ích và dung tích chết của hồ chứa. Phạm vi không gian nghiên cứu bao quát toàn bộ lưu vực sông Nậm Mức từ thượng nguồn đến tuyến công trình thủy điện. Phạm vi thời gian bao gồm chuỗi số liệu quan trắc thực đo lịch sử giai đoạn 1961 - 2012 đối với dòng chảy, giai đoạn 1993 - 2012 đối với bùn cát lơ lửng, cùng các kịch bản biến đổi khí hậu tầm nhìn đến năm 2099 theo công bố của Bộ Tài nguyên và Môi trường. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc bảo vệ tuổi thọ vận hành công trình trên 50 năm, tối ưu hóa quy trình điều tiết hồ chứa và giảm thiểu thiệt hại kinh tế hàng chục tỷ đồng do bồi lắng cơ học.
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 ứng dụng hệ thống lý thuyết thủy văn lưu vực hiện đại kết hợp cơ chế xói mòn - bồi lắng trong môi trường địa hình chia cắt mạnh. Trọng tâm của khung lý thuyết là mô hình SWAT (Soil and Water Assessment Tool) do Tiến sĩ Jeff Arnold thuộc Viện Nghiên cứu Nông nghiệp Hoa Kỳ phát triển. SWAT vận hành dựa trên phương trình cân bằng nước tổng thể, mô tả chi tiết các pha bay hơi, thoát hơi nước, thấm và dòng chảy mặt. Quá trình thấm bề mặt được tính toán thông qua phương pháp số đường cong dòng chảy thực nghiệm kết hợp lý thuyết Green và Ampt.
Đồng thời, quá trình xói mòn và vận chuyển bùn cát trên bề mặt lưu vực được xác định thông qua phương trình mất đất phổ biến cải tiến, tích hợp các yếu tố về năng lượng dòng chảy, độ nhám sườn dốc, độ bền cấu trúc thổ nhưỡng và chỉ số che phủ thảm thực vật. Toàn bộ lưu vực nghiên cứu được chia nhỏ thành các tiểu lưu vực và các đơn vị phản ứng thủy văn đồng nhất. Khung lý thuyết cũng tích hợp báo cáo đánh giá lần thứ năm của Ban liên chính phủ về biến đổi khí hậu, sử dụng hai kịch bản nồng độ khí nhà kính đại diện là kịch bản trung bình và kịch bản phát thải cao để xác định xu thế biến đổi nhiệt ẩm trong tương lai.
Phương pháp nghiên cứu
Nguồn dữ liệu đầu vào của nghiên cứu được xây dựng từ hệ thống thông tin không gian và chuỗi số liệu quan trắc khí tượng thủy văn thực tế. Bản đồ số độ cao địa hình, bản đồ phân loại thổ nhưỡng và bản đồ hiện trạng sử dụng đất tỉnh Điện Biên được số hóa chi tiết, phản ánh diện tích đất nông lâm nghiệp chiếm 79,31% và đất đồi núi dốc chưa sử dụng chiếm 18,41%. Chuỗi số liệu khí tượng thủy văn gồm 52 năm quan trắc lưu lượng dòng chảy (1961 - 2012) và 20 năm số liệu bùn cát lơ lửng (1993 - 2012) tại trạm thủy văn Nậm Mức. Cỡ mẫu nghiên cứu bao gồm toàn bộ chuỗi quan trắc liên tục hàng ngày, được chọn mẫu theo phương pháp toàn thể nhằm đảm bảo tính đại diện thống kê và tính liên tục của chuỗi thời gian thủy văn miền núi.
Phương pháp mô hình hóa toán học với phần mềm SWAT được lựa chọn nhờ khả năng mô phỏng liên tục theo thời gian, tái hiện chi tiết các quá trình vật lý phức tạp trên lưu vực có độ dốc lớn. Để nâng cao độ tin cậy, quy trình hiệu chỉnh và kiểm định bộ thông số được thực hiện thông qua công cụ tự động SWAT-CUP với thuật toán tối ưu hóa đa biến. Tiến trình nghiên cứu được triển khai theo timeline chuẩn hóa: phân tích thống kê chuỗi lịch sử thời kỳ nền (1986 - 2005), hiệu chỉnh mô hình với chuỗi số liệu thực tế, kiểm định độc lập và chạy mô phỏng dự báo dòng chảy, bùn cát theo các mốc thời gian đầu thế kỷ (2016 - 2035), giữa thế kỷ (2046 - 2065) và cuối thế kỷ (2080 - 2099).
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Quá trình mô phỏng thủy văn lưu vực Nậm Mức đã mang lại các kết quả định lượng chuẩn xác về dòng chảy và bùn cát:
Thứ nhất, mô hình SWAT đạt độ tin cậy cao sau khi hiệu chỉnh và kiểm định các chỉ số độ nhám lòng dẫn, hệ số dòng ngầm và dung số đường cong. Chỉ số hiệu quả Nash-Sutcliffe giai đoạn hiệu chỉnh và kiểm định dòng chảy đều đạt trên 0,78, trong khi chỉ số phần trăm sai số thể tích duy trì dưới mức 12%, khẳng định mô hình hoàn toàn phù hợp với điều kiện địa mạo vùng núi Tây Bắc.
Thứ hai, theo kịch bản biến đổi khí hậu, nhiệt độ trung bình năm tại khu vực nghiên cứu có xu thế tăng mạnh từ 0,7°C ở giai đoạn đầu thế kỷ lên 2,3°C ở cuối thế kỷ đối với kịch bản phát thải trung bình, và tăng tới 3,9°C đối với kịch bản phát thải cao so với thời kỳ nền. Lượng mưa năm tại lưu vực biến đổi theo xu hướng gia tăng từ 5% đến trên 15%, đặc biệt lượng mưa cực đoan 1 ngày và 5 ngày lớn nhất có xu thế tăng từ 40% đến 70%.
Thứ ba, dòng chảy mặt gia nhập hồ thủy điện Nậm Mức có sự phân hóa mùa sâu sắc. Lưu lượng mùa lũ (tập trung từ tháng 5 đến tháng 10, chiếm 77% đến 78% tổng lượng nước cả năm) có xu thế gia tăng từ 8,5% đến 18,2% theo các kịch bản tương lai. Ngược lại, dòng chảy mùa kiệt (tháng 11 đến tháng 4 năm sau) suy giảm từ 3,2% đến 7,6%, gây nguy cơ thiếu hụt nước vận hành tổ máy vào các tháng cao điểm mùa khô như tháng 2 và tháng 3.
Thứ tư, tổng lượng bùn cát lơ lửng gia nhập hồ chứa trung bình năm ước tính đạt khoảng 0,408 triệu tấn trong thời kỳ nền, nhưng dưới tác động của mưa lũ cực đoan theo kịch bản phát thải cao, lượng bùn cát mùa lũ dự báo gia tăng từ 16,4% đến 28,5% vào cuối thế kỷ, đẩy nhanh quá trình bồi lắng dung tích hồ chứa.
Thảo luận kết quả
Sự gia tăng đột biến của dòng chảy mùa lũ và lượng bùn cát bắt nguồn từ tác động cộng hưởng giữa biến đổi khí hậu và đặc thù địa hình địa chất. Địa mạo lưu vực sông Nậm Mức có độ dốc sườn núi lớn, cấu trúc địa chất dạng karst đá vôi xen kẹp sa diệp thạch kém ổn định, kết hợp với các trận mưa rào nhiệt đới cường suất cao làm gia tăng sức xói mòn động năng của hạt mưa và dòng tràn mặt đất.
Khi thảo luận kết quả, diễn biến dòng chảy thực tế được minh chứng rõ nét qua các đồ thị thủy văn mô phỏng dạng đường phân bố dòng chảy trung bình tháng và biểu đồ tương quan phân tán giữa số liệu thực đo và mô phỏng. Bảng tổng hợp dòng chảy cho thấy lưu lượng đỉnh lũ trung bình tháng 8 đạt mức 268,8 m³/s, cao gấp hơn 13 lần so với lưu lượng kiệt nhất vào tháng 3 (chỉ đạt 19,8 m³/s). So sánh với các công trình nghiên cứu bùn cát trên bậc thang thủy điện sông Đà như hồ Hòa Bình và Sơn La (nơi hiệu suất giữ bùn cát của các đập dâng dao động từ 73% đến 87%), hồ thủy điện Nậm Mức với dung tích nhỏ hơn sẽ chịu tốc độ bồi lấp tương đối nhanh hơn nếu không có giải pháp xả cát chủ động. Bên cạnh lượng bùn cát lơ lửng mô phỏng được, lượng bùn cát di đáy ước tính chiếm thêm từ 20% đến 40% tổng lượng bùn cát, tạo sức ép lớn lên hệ thống cửa lấy nước và tua-bin thủy lực của nhà máy.
Đề xuất và khuyến nghị
Dựa trên các kết quả dự báo định lượng, bốn nhóm giải pháp đồng bộ được đề xuất nhằm giảm thiểu tác động tiêu cực của biến đổi khí hậu:
- Tăng cường trồng và bảo vệ rừng phòng hộ đầu nguồn: Triển khai các dự án phục hồi thảm thực vật tự nhiên và giao khoán khoanh nuôi tái sinh, nâng tỷ lệ che phủ rừng lưu vực từ 38,67% hiện nay lên mức tối thiểu 45% vào năm 2030, do Ủy ban nhân dân tỉnh Điện Biên phối hợp cùng Sở Nông nghiệp và Phát triển nông thôn thực hiện nhằm giảm thiểu 15% lượng xói mòn đất bề mặt.
- Điều chỉnh quy trình vận hành hồ chứa thích ứng khí hậu: Xây dựng và ban hành quy trình xả lũ kết hợp xả bùn cát đáy định kỳ cho thủy điện Nậm Mức, đặt mục tiêu giảm thiểu 20% dung tích bồi lắng cơ học trong giai đoạn 2025 - 2035, do Ban Quản lý Nhà máy Thủy điện Nậm Mức chủ trì thực hiện dưới sự giám sát của Bộ Công Thương.
- Nâng cấp mạng lưới quan trắc khí tượng thủy văn tự động: Lắp đặt bổ sung từ 3 đến 5 trạm đo mưa và đo bùn cát tự động truyền số liệu thời gian thực trên lưu vực sông Nậm Mức trước năm 2028, do Đài Khí tượng Thủy văn khu vực Tây Bắc thực hiện nhằm nâng cao độ chính xác dự báo lũ lên trên 90%.
- Ứng dụng công trình cơ học bảo vệ bờ dốc và lòng dẫn: Xây dựng hệ thống kè rọ đá, đập chắn cát dạng ngầm tại các nhánh suối có độ dốc cao trước năm 2030, do Sở Tài nguyên và Môi trường tỉnh Điện Biên phối hợp với chính quyền địa phương triển khai nhằm giảm 25% lượng bồi lắng thô đi vào hồ chứa chính.
Đối tượng nên tham khảo luận văn
Nội dung và kết quả nghiên cứu của luận văn mang lại giá trị thực tiễn cho nhiều nhóm chuyên môn:
- Cán bộ quản lý và kỹ sư vận hành nhà máy thủy điện: Vận dụng số liệu mô phỏng dòng chảy và bùn cát để xây dựng phương án nạo vét lòng hồ, lập lịch bảo dưỡng tua-bin và điều tiết xả cát hợp lý cho công trình công suất 44 MW.
- Các nhà nghiên cứu và giảng viên chuyên ngành Thủy văn - Biến đổi khí hậu: Sử dụng bộ thông số mô hình SWAT, quy trình thiết lập dữ liệu không gian DEM và phương pháp hiệu chỉnh SWAT-CUP làm tài liệu tham khảo cho các nghiên cứu thủy văn lưu vực miền núi dốc.
- Cơ quan quản lý nhà nước về tài nguyên nước và môi trường: Các sở ban ngành tại Điện Biên và khu vực Tây Bắc có thể sử dụng dữ liệu dự báo để xây dựng quy hoạch tổng thể tài nguyên nước và kế hoạch hành động ứng phó biến đổi khí hậu tầm nhìn 2050.
- Đơn vị tư vấn thiết kế công trình thủy lợi - thủy điện: Tham khảo phương pháp tính toán bùn cát lơ lửng và tỷ lệ bùn cát di đáy (20% - 40%) để xác định chính xác cao trình ngưỡng lấy nước và dung tích chết của các dự án thủy điện bậc thang.
Câu hỏi thường gặp
-
Sông Nậm Mức có những đặc trưng hình thái và thủy văn nổi bật nào? Sông Nậm Mức là chi lưu cấp 1 bờ phải của sông Đà với chiều dài 165 km, chảy qua địa hình núi cao hiểm trở có độ chia cắt sâu. Chế độ dòng chảy phân hóa thành hai mùa rõ rệt, trong đó mùa lũ tập trung từ tháng 5 đến tháng 10 chiếm gần 78% tổng lượng dòng chảy năm, lưu lượng đỉnh lũ quan trắc lịch sử từng đạt 460 m³/s tại trạm Nậm Mức.
-
Vì sao mô hình SWAT được lựa chọn để mô phỏng lưu vực sông Nậm Mức? Mô hình SWAT sở hữu khả năng phân chia lưu vực thành các tiểu lưu vực và đơn vị phản ứng thủy văn chi tiết dựa trên bản đồ số độ cao DEM, loại đất và thảm phủ. Mô hình tích hợp công cụ tự động hóa SWAT-CUP giúp tối ưu hóa bộ thông số, tính toán đồng thời dòng chảy và lượng bùn cát xói mòn trong điều kiện chuỗi dữ liệu dài trên 50 năm.
-
Nhiệt độ và lượng mưa lưu vực thay đổi ra sao theo các kịch bản tương lai? Theo kịch bản phát thải trung bình và cao, nhiệt độ trung bình tại Điện Biên dự báo tăng từ 0,7°C lên 3,9°C vào cuối thế kỷ 21. Lượng mưa năm gia tăng từ 5% đến trên 15%, đặc biệt các đợt mưa lớn cực đoan có cường suất tăng từ 40% đến 70%, làm gia tăng đáng kể nguy cơ lũ quét và xói mòn sườn núi.
-
Biến đổi khí hậu tác động cụ thể thế nào đến lượng bùn cát vào hồ Nậm Mức? Lượng bùn cát gia nhập hồ có xu hướng tăng mạnh trong các tháng mùa lũ, với mức tăng dự báo từ 16,4% đến 28,5% theo kịch bản phát thải cao. Hiện tượng gia tăng bùn cát lơ lửng kết hợp bùn di đáy làm bồi lấp nhanh dung tích chết của hồ chứa, gây suy giảm tuổi thọ vận hành và mài mòn các thiết bị cơ khí thủy lực.
-
Lượng bùn cát di đáy được xác định bằng phương pháp nào khi thiếu số liệu đo đạc trực tiếp? Trong điều kiện thiếu thiết bị đo đạc chuyên dụng cho bùn cát đáy, nghiên cứu áp dụng quy chuẩn tính toán thủy văn kết hợp công thức thực nghiệm chuyên ngành, xác định lượng bùn cát di đáy tương đương khoảng 20% đến 40% lượng bùn cát lơ lửng quan trắc được tại mặt cắt khống chế trước khi gia nhập lòng hồ.
Kết luận
- Ứng dụng thành công mô hình SWAT trong việc tái hiện quy luật thủy văn và xói mòn bùn cát trên lưu vực sông Nậm Mức với độ chính xác cao.
- Khai thác hiệu quả chuỗi dữ liệu 52 năm lưu lượng và 20 năm bùn cát thực đo, làm cơ sở khoa học vững chắc cho việc hiệu chỉnh bộ thông số mô hình.
- Định lượng rõ nét xu thế phân hóa dòng chảy theo mùa, trong đó dòng chảy và bùn cát mùa lũ dự báo gia tăng tới 28,5% theo kịch bản phát thải cao.
- Cung cấp luận cứ kỹ thuật quan trọng phục vụ quản lý, vận hành an toàn nhà máy thủy điện Nậm Mức công suất 44 MW trong điều kiện khí hậu biến đổi.
- Thiết lập khung giải pháp kết hợp công trình và phi công trình với mục tiêu nâng độ che phủ rừng lên trên 45% trước năm 2030.
Nghiên cứu mở ra hướng đi mới trong việc tích hợp mô hình phân bố không gian với các kịch bản khí hậu hiện đại phục vụ quản trị rủi ro hồ chứa miền núi. Trong giai đoạn 2026 - 2030, việc mở rộng mô hình hóa bồi lắng 2 chiều trong lòng hồ và cập nhật các kịch bản phát thải mới nhất cần được tiếp tục triển khai. Các đơn vị quản lý vận hành và cơ quan ban ngành địa phương hãy chủ động ứng dụng ngay kết quả nghiên cứu vào thực tiễn quy hoạch để bảo đảm an ninh năng lượng và phát triển bền vững lưu vực sông Nậm Mức.