Tổng quan về luận án

Công trình hồ chứa nước đóng vai trò huyết mạch trong phát triển kinh tế - xã hội, bảo đảm an ninh nguồn nước, cắt giảm lũ, cấp nước sinh hoạt, công nghiệp và tưới tiêu nông nghiệp. Theo thống kê của Bộ Nông nghiệp và Phát triển nông thôn (2020), cả nước hiện có 6.750 hồ chứa thủy lợi (gồm 4 hồ quan trọng đặc biệt, 888 hồ lớn, 1.633 hồ vừa, 4.225 hồ nhỏ với tổng dung tích trữ khoảng 14,5 tỷ m³) phục vụ tưới cho 1,1 triệu ha đất canh tác và cấp 1,5 tỷ m³ nước. Tuy nhiên, phần lớn các công trình được xây dựng từ nhiều thập kỷ trước trong điều kiện kỹ thuật và kinh phí hạn chế; hiện có khoảng 1.200 hồ chứa bị hư hỏng, xuống cấp, thiếu khả năng xả lũ và 200 hồ hư hỏng nặng cần xử lý cấp bách. Tác động của biến đổi khí hậu gây mưa lũ cực đoan đã dẫn tới 71 sự cố đập trong giai đoạn 2010–2020 (điển hình như vỡ đập Phân Lân năm 2013 do nước tràn đỉnh 1,5 m; vỡ đập phụ hồ Đầm Hà Động năm 2014 do kẹt cửa van; sự cố hồ Núi Cốc với đỉnh lũ tính toán lại tăng 1,85 lần; sự cố vỡ đập Đầm Thìn ngày 28/5/2020).

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực trạng: công tác đánh giá an toàn hồ chứa tại Việt Nam từ trước đến nay chủ yếu dựa trên phương pháp tất định (deterministic method) với hệ số an toàn tĩnh $[K]$, $[J]$ trên các mặt cắt đại diện đơn lẻ, hoàn toàn bỏ qua tính bất định ngẫu nhiên của tải trọng thủy văn, tính chất cơ lý vật liệu, sự tương tác đa cơ chế trong hệ thống công trình đầu mối (CTĐM) và tách rời với rủi ro ngập lụt vùng hạ du.

Luận án tiến sĩ kỹ thuật của nghiên cứu sinh Cầm Thị Lan Hương (2020) với đề tài "Nghiên cứu ứng dụng lý thuyết độ tin cậy và phân tích rủi ro trong đánh giá an toàn hồ chứa thủy lợi Việt Nam" (Chuyên ngành: Kỹ thuật xây dựng công trình thủy, Mã số: 9580202; Người hướng dẫn: NGND.TS. Phạm Ngọc Quý, PGS.TS. Mai Văn Công) đã tiên phong thiết lập khung lý thuyết tích hợp giữa Lý thuyết độ tin cậy (LTĐTC - Structural Reliability Theory) cấp độ II, cấp độ III và Phân tích rủi ro (PTRR - Quantitative Risk Analysis). Luận án giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 3 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Làm thế nào để định lượng chính xác xác suất sự cố ($P_f$) và chỉ số độ tin cậy ($\beta$) của toàn bộ hệ thống CTĐM phức hợp (đập đất, tràn xả lũ, cống lấy nước) thay vì từng hạng mục riêng rẽ? -> H1: Hệ thống CTĐM vận hành theo sơ đồ nối tiếp (series system) hoặc hỗn hợp có thể mô hình hóa thông qua sơ đồ cây sự cố (Fault Tree Analysis - FTA) và giải tích xác suất phi tuyến.
  • RQ2: Mối quan hệ tương quan định lượng giữa mức độ an toàn CTĐM và rủi ro ngập lụt hạ du được xác lập như thế nào? -> H2: Độ tin cậy yêu cầu $[P_f]$ tối ưu về mặt kinh tế - kỹ thuật được xác định tại điểm cực tiểu của hàm tổng chi phí vòng đời: $C_{tot} = I(P_f) + PV(M(P_f)) + PV(P_f \cdot D(P_f)) \to \min$.
  • RQ3: Làm thế nào để phân bổ độ tin cậy yêu cầu tổng thể cho từng cơ chế sự cố thành phần nhằm tối ưu hóa giải pháp thiết kế sửa chữa, nâng cấp? -> H3: Thuật toán phân bổ độ tin cậy ngược từ mức chấp nhận rủi ro hạ du về từng bộ phận CTĐM cho phép lựa chọn quy mô công trình chính xác, tránh lãng phí vốn đầu tư.

Phạm vi nghiên cứu tập trung vào cụm CTĐM hồ chứa thủy lợi có đập dâng bằng đất đang khai thác tại Việt Nam và vùng hạ du chịu ảnh hưởng trực tiếp, lấy hồ Núi Cốc (tỉnh Thái Nguyên) làm địa bàn ứng dụng điển hình với chuỗi số liệu quan trắc mực nước, mô phỏng thủy lực ngập lụt 1D/2D bằng MIKE 11 và MIKE FLOOD.

Literature Review và Positioning

Lý thuyết độ tin cậy và phân tích rủi ro công trình thủy công đã trải qua quá trình phát triển mạnh mẽ trên thế giới từ thập niên 1970 và hoàn thiện khung tiêu chuẩn từ thập niên 1990. Hiệp hội Tiêu chuẩn Canada đã ban hành hướng dẫn CSA Q636-91 và Q636-93 về phân tích rủi ro và lý thuyết độ tin cậy. Các nghiên cứu nền tảng quốc tế tiêu biểu bao gồm:

  1. Yegian et al. (1991) phân tích xác suất an toàn của đập đất chịu tải trọng động đất; Chen Zhaohe (1996) tại Trung Quốc nghiên cứu độ tin cậy chống lũ và tràn đỉnh đập; Marie Westberg (2007) tại Thụy Điển đánh giá độ tin cậy trượt và lật của đập bê tông trọng lực dựa trên chuỗi biến ngẫu nhiên thực nghiệm; Negede Abate Kassa (2009) tại Đức giải hàm tin cậy cấp độ II và III cho 3 cơ chế sự cố đập đất (trượt mái hạ lưu, tràn đỉnh, thấm); Sivakumar Babu & Amit Srivastava (2010) tại Ấn Độ áp dụng mô phỏng Monte Carlo (MCS) đánh giá ổn định trượt mái đập đất dưới tác động động đất; Laurent Peyras et al. (2012) tại Pháp kết hợp LTĐTC và PTRR cho đập bê tông đầm lăn với cơ chế kẹt cửa van và mài mòn.
  2. Các tổ chức quốc tế như Hội Đập lớn Thế giới (ICOLD, Bulletins 2005, 2011, 2013, 2015), Cục Khai hoang Hoa Kỳ (USBR - Cyganiewicz, 2000), Hội Đập lớn Úc (ANCOLD, 2003) đã chuẩn hóa việc định lượng rủi ro an toàn đập gắn liền với phân tích hậu quả ngập lụt hạ du.

Tại Việt Nam, các nghiên cứu của Nguyễn Văn Mạo (2000, 2014) đặt nền móng lý thuyết ngẫu nhiên cho công trình thủy; Phạm Hồng Cường (2009) phát triển phần mềm DTC2007 cho hệ thống kênh; Mai Văn Công (2010) đánh giá độ tin cậy đê biển Nam Định; Lê Xuân Bảo (2017) nghiên cứu kiểm soát ngập lũ đa thành phần TP. Hồ Chí Minh; Trần Quang Hoài (2018) nghiên cứu hệ thống đê đồng bằng sông Hồng; Nguyễn Lan Hương (2014, 2017) phát triển phần mềm SYPRO2016 đánh giá độ tin cậy đầu mối hồ chứa làm việc nối tiếp. Về phân tích rủi ro hạ du, Dự án Hợp tác Việt Nam - New Zealand (2016–2020) xây dựng sổ tay đánh giá nhanh, các dự án WB8 và ADB đưa ra tiêu chuẩn kiểm tra lũ theo rủi ro hạ du, cùng với Tiêu chuẩn kỹ thuật TCKT 03:2015 và TCVN 9905:2014.

Tuy nhiên, tồn tại 2 trường phái tranh luận và xung đột học thuật lớn:

  • Trường phái tất định truyền thống vs. Trường phái xác suất ngẫu nhiên: Quan điểm tất định cho rằng hệ số an toàn $[K]$ cố định là đủ để đảm bảo công trình an toàn; ngược lại, trường phái xác suất chứng minh rằng hai công trình có cùng hệ số $[K]$ vẫn có thể có xác suất vỡ đập chênh lệch hàng trăm lần do mức độ biến động (phương sai) của tải trọng và sức kháng khác nhau.
  • Trường phái an toàn công trình đơn lẻ vs. Trường phái quản lý rủi ro hệ thống: Nhiều nghiên cứu trước đây chỉ tập trung vào đập chính hoặc chỉ lập bản đồ ngập lụt hạ du một cách độc lập mà không liên kết định lượng giữa xác suất sự cố $P_f$ của CTĐM với đường cong tổn thất kinh tế - xã hội hạ du.

So sánh với nghiên cứu của Negede Abate Kassa (2009) tại Đức (chỉ xét 3 cơ chế đập đất và giả định các biến ngẫu nhiên độc lập) và Laurent Peyras et al. (2012) tại Pháp (chỉ xét đập bê tông và sự cố cửa van), luận án của Cầm Thị Lan Hương đã tạo ra bước định vị học thuật vượt trội: xây dựng hệ thống cây sự cố toàn diện cho toàn bộ cụm CTĐM (đập đất, tràn xả lũ, cống lấy nước ngầm) với 12 cơ chế phá hủy thành phần, xử lý tương quan phi tuyến của các biến ngẫu nhiên và tích hợp mô hình kinh tế tối ưu rủi ro ngập lụt hạ du trong điều kiện đặc thù của Việt Nam.

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

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

Luận án mở rộng và phát triển các lý thuyết nền tảng của Freudenthal, Cornell, Hasofer-Lind và Rackwitz-Fiessler trong cơ học kết cấu công trình thủy thông qua 3 đóng góp lý thuyết căn bản:

  1. Hoàn thiện mô hình trạng thái giới hạn đa cơ chế cho cụm CTĐM: Xây dựng tường minh 12 hàm tin cậy $Z_i = R_i - S_i$ mô phỏng chính xác các quá trình vật lý - cơ học phá hủy bao gồm: nước tràn đỉnh đập ($Z_1$), trượt mái hạ lưu theo phương pháp Bishop cải tiến ($Z_2$), biến hình thấm xói ngầm thân và nền đập ($Z_3, Z_4$), mất ổn định trượt/lật thân tràn xả lũ ($Z_5, Z_6$), khí thực dốc nước ($Z_7$), thấm dọc mang cống ngầm ($Z_8$), cường độ chịu lực thân cống ngầm ($Z_9$) và khí thực sau cửa van cống ($Z_{10}$).
  2. Xác lập mô hình toán tối ưu hóa rủi ro kinh tế toàn diện: Mở rộng lý thuyết cân bằng kinh tế rủi ro của ICOLD và Vrijling, thiết lập hàm mục tiêu chi phí tổng thể $C_{tot}(P_f)$ tích hợp giá trị hiện tại thuần của chi phí đầu tư nâng cấp $I(P_f)$, chi phí vận hành bảo trì $PV(M_{Pf})$ và kỳ vọng rủi ro thiệt hại hạ du $PV(P_f \cdot D)$: $$C_{tot} = I_0 + I_{Pf} + PV(M_{Pf}) + PV(P_f \cdot D) \to \min$$ Phương trình vi phân xác định xác suất sự cố tối ưu $P_{f\text{-opt}}$ thỏa mãn điều kiện: $\frac{d I(P_f)}{d P_f} + \frac{d PV(M_{Pf})}{d P_f} + D(P_f) = 0$.
  3. Chuyển dịch hệ hình (Paradigm Shift): Chuyển đổi căn bản triết lý thiết kế công trình thủy lợi từ kiểm soát an toàn cục bộ tất định sang quản trị rủi ro thích ứng ngẫu nhiên theo hệ thống, gắn liền ngưỡng an toàn công trình với mức chấp nhận rủi ro của cộng đồng hạ du.

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

Khung phân tích của luận án là sự hợp nhất của 3 trụ cột lý thuyết: Lý thuyết độ tin cậy hệ thống + Thủy lực - Thủy văn ngẫu nhiên 1D/2D + Kinh tế học rủi ro công trình.

[Phân tích Độ tin cậy CTĐM (FTA, FORM, MCS)] 
          [Xác suất sự cố hệ thống Pf]
      [Mô phỏng Thủy lực Hạ du (MIKE)]     [Hàm Thiệt hại D(h, v)]
                      [Rủi ro Ngập lụt R = Pf × D]
             [Tối ưu hóa Chi phí Vòng đời Ctot ➔ Min]
          [Độ tin cậy Yêu cầu [Pf], [β] & Phân bổ Nâng cấp]

Cách tiếp cận phân tích giải quyết bài toán qua 3 bước liên hoàn:

  • Bài toán 1: Xác định $P_f$ và $\beta$ hiện trạng của hệ thống CTĐM bằng sơ đồ cây sự cố và giải hàm tin cậy ở Cấp độ II/III.
  • Bài toán 2: Xác định độ tin cậy yêu cầu $[P_f]$ dựa trên đường cong quan hệ giữa chi phí đầu tư nâng cấp $I(P_f)$ và tổn thất rủi ro kinh tế - xã hội hạ du $R(P_f)$.
  • Bài toán 3: Thiết kế nâng cấp kỹ thuật tối ưu thông qua giải thuật phân bổ độ tin cậy ngược từ $[P_f]$ yêu cầu về các cơ chế thành phần xung yếu nhất.

Điều kiện biên của khung phân tích: Áp dụng cho các hồ chứa có đập đất dâng nước, vùng hạ du chịu tác động ngập lụt chủ đạo do một hồ chứa điều tiết/vỡ đập, không xét ảnh hưởng triều cường cực đoan hoặc lũ liên hồ phức tạp chưa kiểm soát.

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

Thiết kế nghiên cứu

Nghiên cứu đứng trên lập trường bản thể luận hiện thực phản biện (critical realism) và phương pháp luận thực chứng định lượng (positivism). Thiết kế nghiên cứu đa tầng bậc tích hợp (multi-level integrated design):

  • Cấp độ vi mô (Micro-level): Phân tích trạng thái giới hạn ứng suất, biến dạng, trường thấm và thủy lực tại từng phân tố của thân đập, tràn, cống.
  • Cấp độ trung mô (Meso-level): Phân tích tương tác hệ thống CTĐM thông qua cấu trúc cây sự cố (cổng logic OR/AND).
  • Cấp độ vĩ mô (Macro-level): Mô phỏng lan truyền sóng lũ vỡ đập/xả lũ khẩn cấp dọc theo lưu vực sông hạ du và định lượng rủi ro kinh tế - xã hội.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý số liệu tuân thủ các chuẩn mực khoa học nghiêm ngặt:

  • Thu thập dữ liệu: Chuỗi số liệu thủy văn quan trắc mực nước hồ cực trị năm ($Z_{mn}$) dài hơn 40 năm; khảo sát địa chất công trình xác định các chỉ tiêu cơ lý đất đắp và nền ($\gamma, c, \varphi, k_t, [J]$); trắc dọc, trắc ngang lòng dẫn hạ du sông Công.
  • Tam giác hóa phương pháp (Methodological Triangulation): Kết hợp phân tích giải tích xấp xỉ bậc một (FORM), mô phỏng số ngẫu nhiên Monte Carlo ($N = 10^6$ mẫu), mô hình phần tử hữu hạn/cân bằng giới hạn ổn định mái dốc GeoSlope 2007 và mô hình thủy lực 2 chiều MIKE 11 / MIKE FLOOD.
  • Kiểm định độ tin cậy và chuẩn hóa phân bố: Sử dụng phần mềm BestFit kiểm định hàm phân bố xác suất của các biến ngẫu nhiên (chuẩn, loga chuẩn, Gumbel, Weibull); kiểm soát sai số tương đối $\varepsilon \le 1%$ và hệ số biến thiên $V_{Pf} \le 0,05$ trong mô phỏng Monte Carlo.

Data và phân tích

Trong trường hợp hồ Núi Cốc, các biến ngẫu nhiên đầu vào được xác lập chuẩn xác:

  • Biến ngẫu nhiên mực nước hồ cực trị năm tuân theo phân bố Log-Normal: $Z_{mn} \sim \text{LN}(\mu = 46,51\text{ m}; \sigma = 0,51\text{ m})$; mực nước lũ thiết kế cập nhật theo chu kỳ lặp là $+48,70$ m (tăng đáng kể so với thiết kế ban đầu $+46,51$ m).
  • Chỉ tiêu cơ lý đất đắp đập chính: Dung trọng tự nhiên $\gamma_w = 1,78 \pm 0,08\text{ g/cm}^3$, lực dính hiệu quả $c' = 0,22 \pm 0,04\text{ kg/cm}^2$, góc ma sát trong $\varphi' = 19^\circ \pm 2,1^\circ$, hệ số thấm $k_t = (1,8 \pm 0,35) \times 10^{-4}\text{ cm/s}$, gradient thấm cho phép $[J] = 0,50 \div 0,65$.
  • Công cụ phần mềm chuyên dụng:
    • OpenFTA: Thiết lập và giải cấu trúc cây sự cố hệ thống với cổng OR cho các cơ chế nối tiếp.
    • VaP 3.0 & BestFit: Giải hàm độ tin cậy phi tuyến bằng thuật toán Rackwitz-Fiessler (Cấp độ II) và Monte Carlo (Cấp độ III).
    • GeoStudio (GeoSlope 2007): Tính toán hệ số an toàn trượt mái ngẫu nhiên kết nối với hàm tin cậy Bishop.
    • MIKE 11 & MIKE FLOOD: Mô phỏng dòng chảy ngập lụt hạ du sông Công qua các huyện Đại Từ, TP. Thái Nguyên, TP. Sông Công, TX. Phổ Yên và huyện Phú Bình ứng với 7 kịch bản xả lũ và vỡ đập.

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 cốt lõi mang tính đột phá khoa học:

  1. Định lượng chính xác hiện trạng rủi ro hồ Núi Cốc: Đánh giá hiện trạng cho thấy hệ thống CTĐM hồ Núi Cốc có xác suất sự cố hàng năm rất cao: $P_f = 1,83 \times 10^{-2}$ (tương ứng chỉ số độ tin cậy $\beta = 2,09$), vượt xa ngưỡng rủi ro cho phép quốc tế ($[P_f] \le 10^{-4} \div 10^{-5}$, $[\beta] \ge 3,71$). Trong đó, cơ chế nước tràn đỉnh đập chiếm tỷ trọng rủi ro chi phối ($P_{f1} = 1,75 \times 10^{-2}$), tiếp theo là cơ chế trượt mái hạ lưu đập chính ($P_{f2} = 8,1 \times 10^{-4}$) và biến hình thấm đặc biệt nền/thân đập ($P_{f3} = 1,2 \times 10^{-4}$).
  2. Mức độ đóng góp rủi ro phi đối xứng: Cơ chế tràn đỉnh và mất ổn định mái hạ lưu chiếm tới 96,2% tổng xác suất sự cố hệ thống, trong khi các sự cố về cống ngầm và tràn xả lũ chỉ chiếm 3,8%. Kết quả này chứng minh rằng nâng cao đỉnh đập và xử lý thoát nước gia cố mái hạ lưu là giải pháp then chốt để hạ thấp rủi ro toàn hệ thống.
  3. Phát hiện nghịch trực giác về tương quan kinh tế - an toàn: Gia tăng đầu tư cho CTĐM không tạo ra mức giảm rủi ro tuyến tính. Khi $P_f$ giảm từ $10^{-2}$ xuống $10^{-4}$, tổn thất rủi ro hạ du giảm đột biến từ 1.850 tỷ VNĐ xuống còn 18,5 tỷ VNĐ; nhưng nếu tiếp tục đầu tư để hạ $P_f$ từ $10^{-5}$ xuống $10^{-7}$, chi phí xây dựng tăng vọt theo hàm số mũ trong khi mức giảm rủi ro biên không đáng kể.
  4. Xác định điểm tối ưu kinh tế $[P_f]_{\text{opt}}$ cho hồ Núi Cốc: Điểm tối ưu tổng chi phí vòng đời đạt được tại $P_{f\text{-opt}} = 4,25 \times 10^{-5}$ (tương ứng chỉ số độ tin cậy yêu cầu $[\beta]{\text{req}} = 3,93$). Tại điểm này, tổng chi phí $C{tot}$ đạt cực tiểu với chi phí đầu tư nâng cấp CTĐM ước tính $I_{Pf} \approx 215$ tỷ VNĐ, giảm thiểu được 98,6% rủi ro ngập lụt hạ du.
  5. Bản đồ hóa phân vùng tổn thất hạ du chi tiết: Mô phỏng thủy lực 2D xác định diện tích ngập lụt hạ du sông Công lên tới 8.420 ha trong kịch bản vỡ đập, đe dọa trực tiếp hơn 120.000 dân cư tại TP. Thái Nguyên, Sông Công và Phổ Yên, với tổng thiệt hại kinh tế ước tính tối đa $D_{\max} = 12.450$ tỷ VNĐ nếu không có phương án nâng cấp và cảnh báo sớm.

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung vào hệ thống lý thuyết xây dựng thủy lợi mô hình phân tích an toàn hệ thống liên kết đầy đủ (Fully Coupled Headwork-Downstream System Reliability Model), giải quyết triệt để khiếm khuyết phân mảnh của các mô hình trước đây.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình tích hợp OpenFTA - VaP - MIKE FLOOD tạo thành bộ công cụ tính toán đồng bộ, có thể chuyển giao và áp dụng rộng rãi cho hơn 6.700 hồ chứa thủy lợi trên toàn lãnh thổ Việt Nam.
  • Về mặt thực tiễn kỹ thuật: Đưa ra phương án sửa chữa nâng cấp hồ Núi Cốc chính xác: tôn cao đỉnh đập thêm 1,2 m (đạt cao trình $+49,70$ m), mở rộng tràn xả lũ thêm 1 khoang có cửa van, gia cố cơ hạ lưu và lắp đặt hệ thống tầng lọc thoát nước ngầm thân đập.
  • Về mặt chính sách quản lý: Cung cấp cơ sở khoa học để Bộ Nông nghiệp và PTNT cùng các địa phương xây dựng bảng thứ tự ưu tiên đầu tư công trong chương trình nâng cấp an toàn đập quốc gia (WB8 và nguồn vốn trung hạn), chuyển đổi quy chuẩn thẩm định an toàn đập từ tiếp cận tĩnh sang tiếp cận ma trận rủi ro động theo Nghị định 114/2018/NĐ-CP.

Limitations và Future Research

Luận án thẳng thắn chỉ ra 4 giới hạn nghiên cứu chính:

  1. Giới hạn chuỗi số liệu quan trắc: Do hệ thống quan trắc tự động tại các hồ chứa vừa và nhỏ ở Việt Nam còn thiếu đồng bộ, chuỗi số liệu đo đạc áp lực thấm, ứng suất và chuyển vị chưa đủ dài, dẫn đến một số biến ngẫu nhiên phải sử dụng giả định phân bố chuẩn hoặc cận chuẩn theo khuyến nghị quốc tế.
  2. Giả định tính độc lập thống kê: Mô hình cây sự cố trong luận án bước đầu giả định các biến cố thành phần làm việc theo sơ đồ nối tiếp độc lập; tương quan đồng thời (cross-correlation) giữa các cơ chế do cùng một nguyên nhân (như động đất kết hợp mưa lũ lịch sử) chưa được mô hình hóa toàn diện bằng hàm Copula.
  3. Định lượng thiệt hại phi kinh tế: Mô hình tổn thất hạ du mới tập trung lượng hóa thiệt hại vật chất và tài sản trực tiếp; thiệt hại về sinh mạng con người, tổn thương tâm lý xã hội và suy thoái hệ sinh thái môi trường mới được quy đổi gần đúng qua chỉ số điểm rủi ro mà chưa có hàm lượng giá kinh tế tài nguyên sinh thái chi tiết.

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

  • Phát triển mô hình độ tin cậy phụ thuộc thời gian (Time-variant Reliability Analysis) xét đến quá trình lão hóa bê tông, xói mòn và suy giảm cường độ kháng cắt của đất đắp theo tuổi thọ công trình.
  • Ứng dụng lý thuyết mờ - ngẫu nhiên (Fuzzy-Random Reliability) và mạng Bayes (Bayesian Networks) để cập nhật xác suất sự cố thời gian thực từ dữ liệu cảm biến IoT quan trắc đập.
  • Nghiên cứu an toàn hồ chứa bậc thang và rủi ro vỡ đập dây chuyền (Cascade Dam Failure) trên các lưu vực sông xuyên biên giới.

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

Nghiên cứu tạo ra tác động sâu rộng trên cả 3 phương diện học thuật, công nghiệp và chính sách:

  • Ảnh hưởng học thuật: Thiết lập chuẩn mực nghiên cứu mới trong lĩnh vực an toàn công trình thủy văn - thủy lợi tại Việt Nam; cung cấp tài liệu tham khảo chuyên khảo cao cấp cho công tác đào tạo sau đại học tại các viện nghiên cứu và trường đại học kỹ thuật công trình.
  • Chuyển đổi ngành và thực tiễn kỹ thuật: Thay đổi phương thức tính toán của các công ty tư vấn thiết kế thủy lợi; tối ưu hóa chi phí đầu tư hàng ngàn tỷ đồng của các dự án sửa chữa, nâng cấp hồ chứa, tránh tình trạng thiết kế thừa quy mô hoặc xử lý chắp vá không đúng cơ chế xung yếu.
  • Tác động chính sách và an sinh xã hội: Đóng góp luận cứ kỹ thuật phục vụ triển khai hiệu quả Luật Thủy lợi và Nghị định số 114/2018/NĐ-CP của Chính phủ; bảo vệ an toàn tính mạng và sinh kế bền vững cho hàng triệu người dân sinh sống tại các vùng hạ lưu đập trên toàn quốc.

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

  • Nghiên cứu sinh và học giả: Tiếp cận phương pháp luận giải hàm tin cậy phi tuyến Cấp độ II/III và quy trình mô hình hóa cây sự cố hệ thống thủy công phức hợp.
  • Kỹ sư tư vấn thiết kế và thẩm tra: Sở hữu công cụ và chỉ dẫn tính toán định lượng để phân bổ độ tin cậy, xác định kích thước hình học đập, tràn, cống tối ưu kinh tế - kỹ thuật.
  • Doanh nghiệp khai thác công trình thủy lợi: Có quy trình khoa học để đánh giá hiện trạng an toàn đập định kỳ, lập kế hoạch bảo trì phòng ngừa và xây dựng phương án ứng phó khẩn cấp (EAP) sát với thực tiễn.
  • Cơ quan quản lý nhà nước (Bộ Nông nghiệp và PTNT, UBND cấp tỉnh): Có công cụ định lượng để xếp hạng mức độ nguy hiểm của 1.200 hồ chứa xuống cấp, từ đó phân bổ ngân sách đầu tư công minh bạch, hiệu quả cao nhất.

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ì?

Đóng góp độc đáo nhất là việc mở rộng lý thuyết cân bằng rủi ro công trình thủy công thông qua việc thiết lập tường minh mô hình toán tối ưu hóa tổng chi phí vòng đời $C_{tot}(P_f) \to \min$, tích hợp đồng thời hàm xác suất sự cố hệ thống đa cơ chế $P_{f,\text{sys}}$ với hàm thiệt hại ngập lụt hạ du $D(h, v)$, loại bỏ hoàn toàn tính chủ quan của phương pháp hệ số an toàn tất định truyền thống.

2. Điểm mới về phương pháp luận so với các nghiên cứu quốc tế trước đây?

So với nghiên cứu của Negede Abate Kassa (2009) chỉ giải 3 cơ chế đập đất riêng rẽ và Laurent Peyras et al. (2012) chỉ phân tích đập bê tông đầm lăn, luận án đã xây dựng mô hình cây sự cố tổng hợp cho toàn bộ cụm CTĐM (đập đất, tràn, cống ngầm) gồm 12 cơ chế phá hủy vật lý liên kết, giải đồng thời bằng phương pháp giải tích FORM tại điểm thiết kế $X^*$ và kiểm chứng bằng mô phỏng Monte Carlo $10^6$ vòng lặp kết hợp mô hình thủy lực 2D MIKE FLOOD.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm tại hồ Núi Cốc?

Mặc dù hồ Núi Cốc vẫn đang vận hành bình thường theo đánh giá tất định thông thường, nhưng phân tích độ tin cậy ngẫu nhiên chỉ ra xác suất sự cố hiện trạng lên tới $P_f = 1,83 \times 10^{-2}$ (thuộc mức rủi ro rất cao), trong đó riêng nguy cơ tràn đỉnh chiếm tới 95,6% xác suất rủi ro do đỉnh lũ tính toán cập nhật theo biến đổi khí hậu đã tăng 1,85 lần so với thiết kế năm 1970.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) không?

Có. Luận án cung cấp chi tiết toàn bộ thuật toán thiết lập 12 hàm tin cậy $Z_i$, cấu trúc bảng ma trận sự cố, quy trình nhập liệu liên kết giữa OpenFTA, VaP 3.0, GeoSlope và MIKE 11/MIKE FLOOD cùng đầy đủ bảng số liệu thống kê biến ngẫu nhiên trong phần Phụ lục, cho phép tái lập quy trình tính toán cho bất kỳ hồ chứa thủy lợi nào.

5. Định hướng chương trình nghiên cứu 10 năm tiếp theo là gì?

Tập trung phát triển hệ thống giám sát cảnh báo sớm thông minh dựa trên bản sao số (Digital Twin) tích hợp mạng nơ-ron nhân tạo (ANN) và mô hình độ tin cậy phụ thuộc thời gian để dự báo xác suất vỡ đập theo thời gian thực (Real-time Dynamic Reliability Assessment) dưới tác động của kịch bản biến đổi khí hậu cực đoan.

Kết luận

Luận án tiến sĩ của tác giả Cầm Thị Lan Hương đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra với 5 đóng góp cốt lõi:

  1. Xác lập hệ thống cơ sở khoa học và phương pháp luận toàn diện ứng dụng Lý thuyết độ tin cậy và Phân tích rủi ro trong đánh giá an toàn hồ chứa thủy lợi tại Việt Nam.
  2. Thiết lập tường minh 12 hàm độ tin cậy mô phỏng các cơ chế phá hủy vật lý của toàn bộ cụm CTĐM (đập đất, tràn xả lũ, cống ngầm).
  3. Xây dựng giải thuật xác định độ tin cậy yêu cầu $[P_f]$ tối ưu kinh tế gắn liền với mức chấp nhận rủi ro ngập lụt hạ du.
  4. Đề xuất quy trình phân bổ độ tin cậy phục vụ thiết kế sửa chữa, nâng cấp tối ưu kỹ thuật - kinh tế.
  5. Ứng dụng thành công vào hồ Núi Cốc, cung cấp các thông số kỹ thuật chuẩn xác phục vụ trực tiếp công tác cải tạo công trình và bảo vệ an toàn cho lưu vực hạ du tỉnh Thái Nguyên.

Công trình đánh dấu bước tiến quan trọng trong việc hiện đại hóa công tác quản trị an toàn đập tại Việt Nam, mở ra hướng ứng dụng sâu rộng các công cụ toán học xác suất tiên tiến vào giải quyết bài toán kỹ thuật công trình thủy lợi phức tạp mang tầm vóc quốc gia và quốc tế.