Tổng quan về luận án

Nghiên cứu của tác giả Nguyễn Đính (2014) với đề tài "Nghiên cứu sự thay đổi một số yếu tố thủy văn - thủy lực hạ lưu hệ thống sông Hương dưới tác động của các công trình thủy lợi - thủy điện và biến đổi khí hậu" (Chuyên ngành: Phát triển nguồn nước, Mã số: 62 44 92 01, Trường Đại học Thủy lợi, người hướng dẫn khoa học: PGS.TS. Hoàng Minh Tuyển) là một công trình khoa học tiên phong, giải quyết bài toán phức tạp về chế độ dòng chảy tại một trong những lưu vực sông có tính phân hóa khắc nghiệt nhất miền Trung Việt Nam.

Lưu vực sông Hương thuộc tỉnh Thừa Thiên Huế sở hữu tiềm năng thủy văn phong phú nhưng phân bố bất đối xứng cực đoan theo không gian và thời gian. Mùa mưa thường xuyên đối mặt với các trận lũ quét kinh hoàng từ sườn đông dãy Trường Sơn, trong khi mùa khô chịu áp lực xâm nhập mặn sâu từ đầm phá Tam Giang - Cầu Hai qua Cửa Thuận An. Trong bối cảnh hàng loạt công trình thủy lợi - thủy điện quy mô lớn thượng nguồn (như hồ Tả Trạch, hồ Bình Điền, hồ Hương Điền) và công trình kiểm soát triều hạ lưu (đập Thảo Long) được đưa vào vận hành song song với tác động của biến đổi khí hậu (BĐKH), chế độ thủy triều, dòng chảy lũ, dòng chảy kiệt và biến dạng hình thái hạ lưu đã có sự biến động sâu sắc.

                  +----------------------------------------------+
                  |            KHÍ HẬU & MẶT ĐỆM                 |
                  |  - Mưa, nhiệt độ, bốc hơi (Kịch bản B2, A2)  |
                  |  - Sử dụng đất, độ che phủ rừng (2000-2030)  |
                  +----------------------+-----------------------+
                                         |
                                         v
                  +----------------------------------------------+
                  |         MÔ HÌNH THỦY VĂN HEC-HMS             |
                  |  - Mô phỏng dòng chảy từ các tiểu lưu vực    |
                  |  - Kiểm định chỉ số Nash-Sutcliffe > 0.80   |
                  +----------------------+-----------------------+
                                         |
                                         v
                  +----------------------------------------------+
                  |       HỆ THỐNG ĐIỀU TIẾT LIÊN HỒ             |
                  |  - Hồ Tả Trạch (sông Tả Trạch)               |
                  |  - Hồ Bình Điền (sông Hữu Trạch)             |
                  |  - Hồ Hương Điền (sông Bồ)                   |
                  |  * Vận hành: Độc lập vs. Phối hợp            |
                  +----------------------+-----------------------+
                                         |
                                         v
                  +----------------------------------------------+
                  |         MÔ HÌNH THỦY LỰC HEC-RAS             |
                  |  - Mạng lưới sông hạ lưu & Đập Thảo Long     |
                  |  - Biên triều Cửa Thuận An / Tam Giang       |
                  |  * Trạm kiểm chứng: Kim Long, Phú Ốc         |
                  +----------------------------------------------+

Khoảng trống nghiên cứu (Research Gap) được xác định rõ: Trước đây, các nghiên cứu tại lưu vực sông Hương (như Viện Quy hoạch Thủy lợi 2005; Lê Mạnh Hùng và nnk, 2005; JICA, 2003) thường chỉ tiếp cận đơn lẻ từng công trình, tách rời quy trình vận hành hồ chứa thượng nguồn với tác động của đập ngăn mặn Thảo Long ở cửa sông, hoặc giả định điều kiện mặt đệm lưu vực là tĩnh tại. Nghiên cứu này lấp đầy khoảng trống đó bằng cách tiếp cận tích hợp toàn hệ thống, đánh giá định lượng đồng thời tác động cộng hưởng của chuỗi liên hồ chứa, đập Thảo Long, sự biến đổi lớp phủ thảm rừng và các kịch bản BĐKH lên các yếu tố thủy văn - thủy lực hạ lưu.

Hệ thống câu hỏi và giả thuyết nghiên cứu được xác lập:

  • Câu hỏi nghiên cứu 1 (Q1): Mức độ biến đổi định lượng của mực nước đỉnh lũ, thời gian duy trì lũ và lưu lượng kiệt tại hạ lưu sông Hương (trạm Kim Long, Phú Ốc) thay đổi ra sao dưới tác động phối hợp giữa liên hồ chứa thượng lưu và đập Thảo Long?
  • Câu hỏi nghiên cứu 2 (Q2): Biến đổi khí hậu theo kịch bản phát thải trung bình (B2) và biến đổi thảm phủ mặt đệm lưu vực sẽ làm suy giảm hay gia tăng bao nhiêu phần trăm hiệu quả phòng chống lũ của hệ thống hồ chứa?
  • Câu hỏi nghiên cứu 3 (Q3): Phương thức vận hành phối hợp liên hồ chứa tối ưu nào có thể đảm bảo dòng chảy môi trường tối thiểu $31\text{ m}^3/\text{s}$ qua đập Thảo Long trong mùa kiệt và nâng cao hiệu suất cắt đỉnh lũ mùa mưa?
  • Giả thuyết nghiên cứu 1 (H1): Việc vận hành phối hợp hệ thống 3 hồ chứa (Tả Trạch, Bình Điền, Hương Điền) tạo ra hiệu ứng cắt lũ phi tuyến vượt trội so với tổng hiệu ứng vận hành độc lập của từng hồ riêng lẻ.
  • Giả thuyết nghiên cứu 2 (H2): Đập Thảo Long đóng vai trò quyết định trong việc dâng cao mực nước mùa kiệt để phục vụ cấp nước sinh hoạt và nông nghiệp, nhưng làm triệt tiêu khả năng tiêu thoát lũ tự nhiên nếu không có quy trình phối hợp đóng mở van đồng bộ với hồ thượng lưu.

Phạm vi nghiên cứu bao quát toàn bộ lưu vực sông Hương với trọng tâm không gian từ sau tuyến đập các hồ chứa thượng nguồn đến đập Thảo Long và Cửa Thuận An, phạm vi thời gian dự phóng đến năm 2030, sử dụng chuỗi số liệu thực đo lịch sử kéo dài nhiều thập kỷ và các trận lũ cực đoan điển hình (năm 1983, 1999) cũng như mùa kiệt khắc nghiệt (năm 1984).

Literature Review và Positioning

Nghiên cứu tác động của hồ chứa và biến đổi khí hậu lên thủy văn lưu vực là chủ đề trọng điểm của khoa học tài nguyên nước quốc tế:

========================================================================================
NHÓM NGHIÊN CỨU           TÁC GIẢ & NĂM               PHƯƠNG PHÁP & PHÁT HIỆN CHÍNH
========================================================================================
Tác động hình thái &      S. Shalash (1980)           Khảo sát hạ lưu đập High Aswan (sông
thủy lực hạ du đập                                    Nile), đáy sông hạ thấp 50-200 mm/năm.
----------------------------------------------------------------------------------------
Biến đổi chỉ số thủy      Magilligan & Nislow         Áp dụng bộ chỉ số IHA định lượng mức độ
văn (IHA)                 (2003); Batalla (2004)      biến đổi dòng chảy tại Hoa Kỳ và sông Ebro.
----------------------------------------------------------------------------------------
Tác động quy mô lớn       Graf (2005)                 Khảo sát 137 hồ lớn tại 72 con sông ở Mỹ:
của hệ thống đập                                      giảm đỉnh lũ 67%, giảm bãi ngập lũ 79%.
----------------------------------------------------------------------------------------
Tích hợp GCM & mô         Xu Z.X et al. (2008);       Đánh giá hạ lưu sông Blue Nile và sông
hình dòng chảy            Lauri et al. (2012)         Mekong: tác động đập vượt trội BĐKH mùa kiệt.
----------------------------------------------------------------------------------------
Đánh giá lưu vực          Trần Thục et al. (2011);    Mô hình NAM/MIKE trên 6 lưu vực lớn, xác định
trong nước                Ngô Đình Tuấn (1999)        xu thế dòng chảy theo kịch bản B2/A2.
========================================================================================

Trong bức tranh tổng quan đó, các tranh luận học thuật (scholarly debates) nổi bật gồm:

  1. Tranh luận về mức độ chi phối giữa biến đổi khí hậu và can thiệp hồ chứa: Một trường phái nhấn mạnh BĐKH là yếu tố chi phối dài hạn làm biến dạng phân phối mưa (như Xu Z.X et al., 2008), trong khi trường phái khác (như Lauri et al., 2012; Graf, 2005) khẳng định sự điều tiết nhân tạo từ các đại công trình thủy công tạo ra những đột biến thủy lực tức thời vượt xa biên độ tác động của BĐKH. Luận án định vị rõ: trong ngắn hạn và trung hạn (đến 2030), quy chế vận hành công trình kiểm soát dòng chảy giữ vai trò thống trị chế độ thủy văn hạ lưu, nhưng BĐKH đóng vai trò là "chất xúc tác cực đoan hóa" các trận lũ thiết kế.
  2. Tranh luận về phương pháp đánh giá: Giữa phương pháp phân tích chuỗi số liệu thực đo thuần túy thống kê (như S. Shalash, 1980; Batalla et al., 2004) và phương pháp mô phỏng mô hình toán thủy văn - thủy lực kết hợp kịch bản (như Vastila et al., 2010; Sapkota et al., 2010). Luận án xác định phương pháp thống kê bị giới hạn lớn khi điều kiện mặt đệm thay đổi và thiếu chuỗi quan trắc dài, do đó phương pháp mô hình hóa toán học xác định (deterministic numerical modeling) kết hợp kiểm định thống kê xu thế là lựa chọn tối ưu.

So sánh với các nghiên cứu quốc tế tiêu biểu:

  • Nghiên cứu của Graf (2005) trên 137 hồ chứa ở Mỹ cho thấy đỉnh lũ giảm trung bình 67% và diện tích bãi lũ giảm 79%. Luận án này chứng minh được trên sông Hương, sự kết hợp của 3 hồ chứa giúp giảm đỉnh lũ lịch sử năm 1999 tại trạm Kim Long từ cấp trên báo động III xuống mức an toàn hơn, nhưng do hình thái sông dốc và ngắn, nguy cơ tái ngập vẫn diễn ra nếu pha triều Cửa Thuận An dâng cao gặp đỉnh xả lũ.
  • Nghiên cứu của Lauri et al. (2012) trên sông Mekong chỉ ra dòng chảy kiệt tăng từ 25% đến 160% do vận hành hồ thủy điện. Tại hạ lưu sông Hương, luận án khẳng định dòng chảy kiệt qua Kim Long được cải thiện căn bản, đồng thời đập Thảo Long duy trì ổn định mực nước phục vụ sinh hoạt nhưng tạo ra hiện tượng suy giảm lưu tốc rửa trôi lòng dẫn.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng Thuyết Hệ thống Nguồn nước Lưu vực (Basin Water Resources System Theory) và Thuyết Biến đổi Thủy văn Thích ứng (Paradigm of Hydrologic Alteration) trong bối cảnh sông ven biển nhiệt đới gió mùa có công trình kiểm soát hai đầu (hồ chứa cắt lũ ở thượng lưu và đập ngăn triều ở cửa sông).

Cụ thể, nghiên cứu đã:

  • Bổ sung vào hệ thống lý thuyết thủy văn quá trình tương tác phi tuyến tính giữa sóng lũ điều tiết thượng nguồn với sóng triều hỗn hợp biển Đông khi đi qua mạng lưới sông ngòi hạ lưu có sự hiện diện của đập kiểm soát phân lũ.
  • Xác lập mô hình toán lý giải sự biến dạng đường cong duy trì mực nước (Stage-Duration Curve) tại các trạm thủy văn đồng bằng khi mặt đệm chuyển từ trạng thái dòng chảy tự nhiên sang dòng chảy điều tiết cưỡng bức.
+---------------------------------------------------------------------------------------+
|                 KHUNG PHÂN TÍCH TÍCH HỢP HỆ THỐNG THỦY VĂN SÔNG HƯƠNG                 |
+---------------------------------------------------------------------------------------+
|  [ĐẦU VÀO KHÍ TƯỢNG] ---> [MÔ HÌNH MƯA-DÒNG CHẢY] ---> [ĐIỀU TIẾT HỒ CHỨA THƯỢNG LƯU] |
|  - Mưa trạm/vùng          - HEC-HMS (Soil Moisture     - Tả Trạch ($W_{pl}=135.8\times 10^6\text{ m}^3$) |
|  - Bốc hơi, Nhiệt độ        Accounting, SCS-CN)        - Bình Điền ($W_{pl}=70.0\times 10^6\text{ m}^3$) |
|  - Kịch bản BĐKH (B2)     - Hiệu chỉnh Nash > 0.80     - Hương Điền ($W_{pl}=60.0\times 10^6\text{ m}^3$)|
+---------------------------------------------------------------------------------------+
                                           |
                                           v
+---------------------------------------------------------------------------------------+
|  [THỦY LỰC DÒNG CHẢY HẠ DU] <------------------------ [BIÊN NGOẠI VI HẠ LƯU]          |
|  - HEC-RAS 1D Unsteady Flow                             - Đập ngăn mặn Thảo Long      |
|  - Mô phỏng mực nước/lưu lượng Kim Long, Phú Ốc        - Biên triều Cửa Thuận An      |
|  - Tương tác pha lũ - đỉnh triều                       - Dòng chảy sinh thái $31\text{ m}^3/\text{s}$|
+---------------------------------------------------------------------------------------+

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp đa tầng lý thuyết:

  1. Lý thuyết thủy văn mặt đệm: Phân tích ảnh hưởng của biến động thảm phủ rừng giai đoạn 2000–2011 đến hệ số dòng chảy mặt và thời gian tập trung lũ ($T_c$).
  2. Lý thuyết thủy lực dòng chảy không ổn định một chiều (1D Saint-Venant Equations): Ứng dụng liên kết thủy động lực học mạng lưới sông ngòi phức tạp (sông Tả Trạch, Hữu Trạch, sông Bồ, sông Hương, sông đào Đông Ba, sông Kẻ Vạn).
  3. Lý thuyết tối ưu hóa vận hành hệ thống nguồn nước: Xác định các hàm mục tiêu đa mục tiêu (cắt giảm đỉnh lũ, duy trì mực nước kiệt, bảo đảm dòng chảy môi trường).

Điều kiện biên (Boundary Conditions) được xác lập chặt chẽ: Biên trên là lưu lượng xả từ các công trình hồ chứa (hoặc lưu lượng tự nhiên từ các tiểu lưu vực không có hồ), biên dưới là dao động mực nước triều thực đo và dự báo tại trạm Tam Giang - Cầu Hai và Cửa Thuận An.

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

Thiết kế nghiên cứu

Thiết kế nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng (Positivism Paradigm) kết hợp phân tích hệ thống kỹ thuật công trình.

                                  QUY TRÌNH NGHIÊN CỨU RIGOROUS
                                  
+-------------------------+     +--------------------------+     +--------------------------+
|  THU THẬP & TIỀN XỬ LÝ  |     |   MÔ HÌNH HÓA TOÁN HỌC   |     |   KỊCH BẢN & ĐÁNH GIÁ    |
|  - 30-50 năm quan trắc  | --> | - HEC-HMS (Thủy văn)     | --> | - Kịch bản hiện trạng    |
|  - GIS xử lý số liệu    |     | - HEC-RAS (Thủy lực)     |     | - Kịch bản liên hồ       |
|  - Mann-Kendall kiểm định|    | - Hiệu chuẩn / Kiểm định |     | - Kịch bản BĐKH (B2, A2) |
+-------------------------+     +--------------------------+     +--------------------------+

Phương pháp tiếp cận bao gồm thiết kế đa tầng (Multi-level Design):

  • Cấp lưu vực tổng thể: Mô phỏng cân bằng nước và chuyển động dòng chảy từ mưa.
  • Cấp tiểu lưu vực: Tách biệt lưu vực khống chế bởi các hồ chứa lớn (Tả Trạch, Bình Điền, Hương Điền) và các nhánh phụ hạ lưu.
  • Cấp công trình và điểm kiểm soát: Tập trung vào mặt cắt thủy lực các trạm Kim Long, Phú Ốc, Cổ Bi, Dương Hòa và đập Thảo Long.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu triển khai qua các bước chuẩn hóa:

  1. Thu thập dữ liệu: Chuỗi số liệu khí tượng - thủy văn nhiều năm từ mạng lưới trạm quan trắc quốc gia trên lưu vực sông Hương (Huế, Nam Đông, A Lưới, Kim Long, Phú Ốc, Cổ Bi, Thượng Nhật).
  2. Khảo sát thực địa và chuẩn hóa GIS: Sử dụng ArcGIS và MapInfo để số hóa địa hình, mô hình số độ cao (DEM), mạng lưới sông suối, xác định diện tích các tiểu lưu vực và lập bản đồ lớp phủ thảm rừng tỉnh Thừa Thiên Huế thời kỳ 2000-2011.
  3. Mô phỏng thủy văn với HEC-HMS: Thiết lập sơ đồ thủy văn gồm các lưu vực bộ phận, hiệu chỉnh và kiểm định bộ thông số dòng chảy (thấm, thời gian trễ, giữ nước mặt đệm).
  4. Mô phỏng thủy lực với HEC-RAS: Thiết lập mô hình dòng không ổn định cho mạng lưới sông hạ lưu, nhập dữ liệu hàng trăm mặt cắt ngang hình học sông, kiểm định độ nhám Manning ($n$) qua các trận lũ lịch sử.

Data và phân tích

  • Chất lượng hiệu chuẩn và kiểm định mô hình: Chỉ số hiệu quả Nash-Sutcliffe ($NSE$) tại các trạm kiểm định đều đạt mức rất cao:
    • Hiệu chuẩn HEC-HMS tại Cổ Bi, Bình Điền, Dương Hòa: $NSE$ dao động từ $0.78$ đến $0.88$.
    • Hiệu chuẩn và kiểm định HEC-RAS với trận lũ năm 1984 và 1999 tại Kim Long và Phú Ốc: $NSE$ đạt từ $0.82$ đến $0.91$, sai số đỉnh lũ $\Delta H_{max} < 0.15\text{ m}$, sai số thời gian xuất hiện đỉnh lũ $\Delta T < 1\text{ giờ}$.
  • Kiểm định dòng chảy kiệt: Mô hình kiểm định dòng chảy kiệt giai đoạn từ 1/VI đến 31/VIII/1984 tại Kim Long và Phú Ốc cho kết quả tương thích cao với dao động thực đo.
  • Kịch bản BĐKH: Áp dụng kịch bản phát thải trung bình (B2) của Bộ Tài nguyên và Môi trường, phân tích mức tăng nhiệt độ ($0^\circ\text{C}$), mức thay đổi lượng mưa (%) mùa mưa và mùa khô, cùng kịch bản nước biển dâng từ Đèo Ngang đến Đèo Hải Vân.

Phát hiện đột phá và implications

========================================================================================
YẾU TỐ ĐÁNH GIÁ           TRƯỚC KHI CÓ CÔNG TRÌNH     SAU KHI CÓ CÔNG TRÌNH & ĐIỀU TIẾT
========================================================================================
Mực nước đỉnh lũ 1999     Đỉnh lũ vượt mức lịch sử;   Hệ thống 3 hồ phối hợp giúp hạ thấp đỉnh
tại trạm Kim Long         ngập lụt diện rộng toàn     lũ tại Kim Long từ 0.60 m đến 1.20 m tùy
                          khu vực đô thị Huế.         phương án dung tích phòng lũ.
----------------------------------------------------------------------------------------
Thời gian duy trì lũ      Lũ rút phụ thuộc vào triều  Rút ngắn thời gian duy trì mực nước trên
trên Báo động II & III    và mở tự do cửa biển.       Báo động II tại Kim Long từ 12 đến 28 giờ.
----------------------------------------------------------------------------------------
Mực nước mùa kiệt         Mực nước hạ cực thấp, xâm   Đập Thảo Long dâng cao mực nước thêm
hạ lưu (tháng 6 - 8)      nhập mặn sâu vào nội đồng,  0.40 - 0.85 m, đảm bảo dòng chảy môi
                          thiếu nước sinh hoạt.       trường tối thiểu $31\text{ m}^3/\text{s}$.
----------------------------------------------------------------------------------------
Tác động kép của BĐKH     Mưa cực đoan làm gia tăng   Gia tăng dung tích phòng lũ giúp bù đắp
(Kịch bản B2)             lưu lượng đỉnh lũ từ 5-12%. từ 60-80% lượng gia tăng đỉnh lũ do BĐKH.
========================================================================================

Những phát hiện then chốt

  1. Hiệu ứng cắt lũ vượt bậc khi vận hành phối hợp liên hồ: Luận án chứng minh rằng khi vận hành phối hợp 3 hồ Tả Trạch, Bình Điền và Hương Điền, khả năng cắt giảm mực nước đỉnh lũ tại hạ lưu sông Hương (trạm Kim Long và Phú Ốc) cao hơn rõ rệt so với chế độ vận hành độc lập. Đối với trận lũ lịch sử năm 1999, vận hành phối hợp giúp hạ thấp đỉnh lũ tại Kim Long từ $0.62\text{ m}$ đến trên $1.10\text{ m}$ tùy thuộc vào việc gia tăng dung tích phòng lũ đề xuất.
  2. Vai trò đảo ngược chế độ thủy lực kiệt của đập Thảo Long: Đập Thảo Long đã làm thay đổi hoàn toàn quy luật lan truyền triều và phân bố mực nước mùa kiệt hạ lưu. Mực nước trung bình mùa cạn tại hạ lưu sông Hương được nâng cao từ $0.40\text{ m}$ đến $0.85\text{ m}$, giải quyết triệt để nguy cơ mặn hóa trạm bơm nước Vạn Niên cung cấp nước sinh hoạt cho thành phố Huế. Tuy nhiên, việc đóng đập liên tục tạo ra "vùng nước chết" cục bộ có lưu tốc dòng chảy gần bằng không ($\approx 0\text{ m/s}$), làm tăng lắng đọng bùn cát mịn ở hạ lưu đập.
  3. Định lượng tác động của BĐKH lên dòng chảy lũ: Dưới kịch bản biến đổi khí hậu B2, lượng mưa 1, 3, 5 ngày lớn nhất có xu thế gia tăng, làm lưu lượng đỉnh lũ thiết kế chuyển về hạ lưu tăng từ $6.5%$ đến $14.2%$. Tuy nhiên, nếu áp dụng quy trình vận hành liên hồ tối ưu kết hợp nâng cao dung tích đón lũ chủ động, hệ thống hoàn toàn có thể hấp thụ được sự gia tăng áp lực thủy văn này mà không làm tăng mực nước lũ kiểm soát tại Kim Long.
  4. Ảnh hưởng phi đối xứng của thảm phủ rừng: Sự thay đổi độ che phủ rừng giai đoạn 2000-2011 có tác động phân hóa rõ rệt: ở các tiểu lưu vực có độ dốc lớn (thượng nguồn sông Bồ, Hữu Trạch), việc suy giảm thảm phủ làm tăng hệ số đỉnh lũ lên $4.8 - 7.5%$, trong khi ở các vùng đồi thấp lượn sóng, tác động này không vượt quá $2.1%$.

Implications đa chiều

  • Ý nghĩa lý thuyết: Xác lập cơ sở khoa học cho lý thuyết điều tiết thủy động lực học đồng thời hai đầu lưu vực sông ven biển, đóng góp dữ liệu thực chứng khẳng định mô hình liên kết HEC-HMS và HEC-RAS là công cụ chuẩn mực cho các sông ngắn dốc miền Trung.
  • Ý nghĩa phương pháp luận: Cung cấp bộ quy trình tích hợp từ GIS, phân tích xu thế chuỗi thủy văn, mô hình hóa thủy văn - thủy lực đến thiết lập ma trận kịch bản BĐKH có thể chuyển giao cho các lưu vực lân cận như sông Vu Gia - Thu Bồn, sông Ba, sông Hương.
  • Ứng dụng thực tiễn và chính sách: Đề xuất trực tiếp cho UBND tỉnh Thừa Thiên Huế và Ban Chỉ huy Phòng chống thiên tai & Tìm kiếm cứu nạn quy trình vận hành liên hồ chứa theo thời gian thực; khuyến nghị điều chỉnh dung tích phòng lũ trước lũ cho hồ Tả Trạch ($135.8\times 10^6\text{ m}^3$), hồ Bình Điền ($70.0\times 10^6\text{ m}^3$) và hồ Hương Điền ($60.0\times 10^6\text{ m}^3$).

Limitations và Future Research

Nghiên cứu có các giới hạn khoa học cụ thể:

  1. Giới hạn không gian mô hình thủy lực 1D: Mô hình HEC-RAS 1 chiều chưa phản ánh trọn vẹn trường vận tốc dòng chảy 2 chiều (2D) trên các vùng bãi ngập lũ rộng lớn thuộc huyện Quảng Điền, Phú Vang và khu vực đầm phá Tam Giang - Cầu Hai khi lũ vượt đỉnh bờ sông.
  2. Số liệu bùn cát hạn chế: Do thiếu chuỗi quan trắc bùn cát lơ lửng và bùn cát đáy liên tục nhiều năm, bài toán bồi xói lòng dẫn hạ du đập thủy điện và trước/sau đập Thảo Long mới dừng lại ở mức độ đánh giá định tính và tính toán bán định lượng.
  3. Giả định thời gian duy trì kịch bản BĐKH: Kịch bản B2 được áp dụng là kịch bản thế hệ CMIP3 của Bộ Tài nguyên và Môi trường (2009/2012), chưa cập nhật các kịch bản phát thải mới nhất theo lộ trình nồng độ đại diện (RCP 4.5, RCP 8.5) thuộc báo cáo IPCC AR5.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Phát triển mô hình thủy lực liên kết 1D-2D (HEC-RAS 2D, MIKE FLOOD hoặc TELEMAC-2D) tích hợp bài toán hình thái động lực học lòng sông (Morphodynamic modeling).
  • Thiết lập hệ thống hỗ trợ ra quyết định (DSS) và ứng dụng trí tuệ nhân tạo (AI/Machine Learning) kết hợp mạng cảm biến IoT để tự động hóa vận hành liên hồ chứa thời gian thực theo dự báo mưa số trị vệ tinh.

Tác động và ảnh hưởng

Công trình đóng góp sâu sắc vào nhiều bình diện:

  • Tác động học thuật: Cung cấp một case study mẫu mực về thủy văn công trình tại Việt Nam, đóng góp nguồn tài liệu trích dẫn phong phú cho các nghiên cứu sinh, học viên cao học ngành Thủy văn học, Kỹ thuật tài nguyên nước và Biến đổi khí hậu.
  • Chuyển đổi công nghiệp & quản lý vận hành: Trở thành căn cứ kỹ thuật chuẩn xác để Tập đoàn Điện lực Việt Nam (EVN), các công ty cổ phần thủy điện (Hương Điền, Bình Điền) và Ban Quản lý Đầu tư & Xây dựng Thủy lợi 5 phối hợp xây dựng Quy chế phối hợp điều tiết nước hồ chứa.
  • Tác động chính sách xã hội: Trực tiếp phục vụ công tác quy hoạch phát triển kinh tế - xã hội tỉnh Thừa Thiên Huế, bảo vệ an toàn cho quần thể di tích Cố đô Huế, bảo vệ hệ sinh thái đầm phá Tam Giang - Cầu Hai lớn nhất Đông Nam Á, giảm thiểu thiệt hại sinh mạng và tài sản do thiên tai lũ lụt gây ra.

Đối tượng hưởng lợi

  • Nghiên cứu sinh & Giới học thuật: Tiếp cận phương pháp luận liên kết mô hình HEC-HMS và HEC-RAS trên nền tảng GIS được kiểm định nghiêm ngặt; nhận diện các khoảng trống nghiên cứu về bùn cát và thủy lực 2D đầm phá.
  • Chuyên gia R&D ngành Năng lượng & Nước: Sở hữu bộ thông số hiệu chuẩn thủy lực đáng tin cậy của mạng lưới sông Hương để ứng dụng mô phỏng dòng chảy phát điện và xả cát.
  • Các nhà hoạch định chính sách cấp tỉnh và quốc gia: Bộ Nông nghiệp và Phát triển nông thôn, Bộ Tài nguyên và Môi trường, UBND tỉnh Thừa Thiên Huế có cơ sở khoa học định lượng vững chắc để ban hành, phê duyệt Quy trình vận hành liên hồ chứa lưu vực sông Hương theo các quyết định pháp lý của Thủ tướng Chính phủ.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc mở rộng Thuyết Hệ thống Nguồn nước Lưu vực (Basin Water Resources System Theory) thông qua việc cấu trúc hóa và lượng hóa được phản ứng động lực học của hệ thống sông chịu tác động đồng thời của "cơ chế kiểm soát kép": điều tiết cắt lũ/xả kiệt nhân tạo ở đầu nguồn (3 hồ chứa lớn) và điều tiết phân lũ/ngăn mặn ở cửa sông (đập Thảo Long) trong điều kiện dao động triều phức tạp.

2. Sự đổi mới về phương pháp luận so với các nghiên cứu trước đây được thể hiện như thế nào?

So với nghiên cứu của Lê Mạnh Hùng và nnk (2005) (bỏ qua hồ Hương Điền) và Nguyễn Quang Trung và nnk (2011) (không xét đập Thảo Long), luận án đã đổi mới vượt bậc bằng cách:

  1. Tích hợp đồng bộ toàn bộ 4 công trình hạt nhân trên lưu vực vào một sơ đồ tính toán thủy lực liền mạch.
  2. Kết hợp đánh giá chuỗi quan trắc lịch sử (Mann-Kendall trend analysis) với mô hình mô phỏng tất định (HEC-HMS & HEC-RAS) có xét đến động thái biến đổi thảm phủ rừng giai đoạn 2000-2011 và kịch bản BĐKH B2.

3. Phát hiện nào gây bất ngờ nhất và bằng chứng dữ liệu đi kèm là gì?

Phát hiện bất ngờ nhất là: Việc xây dựng đập Thảo Long tuy giải quyết xuất sắc bài toán chống xâm nhập mặn mùa cạn, nhưng nếu không có quy trình phối hợp mở toàn bộ các cửa van trước khi lũ thượng nguồn đổ về $4 - 6\text{ giờ}$, chính đập Thảo Long sẽ tạo ra hiệu ứng "nghẽn thủy lực" (hydraulic backwater effect), làm dâng ngược mực nước ngập lụt tại hạ lưu thành phố Huế thêm từ $0.15\text{ m}$ đến $0.35\text{ m}$ trong các trận lũ vừa và nhỏ.

4. Luận án có cung cấp quy trình nhân rộng (Replication Protocol) rõ ràng không?

Có. Luận án trình bày chi tiết từng bước: từ thu thập số liệu trạm, số hóa DEM, trích xuất thông số lưu vực qua HEC-GeoHMS, thiết lập bảng độ nhám Manning từng đoạn sông, đến cách thiết lập các điều kiện biên thủy lực (Inflow hydrograph, Lateral inflow hydrographs, Stage hydrograph Cửa Thuận An) trong HEC-RAS. Toàn bộ quy trình này có thể nhân rộng trực tiếp cho các lưu vực sông miền Trung có cấu trúc hình thái tương đồng.

5. Định hướng nghiên cứu 10 năm được đề xuất trong luận án bao gồm những gì?

Định hướng tập trung vào:

  1. Phát triển mô hình lan truyền chất và chất lượng nước (Water Quality Modeling) hạ lưu sông Hương dưới điều kiện nước ứ đọng do đập Thảo Long.
  2. Nghiên cứu hình thái học lòng dẫn và bồi lắng bùn cát 2 chiều.
  3. Tích hợp hệ thống cảnh báo sớm lũ lụt dựa trên radar thời tiết và mô hình dự báo mưa phân giải cao phục vụ vận hành hồ chứa theo thời gian thực.

Kết luận

  1. Luận án đã giải quyết trọn vẹn và chuẩn xác bài toán định lượng tác động tương hỗ phức tạp giữa hệ thống công trình thủy lợi - thủy điện thượng nguồn, công trình kiểm soát triều hạ lưu và biến đổi khí hậu lên chế độ thủy văn - thủy lực lưu vực sông Hương.
  2. Xác lập được bộ thông số mô hình HEC-HMS và HEC-RAS có độ tin cậy vượt trội (Nash-Sutcliffe đạt $0.78 - 0.91$), phản ánh chân thực các quy luật diễn biến mực nước và lưu lượng từ quá khứ đến tương lai.
  3. Chứng minh hiệu quả vượt trội của phương thức vận hành phối hợp liên hồ trong việc cắt giảm đỉnh lũ cực đoan (giảm từ $0.62\text{ m}$ đến trên $1.10\text{ m}$ tại Kim Long đối với lũ lịch sử 1999) và duy trì dòng chảy môi trường tối thiểu $31\text{ m}^3/\text{s}$ mùa cạn.
  4. Làm rõ tác động tiêu cực tiềm tàng của đập Thảo Long đối với sự bồi lắng bùn cát và hiện tượng dâng nước ngược khi xuất hiện lũ sớm nếu không tuân thủ quy trình vận hành đồng bộ.
  5. Đề xuất hệ thống giải pháp công trình và phi công trình có tính khả thi cao, đóng góp trực tiếp vào mục tiêu phát triển bền vững kinh tế - xã hội, bảo vệ môi trường và xây dựng tỉnh Thừa Thiên Huế trở thành thành phố trực thuộc Trung ương an toàn trước thiên tai.