Tổng quan về luận án
Tính cấp thiết và bối cảnh khoa học của công trình nghiên cứu bắt nguồn từ thực trạng biến đổi khí hậu toàn cầu đang làm gia tăng tần suất và cường độ của các hình thái thời tiết cực đoan tại Việt Nam, đặc biệt là dải đất Bắc Trung Bộ (bao gồm 6 tỉnh: Thanh Hóa, Nghệ An, Hà Tĩnh, Quảng Bình, Quảng Trị, Thừa Thiên - Huế). Khu vực này có địa hình hẹp, dốc, bị chia cắt mạnh bởi các khối núi đâm ngang ra biển như dãy Hoàng Mai, Hồng Lĩnh và dải Trường Sơn Bắc, hình thành hệ thống lưu vực sông suối ngắn, độ dốc lớn, thời gian tập trung nước lũ cực nhanh. Theo số liệu thống kê của Tổng cục Thủy lợi (2014), toàn quốc hiện có 6.831 hồ chứa nước với tổng dung tích 49,88 tỷ $m^3$, trong đó vùng Bắc Trung Bộ tập trung mật độ hồ chứa nhỏ rất cao: hồ có dung tích từ $1 \div 3\text{ triệu }m^3$ chiếm 29,6% và hồ từ $0,2 \div 1\text{ triệu }m^3$ chiếm 32,6% tổng số hồ cùng loại của cả nước. Hầu hết các công trình đầu mối này được xây dựng từ thập niên 1970–1980 trong điều kiện khảo sát, thiết kế và thi công hạn chế, kết cấu thân đập chủ yếu là đập đất đắp đồng chất kết hợp tràn xả lũ tự do không có cửa van điều tiết. Trải qua 30–40 năm vận hành, phần lớn hồ đã xuống cấp trầm trọng, mất an toàn đập do bão hòa nước thân đập, xói ngầm và đặc biệt là nguy cơ nước lũ tràn qua đỉnh đập khi gặp mưa lớn vượt tần suất thiết kế.
+------------------------------------------------------------------------------------+
| BỐI CẢNH VÀ NGUY CƠ HỒ CHỨA BẮC TRUNG BỘ |
+------------------------------------------------------------------------------------+
| • 6.831 hồ chứa toàn quốc (49,88 tỷ m3) -> BTB chiếm 29,6% (1-3M m3) & 32,6% (0,2-1M m3) |
| • Đập đất đắp xây dựng 1970-1980, tràn tự do, thiếu số liệu chuỗi quan trắc dài |
| • 1.004/5.039 đập nhỏ gặp sự cố; 59% (10/17 vụ vỡ đập) do mưa lũ tràn đỉnh đập đất |
+------------------------------------------------------------------------------------+
+------------------------------------------------------------------------------------+
| KHOẢNG TRỐNG NGHIÊN CỨU (RESEARCH GAP) |
+------------------------------------------------------------------------------------+
| 1. Thiếu phương pháp luận phân cấp nguy cơ sự cố định lượng qua chỉ số đặc trưng |
| 2. Chưa xác lập quy luật phân bố mưa thời đoạn 1h trong trận mưa 24h liên tục max |
| 3. Tiêu chuẩn thiết kế cũ (C6-77) đánh giá thấp đỉnh lũ và sức chịu tải tràn xả lũ|
+------------------------------------------------------------------------------------+
Khoảng trống nghiên cứu (Research gap) được xác định rõ rệt: Các nghiên cứu an toàn đập truyền thống tại Việt Nam chủ yếu tiếp cận đơn lẻ theo cơ học đất (ổn định thấm, trượt mái) hoặc đánh giá định tính cảm quan, chưa xây dựng được phương pháp luận khoa học phân cấp định lượng mức độ nguy cơ sự cố (NCSC) do mưa lũ cho tập hợp quy mô lớn các hồ chứa nhỏ. Đồng thời, các quy chuẩn thủy văn hiện hành (như Quy phạm C6-77, TCVN 5060-90) còn phụ thuộc vào chỉ tiêu mưa 1 ngày lịch biểu mà bỏ qua quy luật phân bố thời gian thực của lượng mưa 24 giờ liên tục lớn nhất (24h LTLN) và cường độ mưa đỉnh thời đoạn 1 giờ ($I_{1h-max}$).
Nhằm giải quyết bài toán cấp bách này theo tinh thần Chỉ thị số 21/CT-TTg ngày 14/10/2013 của Thủ tướng Chính phủ về việc tăng cường công tác quản lý và xác định thứ tự ưu tiên sửa chữa, nâng cấp an toàn hồ chứa nước, luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật tài nguyên nước (Mã số: 62 58 02 12) của Nghiên cứu sinh Nguyễn Văn Lợi dưới sự hướng dẫn của PGS. Đoàn Doãn Tuấn và PGS. Nguyễn Văn Hoàng tại Viện Khoa học Thủy lợi Việt Nam (2017) đã giải quyết 2 câu hỏi nghiên cứu cốt lõi:
- Câu hỏi nghiên cứu 1 ($RQ_1$): Quy luật phân bố nội tại theo thời gian của trận mưa 24h LTLN và thông số phân bố cường độ mưa thời đoạn 1 giờ ảnh hưởng như thế nào đến quá trình hình thành dòng chảy đỉnh lũ $Q_{max}$ đổ về hồ chứa nhỏ?
- Câu hỏi nghiên cứu 2 ($RQ_2$): Làm thế nào để xây dựng một khung phương pháp luận phân cấp nguy cơ sự cố do mưa lũ dựa trên các chỉ số vật lý - thủy văn phản ánh đầy đủ phương trình cân bằng nước để sàng lọc nhanh, chính xác hàng trăm hồ chứa phục vụ quy hoạch nâng cấp?
Hai giả thuyết nghiên cứu được kiểm chứng tương ứng:
- $H_1$: Cường độ mưa thời đoạn 1 giờ chuẩn hóa trong trận mưa 24h LTLN tại Bắc Trung Bộ không phân bố đều mà tuân theo hàm phân bố lệch chuẩn (Skew-Normal distribution), trong đó độ lệch chuẩn $\sigma$ càng nhỏ thì đỉnh mưa càng tập trung, làm đỉnh lũ dòng chảy $Q_{max}$ tăng đột biến vượt quá sức chịu tải của tràn xả lũ tự do.
- $H_2$: Nguy cơ sự cố tràn đập do mưa lũ của hồ chứa nhỏ được quyết định đồng thời bởi 3 chỉ số không thứ nguyên tương quan chặt chẽ: tỷ số dung tích hồ trên diện tích lưu vực ($K_V = V_{ho}/F_{lv}$), tỷ số diện tích mặt nước hồ trên diện tích lưu vực ($K_S = F_{ho}/F_{lv}$) và tỷ số tải trọng lưu lượng đỉnh lũ thiết kế trên bề rộng tràn xả lũ ($K_Q = Q_{1%}/B_{tran}$).
Phạm vi nghiên cứu tập trung vào hệ thống hồ chứa thủy lợi dung tích nhỏ chia làm hai nhóm ($1 \div 3\text{ triệu }m^3$ và $0,5 \div 1\text{ triệu }m^3$) kết cấu đập đất, tràn tự do tại 3 tỉnh trọng điểm Nghệ An, Hà Tĩnh và Quảng Trị với chuỗi quan trắc khí tượng thủy văn kéo dài từ năm 1956 đến 2012. Ý nghĩa đột phá của luận án là cung cấp cơ sở khoa học định lượng giúp các nhà quản lý chuyển dịch từ bị động ứng phó sang chủ động phân loại, xếp hạng ưu tiên đầu tư gia cố an toàn hồ đập thích ứng với biến đổi khí hậu.
Literature Review và Positioning
Nghiên cứu an toàn hồ đập và phân tích nguy cơ sự cố do mưa lũ là lĩnh vực học thuật thu hút sự quan tâm sâu sắc của các nhà khoa học thủy lợi quốc tế và trong nước. Tổng quan lịch sử thảm họa hồ đập thế giới cho thấy mưa lũ cực đoan gây tràn đỉnh đập đất luôn là nguyên nhân tàn khốc nhất:
- Thảm họa đập Banqiao và Shimantan trên sông Hoài, tỉnh Hà Nam, Trung Quốc năm 1975 do siêu bão Nina tạo lượng mưa $1.060\text{ mm}$ trong 24 giờ đã gây vỡ liên hoàn 62 hồ chứa bậc thang, làm 230.000 người thiệt mạng, phá hủy hơn 100 km đường sắt Bắc Kinh - Quảng Châu và làm ngập 1 triệu hecta đất canh tác.
- Thảm họa đập Machhu-2 tại bang Gujarat, Ấn Độ ngày 11/8/1979 với lưu lượng xả thiết kế chỉ đạt $5.663\text{ m}^3\text{/s}$ trong khi lưu lượng lũ thực tế do mưa bão cực đoan lên tới $16.300\text{ m}^3\text{/s}$, dẫn đến vỡ đập tràn đỉnh đất cuốn trôi thành phố Morvi và cướp đi sinh mạng của hơn 15.000 người.
- Thảm họa đập Malpasset (Pháp, 1959) làm chết 420 người, đập Vajont (Ý, 1963) gây sóng tràn đỉnh làm 1.910 người thiệt mạng, đập Baldwin Hills (Mỹ, 1963) gây ngập đô thị Los Angeles, hay đập đất Trirlian (Nga, 1994) bị tràn đỉnh do cửa xả tê liệt làm 28 người chết.
Về phương pháp luận đánh giá an toàn, Morris et al. (2004, 2012) tại Anh Quốc đã thiết lập khung đánh giá định lượng rủi ro an toàn hồ chứa (Interim Guide to Quantitative Risk Assessment for UK Reservoirs) gồm 4 giai đoạn: nhận diện mối nguy, ước lượng xác suất, đánh giá hậu quả và kiểm soát rủi ro. Tại Hoa Kỳ, Cục Bảo tồn Tài nguyên Thiên nhiên (NRCS - USDA, 2009) thuộc Chương trình Phục hồi Lưu vực (Watershed Rehabilitation Program) đã triển khai đánh giá an toàn và cải tạo cho hàng chục ngàn hồ đập đất chậm lũ, trong đó chuyển dịch triệt để từ việc sử dụng mưa 6 giờ sang sử dụng chuẩn mưa 24h LTLN và mưa nhiều thời đoạn ngắn.
Về lý thuyết tính toán thủy văn lũ cực trị, trường phái xác suất thống kê thế giới dựa trên định lý cực trị Fisher - Tippett (1928), phát triển thành các dạng hàm phân phối xác suất Gumbel Loại I (1941), Fréchet Loại II (1927), Weibull Loại III (1939) và phân phối Giá trị cực trị tổng quát (GEV) ứng dụng xác định Lượng mưa lớn nhất có thể (PMP) theo Hershfield (1961) và Demetris Koutsoyiannis (1998). Tại Liên Xô (cũ) và Đông Âu, hàm phân phối Gamma 3 thông số Kritxky - Menkel và Pearson Loại III (áp dụng phương pháp 3 điểm của Alekseev) được đưa vào tiêu chuẩn quốc gia để tính toán lũ cực đại. Dmytro Stefanyshyn & Julia Stefanyshyna (2003) khi nghiên cứu lưu lượng lũ sông Dnipro (Ukraine) đã chứng minh tính nhạy cảm và độ bất định lớn giữa các đường tần suất lý luận ở dải xác suất vượt hiếm gặp ($P = 0,01% \div 1%$).
Tại Việt Nam, các công trình của Nguyễn Chiến (2013), Phạm Ngọc Quý et al. (2005, 2013–2015), Vũ Đình Hùng (2007) đã phân tích sâu các nguyên nhân gây mất ổn định đập đất như thấm, xói ngầm thân đập, trượt mái và đề xuất xây dựng tràn sự cố để thoát lũ vượt thiết kế. Tuy nhiên, tồn tại 2 luồng quan điểm đối nghịch trong thiết kế thủy văn:
- Quan điểm truyền thống (Quy phạm C6-77): Dựa trên lượng mưa 1 ngày đo tại thời điểm 19h hàng ngày, giả định mưa phân bố đều hoặc theo mô hình khối đơn giản, xem nhẹ diễn biến mưa cường suất lớn thời đoạn ngắn.
- Quan điểm hiện đại thích ứng BĐKH: Khẳng định mưa 1 ngày lịch biểu làm suy giảm từ 15% đến 30% giá trị lượng mưa đỉnh thực tế so với mưa 24h LTLN trượt liên tục, dẫn đến tính toán thiếu hụt dung tích điều tiết và khẩu độ thoát lũ của tràn xả.
Positioning học thuật của luận án: Luận án định vị là công trình tiên phong tại Việt Nam kết nối đồng bộ giữa phân tích thống kê toán học chuỗi mưa 24h LTLN, thiết lập mô hình phân bố cường độ mưa giờ lệch chuẩn, tích hợp với mô hình số trị thủy văn HEC-HMS và hệ thống chỉ số phân cấp ma trận 3 chiều ($K_V, K_S, K_Q$), giải quyết triệt để bài toán sàng lọc nguy cơ sự cố hồ chứa quy mô lưu vực mà các nghiên cứu trước đây chưa lượng hóa được.
+---------------------------------------------------------------------------------------------------------+
| BẢN ĐỒ VỊ TRÍ HỌC THUẬT (ACADEMIC POSITIONING) |
+---------------------------------------------------------------------------------------------------------+
| Hướng nghiên cứu Quốc tế: |
| • Morris et al. (UK, 2004/2012): Khung QRA 4 giai đoạn định lượng rủi ro tổng thể |
| • NRCS-USDA (Mỹ, 2009): Chuyển đổi chuẩn thiết kế sang mưa 24h LTLN & cường độ thời đoạn ngắn |
| • Gumbel, Pearson III, Kritxky-Menkel, Stefanyshyn (2003): Bất định hàm xác suất lũ cực trị |
+---------------------------------------------------------------------------------------------------------+
+---------------------------------------------------------------------------------------------------------+
| Hướng nghiên cứu Việt Nam trước đây: |
| • Nguyễn Chiến (2013), Vũ Đình Hùng (2007): Phân tích nguyên nhân vỡ đập, kiến nghị sửa đổi C6-77 |
| • Phạm Ngọc Quý et al. (2013-2015): Bộ 5 tiêu chí an toàn đập (Lũ, Thấm, Kết cấu, Địa chất, Quản lý) |
| • Khoảng trống: Chưa có phương pháp luận sàng lọc nhanh NCSC định lượng cho hồ nhỏ qua chỉ số cân bằng nước |
+---------------------------------------------------------------------------------------------------------+
+---------------------------------------------------------------------------------------------------------+
| ĐÓNG GÓP ĐỊNH VỊ CỦA LUẬN ÁN (NGUYỄN VĂN LỢI, 2017): |
| • Xác lập tương quan mưa 1 ngày LN vs mưa 24h LTLN; chuẩn hóa phân bố mưa 1h dạng LỆCH CHUẨN |
| • Phát triển hệ chỉ số KV, KS, KQ thiết lập ma trận phân cấp nguy cơ sự cố 3 chiều |
| • Tích hợp mô hình số trị HEC-HMS định lượng sức chịu tải và yêu cầu xả lũ công trình hồ chứa nhỏ |
+---------------------------------------------------------------------------------------------------------+
Đó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 3 bước tiến lý thuyết then chốt cho chuyên ngành Kỹ thuật tài nguyên nước và Thủy văn công trình:
- Mở rộng lý thuyết phân phối xác suất mưa cực trị thời đoạn ngắn: Khẳng định dạng phân bố của cường độ mưa 1 giờ chuẩn hóa trong đợt mưa 24h LTLN tuân theo quy luật phân bố lệch chuẩn (Skew-Normal distribution), đặc trưng bởi phương sai và hệ số bất đối xứng rõ rệt. Phát hiện này bác bỏ giả định phân bố đối xứng Gauss hay phân bố đồng đều thường dùng trong các bài toán diễn toán lũ kinh nghiệm trước đây.
- Thiết lập mệnh đề tương quan phi tuyến giữa cấu trúc mưa và đỉnh lũ: Luận án chứng minh mệnh đề: "Với cùng một tổng lượng mưa 24h LTLN, thời điểm xuất hiện tâm mưa cường độ 1h lớn nhất và giá trị độ lệch chuẩn $\sigma$ của phân bố mưa đóng vai trò quyết định đến độ dốc đường lũ lên và lưu lượng đỉnh lũ $Q_{max}$ đổ về hồ chứa". Khi độ lệch chuẩn $\sigma$ thu hẹp (mưa tập trung cao độ trong $2 \div 3\text{ giờ}$), đỉnh lũ $Q_{max}$ có thể tăng vọt từ 1,5 đến 2,2 lần so với kịch bản mưa phân tán.
- Hình thành khung lý thuyết phân cấp nguy cơ dựa trên nguyên lý cân bằng nước: Xây dựng hệ thống 3 chỉ số vật lý phản ánh trọn vẹn phương trình vi phân cân bằng nước hồ chứa $\frac{dV}{dt} = Q_{den}(t) - Q_{xa}(t)$:
- Chỉ số lưu vực - dung tích ($K_V$): $K_V = \frac{V_{ho}}{F_{lv}}\text{ }(10^3\text{ m}^3\text{/km}^2\text{ hoặc }m^3\text{/m}^2)$, thể hiện dung tích điều tiết trữ lũ đơn vị của hồ trên một đơn vị diện tích lưu vực hứng nước.
- Chỉ số diện tích mặt nước ($K_S$): $K_S = \frac{F_{ho}}{F_{lv}}\text{ }(%)$, biểu thị tỷ lệ diện tích lòng hồ so với diện tích lưu vực, đại diện cho năng lực dung lượng đệm và dung tích truyền sóng lũ mặt hồ.
- Chỉ số tải trọng xả lũ tràn ($K_Q$): $K_Q = \frac{Q_{1%}}{B_{tran}}\text{ }(\text{m}^3\text{/s}\cdot\text{m})$, biểu thị lưu lượng đỉnh lũ thiết kế tần suất $P = 1%$ tập trung từ lưu vực trên một mét chiều rộng tràn xả lũ tự do.
+------------------------------------------------------------------------------------+
| KHUNG LÝ THUYẾT PHÂN CẤP NGUY CƠ SỰ CỐ (NCSC) |
+------------------------------------------------------------------------------------+
| Phương trình cân bằng nước hồ chứa: dV/dt = Q_den(t) - Q_xa(t) |
+------------------------------------------------------------------------------------+
+---------------------+ +---------------------+ +---------------------+
| CHỈ SỐ KV | | CHỈ SỐ KS | | CHỈ SỐ KQ |
| Kv = V_ho / F_lv | | Ks = F_ho / F_lv | | Kq = Q_1% / B_tran |
| Năng lực tích trữ | | Năng lực đệm mặt hồ | | Áp lực xả qua tràn |
+---------------------+ +---------------------+ +---------------------+
+------------------------------------------------------------------------------------+
| MA TRẬN PHÂN CẤP NGUY CƠ ĐA CHIỀU (3 CẤP ĐỘ) |
| • Nhóm Nguy cơ CAO (I) : Kv thấp, Ks thấp, Kq rất cao |
| • Nhóm Nguy cơ VỪA (II) : Các chỉ số trung bình |
| • Nhóm Nguy cơ THẤP (III) : Kv cao, Ks cao, Kq nhỏ |
+------------------------------------------------------------------------------------+
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp liên hoàn 3 lý thuyết nền tảng: Lý thuyết Thống kê thủy văn cực trị (Extreme Value Theory), Lý thuyết Thủy lực công trình điều tiết lũ (Hydraulics of Spillways & Routing) và Lý thuyết Quản lý rủi ro thiên tai (Disaster Risk Management).
Cách tiếp cận phân tích đa chiều được xây dựng thông qua ma trận kết hợp hai chiều ($K_V \times K_S$) và ba chiều ($K_V \times K_S \times K_Q$). Không gian trạng thái nguy cơ được phân chia thành 3 nhóm cấp độ tường minh:
- Nhóm nguy cơ cao (Cấp I): $K_V$ thấp ($< 100 \times 10^3\text{ m}^3\text{/km}^2$), $K_S$ nhỏ ($< 2%$), $K_Q$ lớn ($> 5\text{ m}^3\text{/s}\cdot\text{m}$). Hồ chứa hứng lượng nước lũ khổng lồ từ lưu vực rộng nhưng dung tích hồ quá nhỏ, bề rộng tràn hẹp, nước lũ nhanh chóng dâng vượt cao trình đỉnh đập gây tràn đập đất.
- Nhóm nguy cơ trung bình (Cấp II): Các chỉ số ở mức trung bình, có nguy cơ khi gặp các đợt mưa lũ liên tiếp kéo dài (lũ chồng lũ).
- Nhóm nguy cơ thấp (Cấp III): $K_V$ lớn ($> 300 \times 10^3\text{ m}^3\text{/km}^2$), $K_S > 5%$, $K_Q$ nhỏ, hồ có năng lực tự điều tiết cắt lũ an toàn.
Điều kiện biên áp dụng (Boundary conditions): Khung phân tích áp dụng tối ưu cho các hồ chứa quy mô vừa và nhỏ (dung tích $0,2 \div 3,0\text{ triệu }m^3$, diện tích lưu vực $F_{lv} \le 50\text{ km}^2$), đập đất đồng chất hoặc đập đất nhiều khối, công trình xả lũ là tràn tự do chảy thẳng hoặc tràn bên không có cửa van điều tiết cơ khí.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Nghiên cứu được thiết kế theo trường phái thực chứng định lượng (Positivism / Quantitative Paradigm), kết hợp phương pháp mô hình hóa số trị và xác thực thực nghiệm từ dữ liệu lịch sử sự cố đập. Thiết kế nghiên cứu đa tầng (Multi-level research design) bao gồm 3 cấp độ:
- Cấp độ 1 (Vùng khí hậu Bắc Trung Bộ): Phân tích thống kê toán học chuỗi số liệu mưa trạm khí tượng thủy văn đại diện.
- Cấp độ 2 (Hệ thống hồ chứa cấp tỉnh): Áp dụng phương pháp luận ma trận chỉ số để phân cấp nguy cơ sự cố cho toàn bộ 627 hồ chứa tại Nghệ An, 341 hồ tại Hà Tĩnh và 181 hồ tại Quảng Trị.
- Cấp độ 3 (Công trình hồ chứa điển hình - Hồ Khe Nu, Nghệ An): Mô phỏng thủy văn - thủy lực chi tiết bằng phần mềm mô hình HEC-HMS nhằm giải mã cơ chế lan truyền sóng lũ và sức chịu tải của tràn xả.
+------------------------------------------------------------------------------------+
| QUY TRÌNH PHƯƠNG PHÁP NGHIÊN CỨU |
+------------------------------------------------------------------------------------+
| BƯỚC 1: Thu thập số liệu KTTV trạm Vinh, Quán Hành, Hương Khê, Đông Hà (1956-2012)|
| Xây dựng đường tần suất P% mưa 1 ngày LN & mưa 24h LTLN |
+------------------------------------------------------------------------------------+
+------------------------------------------------------------------------------------+
| BƯỚC 2: Chuẩn hóa đợt mưa 24h LTLN, thiết lập hàm phân bố lệch chuẩn thời đoạn 1h |
| Xác định tham số vị trí, tỷ lệ và độ lệch chuẩn σ |
+------------------------------------------------------------------------------------+
+------------------------------------------------------------------------------------+
| BƯỚC 3: Tính toán hệ chỉ số Kv, Ks, Kq cho tập hợp hồ chứa nhỏ Nghệ An, Hà Tĩnh, QT|
| Thiết lập ma trận phân cấp mức độ nguy cơ sự cố (NCSC) |
+------------------------------------------------------------------------------------+
+------------------------------------------------------------------------------------+
| BƯỚC 4: Mô phỏng số trị HEC-HMS hồ Khe Nu (6 kịch bản σ của phân bố mưa lệch chuẩn)|
| Định lượng đỉnh lũ Q_max, mực nước lũ lớn nhất và khẩu diện xả tràn tối ưu |
+------------------------------------------------------------------------------------+
Quy trình nghiên cứu rigorous
Quy trình nghiên cứu tuân thủ nghiêm ngặt 4 chuẩn mực khoa học:
- Chiến lược lấy mẫu và thu thập dữ liệu: Thu thập chuỗi dữ liệu mưa ngày và mưa đo mưa tự ghi thời đoạn ngắn tại các trạm khí tượng thủy văn tiêu biểu (Trạm Vinh, Nghi Lộc - Nghệ An; Trạm Hương Khê - Hà Tĩnh; Trạm Đông Hà, Hải Lăng - Quảng Trị) từ năm 1956 đến 2012 (chuỗi số liệu trên 50 năm bảo đảm tính đại diện thống kê). Thu thập hồ sơ thiết kế kỹ thuật, bản đồ lòng hồ, đường cong dung tích $V = f(Z)$, đường cong diện tích mặt nước $F = f(Z)$ và kích thước hình học tràn xả lũ của các hồ chứa thuộc phạm vi nghiên cứu.
- Xử lý thống kê và kiểm định phân phối xác suất: Thiết lập đường tần suất lý luận lượng mưa cực trị theo phân phối Pearson Loại III, Gumbel và Kritxky - Menkel; xác định các tham số thống kê mô tả: lượng mưa trung bình $\bar{X}$, hệ số biến sai $C_v$, hệ số thiên lệch $C_s$. Sử dụng phương pháp 3 điểm của Alekseev để hiệu chỉnh đường tần suất kinh nghiệm (công thức Hazen, Tregođaev).
- Kỹ thuật tam giác đạc (Triangulation): Đối soát chéo giữa kết quả tính toán phân cấp lý thuyết ma trận ($K_V, K_S, K_Q$) với hiện trạng sự cố thực tế của 17 hồ đập bị vỡ hoặc sự cố tràn đỉnh trong lịch sử (như hồ Tây Nguyên vỡ năm 2012, hồ Khe Mơ vỡ năm 2010, đập Đakrông 3 vỡ năm 2012, đập Khe Tuần tràn vỡ cống năm 2013, hồ Hoàng Tân vỡ năm 2013).
- Độ tin cậy và giá trị khoa học: Đánh giá độ nhạy của mô hình thủy văn mưa - dòng chảy với các biến số độ che phủ rừng, nhóm đất thủy văn và hệ số nhám lòng dẫn.
Data và phân tích
Phần mềm mô hình HEC-HMS (Hydrologic Engineering Center - Hydrologic Modeling System) của Công binh Lục quân Hoa Kỳ được sử dụng làm công cụ mô phỏng số trị chính:
- Mô hình tổn thất thấm: Áp dụng phương pháp số đường cong dòng chảy SCS (NRCS Curve Number - CN) dựa trên đặc tính thổ nhưỡng và thảm phủ lưu vực hồ Khe Nu.
- Mô hình chuyển đổi dòng chảy mặt: Sử dụng phương pháp Biểu đồ đơn vị tổng hợp SCS (SCS Unit Hydrograph Method).
- Mô hình diễn toán lũ qua hồ chứa: Ứng dụng phương pháp diễn toán hồ chứa cải biên (Modified Puls Routing Method) giải phương trình vi phân liên tục kết hợp đường cong đặc tính công trình tràn xả lũ $Q_{xa} = m \cdot B_{tran} \cdot \sqrt{2g} \cdot H_o^{3/2}$.
Phân tích độ bền vững (Robustness checks) được thực hiện bằng cách thay đổi diện tích lưu vực giả định từ $0,5\text{ km}^2$ đến $10\text{ km}^2$ và biến thiên 6 kịch bản giá trị độ lệch chuẩn của hàm phân bố mưa lệch chuẩn ($\sigma = 1,5; 2,0; 2,5; 3,0; 3,5; 4,0$) để kiểm tra phản ứng của đỉnh lũ $Q_{den-max}$, mực nước dâng gia cường cực đại ($Z_{max}$) và lưu lượng tràn đỉnh.
Phát hiện đột phá và implications
Những phát hiện then chốt
Nghiên cứu mang lại 5 phát hiện đột phá có cơ sở định lượng thực nghiệm vững chắc:
+------------------------------------------------------------------------------------+
| 5 PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN |
+------------------------------------------------------------------------------------+
| 1. Lượng mưa 24h LTLN lớn hơn mưa 1 ngày lịch biểu từ 15% đến 28% |
| 2. Cường độ mưa 1h trong đợt 24h LTLN tuân theo phân bố LỆCH CHUẨN rõ rệt |
| 3. Tâm mưa lệch chuẩn (σ hẹp) làm đỉnh lũ Q_max tăng từ 35% đến 65% |
| 4. 100% hồ chứa từng bị vỡ đập rơi chính xác vào ô NGUY CƠ CAO của ma trận Kv-Ks-Kq |
| 5. Tồn tại ngưỡng tới hạn của diện tích lưu vực F_lv đối với sức tải tràn tự do |
+------------------------------------------------------------------------------------+
Thứ nhất: Sự chênh lệch có ý nghĩa thống kê giữa lượng mưa 24h LTLN và mưa 1 ngày lớn nhất. Luận án chứng minh lượng mưa 24h liên tục lớn nhất ($W_{24h-LTLN}$) tại tất cả các trạm vùng Bắc Trung Bộ luôn lớn hơn lượng mưa 1 ngày lớn nhất ($W_{1ngay-LN}$) cùng tần suất thiết kế từ 15% đến 28%. Tại trạm Quán Hành (Nghệ An), ứng với tần suất $P = 1%$, mưa 24h LTLN đạt $585\text{ mm}$ trong khi mưa 1 ngày chỉ đạt $470\text{ mm}$; tại Hương Khê (Hà Tĩnh) mưa 24h LTLN đạt $630\text{ mm}$ so với mưa 1 ngày $510\text{ mm}$; tại Đông Hà (Quảng Trị) mưa 24h LTLN đạt $595\text{ mm}$ so với $480\text{ mm}$. Việc sử dụng số liệu mưa 1 ngày truyền thống theo quy chuẩn cũ đã làm thiếu hụt nghiêm trọng tổng lượng dòng chảy lũ thiết kế.
Thứ hai: Quy luật phân bố lệch chuẩn của cường độ mưa thời đoạn 1 giờ. Khi phân tích các trận mưa lũ lịch sử (tháng 10/1991, tháng 10/2010, tháng 9–10/2013), cường độ mưa 1 giờ chuẩn hóa phân bố lệch hẳn về phía sau hoặc phía trước của trận mưa với dạng phân phối Skew-Normal, tạo ra tâm mưa cực đoan có cường suất mưa giờ lớn nhất ($I_{1h-max}$) chiếm tới 25–35% tổng lượng mưa 24 giờ.
Thứ ba: Cơ chế khuếch đại đỉnh lũ do biến thiên độ lệch chuẩn phân bố mưa. Kết quả mô phỏng HEC-HMS trên lưu vực hồ Khe Nu cho thấy: khi tổng lượng mưa 24h không đổi ($W_{24h} = 450\text{ mm}$), nếu độ lệch chuẩn $\sigma$ của phân bố mưa giảm từ 4,0 xuống 1,5 (mưa tập trung cường độ cực đại trong thời gian ngắn), lưu lượng đỉnh lũ đổ về hồ $Q_{den-max}$ tăng từ $34,2\text{ m}^3\text{/s}$ lên $58,6\text{ m}^3\text{/s}$ (tăng 71,3%), làm mực nước hồ vượt cao trình đỉnh đập đất $0,45\text{ m}$, gây tràn đỉnh đập.
"Trong các nguyên nhân hư hỏng đập của hồ chứa thì nguyên nhân hư hỏng do lũ hay tràn đập chiếm tỷ lệ lớn nhất... Trong số 17 lần vỡ đập đã xảy ra ở nước ta, đập đất bị vỡ do nguyên nhân trực tiếp là mưa lũ có 10 đập, chiếm 59%" (Trích dẫn nguồn văn bản luận án).
Thứ tư: Độ chính xác tuyệt đối của Ma trận phân cấp nguy cơ sự cố ($K_V, K_S, K_Q$). Khi định vị các hồ chứa từng xảy ra thảm họa vỡ đập tại miền Trung vào ma trận phân cấp của luận án, 100% các công trình (Đập Tây Nguyên, Khe Mơ, Đồng Đáng, Khe Tuần, Đakrông 3) đều nằm chính xác trong vùng "Nguy cơ sự cố rất cao" (Nhóm I), nơi có $K_V < 120 \times 10^3\text{ m}^3\text{/km}^2$ và $K_Q > 4,5\text{ m}^3\text{/s}\cdot\text{m}$.
Thứ năm: Xác lập ngưỡng giới hạn sức chịu tải công trình theo diện tích lưu vực. Với các hồ nhỏ dung tích $0,5 \div 1,0\text{ triệu }m^3$, đập đất cao dưới 15 m, tràn tự do bề rộng $B_{tran} \le 10\text{ m}$, nếu diện tích lưu vực $F_{lv} > 8\text{ km}^2$ thì nguy cơ tràn đỉnh đập khi xuất hiện mưa 24h LTLN tần suất $P = 1%$ là trên 90% nếu không có tràn sự cố hỗ trợ.
+---------------------------------------------------------------------------------------------------------+
| BẢNG ĐỐI SOÁT KIỂM CHỨNG MA TRẬN NGUY CƠ VỚI CÁC SỰ CỐ VỠ ĐẬP THỰC TẾ |
+---------------------+-------------+--------------+-----------------+----------------+-------------------+
| Tên hồ chứa sự cố | Tỉnh | Năm xảy ra | Chỉ số Kv | Chỉ số Kq | Phân cấp ma trận |
+---------------------+-------------+--------------+-----------------+----------------+-------------------+
| Hồ Quán Hài | Nghệ An | 1978 (vỡ) | Rất thấp | Rất cao | Nguy cơ rất cao |
| Hồ Đồn Húng | Nghệ An | 1978 (vỡ) | Thấp | Cao | Nguy cơ cao |
| Đập Z20 | Hà Tĩnh | 2009 (vỡ) | Cực thấp | Rất cao | Nguy cơ rất cao |
| Hồ Khe Mơ | Hà Tĩnh | 2010 (vỡ) | Thấp (< 100) | Rất cao (> 5) | Nguy cơ rất cao |
| Hồ Tây Nguyên | Nghệ An | 2012 (vỡ) | Thấp | Cao | Nguy cơ cao |
| Đập Đakrông 3 | Quảng Trị | 2012 (vỡ) | Rất thấp | Cực cao | Nguy cơ rất cao |
| Hồ Đồng Đáng | Thanh Hóa | 2013 (vỡ) | Thấp | Cao | Nguy cơ cao |
| Hồ Khe Tuần | Thanh Hóa | 2013 (tràn) | Thấp (< 110) | Cao | Nguy cơ cao |
+---------------------+-------------+--------------+-----------------+----------------+-------------------+
Implications đa chiều
- Về mặt học thuật và lý thuyết: Cung cấp một phương pháp luận toán - thủy văn hoàn chỉnh để phân tích mưa lũ cực trị cho lưu vực nhỏ không có trạm quan trắc dòng chảy; bổ sung quy luật phân bố lệch chuẩn thời đoạn ngắn vào giáo trình giảng dạy kỹ thuật tài nguyên nước.
- Về mặt phương pháp luận: Thiết lập quy trình chuẩn hóa chuỗi mưa 24h LTLN và ma trận sàng lọc nhanh rủi ro thiên tai, có thể nhân rộng chuyển giao cho các vùng địa hình tương đồng như Duyên hải Nam Trung Bộ, Tây Nguyên và Trung du Miền núi phía Bắc.
- Về ứng dụng thực tiễn công trình: Đưa ra công thức xác định khẩu diện tràn xả lũ mở rộng và giải pháp bố trí tràn sự cố (Emergency spillway) kiểu phím piano (PK-Weir) hoặc tràn ngập tự vỡ bằng vật liệu địa phương nhằm bảo vệ an toàn tuyệt đối cho đập chính khi gặp mưa lũ vượt tần suất thiết kế.
- Về chính sách công và quy chuẩn quốc gia: Cung cấp cơ sở khoa học đắt giá để Bộ Nông nghiệp và Phát triển nông thôn, Bộ Khoa học và Công nghệ sửa đổi, thay thế Tiêu chuẩn thiết kế thủy văn công trình thủy lợi (TCVN 8414-2010, TCVN 5060-90) theo hướng bắt buộc áp dụng chuỗi mưa 24h LTLN và kịch bản phân bố mưa giờ cực đoan thích ứng BĐKH.
Limitations và Future Research
Luận án thẳng thắn chỉ ra 4 giới hạn khoa học cần tiếp tục hoàn thiện:
- Giới hạn chuỗi quan trắc mưa tự ghi: Số lượng trạm đo mưa tự ghi đo chi tiết theo từng phút/giờ tại vùng Bắc Trung Bộ còn thưa thớt; phần lớn dữ liệu lịch sử trước năm 1990 là số liệu đo thủ công ngày 2 lần (7h và 19h), dẫn đến việc ngoại suy phân bố mưa thời đoạn cực ngắn ($15\text{ phút}, 30\text{ phút}, 45\text{ phút}$) phải sử dụng phương pháp tương tự khí hậu.
- Phạm vi hình thức kết cấu công trình: Phương pháp luận ma trận chỉ số $K_V, K_S, K_Q$ mới tập trung chủ yếu vào loại hình đập đất đồng chất quy mô nhỏ có tràn xả lũ tự do; chưa bao hàm các hồ chứa lớn có tràn xả sâu kiểm soát bằng cửa van cung hoặc đập bê tông trọng lực.
- Đặc trưng thảm phủ động: Mô hình số trị HEC-HMS trong luận án sử dụng chỉ số số đường cong dòng chảy (CN) cố định theo hiện trạng, chưa mô phỏng động học quá trình suy thoái hoặc phục hồi thảm phủ rừng đầu nguồn theo thời gian nhiều thập kỷ.
Chương trình nghiên cứu tương lai (Future research agenda):
- Mở rộng phân tích quy luật phân bố mưa cực trị ở các thời đoạn siêu ngắn ($15 \div 45\text{ phút}$) ứng dụng công nghệ radar thời tiết Doppler và ảnh vệ tinh phân giải cao.
- Phát triển mô hình tích hợp 2 chiều (2D Hydrodynamic modeling) mô phỏng vết sóng vỡ đập, ngập lụt hạ du và đánh giá rủi ro sinh mạng, thiệt hại kinh tế phục vụ bản đồ số hóa cảnh báo thiên tai.
- Nghiên cứu cơ chế lan truyền xói ngầm phi tuyến trong thân đập đất dưới tác động kết hợp của áp lực cột nước lũ dâng nhanh và sóng gió mặt hồ.
Tác động và ảnh hưởng
Luận án tạo ra tác động sâu rộng trên cả 4 trụ cột học thuật, kinh tế, xã hội và quản lý nhà nước:
- Tác động học thuật: Định hình chuẩn mực mới trong nghiên cứu an toàn hồ đập tại Việt Nam; cung cấp tài liệu tham khảo mẫu mực cho các nghiên cứu sinh, học viên cao học chuyên ngành Thủy văn, Kỹ thuật tài nguyên nước và Xây dựng công trình thủy.
- Tác động ngành và chuyển đổi kỹ thuật: Đóng vai trò là cẩm nang kỹ thuật cho các viện thiết kế thủy lợi, các công ty tư vấn đầu tư xây dựng khi lập dự án sửa chữa, nâng cấp hồ chứa thuộc Dự án Quản lý an toàn đập (WB8 do Ngân hàng Thế giới tài trợ).
- Tác động chính sách vĩ mô: Đóng góp luận cứ kỹ thuật trực tiếp cho Tổng cục Thủy lợi (nay là Cục Thủy lợi - Bộ NN&PTNT) trong việc rà soát, xếp loại danh mục ưu tiên cấp bách hàng ngàn hồ chứa mất an toàn theo Chỉ thị 21/CT-TTg của Thủ tướng Chính phủ.
- Lợi ích xã hội và môi trường: Giúp giảm thiểu tối đa nguy cơ vỡ đập bất ngờ trong mùa mưa bão, bảo vệ tính mạng cho hàng triệu người dân hạ du các tỉnh Bắc Trung Bộ, bảo vệ kết cấu hạ tầng giao thông huyết mạch (Quốc lộ 1A, Đường sắt Bắc - Nam) và ổn định sản xuất nông nghiệp bền vững.
Đối tượng hưởng lợi
+------------------------------------------------------------------------------------+
| BẢN ĐỒ CÁC ĐỐI TƯỢNG HƯỞNG LỢI |
+------------------------------------------------------------------------------------+
| 1. Nhà nghiên cứu & NCS Thủy văn: Thừa hưởng mô hình phân bố mưa lệch chuẩn |
| 2. Kỹ sư tư vấn & Thiết kế thủy lợi: Công thức tính toán thủy lực tràn chuẩn xác |
| 3. Cơ quan Quản lý Nhà nước (Cục Thủy lợi): Ma trận định lượng xếp hạng an toàn đập|
| 4. Đơn vị Quản lý Khai thác Hồ chứa: Quy trình cảnh báo sớm đón lũ cực đoan |
| 5. Cộng đồng dân cư hạ lưu: An toàn tính mạng & sinh kế nông nghiệp bền vững |
+------------------------------------------------------------------------------------+
- Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận phân tích thống kê thủy văn tiên tiến kết hợp mô hình số trị hiện đại; mở ra hướng nghiên cứu mới về phân bố mưa cực đoan thích ứng BĐKH.
- Cơ quan Quản lý Nhà nước (Bộ NN&PTNT, Sở NN&PTNT các tỉnh Bắc Trung Bộ): Sở hữu công cụ sàng lọc, phân nhóm nguy cơ sự cố nhanh chóng, minh bạch và khoa học để phân bổ nguồn vốn ngân sách đầu tư công hiệu quả, tránh lãng phí.
- Đơn vị Tư vấn Thiết kế và Quy hoạch Thủy lợi: Ứng dụng quy trình tính toán mưa 24h LTLN và thông số mưa lệch chuẩn để chuẩn hóa khẩu độ tràn xả lũ, thiết kế bổ sung tràn sự cố đúng kỹ thuật.
- Công ty TNHH MTV Khai thác Công trình Thủy lợi địa phương: Nắm bắt được ngưỡng dung tích và sức chịu tải của từng hồ để xây dựng quy trình vận hành điều tiết đón lũ an toàn, chủ động hạ thấp mực nước hồ trước khi bão đổ bộ.
- Cộng đồng dân cư vùng hạ lưu: Được bảo đảm an toàn tính mạng, nhà cửa, tài sản và diện tích cây trồng nông nghiệp trước hiểm họa vỡ đập bất ngờ.
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 xác lập hàm phân bố lệch chuẩn (Skew-Normal distribution) cho cường độ mưa thời đoạn 1 giờ trong trận mưa 24h LTLN, đồng thời chứng minh mối liên hệ tương quan giữa độ lệch chuẩn $\sigma$ với sự gia tăng đột biến của đỉnh lũ $Q_{max}$. Nghiên cứu đã mở rộng lý thuyết thủy văn cực trị truyền thống (vốn giả định mưa phân bố đều hoặc chuẩn tắc đối xứng) sang lý thuyết phân bố mưa bất đối xứng thực tế thích ứng với hình thái mưa bão ven biển miền Trung.
2. Tính đổi mới về phương pháp luận nghiên cứu so với các công trình quốc tế và trong nước?
So với khung đánh giá rủi ro định tính truyền thống hoặc phương pháp QRA phức tạp đòi hỏi chuỗi số liệu vi mô của Morris et al. (UK), luận án đã sáng tạo ra hệ chỉ số vật lý không thứ nguyên ($K_V, K_S, K_Q$) phản ánh trực tiếp phương trình cân bằng nước hồ chứa. Phương pháp luận này cho phép sàng lọc, phân loại nguy cơ sự cố đồng loạt cho hàng trăm hồ chứa trong điều kiện dữ liệu đo đạc thực địa còn hạn chế, đạt độ tin cậy tương thích tuyệt đối khi kiểm chứng với 100% các vụ vỡ đập lịch sử.
3. Phát hiện bất ngờ nhất (counter-intuitive) có dữ liệu minh chứng từ luận án?
Phát hiện bất ngờ nhất là: "Tổng lượng mưa 24h không phải là yếu tố duy nhất quyết định sự an toàn của hồ chứa nhỏ, mà chính độ lệch chuẩn $\sigma$ của phân bố mưa giờ mới là nhân tố kích hoạt vỡ đập". Cụ thể, cùng một lượng mưa $450\text{ mm}$, nếu mưa phân bố đều trong 24 giờ thì hồ Khe Nu xả lũ an toàn qua tràn tự do; nhưng khi mưa tập trung theo phân bố lệch chuẩn với $\sigma = 1,5$ (tâm mưa đạt cường suất cực đại), lưu lượng đỉnh lũ tăng vọt 71,3%, nước tràn đỉnh đập đất $0,45\text{ m}$ gây nguy cơ vỡ đập hoàn toàn.
4. Luận án có cung cấp quy trình tái lập (replication protocol) đầy đủ không?
Luận án cung cấp quy trình tái lập 8 bước vô cùng chi tiết: từ công thức toán học xác định các chỉ số $K_V, K_S, K_Q$, bảng tham số phân bố lệch chuẩn mưa giờ tại các trạm khí tượng, đến cấu trúc thiết lập thông số lưu vực (CN, thời gian trễ $T_{lag}$, biểu đồ đơn vị SCS) và phương pháp diễn toán dòng chảy trong HEC-HMS, giúp các kỹ sư và nhà nghiên cứu khác có thể áp dụng chính xác cho bất kỳ lưu vực hồ chứa nào.
5. Định hướng chương trình nghiên cứu 10 năm tiếp theo từ nền tảng luận án?
Chương trình 10 năm tới tập trung vào 3 trụ cột: (1) Tích hợp dữ liệu lớn (Big Data) quan trắc mưa tự động thời gian thực và trí tuệ nhân tạo (AI/Machine Learning) để dự báo sớm nguy cơ tràn đỉnh đập trước $6 \div 12\text{ giờ}$; (2) Chuẩn hóa tiêu chuẩn quốc gia về thiết kế tràn sự cố kiểu mới (PK-Weir, Fusegate) cho toàn bộ hệ thống hồ đập nhỏ Việt Nam; (3) Xây dựng hệ thống bản đồ số viễn thám giám sát sức khỏe đập đất thích ứng với các kịch bản BĐKH cực đoan của IPCC.
Kết luận
- Luận án đã giải quyết trọn vẹn và xuất sắc bài toán khoa học - thực tiễn cấp bách về đánh giá nguy cơ sự cố do mưa lũ cho hệ thống hồ chứa nhỏ vùng Bắc Trung Bộ, đóng góp luận cứ trực tiếp thực hiện Chỉ thị số 21/CT-TTg của Thủ tướng Chính phủ.
- Xác lập mối tương quan khoa học giữa mưa 1 ngày lớn nhất và mưa 24h liên tục lớn nhất, chứng minh mưa 1 ngày lịch biểu làm thiếu hụt từ 15% đến 28% lượng mưa tính toán thiết kế công trình xả lũ.
- Khám phá và mô hình hóa thành công quy luật phân bố lệch chuẩn của cường độ mưa thời đoạn 1 giờ trong trận mưa 24h LTLN, làm sáng tỏ cơ chế khuếch đại dòng chảy đỉnh lũ đổ về hồ chứa nhỏ.
- Sáng tạo phương pháp luận phân cấp nguy cơ sự cố độc đáo thông qua bộ 3 chỉ số vật lý cân bằng nước ($K_V, K_S, K_Q$) và ma trận phân cấp 3 chiều, đạt độ chính xác kiểm chứng 100% với các thảm họa vỡ đập thực tế.
- Ứng dụng thành công mô hình số trị HEC-HMS định lượng hóa sức chịu tải của công trình xả lũ hồ Khe Nu và mở rộng cho mạng lưới hồ chứa tại Nghệ An, Hà Tĩnh, Quảng Trị; kiến nghị giải pháp mở rộng tràn và bố trí tràn sự cố nâng cao an toàn đập đất.
- Mở ra các hướng nghiên cứu hiện đại về thủy văn lưu vực nhỏ thời đoạn ngắn thích ứng BĐKH, tạo nền tảng chuyển đổi số trong quản trị rủi ro an toàn đập quốc gia và hội nhập sâu rộng với các tiêu chuẩn an toàn đập quốc tế.
