Tổng quan về luận án

Nghiên cứu cơ sở khoa học điều hành liên hồ chứa chống lũ lưu vực sông Ba là công trình luận án tiến sĩ chuyên ngành Thủy văn học (mã số: 62440224) của nghiên cứu sinh Nguyễn Việt, được thực hiện dưới sự hướng dẫn khoa học của PGS. Nguyễn Hữu Khải và PGS. Nguyễn Tiền Giang tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội (2018). Luận án giải quyết bài toán cốt lõi và cấp bách trong quản lý tài nguyên nước: điều hành thời gian thực hệ thống liên hồ chứa đa mục tiêu trong điều kiện thủy văn bất định cao và mâu thuẫn gay gắt giữa nhiệm vụ phòng chống lũ hạ du với tối ưu hóa hiệu quả phát điện.

                  +------------------------------------------+
                  |     KHÍ TƯỢNG THỦY VĂN BẤT ĐỊNH          |
                  |  (Dự báo 24h: 75-80%; 48-72h: sai số lớn)|
                  +--------------------+---------------------+
                                       |
                                       v
         +------------------------------------------------------------+
         |     BỘ ĐIỀU KHIỂN LOGIC MỜ (FUZZY LOGIC CONTROLLER)        |
         |  - Fuzzification: Q_đến, H_hồ, H_khống chế hạ du           |
         |  - Inference: Mamdani Engine & Rule Base (IF-THEN)         |
         |  - Defuzzification: Trọng tâm (Center of Gravity - COG)    |
         +-----------------------------+------------------------------+
                                       |
                                       v
+-------------------------------------------------------------------------------+
|      MÔ HÌNH THỦY LỰC DIỄN TOÁN TRUYỀN LŨ (MUSKINGUM-CUNGE & MIKE 11)         |
|  - Cụm Thượng lưu (Tây Trường Sơn): Ka Nak (513m->506m), Ayun Hạ (203m->202m) |
|  - Cụm Hạ lưu (Đông Trường Sơn): Krông H'Năng, Sông Hinh, Sông Ba Hạ          |
+--------------------------------------+----------------------------------------+
                                       |
                                       v
+-------------------------------------------------------------------------------+
|                        KẾT QUẢ ĐIỀU TIẾT HẠ DU                                |
|  - Cắt giảm đỉnh lũ Qmax tại Củng Sơn: 20 - 25%                               |
|  - Hạ mực nước đỉnh lũ: Củng Sơn (0.80 - 1.50m), Phú Lâm (0.30 - 0.80m)       |
|  - Bảo toàn cột nước phát điện & chống ngập úng cục bộ không cần thiết        |
+-------------------------------------------------------------------------------+

Bối cảnh nghiên cứu xuất phát từ những bất cập thực tiễn khi vận hành hệ thống hồ chứa theo Quy trình 1757 (Quyết định 1757/QĐ-TTg năm 2010) và Quy trình 1077 (Quyết định 1077/QĐ-TTg năm 2014, sửa đổi bổ sung bởi Quyết định 282/QĐ-TTg năm 2017). Các quy trình hiện hành quy định quy tắc xả cứng nhắc: buộc hạ đồng loạt mực nước đón lũ (MNĐL) tại tất cả các hồ, kể cả khi mưa chỉ xuất hiện cục bộ ở một tiểu lưu vực. Điều này dẫn đến sự suy giảm nghiêm trọng cột nước phát điện và nguy cơ thiếu hụt dung tích tích nước cuối mùa lũ.

Khoảng trống nghiên cứu (research gap) lớn nhất nằm ở sự thiếu hụt phương pháp luận điều hành thích nghi khi thông tin dự báo thủy văn mang tính ngẫu nhiên và sai số cao. Theo ghi nhận trực tiếp từ luận án: "Thời gian dự kiến dự báo quá trình lũ cho khu vực Trung Bộ, Tây Nguyên đối với các sông lớn là 24 giờ với mức đảm bảo dự báo khoảng 75÷80%. Sai số còn lớn đối với dự báo trước 48÷72 giờ." Việc vận hành dựa trên các mô hình tất định truyền thống hoàn toàn bị động trước tính bất định này.

Luận án thiết lập 3 câu hỏi nghiên cứu và 3 giả thuyết khoa học:

  • $RQ_1$: Làm thế nào để mô hình hóa quá trình ra quyết định xả lũ thích nghi của con người dưới điều kiện thông tin dòng chảy đến không chắc chắn? $\to H_1$: Logic mờ (Fuzzy Logic) kết hợp hệ suy luận Mamdani có thể chuyển hóa tri thức chuyên gia thành bộ quy tắc điều khiển định lượng tối ưu.
  • $RQ_2$: Sự kết hợp giữa bộ điều khiển mờ và mô hình diễn toán thủy lực Muskingum-Cunge có phản ánh chính xác tương tác phức tạp của hệ thống hồ bậc thang và song song không? $\to H_2$: Mô hình tích hợp phản ánh trung thực quá trình truyền lũ và biến động mực nước hạ du theo thời gian thực.
  • $RQ_3$: Việc linh hoạt điều chỉnh quy tắc tích/xả đón lũ có nâng cao đồng thời hiệu quả giảm lũ và an ninh năng lượng không? $\to H_3$: Bỏ quy định xả đón lũ đồng loạt bắt buộc sẽ bảo toàn đầu nước phát điện mà vẫn đảm bảo an toàn cắt đỉnh lũ từ 20-25% tại trạm Củng Sơn.

Khung lý thuyết của luận án tích hợp Thuyết tập mờ (Fuzzy Set Theory) của L.A. Zadeh (1965), lý thuyết diễn toán thủy lực Muskingum-Cunge và lý thuyết phân tích hệ thống tài nguyên nước (Loucks & Labadie). Đối tượng và phạm vi nghiên cứu bao gồm lưu vực sông Ba với diện tích 13.900 $\text{km}^2$ cùng hệ thống 5 hồ chứa trọng điểm: Ka Nak – An Khê, Ayun Hạ, Krông H’Năng, Sông Ba Hạ và Sông Hinh. Dữ liệu thực nghiệm bao trùm các chuỗi quan trắc thủy văn lịch sử điển hình qua các trận lũ lớn năm 1981, 1988, 1993, 2009 và 2013 ứng với các tần suất $P = 5%$, $10%$, $20%$.

Literature Review và Positioning

Tổng quan nghiên cứu vận hành hệ thống hồ chứa trên thế giới và Việt Nam cho thấy sự tiến hóa qua bốn trường phái phương pháp luận chính:

Trường phái tối ưu hóa tất định và ngẫu nhiên (Deterministic & Stochastic Optimization): Bắt đầu từ Young (1967) với quy hoạch động Monte-Carlo, tiếp nối bởi Bhaskar & Whitlach (1980), Butcher, Loucks, Karamouz & Houck (1987) với quy hoạch động tất định kết hợp hồi quy tuyến tính nhiều biến (DPR), và quy hoạch động ngẫu nhiên (SDP). Các kỹ thuật này đòi hỏi hàm xác suất dòng chảy chính xác và đối mặt với rào cản "lời nguyền chiều không gian" (curse of dimensionality) khi áp dụng cho hệ thống đa hồ chứa phức tạp.

Trường phái thuật toán tiến hóa (Evolutionary Algorithms): East & Hall (1994), Fahmy, Oliveira & Loucks (1997), Sharif & Wardlaw (2000), Juran Ali Ahmed & Sarma (2005) đã ứng dụng thuật toán di truyền (Genetic Algorithm - GA) để tìm quỹ đạo vận hành tối ưu. Janga Reddy & D. Nagesh Kumar (2006) sử dụng thuật toán tiến hóa đa mục tiêu (MOEA) trên hồ Hirakud (Ấn Độ) để giải bài toán đánh đổi giữa năng lượng và rủi ro lũ. Tuy nhiên, hạn chế cố hữu của thuật toán tiến hóa là tính toán chậm, khó đáp ứng bài toán điều hành lũ theo thời gian thực (real-time).

Trường phái mô phỏng và mô phỏng kết hợp tối ưu (Simulation-Optimization): Các phần mềm mô phỏng tổng thể như HEC-5, HEC-RESSIM (Klipsch et al., 2002), MIKE BASIN, MIKE 11 (DHI) cho phép mô tả chi tiết hình học công trình nhưng thiếu cơ chế ra quyết định tự động linh hoạt. Nghiên cứu của Wurbs (1993) và Labadie (2004) chỉ ra rằng người vận hành thực tế luôn hoài nghi các kết quả tối ưu thuần túy do không phản ánh được yếu tố rủi ro và bất định.

Trường phái Logic mờ (Fuzzy Logic Paradigm): Khởi xướng từ Zadeh (1965), được Savic & Simonovic (1991), Russell & Campbell (1996), Shrestha et al. (1996), Cheng Chuntian (1999), Panigrahi (2001) ứng dụng để mô hình hóa quyết định xả dựa trên ngôn ngữ tự nhiên.

                                  TIẾN TRÌNH TIẾP CẬN LÝ THUYẾT VẬN HÀNH HỒ CHỨA
+-----------------------------+     +-----------------------------+     +-----------------------------+
|    TỐI ƯU HÓA KINH ĐIỂN     |     |      MÔ PHỎNG THỦY LỰC      |     |     ĐIỀU KHIỂN MỜ (FIS)     |
|   (DP, SDP, LP, GA, MOEA)   | --> |    (HEC-RESSIM, MIKE 11)    | --> |  (Mamdani, Zadeh, COG)      |
|  * Hạn chế: "Lời nguyền số  |     |  * Hạn chế: Thử dần thụ     |     |  * Ưu thế: Xử lý bất định,  |
|    chiều", chạy chậm        |     |    động, không tự thích nghi|     |    thời gian thực           |
+-----------------------------+     +-----------------------------+     +-----------------------------+
                                                                                       |
                                                                                       v
                                                                        +-----------------------------+
                                                                        |      VỊ TRÍ LUẬN ÁN         |
                                                                        | Tích hợp FIS + Muskingum-   |
                                                                        | Cunge + Topology Bậc thang/ |
                                                                        | Song song (LVS Sông Ba)     |
                                                                        +-----------------------------+

Về mặt tranh luận học thuật, tồn tại hai quan điểm đối nghịch rõ nét:

  • Quan điểm 1: Coi việc điều hành hồ chứa là bài toán tối ưu toán học chính xác tuyệt đối, nơi mọi biến thủy văn phải được lượng hóa thông qua dự báo tất định dài hạn (3–5 ngày).
  • Quan điểm 2 (Quan điểm của luận án): Cho rằng bản chất dòng chảy lũ là một quá trình ngẫu nhiên phi tuyến có sai số dự báo rất lớn; việc cố định hóa các ngưỡng xả chỉ làm tăng rủi ro ngập lụt giả tạo và tổn thất tài nguyên. Vận hành hồ chứa thực chất là hệ thống kiểm soát thích nghi dựa trên các tín hiệu trạng thái mờ.

So sánh với các nghiên cứu quốc tế điển hình:

  • Mô hình của Dubrovin et al. (2002) tại hồ Päijänne (Phần Lan) áp dụng hệ luật mờ theo mùa cho hồ đơn lẻ, trong khi luận án mở rộng cho hệ thống liên hồ hỗn hợp cả bậc thang và song song gồm 5 hồ lớn.
  • Nghiên cứu của Janga Reddy & Kumar (2006) tại hồ chứa Hirakud trên lưu vực sông Mahanadi (Ấn Độ) sử dụng kỹ thuật Folded Dynamic Programming và MOEA giải quyết bài toán quy hoạch dài hạn, trong khi luận án của Nguyễn Việt giải quyết trực tiếp bài toán điều tiết lũ ngắn hạn theo từng bước thời gian thực (time-step).

Luận án định vị chính xác khoảng trống: Là công trình tiên phong tại Việt Nam tích hợp hoàn chỉnh hệ thống suy luận mờ (FIS) với mô hình thủy lực diễn toán dòng chảy lũ Muskingum-Cunge trên toàn bộ mạng lưới phân bố của lưu vực sông Ba.

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

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

Luận án đóng góp vào việc mở rộng Lý thuyết tập mờ của L.A. Zadeh (1965) và Quy tắc suy diễn Modus Ponens tổng quát hóa (Generalized Modus Ponens) của E.H. Mamdani (1977) vào cơ học chất lỏng ứng dụng và thủy văn công trình.

Tập mờ $A$ được xác định trên không gian tham chiếu $X$ thông qua hàm liên thuộc $\mu_A(x): X \to [0, 1]$: $$A = \left{ \left( x, \mu_A(x) \right) \mid x \in X, \mu_A(x) \in [0, 1] \right}$$

Luận án chuẩn hóa các hàm liên thuộc đầu vào và đầu ra dưới hai dạng hình học chủ đạo: dạng hình thang $\text{trapmf}(x; a, b, c, d)$ và hình tam giác $\text{trimf}(x; a, b, c)$:

$$\mu_{\text{trap}}(x; a, b, c, d) = \max\left( 0, \min\left( \frac{x-a}{b-a}, 1, \frac{d-x}{d-c} \right) \right)$$

$$\mu_{\text{tri}}(x; a, b, c) = \max\left( 0, \min\left( \frac{x-a}{b-a}, \frac{c-x}{c-b} \right) \right)$$

Hệ suy luận mờ đa điều kiện với cấu trúc IF-THEN được xây dựng dựa trên phép hợp thành Max-Min (Zadeh): $$R_i: \text{IF } Q_{\text{đến}} \text{ is } A_i \text{ AND } H_{\text{hồ}} \text{ is } B_i \text{ AND } H_{\text{hạ du}} \text{ is } C_i \text{ THEN } Q_{\text{xả}} \text{ is } D_i$$

Giá trị điều khiển xả rõ ($Q_{\text{xả}}^$) được lượng hóa chính xác thông qua phương pháp giải mờ Trọng tâm (Center of Gravity - COG): $$Q_{\text{xả}}^ = \frac{\int z \cdot \mu_D(z) , dz}{\int \mu_D(z) , dz} \quad \text{hoặc} \quad Q_{\text{xả}}^* = \frac{\sum_{j=1}^n z_j \cdot \mu_D(z_j)}{\sum_{j=1}^n \mu_D(z_j)}$$

          HÀM LIÊN THUỘC (MF)                 HỆ LUẬT MỜ (MAMDANI)                GIẢI MỜ TRỌNG TÂM (COG)
   1.0 +     /-------\                  IF Q_đến = Rất Lớn                     Q_xả* = SUM(z*u) / SUM(u)
       |    /         \                 AND H_hồ = Cao                                   |
   0.5 |   /           \       ===>     AND H_AnKhê = Trung Bình       ===>              v
       |  /             \               THEN Q_xả = Cắt Lũ Tối Đa              [Lưu Lượng Xả Rõ m3/s]
   0.0 +-+---------------+--->
         a   b     c     d

Luận án xác lập 3 mệnh đề lý thuyết (propositions):

  • $P_1$: Quá trình điều tiết lũ của một hệ thống liên hồ chứa đa bậc là một hàm phi tuyến phụ thuộc đồng thời vào trạng thái dung tích tức thời, gradient gia tăng lũ đến và sức chứa thủy lực khả dụng tại các điểm khống chế hạ lưu.
  • $P_2$: Mức độ phụ thuộc $\mu(x)$ của các biến trạng thái thủy văn cho phép triệt tiêu hiện tượng dao động cưỡng bức (hunting effect) tại các cửa van xả tràn so với các quy tắc đóng/mở bậc thang cơ học truyền thống.
  • $P_3$: Cấu trúc không gian lưu vực chi phối trực tiếp đến thứ tự ưu tiên trong hệ luật mờ: cụm hồ bậc thang ưu tiên giải phóng dung tích đón lũ tuần tự từ hạ lưu lên thượng lưu, trong khi cụm hồ song song ưu tiên cân bằng động dung tích phòng lũ tỷ lệ thuận với mô đun dòng chảy đến.

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

Khung phân tích của luận án kết hợp ba trụ cột lý thuyết: (1) Thuyết tập mờ mô phỏng quyết định hành vi, (2) Thuyết diễn toán sóng lũ thủy lực Muskingum-Cunge, và (3) Lý thuyết phân tích cấu trúc hệ thống mạng lưới sông ngòi.

Về cấu trúc mạng lưới, luận án mô hình hóa lưu vực sông Ba thành hai phân hệ thủy lực đặc thù:

  • Phân hệ hồ bậc thang (Cascade Subsystem): Hồ Ka Nak xả nước qua đập dâng An Khê rồi chuyển nước xuống lưu vực sông Côn hoặc trả lại dòng chính sông Ba, kết nối trực tiếp với hồ thủy điện Sông Ba Hạ ở hạ lưu.
  • Phân hệ hồ song song (Parallel Subsystem): Hồ Ayun Hạ (trên sông Ayun), hồ Krông H’Năng (trên sông Krông H’Năng) và hồ Sông Hinh (trên sông Hinh) nằm trên các nhánh sông độc lập cùng đổ dồn lưu lượng về hợp lưu dòng chính tại trạm thủy văn Củng Sơn trước khi đổ ra Phú Lâm.

Biên điều kiện vận hành (Boundary Conditions) được xác lập dựa trên các cấp báo động thủy văn quốc gia tại các trạm khống chế:

  • Trạm An Khê: Báo động I (403,50 m), Báo động II (404,50 m), Báo động III (405,50 m).
  • Trạm Ayun Pa: Báo động I (153,00 m), Báo động II (154,50 m), Báo động III (156,00 m).
  • Trạm Củng Sơn: Báo động I (31,50 m), Báo động II (33,00 m), Báo động III (34,50 m).
  • Trạm Phú Lâm: Báo động I (1,70 m), Báo động II (2,70 m), Báo động III (3,70 m).

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng (Positivism) kết hợp phương pháp luận kỹ thuật tính toán định lượng (Computational Engineering Hydrology). Thiết kế nghiên cứu phân tầng đa cấp độ (Multi-level design):

  • Cấp độ nút (Node level): Mô hình hóa động học hồ chứa đơn lẻ dựa trên phương trình cân bằng liên tục: $\frac{dW}{dt} = Q_{\text{đến}}(t) - Q_{\text{xả}}(t)$.
  • Cấp độ nhánh (Reach level): Diễn toán quá trình biến dạng và dịch chuyển sóng lũ trên các đoạn sông tự nhiên bằng phương pháp sai phân Muskingum-Cunge.
  • Cấp độ toàn mạng lưới (Network level): Bộ điều khiển mờ (FLC) trung tâm đóng vai trò điều phối tổng thể thông tin toàn lưu vực.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu trải qua các bước chuẩn hóa nghiêm ngặt:

  1. Thu thập và số hóa dữ liệu: Đo đạc bình đồ địa hình lòng dẫn, đường đặc tính dung tích hồ $V = f(H)$, đường đặc tính xả của công trình tràn $Q = f(H, a)$, số liệu mưa và mực nước thực đo từ 15 trạm khí tượng thủy văn trên lưu vực sông Ba.
  2. Thiết lập mô hình diễn toán sóng lũ Muskingum: Phương trình dung tích trữ của đoạn sông: $$S = K \cdot O + K \cdot X \cdot (I - O) = K \left[ X \cdot I + (1 - X) O \right]$$ Trong đó $K$ là thời gian truyền sóng lũ; $X$ là hệ số trọng số không thứ nguyên ($0 \le X \le 0,5$); $I$ và $O$ lần lượt là lưu lượng dòng chảy vào và ra đoạn sông. Phương trình sai phân lưu lượng ra ở bước thời gian $t$: $$O_t = C_1 I_t + C_2 I_{t-1} + C_3 O_{t-1}$$ Với các hệ số thỏa mãn điều kiện $C_1 + C_2 + C_3 = 1$: $$C_1 = \frac{\Delta t - 2KX}{2K(1-X) + \Delta t}; \quad C_2 = \frac{\Delta t + 2KX}{2K(1-X) + \Delta t}; \quad C_3 = \frac{2K(1-X) - \Delta t}{2K(1-X) + \Delta t}$$
  3. Hiệu chỉnh và kiểm định: Sử dụng các trận lũ lớn lịch sử năm 1988 và 1993 để hiệu chỉnh bộ thông số $K, X$ cho từng đoạn sông (An Khê – Củng Sơn; Ayun Pa – Củng Sơn; Củng Sơn – Phú Lâm); kiểm định độc lập với trận lũ năm 1981 và 2013.
  4. Triangulation (Tam giác đạc phương pháp): So sánh chéo kết quả tính toán thủy lực giữa mô hình viết trên môi trường MATLAB với mô hình thủy lực 1D MIKE 11 (DHI) và mô hình cân bằng HEC-RESSIM (US Army Corps of Engineers).
+---------------------------------------------------------------------------------------+
|                       MA TRẬN DỮ LIỆU ĐẶC TÍNH 5 HỒ CHỨA CHÍNH                        |
+--------------+-------------+-------------+---------------+-------------+--------------+
| Hồ chứa      | Ka Nak      | Ayun Hạ     | Krông H'Năng  | Sông Ba Hạ  | Sông Hinh    |
+--------------+-------------+-------------+---------------+-------------+--------------+
| MN Chết (m)  | 485,00      | 196,00      | 220,00        | 101,00      | 196,00       |
| MN Dâng BT(m)| 515,00      | 205,00      | 255,00        | 105,00      | 209,00       |
| MNCNTL (m)   | 513,00      | 203,00      | 252,50        | 103,00      | 207,00       |
| MN Đón Lũ (m)| 506,00      | 202,00      | 251,50        | 102,00      | 204,50       |
| V Toàn bộ    | 314,0 tr.m3 | 253,0 tr.m3 | 167,0 tr.m3   | 349,0 tr.m3 | 357,0 tr.m3  |
+--------------+-------------+-------------+---------------+-------------+--------------+

Data và phân tích

Toàn bộ thuật toán điều khiển mờ được lập trình và xử lý trong môi trường MATLAB Fuzzy Logic Toolbox kết hợp các đoạn mã script tự động hóa tính toán thủy văn.

Biến mờ hóa bao gồm:

  • Biến $Q_{\text{đến}}$: Phân thành các tập mờ ngôn ngữ {Rất Nhỏ, Nhỏ, Vừa, Lớn, Rất Lớn, Cực Lớn}.
  • Biến $H_{\text{hồ}}$: Phân thành {Thấp, Trung Bình, Cao, Rất Cao}.
  • Biến $H_{\text{hạ du}}$: Phân theo các mức cảnh giới {Dưới BĐ I, BĐ I - BĐ II, BĐ II - BĐ III, Trên BĐ III}.
  • Biến đầu ra $Q_{\text{xả}}$: Được chia nhỏ thành 7 mức điều tiết lưu lượng từ không xả đến xả cực đại qua tràn sâu và tràn mặt.

Kiểm tra độ nhạy (Robustness checks) được thực hiện với các bước thời gian tính toán $\Delta t = 1\text{h}, 3\text{h}, 6\text{h}$, đảm bảo sự hội tụ số học và tính ổn định của sóng lũ truyền dọc hệ thống sông không phụ thuộc vào bước rời rạc hóa.

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

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

Luận án mang lại 4 phát hiện mang tính bước ngoặt khoa học và thực tiễn:

Thứ nhất, chứng minh sự bất hợp lý nghiêm trọng của Quy trình 1077 khi bắt buộc đồng loạt 5 hồ phải hạ mực nước đón lũ trong mọi hình thế thời tiết. Bằng chứng trích xuất từ luận án chỉ ra: "Thực tế vận hành các hồ chứa trong mùa lũ theo Quy trình 1077 cho đến nay đã bộc lộ một số bất cập, như khi chỉ có mưa lớn sinh lũ ở phần hạ lưu sông Ba hoặc khi các tiểu lưu vực khác trên sông Ba không có mưa mà chỉ có mưa cục bộ tại tiểu lưu vực hồ Krông H’Năng thì việc vận hành xả nước đón lũ của các hồ ở thượng du không có lưu lượng lũ đến lớn như Ka Nak, Ayun Hạ theo quy trình liên hồ hiện hành là không mang lại hiệu quả cho chống lũ đồng thời giảm hiệu quả phát điện." Khi chuyển sang bộ điều khiển mờ thích nghi, hồ Ka Nak và Ayun Hạ chỉ xả hạ mực nước khi thực sự xuất hiện tín hiệu lũ thượng nguồn, tiết kiệm hàng chục triệu mét khối nước phục vụ phát điện.

Thứ hai, định lượng hóa chính xác hiệu quả cắt giảm đỉnh lũ của mô hình mờ. Kết quả mô phỏng điều tiết hệ thống với các trận lũ kiểm định chứng minh: "Khi điều tiết theo quy trình đề xuất, có thể giảm đỉnh lũ $Q_{\max}$ tại Củng Sơn xuống 20-25%, mực nước tại Củng Sơn giảm 0,80-1,50m, còn tại Phú Lâm giảm từ 0,30-0,80m."

                           HIỆU QUẢ CẮT GIẢM MỰC NƯỚC ĐỈNH LŨ
      Trạm Củng Sơn                                     Trạm Phú Lâm
+-------------------------+                       +-------------------------+
| Giảm đỉnh: 0.80 - 1.50 m|                       | Giảm đỉnh: 0.30 - 0.80 m|
| Giảm Qmax: 20 - 25%     |                       | Khống chế < BĐ III      |
+-------------------------+                       +-------------------------+

Thứ ba, phát hiện ngưỡng giới hạn vật lý của dung tích phòng lũ lưu vực. Luận án chỉ rõ quy tắc biên: "Không thể khống chế mực nước tại Phú Lâm xuống dưới BĐ III đối với lũ vượt tần suất $P=10%$, chỉ có thể khống chế đối với lũ nhỏ hơn mức $P=20%$ (tương đương với lũ năm 1986, 1988) trở xuống." Phát hiện này phản bác các kỳ vọng phi thực tế của chính quyền địa phương về việc giao phó hoàn toàn nhiệm vụ chống ngập triệt để cho hệ thống hồ chứa khi gặp các trận đại hồng thủy.

Thứ tư, giải quyết mâu thuẫn giữa dung tích phòng lũ và tích nước cuối vụ. Trích dẫn nguyên văn điều khoản hạn chế trong văn bản cũ: "Việc thực hiện theo quy định tại khoản 3 Điều 5 của Quy trình 'Trong thời kỳ mùa lũ, khi chưa tham gia vận hành giảm lũ cho hạ du, mực nước các hồ chứa không được vượt mực nước cao nhất trước lũ...' cũng làm giảm đầu nước phát điện và tiềm ẩn nguy cơ không tích được đủ nước vào cuối mùa lũ." Mô hình mờ cho phép duy trì mực nước cao linh hoạt gần MNDBT trong giai đoạn cảnh báo sớm mà vẫn đảm bảo thời gian hạ thấp an toàn trước khi đỉnh lũ thực tế xuất hiện tại tuyến công trình.

Implications đa chiều

Về mặt lý thuyết: Mở ra hướng tiếp cận Soft Computing trong Thủy văn học công trình tại Việt Nam, chứng minh tính ưu việt của Logic mờ so với các mô hình toán kinh điển trong bối cảnh dữ liệu đầu vào có độ bất định cao.

Về mặt phương pháp luận: Thiết lập quy trình chuẩn tích hợp công cụ trí tuệ nhân tạo (FLC) với các phần mềm thủy lực truyền thống (Muskingum-Cunge, MIKE 11, HEC-RESSIM), có khả năng chuyển giao áp dụng cho các lưu vực sông phức tạp khác như sông Vu Gia – Thu Bồn, sông Cả, sông Mã, sông Đồng Nai.

Về mặt chính sách và thực tiễn quản lý: Cung cấp luận cứ khoa học trực tiếp cho Bộ Tài nguyên và Môi trường, Bộ Công Thương trong việc rà soát, điều chỉnh Quy trình vận hành liên hồ chứa lưu vực sông Ba (thay thế hoặc bổ sung cho Quyết định 1077/QĐ-TTg và Quyết định 282/QĐ-TTg).

Limitations và Future Research

Luận án thừa nhận 4 giới hạn khoa học khách quan:

  1. Mô hình diễn toán thủy lực nhánh sông mới áp dụng phương pháp thủy văn Muskingum-Cunge 1 chiều đơn giản hóa, chưa phản ánh đầy đủ hiệu ứng nước dâng do thủy triều và ảnh hưởng của sóng cồn bão tại vùng cửa sông Đà Rằng (Phú Lâm).
  2. Hệ quy tắc mờ (Rule base) được xây dựng chủ yếu dựa trên tri thức chuyên gia chuyên ngành kết hợp số liệu lịch sử tĩnh; chưa tích hợp thuật toán tự học thích nghi thời gian thực (Adaptive Neuro-Fuzzy Inference System - ANFIS).
  3. Độ tin cậy của mô hình phụ thuộc mật thiết vào mạng lưới quan trắc mưa tự động và mật độ trạm đo trên các tiểu lưu vực vùng núi cao hẻo lánh phía Tây Trường Sơn.

Hướng nghiên cứu tiếp nối trong giai đoạn 5–10 năm tới bao gồm:

  • Mở rộng phát triển hệ thống điều hành mờ tự thích nghi trực tuyến (Real-time Online ANFIS) kết nối API trực tiếp với dữ liệu radar thời tiết và ảnh mây vệ tinh phân giải cao.
  • Nâng cấp mô hình thủy lực 1D lên mô hình ngập lụt 2D thủy động lực học hoàn chỉnh (2D Saint-Venant) để mô phỏng chính xác mức độ ngập lụt chi tiết cho các vùng đô thị hạ du như Thành phố Tuy Hòa và thị xã An Khê.
  • Nghiên cứu cơ chế vận hành tối ưu đa mục tiêu có tính đến tác động cực đoan của Biến đổi khí hậu toàn cầu theo các kịch bản phát thải mới nhất của IPCC.

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

Về mặt học thuật: Luận án tạo ra bước đột phá về phương pháp luận, mở ra dòng nghiên cứu ứng dụng Trí tuệ tính toán (Computational Intelligence) trong quản lý rủi ro thiên tai thủy văn tại các trường đại học và viện nghiên cứu chuyên ngành.

Về mặt công nghiệp và năng lượng: Đóng góp giải pháp kỹ thuật trực tiếp cho Tập đoàn Điện lực Việt Nam (EVN), Công ty Cổ phần Sông Ba (chủ quản NMTĐ Krông H’Năng) và các công ty thủy điện trên lưu vực sông Ba nhằm tối ưu hóa sản lượng điện năng, giảm thiểu lượng nước xả thừa qua tràn, gia tăng hàng trăm tỷ đồng giá trị kinh tế mỗi năm từ nguồn thủy năng.

Về mặt xã hội và an sinh: Giúp chính quyền hai tỉnh Phú Yên và Gia Lai chủ động chuyển dịch từ ứng phó thiên tai thụ động sang quản trị rủi ro thiên tai thích ứng, giảm thiểu thiệt hại về người và tài sản cho hàng triệu người dân hạ du trong mùa mưa bão.

                                  MA TRẬN TÁC ĐỘNG CỦA CÔNG TRÌNH
+-----------------------------------------------------------------------------------------------+
| HỌC THUẬT & ĐÀO TẠO    | Đặt nền móng lý thuyết Logic mờ ứng dụng trong Thủy văn học VN       |
+------------------------+----------------------------------------------------------------------+
| KINH TẾ & NĂNG LƯỢNG   | Tăng sản lượng điện, bảo toàn đầu nước, giảm xả thừa cho EVN         |
+------------------------+----------------------------------------------------------------------+
| AN SINH XÃ HỘI         | Cắt 20-25% đỉnh lũ, giảm ngập lụt cho hạ du tỉnh Phú Yên và Gia Lai  |
+------------------------+----------------------------------------------------------------------+
| CHÍNH SÁCH NHÀ NƯỚC    | Luận cứ khoa học sửa đổi Quy trình vận hành liên hồ chứa Quốc gia    |
+-----------------------------------------------------------------------------------------------+

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

  • Nghiên cứu sinh và Giảng viên chuyên ngành Thủy văn, Tài nguyên nước: Kế thừa khung phương pháp luận kết hợp FIS và thủy lực diễn toán dòng chảy lũ.
  • Kỹ sư điều độ hệ thống điện và Thủy điện: Nắm bắt công cụ ra quyết định xả lũ thích nghi khoa học, chính xác, tự động hóa cao.
  • Cơ quan quản lý Nhà nước (Cục Quản lý Tài nguyên nước, Ban Chỉ đạo Quốc gia về Phòng chống thiên tai): Sở hữu căn cứ thực nghiệm để hoàn thiện thể chế và quy chuẩn điều hành liên hồ.
  • Cộng đồng dân cư lưu vực sông Ba: Hưởng lợi trực tiếp từ việc giảm thiểu tần suất và biên độ ngập lụt cực đoan.

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 độc đáo nhất là việc mở rộng Thuyết tập mờ của L.A. Zadeh (1965) và Quy tắc suy luận Mamdani vào bài toán điều phối liên hồ chứa hỗn hợp bậc thang – song song. Thay vì xử lý các bài toán hồ đơn lẻ tĩnh, luận án đã toán học hóa mối tương tác động giữa các biến mờ đa chiều ($Q_{\text{đến}}, H_{\text{hồ}}, H_{\text{hạ du}}$) thông qua thuật toán giải mờ COG, tạo nền tảng cho việc điều hành theo thời gian thực thích nghi với sai số dự báo thủy văn.

  2. Tính đổi mới về mặt phương pháp luận so với các nghiên cứu tiền nhiệm? So với các nghiên cứu tối ưu kinh điển như Young (1967) hay Janga Reddy & Kumar (2006) vốn phụ thuộc vào chuỗi dữ liệu tất định và tốn nhiều thời gian tính toán, luận án đã tích hợp bộ điều khiển mờ trực tiếp với mô hình diễn toán Muskingum-Cunge. Phương pháp này cho phép phản hồi tức thời trong vài giây tại mỗi bước thời gian thực mà không đòi hỏi dự báo chính xác tuyệt đối lưu lượng lũ đến trước 3–5 ngày.

  3. Phát hiện bất ngờ nhất từ thực nghiệm dữ liệu là gì? Phát hiện bất ngờ nhất là việc tuân thủ cứng nhắc quy định duy trì hồ ở "Mực nước cao nhất trước lũ" theo Quy trình 1077 thực chất không mang lại lợi ích cắt lũ thực tế khi mưa chỉ xảy ra cục bộ ở hạ lưu hoặc tiểu lưu vực biệt lập, ngược lại còn gây lãng phí nghiêm trọng nguồn năng lượng thủy điện và dẫn đến nguy cơ không thể tích đủ nước vào cuối mùa lũ.

  4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không? Luận án cung cấp toàn bộ quy trình tái lập chi tiết bao gồm: ma trận phân chia tập mờ, bảng 10 luật mờ điều khiển trung tâm, phương trình giải tích giải mờ COG, các tham số sai phân $C_1, C_2, C_3$ của phương pháp Muskingum-Cunge và sơ đồ thuật toán xử lý dữ liệu trong MATLAB.

  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào? Lộ trình 10 năm tập trung vào 3 trọng tâm: (1) Tích hợp học máy tự động (ANFIS, Deep Learning) vào cơ chế sinh luật mờ; (2) Mở rộng mạng lưới dự báo lũ số trị thời gian thực liên kết mô hình khí tượng phân giải cao WRF; (3) Tích hợp bài toán điều hành lũ với bài toán xả mặn và cân bằng sinh thái hạ lưu sông trong bối cảnh biến đổi khí hậu.

Kết luận

Luận án của nghiên cứu sinh Nguyễn Việt xác lập 5 đóng góp học thuật và thực tiễn cốt lõi:

  1. Xác lập thành công cơ sở khoa học và phương pháp luận ứng dụng Logic mờ trong bài toán điều hành liên hồ chứa chống lũ tại Việt Nam.
  2. Thiết lập mô hình toán tích hợp FLC và diễn toán Muskingum-Cunge, giải quyết triệt để bài toán điều khiển hệ thống liên hồ phức hợp (bậc thang Ka Nak – An Khê – Sông Ba Hạ kết hợp song song Ayun Hạ – Krông H’Năng – Sông Hinh).
  3. Chứng minh và lượng hóa năng lực cắt giảm đỉnh lũ hạ du đạt 20–25% tại trạm Củng Sơn, hạ mực nước ngập lụt từ 0,80–1,50 m tại Củng Sơn và 0,30–0,80 m tại Phú Lâm.
  4. Phát hiện các bất cập căn bản của quy trình vận hành liên hồ hiện hành (Quy trình 1077), đưa ra giải pháp bảo toàn cột nước phát điện và dung tích tích nước cuối vụ mà vẫn đảm bảo an toàn tuyệt đối cho công trình và hạ du.
  5. Định hình mô hình quản trị thích ứng trong Thủy văn học công trình, mở đường cho việc xây dựng các hệ thống điều hành liên hồ chứa thông minh theo thời gian thực trên toàn quốc.