Tổng quan nghiên cứu

Lưu vực sông Vu Gia – Thu Bồn có diện tích tự nhiên 10.350 km2, nằm trên địa bàn tỉnh Quảng Nam, thành phố Đà Nẵng và một phần tỉnh Kon Tum. Đây là khu vực có mạng lưới sông ngòi dày đặc, địa hình đồi núi dốc và chịu tác động mạnh mẽ của các hình thái thời tiết cực đoan. Chế độ dòng chảy trên lưu vực phân hóa sâu sắc theo mùa, trong đó lượng dòng chảy mùa lũ chiếm từ 60% đến 70% tổng lượng dòng chảy cả năm, tập trung cao điểm vào các tháng 10 và tháng 11 với lượng nước chiếm 25% đến 30% lượng dòng chảy năm. Các trận mưa lũ lịch sử đã gây ngập lụt nghiêm trọng vùng đồng bằng hạ du từ 2 m đến hơn 4 m, làm tê liệt hạ tầng giao thông và đe dọa trực tiếp đến đời sống của hơn 200.000 hộ dân cùng 10.000 ha đất sản xuất nông nghiệp.

Mặc dù quy hoạch phát triển thủy điện ban đầu xác định tổng dung tích hữu ích là 2.857 triệu m3 với dung tích phòng lũ 1.187 triệu m3 (chiếm 42%), nhưng thực tế triển khai các công trình ưu tiên phát điện khiến dung tích phòng lũ thực tế bị thu hẹp chỉ còn 146,44 triệu m3, chiếm khoảng 10% tổng dung tích hiệu dụng 1.453,24 triệu m3. Sự thiếu hụt dung tích phòng lũ cùng với các bất cập trong quy trình vận hành đơn hồ làm gia tăng xung đột giữa an toàn hạ du và hiệu quả kinh tế phát điện. Trước thực trạng đó, nghiên cứu được thực hiện nhằm xây dựng phương án dự báo lũ thời gian thực trước 24 giờ phục vụ vận hành tối ưu hệ thống hồ chứa thượng lưu, nâng cao năng lực cắt giảm lũ từ 15% đến 25% cho các trọng điểm hạ du như Ái Nghĩa và Câu Lâu.

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 lý thuyết mô hình hóa thủy văn mưa – dòng chảy tập trung và lý thuyết thủy lực dòng chảy một chiều không ổn định trong môi trường sông ngòi tự nhiên. Hai mô hình cốt lõi được tích hợp bao gồm:

  • Mô hình thủy văn MIKE-NAM: Mô phỏng quá trình hình thành dòng chảy từ mưa trên lưu vực đồi núi thông qua cơ chế cân bằng ẩm liên tục giữa 4 tầng trữ nước (bề mặt phủ, vùng rễ cây, tầng ngầm nông và tầng ngầm sâu).
  • Mô hình thủy lực 1 chiều MIKE 11: Giải hệ phương trình Saint-Venant phi tuyến tính bằng phương pháp sai phân hữu hạn ẩn để diễn toán sóng lũ trên mạng lưới sông phức tạp, có xét đến tác động lan truyền triều từ biển với biên độ dao động từ 80 cm đến 160 cm.

Các khái niệm chuyên ngành then chốt bao gồm: cường suất lũ (tốc độ dâng mực nước đạt từ 70 cm/giờ đến 140 cm/giờ ở thượng lưu), thời gian truyền lũ (khoảng 3,8 giờ cho đoạn 19 km từ Sơn Tân đến Nông Sơn và 15 giờ cho đoạn 23 km từ Ái Nghĩa đến Cẩm Lệ), và dung tích đón lũ điều tiết phòng cắt lũ hạ du.

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

Nguồn dữ liệu và cỡ mẫu nghiên cứu được thu thập từ chuỗi quan trắc khí tượng thủy văn nhiều năm giai đoạn 1964 – 2011 tại 11 trạm khí tượng đo mưa (Thành Mỹ, Hội An, Ái Nghĩa, Nông Sơn, Khâm Đức, Trà My...) và hệ thống trạm thủy văn cấp I, cấp II trên toàn lưu vực. Cỡ mẫu phân tích chuyên sâu lựa chọn 5 trận lũ lớn điển hình (các trận lũ tháng 11/1998, 11/1999, 10/2008, 9/2009 và 11/2011) với hàng nghìn điểm đo mực nước và lưu lượng theo từng giờ.

Phương pháp chọn mẫu có chủ đích được áp dụng nhằm lựa chọn những trận lũ cực đoan có đỉnh lũ vượt báo động 3, đồng thời có đầy đủ chuỗi số liệu mưa số trị và thông số vận hành hồ chứa. Lý do lựa chọn tổ hợp phương pháp phân tích MIKE-NAM và MIKE 11 là bởi khả năng liên kết động giữa dòng chảy từ các sườn dốc thượng nguồn với mạng lưới kênh rạch chằng chịt hạ lưu. Timeline nghiên cứu được triển khai đồng bộ từ khâu chuẩn hóa dữ liệu, hiệu chỉnh bộ thông số đến kiểm định và chạy dự báo thử nghiệm trong giai đoạn 2013 – 2015.

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

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

Thứ nhất, nghiên cứu làm rõ tính chất cực đoan của chế độ mưa sinh lũ trên lưu vực. Trận lũ tháng 11/1999 ghi nhận lượng mưa ngày kỷ lục lên tới 666,6 mm tại Hội An và 621,6 mm tại Thành Mỹ. Lưu lượng đỉnh lũ lớn nhất đạt 10.600 m3/s tại Nông Sơn và 7.000 m3/s tại Thành Mỹ, chỉ đứng sau trận lũ lịch sử năm 1964 với lưu lượng ước tính 18.000 m3/s. Tốc độ truyền lũ ở thượng lưu đạt từ 5 km/h đến 7 km/h, nhưng khi lan truyền về hạ du từ Ái Nghĩa đến Cẩm Lệ thì tốc độ giảm hơn 50% do lòng sông mở rộng và chịu ảnh hưởng của thủy triều.

Thứ hai, kết quả hiệu chỉnh và kiểm định mô hình MIKE-NAM cho các tiểu lưu vực thượng lưu đạt độ tin cậy rất cao. Hệ số hiệu quả Nash-Sutcliffe tại hai trạm khống chế chính là Nông Sơn và Thành Mỹ đạt từ 0,85 đến 0,92; sai số lưu lượng đỉnh lũ dao động trong khoảng 5% đến 8%, và sai số tổng lượng dòng chảy toàn trận lũ luôn dưới 7%.

Thứ ba, mô hình thủy lực MIKE 11 mô phỏng mực nước ngập lụt tại các trạm khống chế hạ du Ái Nghĩa và Câu Lâu cho kết quả sát thực tế với sai số đỉnh lũ chỉ từ 0,08 m đến 0,25 m, hệ số tương quan đạt trên 0,90 trong các đợt lũ kiểm định 2008, 2009 và 2011.

Thứ tư, phương án dự báo dòng chảy đến hồ trước 24 giờ kết hợp xử lý mưa số trị cho phép 5 hồ chứa lớn (A Vương, Sông Tranh 2, Đak Mi 4, Sông Bung 2, Sông Bung 4) chủ động xả đón lũ từ 50 đến 80 triệu m3, giúp giảm mực nước ngập lụt hạ du từ 0,11 m đến 0,61 m tại Hội Khách và Ái Nghĩa.

Thảo luận kết quả

Nguyên nhân chính dẫn đến ngập lụt hạ lưu là sự kết hợp đồng thời giữa mưa cường suất lớn ở vùng núi cao và dòng chảy phân lưu phức tạp qua sông Quảng Huế. Khi không có công cụ dự báo sớm, các nhà máy thủy điện thường phải xả tràn bị động khi hồ đã đầy, tạo đỉnh lũ nhân tạo gây sốc cho hạ du với biên độ ngập tăng thêm từ 1,5 m đến 4 m. So với các lưu vực sông Hồng hay sông Ba, lưu vực Vu Gia – Thu Bồn có độ dốc địa hình lớn hơn và thời gian tập trung nước ngắn hơn từ 30% đến 40%, đòi hỏi mô hình dự báo phải có độ trễ tính toán cực thấp.

Để trực quan hóa toàn bộ phát hiện này, các kết quả nghiên cứu có thể được thể hiện qua biểu đồ đường quá trình lưu lượng thực đo so với dự báo theo từng giờ tại Nông Sơn và Thành Mỹ nhằm phản ánh chính xác hình thái lũ đơn hoặc lũ kép. Đồng thời, bảng tổng hợp so sánh sai số mực nước cực đại và bảng cân bằng dung tích điều tiết giữa 5 hồ chứa sẽ minh chứng rõ ràng năng lực cắt giảm lũ của phương án đề xuất.

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

Thứ nhất, lắp đặt và nâng cấp mạng lưới quan trắc khí tượng thủy văn tự động. Cần bổ sung 15 trạm đo mưa chuyên dùng và 5 trạm quan trắc lưu lượng tự động trên các tiểu lưu vực nhánh thượng nguồn nhằm giảm sai số đo đạc mưa xuống dưới 10%. Giải pháp này do Trung tâm Dự báo Khí tượng Thủy văn Quốc gia chủ trì phối hợp cùng Ủy ban nhân dân tỉnh Quảng Nam và thành phố Đà Nẵng triển khai trong giai đoạn 2026 – 2028.

Thứ hai, thiết lập hệ thống phần mềm hỗ trợ ra quyết định vận hành liên hồ chứa theo thời gian thực. Tự động hóa việc tiếp nhận dữ liệu mưa số trị viễn thám cập nhật mỗi 3 giờ để chạy dự báo dòng chảy về hồ trước 24 đến 48 giờ với độ chính xác đạt trên 85%. Tập đoàn Điện lực Việt Nam và các công ty thủy điện trên lưu vực chịu trách nhiệm vận hành trong giai đoạn 2027 – 2029.

Thứ ba, điều chỉnh linh hoạt quy trình vận hành liên hồ chứa mùa lũ. Cho phép chủ động hạ thấp mực nước hồ để tạo dung tích đón lũ tăng thêm từ 15% đến 20% khi xuất hiện bản tin dự báo mưa diện rộng trên 300 mm/24h, bảo đảm mục tiêu cắt giảm đỉnh lũ hạ du từ 0,3 m đến 0,5 m. Bộ Nông nghiệp và Phát triển nông thôn chủ trì phối hợp với Bộ Công Thương ban hành quy trình sửa đổi trong giai đoạn 2026 – 2027.

Thứ tư, nạo vét và chỉnh trị lòng dẫn phân lưu sông Quảng Huế. Gia cố bờ kè và nạo vét luồng tiêu thoát lũ nhằm nâng cao khả năng thoát nước thêm 20% đến 25%, giảm thiểu nguy cơ ngập úng kéo dài cho hơn 200.000 hộ dân hạ du. Ban Chỉ huy Phòng chống thiên tai các địa phương tổ chức thi công liên tục trong giai đoạn 2026 – 2030.

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

  • Cơ quan quản lý nhà nước và Ban Chỉ huy Phòng chống thiên tai: Sử dụng luận văn làm cơ sở khoa học để thẩm định, xây dựng và phê duyệt quy trình vận hành liên hồ chứa liên tỉnh, nâng cao hiệu quả phòng tránh ngập lụt vùng ven biển miền Trung.
  • Các công ty quản lý và vận hành nhà máy thủy điện: Ứng dụng thuật toán dự báo dòng chảy lũ đến hồ trước 24 giờ để xây dựng kế hoạch xả đón lũ và tích nước tối ưu, giúp gia tăng từ 5% đến 10% sản lượng điện cuối vụ mà vẫn bảo đảm an toàn đập.
  • Kỹ sư thủy văn, chuyên gia tài nguyên nước và môi trường: Tham khảo phương pháp hiệu chỉnh bộ thông số mô hình MIKE-NAM và liên kết thủy lực MIKE 11 trong điều kiện thiếu dữ liệu địa hình chi tiết trên các lưu vực sông dốc ngắn.
  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Thủy văn: Khai thác tài liệu như một công trình mẫu mực về ứng dụng công nghệ mô hình số trị, xử lý dữ liệu khí tượng viễn thám và phân tích diễn toán sóng lũ thời gian thực.

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

Tại sao công tác dự báo lũ trên lưu vực sông Vu Gia – Thu Bồn lại gặp nhiều thách thức? Địa hình lưu vực có độ dốc rất lớn với độ chênh cao đạt từ 400 m đến 600 m trên mỗi kilômét chiều dài sông thượng nguồn, khiến lũ tập trung rất nhanh và thời gian truyền lũ về hạ du chỉ dưới 10 giờ. Thêm vào đó, mạng lưới trạm đo mưa vùng núi cao còn thưa thớt, gây khó khăn cho việc nắm bắt kịp thời các trận mưa đặc biệt lớn vượt 600 mm/ngày.

Mô hình liên kết MIKE-NAM và MIKE 11 mang lại ưu điểm gì vượt trội? Mô hình MIKE-NAM có khả năng mô phỏng nhanh quá trình mưa - dòng chảy trên các sườn dốc để xác định chính xác lưu lượng nước đổ về các hồ chứa. Sau đó, mô hình thủy lực MIKE 11 tiếp nhận dòng chảy này để tính toán chi tiết mực nước ngập lụt tại hạ du Ái Nghĩa và Câu Lâu với sai số mực nước chỉ từ 0,08 m đến 0,25 m.

Dung tích phòng lũ thực tế của các hồ chứa hiện nay đáp ứng ra sao so với quy hoạch? Theo quy hoạch ban đầu năm 2000, tổng dung tích phòng lũ dự kiến đạt 1.187 triệu m3, chiếm 42% dung tích hữu ích. Tuy nhiên, hiện tại chỉ có 5 hồ tham gia cắt lũ với tổng dung tích phòng lũ thực tế là 146,44 triệu m3, chiếm 10% tổng dung tích hiệu dụng, đòi hỏi phải có công cụ dự báo sớm 24 giờ để điều tiết linh hoạt.

Dữ liệu mưa dự báo số trị được xử lý như thế nào trước khi đưa vào mô hình? Số liệu mưa từ các mô hình dự báo thời tiết số trị được phân tích, hiệu chỉnh trọng số không gian qua phương pháp đa giác Thiessen và so sánh tương quan với chuỗi mưa thực đo tại 11 trạm khí tượng. Quy trình này giúp loại bỏ sai số cục bộ, cung cấp dữ liệu đầu vào chuẩn xác cho mô hình thủy văn mô phỏng dòng chảy.

Hiệu quả cắt giảm lũ hạ du khi vận hành hồ chứa theo dự báo sớm đạt mức nào? Dự báo dòng chảy về hồ trước 24 giờ giúp các chủ đập hạ thấp mực nước hồ đón lũ đúng thời điểm. Thử nghiệm trên các trận lũ năm 2008 và 2009 cho thấy phương án này làm giảm mực nước đỉnh lũ hạ lưu tại trạm Hội Khách tới 0,61 m và tại Ái Nghĩa tới 0,11 m so với phương thức vận hành thụ động trước đây.

Kết luận

  • Luận văn hệ thống hóa toàn diện các đặc trưng khí tượng, thủy văn và quy luật hình thành lũ cực đoan trên lưu vực 10.350 km2 của sông Vu Gia – Thu Bồn.
  • Thiết lập thành công chuỗi mô hình kết hợp MIKE-NAM và MIKE 11 với hệ số Nash-Sutcliffe đạt từ 0,85 đến 0,92, mô phỏng chuẩn xác dòng chảy lũ đến hồ và diễn toán ngập lụt hạ du.
  • Đề xuất quy trình dự báo lũ sớm trước 24 giờ từ số liệu mưa số trị, giải quyết thỏa đáng bài toán cân bằng giữa phát điện thương mại và cắt giảm lũ bảo vệ hạ du.
  • Chứng minh hiệu quả giảm mực nước đỉnh lũ hạ lưu từ 0,11 m đến 0,61 m thông qua điều tiết tối ưu 5 hồ chứa lớn gồm A Vương, Sông Tranh 2, Đak Mi 4, Sông Bung 2 và Sông Bung 4.
  • Cung cấp luận cứ khoa học thực tiễn vững chắc phục vụ các bộ, ngành và địa phương hoàn thiện quy chế vận hành liên hồ chứa theo thời gian thực.

Luận văn là tài liệu tham khảo giá trị cho công tác quản lý rủi ro thiên tai và quy hoạch tài nguyên nước. Trong giai đoạn 2026 – 2030, việc mở rộng tích hợp trí tuệ nhân tạo và dữ liệu radar thời tiết sẽ tiếp tục tối ưu hóa độ chính xác của mô hình. Các cơ quan quản lý và đơn vị vận hành hãy nghiên cứu, áp dụng ngay giải pháp công nghệ này để nâng cao an toàn cho cộng đồng và thúc đẩy phát triển bền vững lưu vực sông.