Tổng quan nghiên cứu

Lưu vực sông Đà là phụ lưu lớn nhất trong hệ thống sông Hồng với tổng diện tích tự nhiên đạt 52.900 km², trong đó phần lưu vực nằm trên lãnh thổ Việt Nam trải dài 570 km qua các tỉnh Tây Bắc. Mặc dù chỉ chiếm 37% tổng diện tích lưu vực sông Hồng tính đến trạm thủy văn Sơn Tây, nhưng lượng dòng chảy của sông Đà đóng góp tới 48% tổng lượng nước toàn hệ thống, tương đương 55,6 tỷ m³ dòng chảy năm tại vị trí tuyến đập Hòa Bình. Chế độ thủy văn tại đây biến động vô cùng phức tạp, trong đó lượng dòng chảy tập trung chủ yếu vào mùa lũ từ tháng 5 đến tháng 10, chiếm từ 69% đến 78% tổng lượng dòng chảy cả năm.

Vấn đề cốt lõi của công tác quản lý tài nguyên nước trên lưu vực là việc giải quyết mâu thuẫn giữa nhiệm vụ tích trữ nước phát điện và nhiệm vụ phòng chống lũ bảo vệ đồng bằng Bắc Bộ. Hệ thống bậc thang thủy điện trên sông Đà bao gồm các công trình trọng điểm quốc gia như hồ Hòa Bình với công suất 1.920 MW, dung tích toàn bộ 9.450 triệu m³ và hồ Sơn La với công suất 2.400 MW, dung tích 9.260 triệu m³. Để đảm bảo an toàn đập và kiểm soát dung tích phòng lũ 3.000 triệu m³ tại hồ Hòa Bình, việc nâng cao độ chính xác trong dự báo dòng chảy lũ ngắn hạn (24h đến 48h) và mở rộng (72h đến 120h) là yêu cầu cấp thiết.

Mục tiêu cụ thể của nghiên cứu là thiết lập chuỗi mô hình toán học tích hợp giữa mô hình thủy văn mưa - dòng chảy HEC-HMS và mô hình điều tiết liên hồ chứa HEC-ResSim, nhằm dự báo chính xác diễn biến lưu lượng lũ đổ về các hồ chứa lớn. Nghiên cứu mang ý nghĩa chiến lược trong việc cung cấp công cụ hỗ trợ ra quyết định thời gian thực, giảm thiểu rủi ro ngập lụt hạ du và tối ưu hóa hiệu quả kinh tế năng lượng.

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 phương trình cân bằng nước tổng quát trong thời kỳ lũ:

$$Y = X - Z - P$$

Trong đó $Y$ biểu thị độ sâu dòng chảy sinh ra từ lượng mưa $X$, sau khi đã trừ đi tổn thất bốc hơi $Z$ và tổn thất do thấm $P$. Để mô phỏng quá trình tổn thất, nghiên cứu áp dụng phương pháp tốc độ thấm ban đầu và thấm ổn định (Initial and Constant Rate), phản ánh trạng thái bão hòa ẩm của tầng đất mặt đệm lưu vực trước mỗi đợt mưa lớn.

Quá trình chuyển đổi mưa hiệu quả thành đường quá trình dòng chảy mặt được mô hình hóa bằng lý thuyết đường lũ đơn vị tổng hợp Snyder (Snyder Unit Hydrograph). Mô hình này xác định các đặc trưng then chốt gồm thời gian trễ lưu vực ($t_p$) và mô đun lưu lượng đỉnh lũ ($q_p$) thông qua mối quan hệ hình thái lòng dẫn và khoảng cách trọng tâm lưu vực ($L \cdot L_c$). Đối với quá trình diễn toán dòng chảy trong sông và điều tiết hồ chứa, nghiên cứu kết hợp phương pháp lượng trữ thủy văn Muskingum dựa trên dạng rút gọn của hệ phương trình Saint-Venant, cùng thuật toán mô phỏng vận hành liên hồ chứa đa mục tiêu trong HEC-ResSim theo ràng buộc mực nước dâng bình thường và mực nước gia cường.

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

Nguồn số liệu nghiên cứu được thu thập đồng bộ từ mạng lưới quan trắc khí tượng thủy văn trên toàn lưu vực sông Đà, bao gồm số liệu mưa thực đo từ hàng chục trạm đo bề mặt và số liệu lưu lượng tại 6 trạm thủy văn cơ bản: Mường Tè, Nậm Giàng, Mường Lay, Tà Gia, Nậm Mức và Sơn La. Bộ dữ liệu chuỗi thời gian được phân chia thành hai giai đoạn rõ ràng: chuỗi số liệu các trận lũ năm 2011 và 2013 được sử dụng cho công tác hiệu chỉnh thông số mô hình; số liệu trận lũ năm 2014 phục vụ công tác kiểm định độc lập.

Phương pháp mô phỏng số được lựa chọn thay thế cho các phương pháp phân tích thống kê cổ điển nhờ khả năng phản ánh cơ chế vật lý của từng tiểu lưu vực bộ phận. Địa hình lưu vực được số hóa tự động thông qua công cụ HEC-GeoHMS từ mô hình số hóa độ cao DEM. Việc tích hợp HEC-HMS và HEC-ResSim cho phép tái lập chính xác các kịch bản vận hành xả lũ liên hoàn theo thời gian thực tại các bậc thang thủy điện Lai Châu, Bản Chát, Sơn La và Hòa Bình.

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

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

Quá trình hiệu chỉnh mô hình HEC-HMS cho các trận lũ năm 2011 và 2013 tại 6 trạm thủy văn trọng điểm cho thấy mức độ tương thích cao giữa đường quá trình lưu lượng tính toán và thực đo. Hệ số hiệu quả Nash-Sutcliffe tại các trạm kiểm soát dòng chính như Mường Tè, Mường Lay và Sơn La đều đạt giá trị ấn tượng từ 0,82 đến 0,91, vượt xa ngưỡng yêu cầu của quy chuẩn dự báo thủy văn.

Sai số lưu lượng đỉnh lũ ($\Delta Q_{max}$) trong giai đoạn kiểm định năm 2014 được kiểm soát chặt chẽ dưới 8,5%, và sai số thời gian xuất hiện đỉnh lũ ($\Delta t_{max}$) dao động trong khoảng từ 1 đến 3 giờ. Thử nghiệm dự báo tác nghiệp trận lũ năm 2015 đổ về hồ thủy điện Sơn La và Hòa Bình chứng minh hiệu quả vượt trội của mô hình ở các thời đoạn dự báo 24h, 48h và 72h. Tại thời đoạn 24h, sai số lưu lượng dòng chảy đến hồ Sơn La đạt mức dưới 6,2%, đảm bảo tiêu chuẩn chất lượng phương án dự báo loại tốt theo Quy phạm dự báo thủy văn hiện hành (QP 94/TCN-91).

Thảo luận kết quả

Khi biểu diễn trên đồ thị thủy văn, các đường quá trình lưu lượng dự báo tái hiện chính xác hình thái đỉnh lũ dạng dốc đứng và chu kỳ lên xuống đặc trưng của vùng núi cao Tây Bắc. Kết quả tính toán chỉ ra rằng các nhánh sông sườn phía tây dãy Hoàng Liên Sơn như Nậm Na, Nậm Mu tạo ra mô đun dòng chảy đỉnh lũ rất cao, đạt trên 2.000 l/s/km², tạo nên áp lực dòng chảy gia nhập cực lớn vào lòng hồ Sơn La trong thời gian ngắn.

Dữ liệu mô phỏng từ HEC-ResSim thể hiện rõ rệt tác động điều tiết cắt lũ của hồ Bản Chát và Sơn La đối với lưu lượng xả về hạ lưu hồ Hòa Bình. Sự kết hợp giữa mô hình thủy văn và mô hình điều tiết giúp triệt tiêu hiện tượng trễ sóng lũ và sai lệch biên dòng chảy mà các mô hình hồi quy truyền thống thường mắc phải. Nghiên cứu chứng minh rằng khi kết hợp cập nhật dòng chảy thực đo tại thời điểm phát báo, độ tin cậy của bản tin dự báo 48h được cải thiện thêm 18,4% so với việc chỉ sử dụng dữ liệu mưa dự báo đơn thuần.

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

Nhằm nâng cao năng lực ứng phó thiên tai và hiện đại hóa công tác vận hành hệ thống hồ chứa lưu vực sông Đà, nghiên cứu đề xuất 4 nhóm giải pháp trọng tâm:

  1. Nâng cấp mạng lưới quan trắc khí tượng thủy văn tự động: Cần bổ sung tối thiểu 20 trạm đo mưa tự động và 5 trạm đo mực nước tự động tại các khu vực sườn núi hiểm trở thuộc thượng nguồn sông Nậm Na và Nậm Pô, nhằm tăng mật độ phủ dữ liệu mặt đệm lên 95% trước năm 2027, do Trung tâm Dự báo Khí tượng Thủy văn Quốc gia phối hợp với các địa phương triển khai.
  2. Tích hợp mô hình dự báo mưa số trị độ phân giải cao: Tiến hành liên kết trực tiếp dữ liệu mưa dự báo số trị từ mô hình WRF và HRM vào mô hình HEC-HMS để nâng thời gian dự kiến dự báo dòng chảy lũ từ 48h lên 120h với độ tin cậy đạt trên 80% trong giai đoạn 2026-2028, dưới sự chủ trì của Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu.
  3. Chuẩn hóa cơ sở dữ liệu số hóa địa hình lòng dẫn và hồ chứa: Cập nhật định kỳ đường đặc tính dung tích lòng hồ $V(Z)$ và mô hình mặt cắt ngang lòng sông chi tiết trên toàn tuyến 570 km hạ lưu các công trình đập để giảm sai số tính toán thủy lực xuống dưới 5% trước năm 2029, do Tập đoàn Điện lực Việt Nam (EVN) làm đầu mối thực hiện.
  4. Tự động hóa hệ sinh thái hỗ trợ điều hành liên hồ chứa thời gian thực: Thiết lập nền tảng tính toán tự động kết nối mô hình HEC-ResSim với hệ thống giám sát SCADA của các nhà máy thủy điện Hòa Bình, Sơn La, Lai Châu, Bản Chát, nhằm tối ưu hóa việc phân bổ dung tích phòng lũ 3.000 triệu m³ theo đúng Quyết định 1662/QĐ-TTg, do Ban Chỉ đạo Quốc gia về Phòng chống thiên tai điều phối hoàn thành vào năm 2027.

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

Tài liệu luận văn mang lại giá trị học thuật và thực tiễn phong phú cho các nhóm độc giả chuyên ngành:

  • Kỹ sư vận hành công trình thủy điện và cán bộ điều độ hệ thống điện: Cung cấp tài liệu tham khảo chuyên sâu để thiết lập các kịch bản vận hành xả lũ khẩn cấp, cân đối lưu lượng phát điện cực đại qua tuabin và quản lý an toàn dung tích hồ chứa.
  • Cơ quan quản lý nhà nước về phòng chống thiên tai và tài nguyên nước: Hỗ trợ khung công cụ khoa học để xây dựng bản đồ cảnh báo ngập lụt, phân vùng rủi ro thiên tai và ban hành các lệnh vận hành liên hồ chứa chính xác trong mùa mưa bão.
  • Nhà nghiên cứu, giảng viên và sinh viên ngành Kỹ thuật Thủy văn, Thủy lực: Cung cấp phương pháp luận chuẩn mực về việc ghép nối bộ phần mềm HEC (HEC-GeoHMS, HEC-HMS, HEC-ResSim, HEC-DSS) trên các lưu vực sông lớn có cấu trúc phân nhánh phức tạp.
  • Các tổ chức tư vấn thiết kế công trình thủy lợi, thủy điện: Khai thác bộ thông số thủy văn - địa hình đã được kiểm chứng thực tế để ứng dụng trong công tác tính toán lũ thiết kế, lũ kiểm tra và đánh giá an toàn đập.

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

Tại sao mô hình HEC-HMS lại được chọn để tính toán mưa - dòng chảy trên lưu vực sông Đà? Mô hình HEC-HMS có khả năng phân chia lưu vực rộng lớn 52.900 km² thành nhiều tiểu lưu vực hình nhánh cây và áp dụng linh hoạt phương pháp đường đơn vị Snyder kết hợp tổn thất Initial and Constant, phù hợp tối ưu với địa hình chia cắt mạnh tại vùng núi phía Bắc Việt Nam.

Vai trò chính của mô hình HEC-ResSim trong nghiên cứu này là gì? HEC-ResSim đóng vai trò mô phỏng chính xác quy trình vận hành liên hồ chứa bậc thang trên lưu vực sông Đà gồm Lai Châu, Bản Chát, Sơn La và Hòa Bình, xử lý các ràng buộc phức tạp về dung tích hữu ích, mực nước dâng bình thường và quy định xả lũ khẩn cấp.

Chất lượng dự báo dòng chảy lũ 24h và 48h đến các hồ chứa đạt mức độ tin cậy ra sao? Kết quả kiểm định cho thấy hệ số Nash-Sutcliffe đạt trên 0,85 và sai số lưu lượng đỉnh lũ duy trì dưới mức 10%, đáp ứng đầy đủ tiêu chuẩn dự báo thủy văn nghiệp vụ loại tốt theo quy chuẩn kỹ thuật quốc gia.

Quy trình vận hành liên hồ chứa theo Quyết định 1662/QĐ-TTg được tích hợp như thế nào? Các quy định về duy trì mực nước đón lũ, thời điểm xả lũ qua cửa xả đáy/mặt và giới hạn lưu lượng xả của hồ Sơn La, Hòa Bình theo Quyết định 1662/QĐ-TTg được số hóa trực tiếp thành các quy tắc vận hành (rules) trong mô hình HEC-ResSim.

Yếu tố nào gây ảnh hưởng lớn nhất đến sai số của mô hình dự báo mở rộng đến 120h? Sai số dự báo hạn dài từ 72h đến 120h chủ yếu phụ thuộc vào chất lượng và độ phân giải không gian của trường mưa dự báo số trị đầu vào, cũng như sự biến động đột ngột của các tổ hợp hình thế synop gây mưa lớn thượng nguồn.

Kết luận

  • Nghiên cứu đã ứng dụng thành công chuỗi mô hình toán học HEC-HMS và HEC-ResSim để mô phỏng và dự báo toàn diện quá trình hình thành dòng chảy lũ trên lưu vực sông Đà với diện tích 52.900 km².
  • Xác lập được bộ thông số thủy văn tối ưu cho các tiểu lưu vực thông qua số liệu thực đo các trận lũ lớn năm 2011, 2013 và kiểm định thành công với số liệu lũ năm 2014.
  • Đánh giá và lượng hóa chính xác diễn biến dòng chảy lũ gia nhập vào hệ thống bậc thang hồ chứa lớn gồm Lai Châu, Bản Chát, Sơn La và Hòa Bình ở các thời đoạn dự báo từ 24h đến 120h.
  • Đóng góp cơ sở khoa học vững chắc giúp nâng cao hiệu quả vận hành liên hồ chứa đa mục tiêu theo Quyết định 1662/QĐ-TTg, đảm bảo an toàn đập và giảm thiểu nguy cơ ngập lụt cho đồng bằng sông Hồng.
  • Lộ trình phát triển tiếp theo giai đoạn 2026-2030 tập trung vào việc tự động hóa cập nhật số liệu vệ tinh viễn thám, radar và tích hợp trí tuệ nhân tạo để tối ưu hóa thời gian dự báo nghiệp vụ.

Các cơ quan quản lý và đơn vị vận hành nhà máy thủy điện có thể trực tiếp ứng dụng kết quả và bộ thông số của luận văn nhằm phục vụ hiệu quả công tác vận hành an toàn hệ thống hồ chứa trong các mùa lũ tiếp theo.