Tổng quan về luận án
Bối cảnh biến đổi khí hậu và sự gia tăng của các hình thái thời tiết cực đoan tại miền Trung Việt Nam đã đặt ra thách thức nghiêm trọng cho công tác quản lý tài nguyên nước và phòng chống thiên tai. Lưu vực sông Vu Gia – Thu Bồn với diện tích tự nhiên $10.350\text{ km}^2$, địa hình sườn Đông dãy Trường Sơn dốc và chia cắt mạnh, là khu vực thường xuyên hứng chịu các trận lũ đặc biệt lớn có cường suất lũ cao, lũ lên nhanh và rút nhanh. Trong khi đó, hệ thống các công trình thủy điện bậc thang thượng nguồn (A Vương, Đăk Mi 4, Sông Tranh 2, Sông Bung 2, Sông Bung 4) được xây dựng với mục tiêu chính là phát điện, khiến dung tích phòng lũ dành cho hạ du bị suy giảm nghiêm trọng so với quy hoạch ban đầu—tổng dung tích phòng lũ thực tế chỉ đạt $146,44\text{ triệu m}^3$, chiếm vỏn vẹn khoảng $10%$ tổng dung tích hữu ích ($1.453,24\text{ triệu m}^3$), thấp hơn rất nhiều so với con số $1.187\text{ triệu m}^3$ (chiếm $42%$ dung tích hữu ích) được phê duyệt trong Quyết định số 1336/QĐ-TTg.
Thực tiễn vận hành trong các mùa lũ lịch sử (như trận lũ năm 2009 trên sông Vu Gia – Thu Bồn và năm 2010 trên sông Ba) đã bộc lộ lỗ hổng nguy hiểm: việc áp dụng các "quy trình vận hành cứng" thiếu công cụ dự báo và cảnh báo thời gian thực đã dẫn đến tình trạng "quá trình xả gây 'sốc' cho hạ du do lưu lượng xả tăng đột biến trong một thời gian ngắn" hoặc gây ra hiện tượng "lũ chồng lên lũ".
Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ việc các quy trình liên hồ chứa hiện hành chủ yếu dựa vào dự báo dòng chảy ngắn hạn (6–24 giờ) hoặc các quy tắc vận hành tĩnh, chưa tích hợp được mô hình dự báo mưa số trị hạn vừa (3–5 ngày) vào việc ra quyết định điều tiết hồ chứa theo thời gian thực cho các lưu vực sông dốc miền Trung.
Luận án của NCS. Tô Thúy Nga (2014) dưới sự hướng dẫn của GS. Nguyễn Thế Hùng và GS. Hà Văn Khối tại Trường Đại học Bách khoa – Đại học Đà Nẵng đã giải quyết trực tiếp khoảng trống này thông qua việc thiết lập bài toán điều tiết tối ưu – mô phỏng động. Nghiên cứu tập trung giải quyết ba câu hỏi nghiên cứu (Research Questions):
- RQ1: Làm thế nào để xây dựng một khung toán học vận hành thời gian thực tích hợp dự báo khí tượng hạn vừa (3–5 ngày) nhằm giải quyết mâu thuẫn giữa an toàn tích nước phát điện và cắt giảm đỉnh lũ hạ du?
- RQ2: Cấu trúc liên kết mô hình mưa – dòng chảy, điều tiết hồ chứa thượng lưu và diễn toán thủy lực một chiều vùng đồng bằng hạ lưu cần được thiết lập như thế nào để đảm bảo tính ổn định số và tốc độ tính toán trong điều kiện tác nghiệp?
- RQ3: Mức độ điều chỉnh dung tích đón lũ và quỹ đạo xả lũ tối ưu nào cho phép giảm thiểu nguy cơ ngập lụt tại các điểm kiểm soát khống chế (Hội Khách, Ái Nghĩa, Giao Thủy, Câu Lâu) mà tổn thất điện năng không đáng kể?
Ba giả thuyết nghiên cứu (Hypotheses) tương ứng được kiểm chứng:
- H1: Kéo dài thời gian dự kiến dự báo lũ lên 3–5 ngày bằng mô hình số trị quy mô vừa (HRM, ETA) cho phép chủ động hạ thấp mực nước đón lũ của các hồ chứa trước khi đỉnh lũ xuất hiện mà không làm suy giảm khả năng tích đầy hồ cuối vụ.
- H2: Mô hình tích hợp MOPHONG-LU kết nối MIKE 11 phản ánh chính xác tương tác thủy lực phức tạp giữa lưu vực sông dốc thượng nguồn và hệ thống phân lưu hạ lưu (sông Quảng Huế, sông Vĩnh Điện).
- H3: Tồn tại giải pháp thỏa hiệp tối ưu (Pareto-optimal trade-off) cho phép tăng dung tích phòng lũ hữu hiệu thêm $175,71 - 199,20\text{ triệu m}^3$, giúp giảm mực nước ngập lụt hạ du từ $0,10 - 0,60\text{ m}$ với mức độ suy giảm sản lượng điện dưới ngưỡng $3%$.
Literature Review và Positioning
Nghiên cứu về điều tiết hệ thống hồ chứa trên thế giới và tại Việt Nam phát triển qua ba trường phái phương pháp luận chính:
- Phương pháp tối ưu hóa thuần túy (Pure Optimization): Ứng dụng quy hoạch tuyến tính (LP), quy hoạch phi tuyến (NLP), quy hoạch động ngẫu nhiên (SDP) và thuật toán tiến hóa (GA, PSO). Tiêu biểu như El-Said Mohamed Said Ahmed và Larry W. Mays (1992) áp dụng điều khiển tối ưu rời rạc cho hệ thống hồ chứa Lake Travis (Texas); Aleix Serrat-Capdevila và Juan B. Valdes (2007) ứng dụng SDP cho lưu vực quốc tế Rio Grande; Li-Chiu Chang (2008) và Chih-Chiang Wei, Nien-Sheng Hsu (2009) áp dụng thuật toán di truyền (GA) và cây quyết định cho hệ thống hồ chứa Shihmen (Đài Loan). Hạn chế của trường phái này là số lượng ràng buộc phi tuyến thủy lực quá lớn, thời gian hội tụ chậm, khó áp dụng cho bài toán điều hành tác nghiệp thời gian thực với mạng lưới sông phức tạp.
- Phương pháp mô phỏng thuần túy (Simulation Approach): Tập trung vào việc giải các phương trình thủy văn, thủy lực với các kịch bản định trước. Chung-Tian Cheng và K. Chau (2004) phát triển hệ thống SOPs cho lưu vực sông Dương Tử và Liaohe; Xiang-Yang Li và cộng sự (2006) xây dựng hệ thống WFFS trên nền tảng Web cho lưu vực Shuangpai; Giha Lee và cộng sự (2011) ứng dụng mô hình OHyMoS cho đập Daechung (Hàn Quốc); Wan Xin-yu và cộng sự (2012) mô phỏng hệ thống kiểm soát lũ Zhengyanghuan. Phương pháp này mô tả chi tiết dòng chảy nhưng thiếu cơ chế tìm kiếm giải pháp tối ưu tự động giữa các mục tiêu xung đột.
- Phương pháp tích hợp Mô phỏng – Tối ưu hóa (Simulation-Optimization Coupling): Trường phái hiện đại kết hợp tính chính xác của mô hình thủy lực với năng lực tối ưu đa mục tiêu (Pareto). Điển hình là nghiên cứu của Long Le Ngo, Henrik Madsen, Dan Rosbjerg (2007) và Bertrand Richaud và cộng sự (2011) cho hồ chứa Hòa Bình trên lưu vực sông Hồng; Bahram Malekmohammadi, Banafsheh Zahraie, Reza Kerachian (2009, 2011) kết hợp GA, GIS và HEC-RAS cho hệ thống sông Bakhtiari – Dez (Iran) và sông Han (Hàn Quốc); Salehi Neyshabouri (2011) tích hợp MIKE 11 và NSGA-II.
Tại Việt Nam, các nghiên cứu vận hành thời gian thực trước đây tập trung chủ yếu ở lưu vực sông Hồng (Nguyễn Lan Châu & Nguyễn Quốc Anh, 2005; Ngô Huy Cẩn và cs., 2006; Đoàn Xuân Thủy, Hà Ngọc Hiến, Nguyễn Văn Điệp, 2006; Hà Ngọc Hiến và cs., 2010; Nguyễn Lan Châu & Bùi Đình Lập, 2011). Đối với các lưu vực sông miền Trung, phần lớn chỉ dừng lại ở việc đánh giá quy hoạch hoặc thiết lập quy trình vận hành liên hồ tĩnh (Hoàng Minh Tuyển, 2010 cho sông Ba; Ngô Lê Long, 2011 cho sông Srêpôk; Hoàng Thanh Tùng và cs., 2011 cho sông Cả; Đặng Thanh Mai và cs., 2013).
TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT VẬN HÀNH HỒ CHỨA
Tối ưu hóa thuần túy (Mays 1992, Valdes 2007) Mô phỏng thuần túy (Chau 2004, Lee 2011)
\ /
\ /
---> Mô phỏng - Tối ưu hóa (Madsen 2007, Malekmohammadi 2009) <---
|
Khoảng trống: Thiếu khung điều tiết thời gian thực tích hợp dự báo khí tượng 3-5 ngày
cho lưu vực dốc, lũ tập trung nhanh, đa hồ chứa miền Trung
|
v
LUẬN ÁN TÔ THÚY NGA (2014): MOPHONG-LU + MIKE 11
Luận án định vị tại giao điểm của trường phái Mô phỏng – Tối ưu hóa thỏa hiệp: Luận án không tìm kiếm nghiệm tối ưu hình thức khó khả thi trong điều kiện dữ liệu thực tế, mà phát triển mô hình tính toán hỗ trợ ra quyết định thời gian thực, liên kết trực tiếp kết quả dự báo mưa số trị 3–5 ngày với mô hình phân tích kịch bản dòng chảy nhằm xác định phương án điều tiết "hợp lý", giải quyết triệt để nghịch lý giữa an toàn xả lũ hạ du và hiệu quả kinh tế phát điện.
Đó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 Lý thuyết Điều khiển Thích nghi Hệ thống Tài nguyên nước (Adaptive Water Resources Systems Theory) và Lý thuyết Vận hành Đa mục tiêu (Multi-Objective Reservoir Operation Theory) qua việc định hình toán học cho hàm ra quyết định thời gian thực:
$$\mathbf{U}_t = \mathcal{F}\Big(\mathbf{S}t, \widehat{\mathbf{I}}{t \to t+k}, \mathbf{\Omega}_t\Big)$$
Trong đó:
- $\mathbf{U}t$: Vectơ quyết định xả lũ của hệ thống $N$ hồ chứa tại thời điểm $t$, $\mathbf{U}t = [Q{x,1}(t), Q{x,2}(t), \dots, Q_{x,N}(t)]^T$.
- $\mathbf{S}t$: Vectơ trạng thái dung tích/mực nước hồ chứa tại thời điểm bắt đầu bước tính, $\mathbf{S}t = [Z{1}(t), Z{2}(t), \dots, Z_{N}(t)]^T$.
- $\widehat{\mathbf{I}}_{t \to t+k}$: Ma trận dự báo dòng chảy đến các nút hồ chứa và các nhập lưu trung gian trong khoảng thời gian dự kiến $k = 3 \div 5\text{ ngày}$ ($72 \div 120\text{ giờ}$), suy biến từ trường mưa số trị dự báo.
- $\mathbf{\Omega}t$: Tập hợp các ràng buộc an toàn công trình, giới hạn xả không gây "sốc" thủy lực hạ du ($\partial Q_x / \partial t \le \Delta Q{\text{allow}}$) và ngưỡng báo động lũ tại các trạm thủy văn hạ lưu (Ái Nghĩa, Giao Thủy, Câu Lâu).
Nghiên cứu mang lại bước chuyển dịch mô thức (paradigm shift): Chuyển từ quan điểm "điều tiết dựa trên dòng chảy đã xuất hiện hoặc dự báo cực ngắn (6–12h)" sang "điều tiết chủ động dựa trên cửa sổ dự báo khí tượng mở rộng (72–120h)". Sự chuyển dịch này phá vỡ tiên đề vận hành tĩnh vốn xem mực nước đón lũ là một hằng số cố định suốt mùa lũ.
KHUNG KHÁI NIỆM ĐIỀU TIẾT THỜI GIAN THỰC
+---------------------------------------------------------------------------------+
| DỮ LIỆU KHÍ TƯỢNG SỐ TRỊ |
| Mô hình dự báo thời tiết quy mô vừa: HRM, ETA, BOLAM (GFS) |
| Thời gian dự kiến: 3 - 5 ngày (72h - 120h) |
+---------------------------------------------------------------------------------+
|
v
+---------------------------------------------------------------------------------+
| KHU VỰC THƯỢNG LƯU: MOPHONG-LU |
| 1. Mô-đun Mưa - Dòng chảy (NAM / TANK): Chuyển đổi mưa dự báo thành dòng chảy |
| 2. Mô-đun Diễn toán Lũ trong Sông (Muskingum / Muskingum-Cunge) |
| 3. Mô-đun Vận hành 5 Hồ Chứa: A Vương, Đăk Mi 4, S. Tranh 2, S. Bung 2 & 4 |
+---------------------------------------------------------------------------------+
|
Biên lưu lượng xả & thủy văn nhập lưu thượng nguồn
|
v
+---------------------------------------------------------------------------------+
| KHU VỰC HẠ LƯU: MIKE 11 HD |
| - Thủy lực 1 chiều Saint-Venant (Hội Khách & Nông Sơn -> Cửa Đại / Đà Nẵng) |
| - Tương tác trao đổi dòng qua sông Quảng Huế & sông Vĩnh Điện |
| - Kiểm soát mực nước ngập lụt: Ái Nghĩa, Giao Thủy, Câu Lâu, Hội An |
+---------------------------------------------------------------------------------+
|
v
+---------------------------------------------------------------------------------+
| RA QUYẾT ĐỊNH ĐIỀU TIẾT THỎA HIỆP PARETO |
| Tối ưu hóa: Maximize (Giảm đỉnh lũ hạ du) & Minimize (Tổn thất điện) |
| Cập nhật chu kỳ Delta t = 6h, 12h, 24h theo thời gian thực |
+---------------------------------------------------------------------------------+
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp ba lý thuyết chuyên ngành:
- Lý thuyết Thủy văn Mô hình Khái niệm (Conceptual Hydrological Theory): Vận dụng cấu trúc bể chứa phi tuyến của mô hình NAM/TANK để mô phỏng cơ chế bão hòa độ ẩm đất, dòng chảy mặt, dòng chảy sát mặt và dòng chảy ngầm trên $10.350\text{ km}^2$ lưu vực.
- Lý thuyết Thủy động lực học và Diễn toán Thủy lực (Hydrodynamic & Hydraulic Routing Theory): Kết hợp phương pháp biến đổi sóng động học/sóng khuếch tán Muskingum-Cunge cho mạng lưới sông miền núi dốc và hệ phương trình vi phân đạo hàm riêng Saint-Venant đầy đủ 1 chiều trong MIKE 11 cho mạng sông đồng bằng hạ du chịu ảnh hưởng triều.
- Lý thuyết Tối ưu hóa Đa mục tiêu Thỏa hiệp Pareto (Pareto Compromise Decision Analysis): Xác lập đường cong đánh đổi biên giữa hai hàm mục tiêu xung đột:
$$\min f_1 = \max_{t} H_{\text{downstream}}(t)$$
$$\min f_2 = \Delta E = \sum_{i=1}^{N} \int (Q_{\text{design}, i} - Q_{\text{actual}, i}) \cdot H_{\text{head}, i} \cdot \eta_i , dt$$
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 phân tích hệ thống kỹ thuật công trình (Systems Engineering Approach). Thiết kế nghiên cứu đa tầng bậc (Multi-level Design) bao gồm:
- Tầng vĩ mô (Lưu vực): Phân tích tương tác khí quyển – thủy văn trên toàn lưu vực Vu Gia – Thu Bồn.
- Tầng trung mô (Nút công trình): Cân bằng nước và động lực học điều tiết tại 5 hồ chứa bậc thang.
- Tầng vi mô (Mạng lưới sông và bãi ngập lũ): Thủy lực chi tiết tại các vị trí phân lưu, hợp lưu phức tạp (sông Quảng Huế tách nước từ Vu Gia sang Thu Bồn, sông Vĩnh Điện chuyển ngược dòng nước).
Quy trình nghiên cứu rigorous
Quy trình nghiên cứu được chuẩn hóa qua 5 giai đoạn liên hoàn:
+---------------------------------------------------------------------------------+
| GIAI ĐOẠN 1: THU THẬP VÀ XỬ LÝ DỮ LIỆU ĐA NGUỒN |
| - Khí tượng thủy văn: Chuỗi số liệu 1976 - 2010 từ 19 trạm đo mưa, 8 trạm thủy văn|
| - Địa hình: Mặt cắt ngang lòng dẫn, trắc dọc sông Vu Gia, Thu Bồn, Quảng Huế |
| - Công trình: Đường đặc tính V(Z), F(Z), đường đặc tính xả lũ tràn và tổ máy |
+---------------------------------------------------------------------------------+
|
v
+---------------------------------------------------------------------------------+
| GIAI ĐOẠN 2: THIẾT LẬP VÀ HIỆU CHỈNH HỆ THỐNG MÔ HÌNH |
| - Hiệu chuẩn (Calibration) và kiểm định (Validation) mô hình NAM cho thượng lưu |
| - Hiệu chỉnh hệ số nhám Manning (n) mô hình MIKE 11 HD hạ lưu |
| - Tiêu chuẩn đánh giá: Chỉ số Nash-Sutcliffe (NSE >= 0.80), Sai số đỉnh lũ (<=10%)|
+---------------------------------------------------------------------------------+
|
v
+---------------------------------------------------------------------------------+
| GIAI ĐOẠN 3: XÂY DỰNG PHẦN MỀM TÍCH HỢP MOPHONG-LU |
| - Lập trình mô-đun kết nối tự động: Nhận mưa dự báo ETA/HRM (3-5 ngày) |
| - Tính toán liên tục quá trình dòng chảy đến hồ và diễn toán qua hồ chứa |
| - Kết xuất biên thủy lực tự động chuyển tiếp sang mô hình MIKE 11 |
+---------------------------------------------------------------------------------+
|
v
+---------------------------------------------------------------------------------+
| GIAI ĐOẠN 4: THỬ NGHIỆM VÀ ĐÁNH GIÁ CÁC KỊCH BẢN ĐIỀU TIẾT |
| - Tái hiện các trận lũ lịch sử: 1998, 1999, 2007, 2009 |
| - Mô phỏng các kịch bản mực nước đón lũ (Z_đl) và quy tắc đóng mở cửa van |
| - So sánh: Vận hành thực tế vs Quy trình 1879 vs Quy tắc đề xuất |
+---------------------------------------------------------------------------------+
|
v
+---------------------------------------------------------------------------------+
| GIAI ĐOẠN 5: XÂY DỰNG BỘ QUY TẮC ĐIỀU HÀNH THỜI GIAN THỰC |
| - Đề xuất quy trình vận hành liên hồ thời gian thực hoàn chỉnh |
| - Xây dựng giải pháp thỏa hiệp Pareto: Tối ưu dung tích lũ vs Sản lượng điện |
+---------------------------------------------------------------------------------+
Data và phân tích
Số liệu khí tượng thủy văn được thu thập xuyên suốt giai đoạn 1976–2010 với mức độ chi tiết cao:
Đặc trưng lượng mưa thời đoạn lớn nhất tại các trạm đo đại diện:
- Trạm Ái Nghĩa: $X_{1\text{ ngày}} = 501\text{ mm}$ (03/11/1999); $X_{3\text{ ngày}} = 1.073\text{ mm}$; $X_{5\text{ ngày}} = 1.230\text{ mm}$; $X_{7\text{ ngày}} = 1.306\text{ mm}$.
- Trạm Trao: $X_{1\text{ ngày}} = 482\text{ mm}$ (17/11/1980); $X_{3\text{ ngày}} = 1.124\text{ mm}$; $X_{5\text{ ngày}} = 1.488\text{ mm}$; $X_{7\text{ ngày}} = 1.677\text{ mm}$.
- Trạm Tiên Phước: $X_{1\text{ ngày}} = 534,4\text{ mm}$ (04/12/1999); $X_{3\text{ ngày}} = 1.107\text{ mm}$; $X_{5\text{ ngày}} = 1.331\text{ mm}$; $X_{7\text{ ngày}} = 1.440\text{ mm}$.
- Trạm Thành Mỹ: $X_{1\text{ ngày}} = 621,9\text{ mm}$ (02/11/1999); $X_{3\text{ ngày}} = 1.070\text{ mm}$; $X_{5\text{ ngày}} = 1.213\text{ mm}$.
- Trạm Hội An: $X_{1\text{ ngày}} = 666,6\text{ mm}$ (03/12/1999); $X_{3\text{ ngày}} = 1.140\text{ mm}$; $X_{5\text{ ngày}} = 1.251\text{ mm}$; $X_{7\text{ ngày}} = 1.315\text{ mm}$.
Đặc trưng lưu lượng và đỉnh lũ lịch sử quan trắc:
- Sông Thu Bồn tại Nông Sơn ($F_{\text{lv}} = 3.150\text{ km}^2$): Lưu lượng đỉnh lũ lớn nhất đo được đạt $Q_{\max} = 10.600\text{ m}^3/\text{s}$ (xuất hiện vào các ngày 20/11/1998, 04/12/1999 và 11/11/2007), ứng với mô số dòng chảy lũ $M_{\max} = 3,37\text{ m}^3/\text{s}\cdot\text{km}^2$, chu kỳ lặp lại xấp xỉ tần suất $P = 5,8%$. Mức độ biến động dòng chảy cực lớn: $Q_{\max}$ năm 2007 gấp $3,8\text{ lần}$ năm lũ nhỏ nhất ($1982: 2.800\text{ m}^3/\text{s}$).
- Sông Vu Gia tại Thành Mỹ ($F_{\text{lv}} = 1.850\text{ km}^2$): Đỉnh lũ lịch sử ghi nhận $Q_{\max} = 7.230\text{ m}^3/\text{s}$ (29/09/2009, tần suất $P = 3,3%$) và $7.000\text{ m}^3/\text{s}$ (20/11/1998), ứng với mô số đỉnh lũ $M_{\max} = 3,90\text{ m}^3/\text{s}\cdot\text{km}^2$. Giá trị đỉnh lũ năm 2009 gấp $5,6\text{ lần}$ so với năm lũ nhỏ nhất ($1978: 1.290\text{ m}^3/\text{s}$).
Thời gian và tốc độ truyền lũ trên mạng lưới sông:
- Đoạn Thành Mỹ – Ái Nghĩa (sông Vu Gia, chiều dài $40,1\text{ km}$): Thời gian truyền lũ trung bình $4\text{ giờ}$ (nhanh nhất $3\text{ giờ}$, chậm nhất $7\text{ giờ}$), tốc độ truyền lũ trung bình $10,1\text{ km/h}$.
- Đoạn Hiệp Đức – Nông Sơn (sông Thu Bồn, $23\text{ km}$): Thời gian truyền lũ trung bình $4\text{ giờ}$, tốc độ $5,8\text{ km/h}$.
- Đoạn Nông Sơn – Giao Thủy (sông Thu Bồn, $27\text{ km}$): Thời gian truyền lũ trung bình $4\text{ giờ}$, tốc độ $6,8\text{ km/h}$.
- Đoạn Giao Thủy – Câu Lâu (sông Thu Bồn, $26\text{ km}$): Thời gian truyền lũ trung bình $6\text{ giờ}$, tốc độ $4,3\text{ km/h}$.
Thông số kỹ thuật 5 công trình hồ chứa nghiên cứu:
- Thủy điện Sông Bung 2: $F_{\text{lv}} = 324\text{ km}^2$; $\text{MNDBT} = 605\text{ m}$; $\text{MNC} = 565\text{ m}$; $W_{\text{tb}} = 94,3 \times 10^6\text{ m}^3$; $W_{\text{hi}} = 73,9 \times 10^6\text{ m}^3$; $W_{\text{fl}} = 7,19 \times 10^6\text{ m}^3$; $N_{\text{lm}} = 100\text{ MW}$; $Z_{\text{đl}}$ đề xuất $= 600\text{ m}$ ($W_{\text{fl}}$ mở rộng $= 13,45 \times 10^6\text{ m}^3$).
- Thủy điện Sông Bung 4: $F_{\text{lv}} = 1.474\text{ km}^2$; $\text{MNDBT} = 222,5\text{ m}$; $\text{MNC} = 205\text{ m}$; $W_{\text{tb}} = 510,8 \times 10^6\text{ m}^3$; $W_{\text{hi}} = 234,0 \times 10^6\text{ m}^3$; $W_{\text{fl}} = 20,14 \times 10^6\text{ m}^3$; $N_{\text{lm}} = 156\text{ MW}$; $Z_{\text{đl}}$ đề xuất $= 218\text{ m}$ ($W_{\text{fl}}$ mở rộng $= 68,0 \times 10^6\text{ m}^3$).
- Thủy điện A Vương: $F_{\text{lv}} = 682\text{ km}^2$; $\text{MNDBT} = 380\text{ m}$; $\text{MNC} = 340\text{ m}$; $W_{\text{tb}} = 343,55 \times 10^6\text{ m}^3$; $W_{\text{hi}} = 266,48 \times 10^6\text{ m}^3$; $W_{\text{fl}}$ thiết kế $= 14,25 \times 10^6\text{ m}^3$; $N_{\text{lm}} = 210\text{ MW}$; $Z_{\text{đl}}$ quy trình $= 376\text{ m}$ ($W_{\text{fl}}$ điều hành $= 43,92 \times 10^6\text{ m}^3$).
- Thủy điện Đăk Mi 4 (Bậc 1): $F_{\text{lv}} = 1.125\text{ km}^2$; $\text{MNDBT} = 258\text{ m}$; $\text{MNC} = 240\text{ m}$; $W_{\text{tb}} = 310,32 \times 10^6\text{ m}^3$; $W_{\text{hi}} = 158,03 \times 10^6\text{ m}^3$; $W_{\text{fl}}$ thiết kế $= 2,2 \times 10^6\text{ m}^3$; $N_{\text{lm}} = 148\text{ MW}$; $Z_{\text{đl}}$ quy trình $= 255\text{ m}$ ($W_{\text{fl}}$ điều hành $= 49,94 \times 10^6\text{ m}^3$).
- Thủy điện Sông Tranh 2: $F_{\text{lv}} = 1.100\text{ km}^2$; $\text{MNDBT} = 175\text{ m}$; $\text{MNC} = 140\text{ m}$; $W_{\text{tb}} = 733,4 \times 10^6\text{ m}^3$; $W_{\text{hi}} = 521,17 \times 10^6\text{ m}^3$; $W_{\text{fl}}$ thiết kế $= 102,66 \times 10^6\text{ m}^3$; $N_{\text{lm}} = 190\text{ MW}$; $Z_{\text{đl}}$ quy trình $= 172\text{ m}$ ($W_{\text{fl}}$ điều hành $= 81,85 \times 10^6\text{ m}^3$).
Phát hiện đột phá và implications
Những phát hiện then chốt
- Khắc phục hoàn toàn hiện tượng sốc thủy lực hạ du: Nghiên cứu chứng minh rằng quy luật vận hành tối ưu không phải là tích nước ngay từ đầu thời điểm lũ lên (khiến hồ nhanh đầy và buộc phải xả tràn ồ ạt đúng lúc đỉnh lũ ngoại biên xuất hiện), mà là: Chủ động xả đón lũ ở lưu lượng an toàn không gây ngập hạ lưu trong giai đoạn đầu lũ lên, tập trung tích nước cắt đỉnh lũ ở cuối thời kỳ lũ lên và chuyển sang tích nước đầy hồ ở thời kỳ lũ xuống.
- Hiệu quả giảm ngập lụt hạ du có định lượng rõ rệt: Tái mô phỏng trận lũ năm 2009 cho thấy nếu vận hành theo phương thức đề xuất của luận án:
- Mực nước đỉnh lũ tại trạm Hội Khách giảm $0,606\text{ m}$.
- Mực nước tại trạm Ái Nghĩa giảm $0,109\text{ m}$.
- Mực nước tại trạm Giao Thủy giảm $0,114\text{ m}$.
- Mực nước tại trạm Câu Lâu giảm $0,012\text{ m}$.
- Mở rộng dung tích phòng lũ khả thi với tổn thất điện năng tối thiểu: Đề xuất điều chỉnh mực nước đón lũ ($Z_{\text{đl}}$) của 5 hồ chứa tạo thêm tổng dung tích phòng lũ thực tế đạt $199,20\text{ triệu m}^3$ (thay vì $146,44\text{ triệu m}^3$ theo thiết kế gốc hoặc $175,71\text{ triệu m}^3$ theo quy trình tạm thời). Nhờ có thông tin dự báo mưa 3–5 ngày, các hồ vẫn đảm bảo xác suất tích đầy về mực nước dâng bình thường vào cuối trận lũ/cuối mùa lũ, tổn thất điện năng ước tính dưới $2,5%$.
- Giải quyết nghịch lý phân chia trách nhiệm dự báo – vận hành: Luận án chỉ ra sự bất cập sâu sắc trong cơ chế phối hợp cũ (khi Trung tâm DBKTTV chỉ dự báo mực nước hạ lưu 24h sau khi hồ xả, còn chủ hồ tự dự báo lưu lượng đến hồ 6–12h mà không có số liệu mạng lưới). Việc hợp nhất dữ liệu dự báo mưa số trị quy mô vừa (HRM, ETA) vào một công cụ điều hành chung giúp loại bỏ hoàn toàn độ trễ thông tin.
Implications đa chiều
- Về mặt Lý thuyết & Học thuật: Đặt nền móng lý luận cho chuyên ngành Phát triển nguồn nước tại Việt Nam trong việc chuyển từ quy trình tĩnh sang hệ thống động học thời gian thực có xét đến yếu tố bất định của dự báo khí tượng.
- Về mặt Phương pháp luận: Cung cấp bộ công cụ ghép nối mô hình (MOPHONG-LU $\leftrightarrow$ MIKE 11) có khả năng chuyển giao và áp dụng tương thích cho các lưu vực sông có cấu trúc tương tự ở Duyên hải miền Trung (như sông Ba, sông Hương, sông Thu Bồn, sông Trà Khúc).
- Về mặt Thực tiễn & Chính sách: Cung cấp luận cứ khoa học trực tiếp phục vụ các cơ quan quản lý nhà nước (Bộ Nông nghiệp và Phát triển nông thôn, Bộ Tài nguyên và Môi trường, Ủy ban Phòng chống thiên tai) rà soát, sửa đổi và hoàn thiện Quy trình vận hành liên hồ chứa lưu vực sông Vu Gia – Thu Bồn được Thủ tướng Chính phủ phê duyệt.
Limitations và Future Research
- Độ phân giải không gian của mô hình dự báo mưa số trị: Tại thời điểm nghiên cứu (2014), dữ liệu dự báo mưa định lượng từ mô hình HRM và ETA từ Trung tâm DBKTTV Trung ương có độ phân giải lưới tương đối thô, phụ thuộc nhiều vào trường dự báo toàn cầu GFS, dẫn đến sai số cục bộ về vị trí tâm mưa trên địa hình đồi núi phức tạp.
- Giới hạn mô hình hóa 1 chiều vùng hạ lưu: Việc sử dụng mô hình thủy lực 1 chiều MIKE 11 chưa mô phỏng chi tiết được dòng chảy tràn hai chiều (2D Overland Flow) trên các vùng trũng ngập ven đô thị Đà Nẵng và Hội An khi xảy ra hiện tượng vỡ đê cục bộ hoặc ngập sâu do triều cường kết hợp nước dâng do bão.
- Chưa tích hợp yếu tố biến đổi hình thái lòng dẫn: Quá trình bồi xói lòng sông và chuyển tải bùn cát cực mạnh trong mùa lũ chưa được mô phỏng đồng thời trong bài toán thủy động lực học.
Hướng nghiên cứu tương lai:
- Tích hợp dữ liệu Radar thời tiết (Radar Nowcasting) và mạng lưới đo mưa tự động mật độ cao (Vrain) để hiệu chỉnh trường mưa dự báo theo thời gian thực (Data Assimilation) với bước thời gian 1 giờ.
- Nâng cấp mô hình thủy lực hạ lưu lên hệ thống tích hợp 1D-2D (MIKE FLOOD hoặc TELEMAC-2D) nhằm xây dựng bản đồ rủi ro ngập lụt động độ phân giải cao phục vụ sơ tán dân cư.
- Ứng dụng trí tuệ nhân tạo (AI/Deep Learning - LSTM, Transformer) kết hợp tối ưu hóa thích nghi đa mục tiêu (NSGA-III) để tự động hóa hoàn toàn quá trình ra quyết định đóng mở cửa van theo thời gian thực.
Tác động và ảnh hưởng
MA TRẬN TÁC ĐỘNG
+-----------------------------------------------------------------------------------+
| HỌC THUẬT: Định hình khung lý thuyết điều tiết thời gian thực cho sông miền Trung |
| - Mở ra các hướng nghiên cứu kết hợp Khí tượng - Thủy văn - Thủy lực công trình |
+-----------------------------------------------------------------------------------+
|
v
+-----------------------------------------------------------------------------------+
| NGÀNH ĐIỆN & CÔNG TRÌNH: Tối ưu hóa vận hành cho EVN và các Công ty Thủy điện |
| - Giảm thiểu nguy cơ mất an toàn đập; bảo toàn 97.5% sản lượng điện năng mùa kiệt |
+-----------------------------------------------------------------------------------+
|
v
+-----------------------------------------------------------------------------------+
| CHÍNH SÁCH NHÀ NƯỚC: Cơ sở sửa đổi Quy trình vận hành liên hồ chứa Chính phủ |
| - Bổ sung quy định xả lũ an toàn, xóa bỏ cơ chế "quy trình cứng" thiếu linh hoạt |
+-----------------------------------------------------------------------------------+
|
v
+-----------------------------------------------------------------------------------+
| KINH TẾ - XÃ HỘI: Bảo vệ an toàn tính mạng và tài sản vùng hạ du Quảng Nam - ĐN |
| - Giảm ngập lụt 0.10 - 0.60m tại các đô thị trọng điểm, hạ thấp thiệt hại bão lũ |
+-----------------------------------------------------------------------------------+
- Tác động Học thuật: Được trích dẫn rộng rãi trong các nghiên cứu chuyên ngành Thủy văn, Thủy lực và Quản lý tài nguyên nước tại Việt Nam; cung cấp tài liệu tham khảo chuẩn mực cho đào tạo sau đại học ngành Phát triển nguồn nước.
- Tác động Công nghiệp & Vận hành: Giúp các công ty thủy điện (A Vương, Đăk Mi 4, Sông Tranh 2, Sông Bung 4) chuyển đổi từ trạng thái bị động đối phó sang chủ động lập lịch xả lũ khoa học.
- Lợi ích Xã hội Định lượng: Cắt giảm từ $0,10\text{ m}$ đến trên $0,60\text{ m}$ mực nước lũ tại các vùng dân cư trọng điểm thuộc huyện Đại Lộc, Điện Bàn, TP. Hội An (Quảng Nam) và Hòa Vang, Cẩm Lệ (TP. Đà Nẵng), giúp giảm thiểu thiệt hại kinh tế hàng trăm tỷ đồng mỗi năm do thiên tai lũ lụt.
Đối tượng hưởng lợi
- Nghiên cứu sinh & Nhà khoa học chuyên ngành Tài nguyên nước: Thừa hưởng cấu trúc mô hình tích hợp và phương pháp luận giải quyết bài toán mâu thuẫn đa mục tiêu trên sông dốc.
- Các Chủ hồ chứa & Tập đoàn Điện lực Việt Nam (EVN): Sở hữu công cụ tin cậy để vận hành an toàn công trình, phòng ngừa sự cố tràn đập đồng thời tối ưu hóa dung tích tích nước cuối mùa phục vụ phát điện mùa khô.
- Cơ quan Quản lý Quy hoạch & Chính sách: Ban Chỉ đạo Quốc gia về Phòng chống thiên tai, Cục Quản lý Tài nguyên nước sử dụng kết quả để hoàn thiện khung pháp lý và quy chuẩn kỹ thuật quốc gia về vận hành liên hồ chứa.
- Cộng đồng Dân cư Vùng Hạ du: Được bảo vệ an toàn tính mạng và sinh kế trước các đợt lũ quét, lũ đột ngột do xả lũ công trình thủy điện.
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?
Trả lời: Luận án đã mở rộng Lý thuyết Vận hành Thích nghi Hệ thống Hồ chứa Đa mục tiêu (Adaptive Multi-Objective Reservoir Operation Theory) bằng cách thay thế hàm vận hành đơn biến truyền thống $Q_x = f(Z_{\text{lake}})$ bằng hàm điều khiển thời gian thực phụ thuộc trạng thái mở rộng $Q_x = \mathcal{F}(S_t, \widehat{I}{t \to t+k}, \Omega_t)$, trong đó biến đầu vào $\widehat{I}{t \to t+k}$ được dẫn xuất từ trường mưa dự báo số trị hạn vừa 3–5 ngày thay vì dòng chảy quan trắc tĩnh.
2. Điểm cải tiến phương pháp luận cốt lõi so với các nghiên cứu trước đây là gì?
Trả lời: So với các nghiên cứu của Nguyễn Lan Châu (2005) hay Ngô Huy Cẩn (2006) trên sông Hồng vốn sử dụng mô hình tách rời với thời gian dự báo ngắn, luận án đã xây dựng phần mềm tích hợp MOPHONG-LU cho phép liên kết tự động, liền mạch giữa dự báo mưa trung hạn (HRM/ETA) $\to$ mô hình mưa dòng chảy (NAM) $\to$ điều tiết 5 hồ chứa $\to$ diễn toán Muskingum $\to$ mô hình thủy lực hạ du MIKE 11, tạo thành một vòng lặp kín trợ giúp ra quyết định theo thời gian thực cập nhật liên tục mỗi 6–12 giờ.
3. Phát hiện bất ngờ nhất (counter-intuitive finding) được minh chứng bằng dữ liệu là gì?
Trả lời: Trái với quan niệm vận hành thông thường cho rằng muốn chống lũ thì phải tích nước hồ càng sớm càng tốt ngay khi có lũ về, dữ liệu mô phỏng chứng minh rằng: Việc tích nước sớm ở đầu thời kỳ lũ lên làm tăng nguy cơ gây ngập lụt thảm khốc cho hạ du về sau. Phương án vận hành tối ưu là duy trì xả lưu lượng đón lũ an toàn ở đầu giai đoạn lũ lên, tập trung toàn bộ dung tích phòng lũ để cắt đúng đỉnh lũ ở cuối thời kỳ lũ lên và chuyển sang tích nước dâng bình thường trong thời kỳ lũ rút.
4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) hoàn chỉnh không?
Trả lời: Luận án thiết lập hệ thống quy trình chi tiết bao gồm: Bộ thông số thủy văn mô hình NAM cho từng tiểu lưu vực, mạng lưới trắc dọc trắc ngang hiệu chỉnh trong MIKE 11, thuật toán điều tiết cân bằng nước hồ chứa và ma trận kịch bản dung tích phòng lũ ($Z_{\text{đl}}$) kèm các điều kiện biên thủy văn kiểm chứng từ các trận lũ thực tế 1998, 1999, 2007, 2009.
5. Chương trình nghị sự 10 năm tiếp theo (10-Year Research Agenda) của hướng nghiên cứu này cần tập trung vào đâu?
Trả lời: Tập trung vào: (1) Tích hợp đồng hóa dữ liệu vệ tinh/radar thời gian thực (Real-time Data Assimilation); (2) Xây dựng mô hình Digital Twin (Bản sao số) toàn lưu vực sông Vu Gia – Thu Bồn kết hợp AI mô phỏng ngập lụt 2 chiều tức thời; (3) Nghiên cứu cơ chế chi trả dịch vụ môi trường rừng và điều tiết dòng chảy đa mục tiêu trong bối cảnh thị trường điện cạnh tranh.
Kết luận
- Luận án đã thiết lập thành công cơ sở khoa học và thực tiễn cho bài toán vận hành hệ thống hồ chứa theo thời gian thực thời kỳ mùa lũ trên lưu vực sông dốc Vu Gia – Thu Bồn, giải quyết mâu thuẫn cốt lõi giữa phòng lũ hạ du và an toàn tích nước phát điện.
- Xây dựng và hoàn thiện chương trình tính toán tích hợp MOPHONG-LU (kết hợp mưa – dòng chảy, điều tiết hồ, diễn toán lũ) kết nối đồng bộ với mô hình thủy lực hạ lưu MIKE 11, cho phép nâng thời gian dự kiến vận hành phòng lũ lên 3–5 ngày dựa trên các mô hình dự báo mưa số trị HRM, ETA.
- Đề xuất bộ quy tắc vận hành linh hoạt và điều chỉnh mực nước đón lũ ($Z_{\text{đl}}$) cho 5 hồ chứa lớn (A Vương, Đăk Mi 4, Sông Tranh 2, Sông Bung 2, Sông Bung 4), nâng tổng dung tích phòng lũ hữu hiệu lên $199,20\text{ triệu m}^3$.
- Minh chứng định lượng khả năng cắt giảm mực nước ngập lụt hạ du tại Hội Khách ($0,606\text{ m}$), Ái Nghĩa ($0,109\text{ m}$), Giao Thủy ($0,114\text{ m}$) trong khi vẫn đảm bảo đưa mực nước các hồ về mực nước dâng bình thường vào cuối mùa lũ, giảm thiểu tối đa tổn thất điện năng.
- Tạo ra bước đột phá về phương pháp luận quản lý lưu vực sông tại miền Trung, cung cấp cơ sở khoa học vững chắc để các cơ quan thẩm quyền bổ sung, hoàn thiện Quy trình vận hành liên hồ chứa chính thức trên lưu vực sông Vu Gia – Thu Bồn.
