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 lưu vực đạt 52.900 km² và chiều dài dòng chính khoảng 660 km (đoạn chảy qua lãnh thổ Việt Nam dài 570 km). Mặc dù diện tích chỉ chiếm 37% lưu vực sông Hồng tính đến trạm Sơn Tây, sông Đà lại đóng góp tới 48% tổng lượng dòng chảy toàn hệ thống, tương đương 55,6 tỷ m³ nước mỗi năm tại trạm Hòa Bình và 41,0 tỷ m³ tại trạm Sơn La. Dòng chảy tập trung chủ yếu vào mùa lũ kéo dài từ tháng 6 đến tháng 10, chiếm từ 69% đến 78,5% tổng lượng dòng chảy cả năm, với những đợt mưa lớn có thể gây ra lũ đặc biệt lớn vượt mức 14.500 m³/s.

Vấn đề nghiên cứu trọng tâm xuất phát từ mâu thuẫn vận hành thực tế giữa các mục tiêu tại hệ thống bậc thang thủy điện lớn nhất cả nước gồm Lai Châu, Bản Chát (220 MW), Sơn La (2.400 MW) và Hòa Bình (1.920 MW). Việc đảm bảo an toàn đập, cắt giảm lũ cho vùng đồng bằng Bắc Bộ nhưng đồng thời phải tích trữ nước tối ưu để phát điện và cấp nước mùa kiệt đòi hỏi các bản tin dự báo dòng chảy đến hồ phải đạt độ tin cậy cao trong khung thời gian từ 24 giờ đến 120 giờ.

Mục tiêu cụ thể của luận văn là nghiên cứu, thiết lập và tích hợp bộ công cụ mô hình toán HEC-HMS và HEC-ResSim nhằm mô phỏng chính xác quá trình mưa - dòng chảy và diễn toán điều tiết lũ trên lưu vực sông Đà đến tuyến đập Hòa Bình. Nghiên cứu tập trung giải quyết bài toán dự báo ngắn hạn và trung hạn phục vụ Quy trình vận hành liên hồ chứa theo Quyết định số 1662/QĐ-TTg, hướng tới nâng cao chỉ số Nash-Sutcliffe lên trên 0,80 và giảm sai số đỉnh lũ xuống dưới 12%, đáp ứng đầy đủ tiêu chuẩn nghiệp vụ QP 94/TCN-91.

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 hai hệ thống lý thuyết nền tảng trong thủy văn học và thủy lực công trình:

Thứ nhất, lý thuyết cân bằng nước và mô hình hóa mưa - dòng chảy mặt thông qua phương trình cân bằng ẩm: lớp dòng chảy hình thành bằng lượng mưa trừ đi tổn thất bốc hơi và tổn thất thấm. Tổn thất thấm được xác định bằng phương pháp thấm ban đầu và thấm ổn định kết hợp hệ số thấm không đổi và lượng trữ ban đầu của đất. Quá trình chuyển đổi mưa hiệu quả thành đường quá trình dòng chảy áp dụng lý thuyết Đường lũ đơn vị tổng hợp Snyder (1938), xác định lưu lượng đỉnh đơn vị và thời gian trễ lưu vực dựa trên các hệ số hình thái lòng dẫn. Diễn toán dòng chảy lũ dọc theo các đoạn sông chính sử dụng phương pháp Muskingum dựa trên mối quan hệ trữ lượng hình nêm và hình lăng trụ với trọng số hiệu chỉnh dao động từ 0 đến 0,3.

Thứ hai, lý thuyết điều tiết và phân phối dòng chảy liên hồ chứa được mô hình hóa trong HEC-ResSim. Cơ chế vận hành mô phỏng dựa trên phương trình cân bằng dung tích liên tục theo từng bước thời gian tính toán, tích hợp các đường đặc tính dung tích - mực nước, quy tắc đóng mở cửa xả đáy, cửa xả mặt và công suất lưu lượng qua tua-bin phát điện (tối đa 2.400 m³/s tại hồ Hòa Bình và 3.438 m³/s tại hồ Sơn La).

Các khái niệm chuyên ngành cốt lõi bao gồm: mô đun dòng chảy đỉnh lũ (đạt tới 2.147 l/s/km² tại Nậm Pô), dung tích phòng lũ (3.000 triệu m³ tại hồ Hòa Bình), mực nước dâng bình thường (115 m tại Hòa Bình, 215 m tại Sơn La, 475 m tại Bản Chát) và thời gian truyền sóng lũ.

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

Nguồn dữ liệu của luận văn được thu thập từ Trung tâm Dự báo Khí tượng Thủy văn Quốc gia và Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu. Dữ liệu bao gồm chuỗi lượng mưa quan trắc tại các trạm khí tượng trên toàn lưu vực, dữ liệu mực nước và lưu lượng thực đo 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, Sơn La), cùng dữ liệu nhật ký xả tràn và phát điện của các hồ chứa. Dữ liệu địa hình được xử lý từ mô hình số độ cao DEM qua công cụ HEC-GeoHMS để phân chia 6 tiểu lưu vực thành phần.

Cỡ mẫu nghiên cứu gồm 5 tập chuỗi số liệu trận lũ lớn điển hình: giai đoạn hiệu chỉnh mô hình sử dụng chuỗi lũ năm 2011 và 2013; giai đoạn kiểm định sử dụng trận lũ năm 2014; và giai đoạn thử nghiệm dự báo nghiệp vụ áp dụng cho trận lũ mùa mưa năm 2015 ở 5 mốc thời gian từ 24h, 48h, 72h, 96h đến 120h. Phương pháp chọn mẫu là chọn mẫu định lượng có chủ đích, ưu tiên các đợt mưa diện rộng kéo dài từ 2 đến 4 ngày với tổng lượng mưa vượt trên 200 mm để kiểm tra độ tin cậy trong các tình huống cực đoan.

Lý do lựa chọn phương pháp phân tích tích hợp HEC-HMS và HEC-ResSim: Mô hình HEC-HMS có cấu trúc dạng nhánh cây linh hoạt, thích hợp với địa hình phân cắt mạnh có độ cao trung bình 1.130 m của sông Đà. Mô hình HEC-ResSim cho phép lập trình chi tiết các ràng buộc vận hành phức tạp theo Quyết định 1662/QĐ-TTg, khắc phục hạn chế của các phương pháp thống kê truyền thống vốn chỉ đạt mức bảo đảm khoảng 65% đến 70%.

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

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

Quá trình tính toán và kiểm nghiệm mô hình đã mang lại bốn phát hiện khoa học quan trọng:

Thứ nhất, kết quả hiệu chỉnh mô hình HEC-HMS trên chuỗi số liệu lũ năm 2011 và 2013 tại 6 tiểu lưu vực đạt độ chính xác cao. Chỉ số hiệu quả Nash-Sutcliffe tại các trạm kiểm soát dòng chảy đạt từ 0,82 đến 0,91, sai số lưu lượng đỉnh lũ dao động trong biên độ cho phép từ 4,2% đến 11,8%, và sai số thời gian xuất hiện đỉnh lũ không vượt quá 3 giờ.

Thứ hai, trong giai đoạn kiểm định năm 2014, mô hình duy trì tính ổn định vượt trội. Hệ số Nash tại trạm Mường Tè đạt 0,86 và tại tuyến Sơn La đạt 0,89; sai số tổng lượng dòng chảy toàn đợt lũ được khống chế dưới 8,5%, chứng minh tính phù hợp của bộ thông số vật lý lưu vực.

Thứ ba, thử nghiệm dự báo dòng chảy lũ năm 2015 đến các hồ Lai Châu, Bản Chát, Sơn La và Hòa Bình khẳng định tính khả thi cao trong thực tế. Đối với thời hạn dự báo ngắn 24 giờ và 48 giờ, mức đảm bảo dự báo đạt trên 85%, sai số đỉnh lũ dưới 10%. Khi kéo dài thời hạn dự báo lên 72 giờ và 120 giờ, chỉ số Nash giảm dần từ 0,88 xuống 0,73 nhưng vẫn đáp ứng đầy đủ tiêu chuẩn nghiệp vụ dự báo thủy văn cấp 1 theo quy chuẩn QP 94/TCN-91.

Thứ tư, kết quả điều tiết của HEC-ResSim chứng minh năng lực cắt lũ vượt trội của hệ thống hồ chứa bậc thang. Khi xuất hiện đỉnh lũ trên 14.000 m³/s ở thượng nguồn, việc phối hợp vận hành dung tích phòng lũ giữa hồ Sơn La (dung tích hữu ích 6.504 triệu m³) và hồ Hòa Bình (dung tích phòng lũ 3.000 triệu m³) giúp làm giảm từ 28% đến 35% lưu lượng đỉnh lũ đổ về hạ du sông Hồng.

Thảo luận kết quả

Độ chính xác cao của kết quả mô phỏng bắt nguồn từ việc tính toán chi tiết các đặc trưng hình thái lưu vực sông Đà. Mật độ sông suối dày từ 1,5 đến 1,78 km/km² ở sườn tây Hoàng Liên Sơn tạo ra mô đun đỉnh lũ cục bộ rất lớn (trạm Nậm Chiến đạt 2.048 l/s/km², Suối Sập đạt 901 l/s/km²). Nhờ cơ chế chia nhỏ lưu vực trong HEC-HMS, dòng chảy gia nhập từ các phụ lưu không có trạm đo đã được tập hợp đầy đủ trước khi đổ vào dòng chính.

Về mặt trình bày dữ liệu, các kết quả mô phỏng được biểu diễn trực quan qua hệ thống biểu đồ so sánh đường quá trình lưu lượng thực đo và tính toán theo bước thời gian từ 1 giờ đến 6 giờ. Biểu đồ đường quá trình thể hiện sự trùng khớp chặt chẽ ở cả sườn lũ lên, đỉnh lũ và sườn lũ rút. Đồng thời, bảng thống kê phân phối sai số xác suất theo tiêu chuẩn QP 94/TCN-91 xác nhận tỷ lệ bản tin dự báo đạt mức tốt và xuất sắc chiếm trên 82% tổng số đợt phát báo hạn ngắn.

So với các nghiên cứu trước đây áp dụng mô hình TANK hoặc phương pháp hồi quy thống kê vốn chỉ đạt mức bảo đảm trung bình khoảng 65% và sai số đỉnh lũ thường vượt 20%, phương pháp tích hợp HEC-HMS và HEC-ResSim thể hiện sự vượt trội rõ rệt. Phương pháp này giải quyết triệt để sự tương tác thủy lực phức tạp giữa các hồ chứa nối tiếp trên cùng một lưu vực.

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

Dựa trên kết quả nghiên cứu, luận văn đưa ra 4 giải pháp trọng tâm nhằm nâng cao năng lực dự báo và quản lý thiên tai:

  1. Hiện đại hóa mạng lưới trạm quan trắc tự động: Lắp đặt bổ sung từ 12 đến 15 trạm đo mưa tự động và 4 trạm đo lưu lượng tại các khu vực thượng nguồn biên giới giáp Trung Quốc thuộc tỉnh Lai Châu và Điện Biên. Mục tiêu giảm sai số ước lượng mưa bình quân lưu vực xuống dưới 7%, hoàn thành trước tháng 5/2027. Đơn vị chủ trì: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia.

  2. Tích hợp mô hình dự báo mưa số trị độ phân giải cao: Ghép nối dữ liệu mưa thời gian thực từ các mô hình khí tượng số trị phân giải cao (WRF, HRM) vào mô hình HEC-HMS. Mục tiêu nâng cao độ chính xác của bản tin dự báo lũ hạn 72h đến 120h, duy trì hệ số Nash trên 0,78 trong giai đoạn 2026 - 2028. Đơn vị thực hiện: Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu.

  3. Chuẩn hóa quy trình điều hành liên hồ chứa thời gian thực: Cập nhật các bảng tra lưu lượng xả và quy tắc khống chế mực nước trước lũ (104 m tại Hòa Bình và 205 m tại Sơn La) vào mô hình HEC-ResSim. Hướng tới cắt giảm tối thiểu 30% lưu lượng đỉnh lũ khi xuất hiện lũ lớn trên 14.000 m³/s, triển khai từ mùa lũ năm 2027. Đơn vị thực hiện: Ban Chỉ đạo Quốc gia về Phòng chống Thiên tai phối hợp với Tập đoàn Điện lực Việt Nam (EVN).

  4. Đào tạo và chuyển giao công nghệ mô hình toán: Tổ chức định kỳ 4 khóa tập huấn chuyên sâu hàng năm về hiệu chỉnh thông số mô hình HEC-HMS và điều hành HEC-ResSim cho hơn 50 kỹ sư thủy văn tại các nhà máy thủy điện Hòa Bình, Sơn La, Lai Châu và Bản Chát, bắt đầu từ quý IV/2026. Đơn vị chủ trì: Trường Đại học Thủy lợi.

Đố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 khoa học và thực tiễn giá trị cho 4 nhóm đối tượng chính:

  1. Kỹ sư và cán bộ vận hành nhà máy thủy điện: Đội ngũ kỹ thuật tại Thủy điện Hòa Bình (1.920 MW), Sơn La (2.400 MW), Lai Châu (1.200 MW) và Bản Chát (220 MW) có thể áp dụng trực tiếp quy trình HEC-ResSim để lập kịch bản xả lũ và tối ưu hóa sản lượng điện hàng năm đạt trên 10.227 triệu kWh tại Sơn La và 8.160 triệu kWh tại Hòa Bình.

  2. Dự báo viên khí tượng thủy văn: Cán bộ dự báo tại Trung tâm Dự báo Khí tượng Thủy văn Quốc gia và Đài Khí tượng Thủy văn khu vực Tây Bắc sử dụng bộ thông số mô hình HEC-HMS đã được kiểm định (Nash đạt 0,82 - 0,91) để phục vụ công tác phát báo định kỳ từ 1 đến 5 ngày.

  3. Cơ quan quản lý nhà nước về phòng chống thiên tai: Ban Chỉ đạo Quốc gia về Phòng chống Thiên tai và Cục Quản lý Đê điều sử dụng các phân tích điều tiết dung tích phòng lũ 3.000 triệu m³ để hoàn thiện phương án ứng phó khẩn cấp và bảo vệ an toàn đồng bằng sông Hồng.

  4. Giảng viên, học viên cao học và sinh viên chuyên ngành: Các cơ sở nghiên cứu và đào tạo ngành Thủy văn, Kỹ thuật Tài nguyên nước sử dụng công trình làm tài liệu tham khảo tiêu biểu về ứng dụng mô hình toán thủy văn - thủy lực trong lưu vực sông lớn có điều tiết hồ chứa.

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

  1. Tại sao cần kết hợp đồng thời hai mô hình HEC-HMS và HEC-ResSim trên sông Đà?
    Lưu vực sông Đà có hệ thống bậc thang thủy điện lớn kiểm soát hơn 55,6 tỷ m³ dòng chảy hàng năm. Mô hình HEC-HMS đóng vai trò tính toán dòng chảy sinh ra từ mưa trên các lưu vực tự nhiên, trong khi HEC-ResSim tiếp nhận dòng chảy này để mô phỏng chính xác việc điều tiết qua các cửa xả của 4 hồ chứa theo đúng Quyết định 1662/QĐ-TTg.

  2. Chỉ số Nash-Sutcliffe đạt trên 0,85 phản ánh điều gì về độ tin cậy của mô hình?
    Chỉ số Nash đạt từ 0,82 đến 0,91 trong giai đoạn hiệu chỉnh (2011, 2013) và trên 0,85 khi kiểm định (2014) cho thấy đường quá trình dòng chảy tính toán gần như trùng khít với số liệu thực đo. Theo quy chuẩn thủy văn quốc tế, chỉ số vượt mức 0,75 được xếp loại rất tốt, khẳng định mô hình đủ tiêu chuẩn phục vụ dự báo tác nghiệp.

  3. Mức độ chính xác của dự báo hạn 24h và hạn 120h khác nhau như thế nào theo tiêu chuẩn QP 94/TCN-91?
    Kết quả thử nghiệm năm 2015 chỉ ra rằng dự báo hạn ngắn 24h - 48h đạt mức đảm bảo trên 85% với sai số đỉnh lũ dưới 10%. Đối với hạn dài 120h, mức đảm bảo đạt khoảng 72% đến 75%. Độ chính xác được duy trì nhờ cơ chế cập nhật liên tục lưu lượng thực đo tại các trạm biên như Mường Tè và Nậm Giàng.

  4. Dung tích phòng lũ của hồ Sơn La và Hòa Bình đóng vai trò gì trong việc giảm nhẹ thiên tai?
    Hồ Sơn La và Hòa Bình có tổng dung tích phòng lũ chuyên dụng lớn, riêng hồ Hòa Bình là 3.000 triệu m³. Khi xuất hiện các trận lũ thượng nguồn vượt 14.000 m³/s, việc phối hợp đóng mở cửa xả giúp cắt giảm từ 28% đến 35% lưu lượng đỉnh lũ đổ về hạ lưu, bảo vệ vững chắc đê điều vùng đồng bằng Bắc Bộ.

  5. Làm thế nào để xử lý số liệu mưa tại các tiểu lưu vực thiếu trạm đo ở vùng núi cao?
    Mô hình sử dụng phương pháp đa giác Thiessen và chia nhỏ 6 tiểu lưu vực kết hợp gán trọng số trạm mưa đại diện. Tại các khu vực sườn đón gió có lượng mưa lớn từ 300 mm đến 500 mm, việc chuẩn hóa hệ số dòng chảy và tốc độ thấm giúp mô hình tái hiện chính xác lượng dòng chảy mặt thực tế.

Kết luận

Luận văn thạc sĩ của tác giả Nguyễn Thị Thu Bình đã giải quyết thành công bài toán cấp bách về dự báo dòng chảy lũ phục vụ điều hành liên hồ chứa đa mục tiêu trên lưu vực sông Đà. Năm đóng góp nổi bật của công trình bao gồm:

  • Thiết lập thành công mô hình HEC-HMS cho 6 tiểu lưu vực sông Đà với độ chính xác cao (chỉ số Nash đạt 0,82 - 0,91, sai số đỉnh lũ dưới 12%).
  • Xây dựng mô hình điều tiết HEC-ResSim cho 4 hồ chứa lớn (Lai Châu, Bản Chát, Sơn La, Hòa Bình) bám sát các ràng buộc của Quyết định 1662/QĐ-TTg.
  • Thử nghiệm dự báo lũ thời gian thực hạn 24 giờ đến 120 giờ cho trận lũ năm 2015, đạt tiêu chuẩn nghiệp vụ QP 94/TCN-91 với mức đảm bảo trên 80%.
  • Định lượng hiệu quả cắt giảm từ 28% đến 35% lưu lượng đỉnh lũ cho vùng đồng bằng hạ du thông qua phối hợp vận hành dung tích liên hồ chứa.
  • Cung cấp cơ sở khoa học hoàn chỉnh để chuyển giao công nghệ dự báo số trị cho các đơn vị quản lý vận hành từ mùa lũ năm 2026.

Để đưa kết quả nghiên cứu vào ứng dụng thực tế trên diện rộng, các đơn vị quản lý cần khẩn trương thực hiện lộ trình tích hợp hệ thống mô hình vào quy trình phát báo tự động trong giai đoạn 2026 - 2027. Hãy liên hệ với nhóm nghiên cứu tại Trường Đại học Thủy lợi để tiếp nhận chuyển giao công nghệ và khai thác tối đa hiệu quả kiểm soát lũ thông minh trên lưu vực sông Đà.