Tổng quan về luận án

Bối cảnh khoa học của nghiên cứu xuất phát từ thách thức ngập lụt nghiêm trọng tại vùng hạ du hệ thống sông Đồng Nai – Sài Gòn (ĐN-SG), nơi tập trung trung tâm kinh tế đầu tàu Thành phố Hồ Chí Minh (TP.HCM) với quy mô dân số trên 8 triệu người, đóng góp tới 20,5% tổng GDP cả nước và có tốc độ mở rộng đô thị đạt 4%/năm. Về mặt địa hình, gần 70% diện tích TP.HCM nằm ở cao độ thấp dưới +2,0 m, tạo thành vùng trũng chịu tác động đồng thời của chế độ bán nhật triều không đều từ Biển Đông, lũ thượng nguồn sông Đồng Nai, sông Sài Gòn (qua hồ Trị An và Dầu Tiếng), lũ sông Mê Kông truyền qua sông Vàm Cỏ, cùng hiện tượng mưa cực đoan đô thị gia tăng. Dữ liệu quan trắc tại trạm thủy văn Phú An cho thấy mực nước đỉnh triều lịch sử đã đạt mức +1,68 m vào ngày 18/10/2016, vượt xa mức báo động III (+1,50 m). Trong giai đoạn 2011–2014, số lần triều vượt báo động III xuất hiện 76 lần (tăng 506,67% so với thời kỳ 2006–2010). Đồng thời, kịch bản biến đổi khí hậu và nước biển dâng dự báo mực nước triều tại khu vực sẽ dâng thêm 23–27 cm vào năm 2050 và 59–75 cm vào năm 2100 (đỉnh triều tương ứng đạt 1,76 m và 2,23 m).

Khoảng trống nghiên cứu (research gap) cốt lõi được nhận diện là sự bất cập mang tính hệ thống trong các tiêu chuẩn thiết kế hiện hành tại Việt Nam (TCVN 9902:2016, TCVN 9901:2014, QCVN 04-05:2012/BNNPTNT). Phương pháp truyền thống xác định Mức bảo đảm an toàn (MBĐAT) hay Tiêu chuẩn an toàn (TCAT) chống ngập thuần túy dựa vào cấp công trình (Cấp I tương ứng tần suất thiết kế $1/250$ năm theo Quy hoạch thủy lợi chống ngập úng QH1547) thông qua các chỉ tiêu định tính hoặc gián tiếp như quy mô diện tích và mật độ dân số hành chính. Cách tiếp cận này bộc lộ hai nhược điểm nghiêm trọng:

  1. Áp đặt một giá trị TCAT đồng nhất ($1/250$ năm) cho toàn bộ diện tích bảo vệ gần $2.000\text{ km}^2$, hoàn toàn tách rời quan hệ hàm số giữa quy mô đầu tư công trình với giá trị tài sản kinh tế - xã hội cụ thể của từng tiểu vùng bảo vệ;
  2. Đánh giá an toàn công trình theo phương pháp tiền định (deterministic) với hệ số an toàn cục bộ ($SF \ge [SF]$) trên từng mặt cắt đại diện, bỏ qua tính ngẫu nhiên của điều kiện biên tải trọng thủy lực - địa kỹ thuật và chưa xem xét tương tác tổng thể của hệ thống đa thành phần.

Luận án của nghiên cứu sinh Lê Xuân Bảo (Trường Đại học Thủy Lợi, 2017) dưới sự hướng dẫn của GS. Franz Nestmann và PGS. Mai Văn Công đã giải quyết triệt để bài toán này thông qua hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để mô hình hóa xác suất sự cố của các công trình kiểm soát ngập lụt (KSNL) đa thành phần chịu tổ hợp tải trọng ngẫu nhiên của triều cường, lũ và nền đất yếu phức tạp?
  • Câu hỏi nghiên cứu 2 (RQ2): Mối quan hệ tương tác giữa hàm chi phí đầu tư xây dựng/nâng cấp và hàm thiệt hại kinh tế tiềm tàng do ngập lụt được thiết lập ra sao để xác định MBĐAT tối ưu?
  • Câu hỏi nghiên cứu 3 (RQ3): Giá trị MBĐAT kinh tế tối ưu cho từng vùng kiểm soát ngập lụt cụ thể tại TP.HCM (điển hình là vùng IA1-3) có khác biệt như thế nào so với tần suất thiết kế $1/250$ năm theo quy hoạch hiện hành?
  • Giả thuyết 1 (H1): Việc chuyển đổi từ phương pháp hệ số an toàn truyền thống sang phương pháp phân tích độ tin cậy kết cấu bậc một (FORM) và mô phỏng Monte Carlo (MCS) cho phép lượng hóa chính xác xác suất sự cố tích lũy ($P_f$) của từng cơ chế phá hoại và toàn hệ thống.
  • Giả thuyết 2 (H2): Tổng chi phí vòng đời của hệ thống KSNL ($C_{tot}$) tồn tại một điểm cực tiểu cục bộ tương ứng với mức rủi ro kinh tế chấp nhận được, đóng vai trò là MBĐAT kinh tế tối ưu.
  • Giả thuyết 3 (H3): Mức bảo đảm an toàn tối ưu theo quan điểm phân tích rủi ro có tính phân hóa không gian theo mục đích sử dụng đất và giá trị tài sản, giúp tiết kiệm đáng kể nguồn vốn đầu tư công mà vẫn bảo đảm độ tin cậy mục tiêu.

Phạm vi nghiên cứu tập trung vào hệ thống công trình KSNL khép kín bảo vệ phía Nam vùng IA1 (tiểu vùng IA1-3 thuộc Quy hoạch QH1547) tại TP.HCM, bao gồm cống ngăn triều dạng trụ đỡ lớn (cống Phú Định), tuyến đê bao sông Cần Giuộc và tuyến kè tường đứng sông Sài Gòn.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về kiểm soát ngập lụt và an toàn công trình phân hóa thành ba dòng lý thuyết chính:

  1. Trường phái an toàn tiền định và hệ số cục bộ: Khởi nguồn từ cơ học kết cấu cổ điển và tiêu chuẩn xây dựng thế kỷ 20, dòng nghiên cứu này xác định độ an toàn qua tỷ số giữa sức kháng danh định và tải trọng thiết kế ($SF = R/S$). Cách tiếp cận này được quy chuẩn hóa trong các cẩm nang thiết kế của USACE (Công binh Lục quân Hoa Kỳ) trước năm 1990 và các tiêu chuẩn thủy lợi Việt Nam (TCVN 9901:2014, TCVN 9902:2016). Điểm nghẽn là không phản ánh được độ bất định ngẫu nhiên (aleatory uncertainty) và độ bất định nhận thức (epistemic uncertainty) của các thông số đất nền và khí tượng thủy văn.
  2. Trường phái Lý thuyết độ tin cậy kết cấu (Structural Reliability Theory): Phát triển ban đầu bởi Mayer (1920), Freudenthal (1956), Cornell (1969) và hoàn thiện đột phá với chỉ số độ tin cậy hình học bất biến của Hasofer & Lind (1974). Lý thuyết này mô tả biên trạng thái giới hạn $Z = R - S = 0$ và tính xác suất sự cố $P_f = \Phi(-\beta)$ trong không gian biến chuẩn hóa. Tại Việt Nam, phương pháp độ tin cậy đã được GS. Phạm Khắc Hùng nghiên cứu cho công trình biển ngoài khơi, GS. Nguyễn Văn Mạo (2001–2003) ứng dụng sơ bộ cho hồ đập thủy lợi, và GS. Ngô Trí Viềng (2008–2010) áp dụng phân tích tải trọng sóng trên đê biển. Tuy nhiên, các nghiên cứu này chỉ dừng lại ở cấu kiện đơn lẻ, chưa giải quyết bài toán hệ thống đa thành phần tương tác.
  3. Trường phái Quản lý tổng hợp rủi ro ngập lụt (Integrated Flood Risk Management): Tiêu biểu là công trình kinh điển của J.K. Vrijling (1998) tại Đại học Công nghệ Delft (Hà Lan) về tiêu chuẩn an toàn công trình chống lũ ("Safety Standards of Flood Defenses"), tích hợp xác suất sự cố công trình với hàm thiệt hại kinh tế xã hội $Risk = P_f \times Damage$. Dự án FLOODSite của Liên minh Châu Âu (2009) và chương trình nghiên cứu XtremRisk của Đức (2010) tại cửa sông Elbe (Hamburg) đã thiết lập khung đánh giá rủi ro ngập lụt tích hợp mô hình thủy lực và độ tin cậy công trình.
       [TRƯỜNG PHÁI TIỀN ĐỊNH CỔ ĐIỂN]
   (USACE pre-1990, TCVN 9901:2014, TCVN 9902:2016)
         - Hệ số an toàn tĩnh: SF = R/S >= [SF]
         - Tần suất thiết kế danh định đồng nhất (1/250 năm)
   [TRƯỜNG PHÁI ĐỘ TIN CẬY KẾT CẤU (RELIABILITY)]
      (Mayer 1920, Cornell 1969, Hasofer & Lind 1974)
         - Hàm trạng thái giới hạn: Z = R - S = 0
         - Chỉ số độ tin cậy β và xác suất sự cố Pf = Φ(-β)
   [TRƯỜNG PHÁI QUẢN LÝ RỦI RO TÍCH HỢP (IFRM)]
  (Vrijling 1998, FLOODSite EU 2009, XtremRisk Đức 2010)
         - Mô hình tối ưu hóa rủi ro kinh tế: Risk = Pf × Damage
         - Tiêu chuẩn đê biển Hà Lan (1/4.000 - 1/10.000 năm)
                               hạ du sông triều nhiệt đới tại Đông Nam Á)
     VỊ TRÍ TIÊN PHONG CỦA LUẬN ÁN (LÊ XUÂN BẢO, 2017)
     - Tích hợp FORM/MCS + Phân tích Cây sự cố (FTA)
     - Khảo sát vi mô 1.139 hộ & GIS mô hình hóa thiệt hại
     - Tối ưu hóa kinh tế xác định MBĐAT cho vùng IA1-3 TP.HCM

Cuộc tranh luận khoa học lớn nhất diễn ra giữa quan điểm "tuyệt đối hóa an toàn kỹ thuật" (thiết kế theo mực nước cực hạn lịch sử cộng biên tự do) và quan điểm "tối ưu hóa rủi ro kinh tế" (chấp nhận một xác suất sự cố có kiểm soát để đạt tổng chi phí xã hội thấp nhất). Luận án định vị nghiên cứu tại điểm giao thoa giữa Lý thuyết độ tin cậy kết cấu bậc một (FORM Cấp độ II), mô phỏng Monte Carlo (Cấp độ III) và mô hình hóa không gian thiệt hại kinh tế ngập lụt vi mô.

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

  • Trường hợp Hà Lan (Delta Works / Vrijling, 1998): Hà Lan áp dụng tiêu chuẩn bảo vệ đê biển cực cao ($1/4.000$ đến $1/10.000$ năm) nhờ điều kiện kinh tế phát triển và địa hình phần lớn dưới mực nước biển. Tuy nhiên, mô hình này giả định hệ thống đê đồng nhất và địa chất tương đối ổn định. Luận án kế thừa nguyên lý tối ưu hóa kinh tế của Vrijling nhưng cải tiến cấu trúc cây sự cố (Fault Tree Analysis - FTA) để phù hợp với hệ thống công trình đa dạng (cống, đê, kè) trên nền bùn sét yếu của đồng bằng ngập triều TP.HCM.
  • Trường hợp Thành phố New Orleans, Hoa Kỳ (USACE post-Katrina, 2008): Sau thảm họa bão Katrina (2005), Mỹ chuyển đổi hệ thống kiểm soát ngập sang phương pháp độ tin cậy với chỉ số mục tiêu $[\beta] = 4,0$ (tần suất vượt $1/100$ năm có xét đến rủi ro nhân mạng). Nghiên cứu của Lê Xuân Bảo mở rộng bài toán này cho đô thị đang phát triển, nơi dữ liệu vi mô về thiệt hại được thu thập thực địa kết hợp hàm chiều sâu ngập lụt phi tuyến cho 1.139 hộ dân và cơ sở sản xuất kinh doanh tại 17 quận/huyện TP.HCM.

Đó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 độ tin cậy kết cấu của Hasofer & Lind (1974) và lý thuyết rủi ro kinh tế của J.K. Vrijling (1998) trong các khía cạnh:

  1. Thiết lập mô hình trạng thái giới hạn phi tuyến đa biến cho hệ thống kiểm soát ngập hỗn hợp triều - lũ: Luận án mở rộng phương trình trạng thái giới hạn $Z = R - S$ cho 17 cơ chế sự cố cơ bản thuộc 3 loại hình công trình thủy công điển hình: cống ngăn triều trụ đỡ, đê bao đất, và kè sông tường đứng.
  2. Cơ chế sự cố và hàm tin cậy tương ứng:
    • Cơ chế tràn đỉnh: $$Z_1 = Z_{\text{công trình}} - (MN_{\max} + \Delta h_{\text{gió}} + R_{\text{sóng}})$$
    • Cơ chế mất ổn định nền và cọc chịu lực thẳng đứng (cống và kè): $$Z_2 = R_{\text{mũi}} + R_{\text{bên}} - N_t$$
    • Cơ chế mất ổn định tải trọng ngang của cọc: $$Z_3 = Q_c - Q_{\text{tác dụng}}$$
    • Cơ chế mất ổn định thấm và xói ngầm đáy móng đê/cống (theo tiêu chí Bligh và Darcy): $$Z_4 = J_0 - \frac{\Delta H}{L_{\text{thấm}}}$$
    • Cơ chế trượt sâu mái dốc đê bao: $$Z_5 = SF_{\text{tính toán}}(\phi, c, \gamma, \text{mực nước}) - [SF]$$
    • Cơ chế lún cố kết vượt ngưỡng cho phép: $$Z_6 = [S] - \sum \frac{e_{1i} - e_{2i}}{1 + e_{1i}} h_i$$
  3. Mô hình toán tối ưu hóa chi phí rủi ro vòng đời hệ thống: Luận án thiết lập phương trình hàm mục tiêu tổng chi phí quy về giá trị hiện tại ($PV$): $$C_{tot}(P_f) = I(P_f) + PV\big(M(P_f)\big) + PV\big(P_f \cdot D\big) \longrightarrow \min$$ Trong đó:
    • $I(P_f)$ là hàm chi phí đầu tư xây dựng/nâng cấp hệ thống tương ứng với xác suất sự cố $P_f$;
    • $PV(M)$ là hiện giá chi phí duy tu, bảo dưỡng định kỳ trong suốt tuổi thọ công trình $T$ với tỷ suất chiết khấu $r$: $$PV(M) = \sum_{t=1}^{T} \frac{E(M)}{(1 + r)^t}$$
    • $PV(P_f \cdot D)$ là hiện giá kỳ vọng rủi ro thiệt hại kinh tế do ngập lụt khi hệ thống gặp sự cố: $$PV(P_f \cdot D) = \sum_{t=1}^{T} \frac{P_f \cdot E(D)}{(1 + r)^t}$$
 Chi phí ($)
                      Pf,tối ưu
               (MBĐAT Kinh tế tối ưu)

Sự chuyển dịch mô thức (paradigm shift) thể hiện ở việc thay thế hoàn toàn quan điểm an toàn cục bộ tĩnh bằng quan điểm an toàn hệ thống ngẫu nhiên, chứng minh rằng điểm tối ưu kinh tế không nằm ở mức an toàn tuyệt đối ($P_f \to 0$) mà nằm tại điểm cân bằng biên giữa chi phí cận biên xây dựng công trình và rủi ro cận biên giảm thiểu được.

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

Khung phân tích của luận án tích hợp liên ngành giữa ba trụ cột lý thuyết: Cơ học đất - Kết cấu công trình thủy, Thủy lực - Thủy văn sông triều và Kinh tế học không gian ngập lụt.

Hệ thống KSNL được phân rã thành sơ đồ Cây sự cố (Fault Tree Analysis - FTA) với các cổng logic OR và AND:

  • Toàn hệ thống IA1-3 bao gồm 3 thành phần chính mắc nối tiếp: $$P_{f,sys} = 1 - \prod_{i=1}^{m} (1 - P_{f,comp_i}) = 1 - (1 - P_{f,\text{cống}})(1 - P_{f,\text{đê}})(1 - P_{f,\text{kè}})$$
  • Mỗi công trình thành phần lại là một hệ thống nối tiếp của các hạng mục kết cấu và các cơ chế phá hoại: $$P_{f,comp} = 1 - \prod_{j=1}^{k} (1 - P_{f,mechanism_j})$$

Điều kiện biên áp dụng (boundary conditions) được xác định rõ ràng: vùng châu thổ hạ du chịu tác động của triều cường bán nhật triều không đều, đất nền thuộc trầm tích Holocene sét yếu có chiều dày $15\text{ m} \div 30\text{ m}$, hệ số nén lún cao ($e_0 > 1,2$), và khu vực bảo vệ có mức độ đô thị hóa nhanh chóng.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên nền tảng triết học thực chứng (positivism) kết hợp hiện thực phê phán (critical realism), sử dụng thiết kế phương pháp hỗn hợp định lượng bậc cao (advanced quantitative mixed-method design). Thiết kế nghiên cứu đa cấp độ (multi-level research architecture) được triển khai qua 4 tầng phân tích:

  1. Tầng 1 - Cơ chế vật lý cơ sở: Mô hình hóa toán học 17 hàm trạng thái giới hạn cơ học - thủy lực.
  2. Tầng 2 - Độ tin cậy thành phần và hệ thống: Ứng dụng giải thuật FORM (Cấp độ II) và MCS (Cấp độ III) để giải xác suất sự cố $P_f$ và chỉ số $\beta$.
  3. Tầng 3 - Mô phỏng thủy lực và ngập lụt: Sử dụng bộ mô hình MIKE 11 / MIKE 21 để mô phỏng lan truyền sóng lũ, triều dâng và tính toán bản đồ ngập lụt số hóa theo các kịch bản tần suất ($p = 1%, 2%, 5%, 10%, 20%$).
  4. Tầng 4 - Mô hình hóa kinh tế không gian: Tích hợp phần mềm ArcGIS, phân lớp bản đồ hiện trạng và quy hoạch sử dụng đất (đến năm 2025 và 2100), áp dụng hàm thiệt hại vi mô để tính toán đường cong rủi ro kinh tế.
       QUY TRÌNH PHÂN TÍCH ĐỘ TIN CẬY VÀ TỐI ƯU HÓA RỦI RO HỆ THỐNG

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu kiểm soát sai số nghiêm ngặt qua việc phân tích thống kê và tam giác đạc dữ liệu (data triangulation):

  • Phân tích phân bố ngẫu nhiên: Thu thập chuỗi số liệu quan trắc mực nước đỉnh triều từ năm 1983 đến 2014 tại các trạm Phú An, Nhà Bè; số liệu mưa 40 năm (1962–2014) tại trạm Tân Sơn Nhất; số liệu địa chất công trình gồm 45 lỗ khoan tại khu vực quận 8, Bình Chánh, Bình Tân.
  • Xử lý thống kê biến ngẫu nhiên: Sử dụng phần mềm BestFit để kiểm định phân phối xác suất (Goodness-of-Fit tests: Chi-Square, Kolmogorov-Smirnov, Anderson-Darling). Các biến tải trọng thủy lực ($MN_{\max}$, vận tốc gió $v_w$, góc ma sát trong $\phi$, lực dính $c$, trọng lượng thể tích $\gamma$, mô đun biến dạng $E$) được xác định chính xác giá trị kỳ vọng $\mu$, độ lệch chuẩn $\sigma$, hệ số biến thiên $V_x = \sigma/\mu$, và gán các hàm phân bố xác suất chuẩn (Normal), Lognormal, Gumbel hoặc Weibull.
  • Dữ liệu điều tra vi mô thiệt hại ngập lụt: Sử dụng tập dữ liệu khảo sát thực địa toàn diện gồm 1.139 phiếu điều tra phủ khắp 17 quận/huyện TP.HCM của nhóm nghiên cứu Đại học Bách khoa Delft và Đại học Thủy Lợi (2012), kết hợp 126 phiếu điều tra doanh nghiệp/hộ sản xuất kinh doanh và 291 phiếu hộ dân cư tại các quận 5, 6, 8, 11, Bình Tân, Bình Thạnh, Tân Bình, Tân Phú (2014).
  • Hệ thống công cụ phần mềm chuyên dụng:
    • BestFit: Xây dựng hàm mật độ xác suất (PDF) và hàm phân bố tích lũy (CDF) của các biến đầu vào.
    • VAP (Variable Analysis Package): Tính toán chỉ số độ tin cậy $\beta$, xác suất sự cố $P_f$, tọa độ điểm thiết kế (Most Probable Point - MPP) và hệ số ảnh hưởng độ nhạy $\alpha_i^2$ theo thuật giải FORM.
    • OpenFTA: Xây dựng và tính toán sơ đồ cây sự cố toàn hệ thống.
    • PLAXIS 2D / GeoSlope (SLOPE/W, SEEP/W): Tính toán các hệ số ổn định trượt mái, trường thấm và biến dạng địa kỹ thuật làm cơ sở hiệu chỉnh hàm trạng thái giới hạn.
    • ArcGIS: Phân tích không gian, chồng xếp bản đồ ngập lụt lên bản đồ hiện trạng sử dụng đất để xuất ma trận thiệt hại $D_i(h)$.

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

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

  1. Sự chi phối của các cơ chế địa kỹ thuật đối với xác suất sự cố cống ngăn triều: Kết quả phân tích độ tin cậy cống Phú Định cho thấy xác suất sự cố của cống không bị chi phối bởi cơ chế tràn đỉnh ($P_{f,\text{tràn}} = 2,15 \times 10^{-4}$, $\beta = 3,52$) mà bị chi phối mạnh mẽ bởi cơ chế mất khả năng chịu tải trọng ngang của cọc và mất ổn định thấm đẩy trồi đáy móng ($P_{f,\text{ngang}} = 1,42 \times 10^{-2}$, $\beta = 2,19$). Mực nước triều cao và chênh lệch áp lực cột nước trong mùa kiệt tạo ra gradient thấm cục bộ vượt ngưỡng tới hạn tại chân công trình.
  2. Nguy cơ tiềm ẩn từ lún cố kết và xói ngầm thân đê bao Cần Giuộc: Đối với đê bao Cần Giuộc, cơ chế trượt sâu mái đê phía đồng và sự cố lún thân đê vượt trị số cho phép $[S]$ chiếm tỷ trọng áp đảo trong cây sự cố ($P_{f,\text{lún}} = 3,68 \times 10^{-2}$, $\beta = 1,79$). Biến ngẫu nhiên có độ nhạy lớn nhất là lực dính $c$ ($\alpha^2 = 48,2%$) và góc ma sát trong $\phi$ ($\alpha^2 = 35,6%$) của lớp bùn sét yếu bão hòa nước.
  3. Hiệu ứng hệ thống làm suy giảm nghiêm trọng độ an toàn thực tế: Khi phân tích độc lập từng mặt cắt, các công trình đều đạt hệ số an toàn quy phạm ($SF \ge 1,2 \div 1,3$). Tuy nhiên, khi tích hợp toàn bộ hệ thống IA1-3 qua sơ đồ nối tiếp trên cây sự cố (gồm 1 cống lớn, $12,5\text{ km}$ đê bao và $4,8\text{ km}$ kè), xác suất sự cố tổng thể của hệ thống tăng vọt lên mức $P_{f,sys} = 6,85 \times 10^{-2}$ (tương đương chỉ số tin cậy hệ thống $\beta_{sys} = 1,49$, chu kỳ lặp sự cố chỉ khoảng 15 năm). Điều này chứng minh phương pháp thiết kế truyền thống đã đánh giá quá lạc quan về mức độ an toàn thực tế của hệ thống.
  4. Xác lập MBĐAT kinh tế tối ưu có sự phân hóa rõ rệt theo kịch bản quy hoạch: Đường cong tổng chi phí $C_{tot}(P_f)$ cho thấy:
    • Kịch bản 1 (Hiện trạng sử dụng đất 2010): MBĐAT tối ưu đạt được ở tần suất sự cố $1/50$ năm ($P_f = 0,02$).
    • Kịch bản 2 (Quy hoạch sử dụng đất đến 2025): Giá trị tài sản đô thị tăng mạnh làm rủi ro $PV(P_f \cdot D)$ tăng vọt, đẩy MBĐAT tối ưu lên mức $1/100$ năm ($P_f = 0,01$).
    • Kịch bản 3 (Quy hoạch dài hạn 2100 có xét đến nước biển dâng): MBĐAT tối ưu đạt mức $1/500$ năm ($P_f = 0,002$). Phát hiện này bác bỏ tính hợp lý của việc áp đặt cứng nhắc một mức TCAT $1/250$ năm cho mọi giai đoạn phát triển như quy định tại QH1547.
Hạng mục / Kịch bản Xác suất sự cố ($P_f$) Chỉ số độ tin cậy ($\beta$) Cơ chế sự cố chi phối chính Tần suất tối ưu ($T$)
Cống Phú Định $1,42 \times 10^{-2}$ $2,19$ Tải trọng ngang của cọc & Thấm đẩy trồi $1/70$ năm
Đê bao Cần Giuộc $3,68 \times 10^{-2}$ $1,79$ Lún cố kết vượt ngưỡng & Trượt mái dốc $1/27$ năm
Kè sông Sài Gòn $2,84 \times 10^{-2}$ $1,90$ Trượt sâu & Mất sức chịu tải đứng cọc $1/35$ năm
Toàn hệ thống IA1-3 $6,85 \times 10^{-2}$ $1,49$ Tích hợp nối tiếp các thành phần $1/15$ năm
Tối ưu KT - Hiện trạng 2010 $2,00 \times 10^{-2}$ $2,05$ Tối ưu hóa hàm chi phí $C_{tot}$ $1/50$ năm
Tối ưu KT - Quy hoạch 2025 $1,00 \times 10^{-2}$ $2,33$ Tối ưu hóa hàm chi phí $C_{tot}$ $1/100$ năm
Tối ưu KT - Quy hoạch 2100 $2,00 \times 10^{-3}$ $2,88$ Tối ưu hóa hàm chi phí $C_{tot}$ $1/500$ năm

Implications đa chiều

  • Đóng góp lý luận: Bổ sung phương pháp luận hoàn chỉnh về phân tích rủi ro và độ tin cậy cho công trình thủy lợi vùng hạ du sông triều vào hệ thống học thuật Việt Nam.
  • Đổi mới phương pháp luận: Cung cấp quy trình tích hợp khép kín từ mô hình hóa biến ngẫu nhiên, phân tích độ nhạy FORM/MCS, mô phỏng cây sự cố FTA đến phân tích thiệt hại kinh tế vi mô bằng GIS, có thể chuyển giao áp dụng cho toàn bộ các lưu vực sông khác tại Đồng bằng sông Cửu Long và Đồng bằng sông Hồng.
  • Ứng dụng thực tiễn: Giúp các cơ quan quản lý (Trung tâm Điều hành Chương trình Chống ngập nước TP.HCM - TTCN, Sở Nông nghiệp và PTNT) xác định chính xác các điểm xung yếu cục bộ (ví dụ: cần tăng cường cọc chịu lực ngang cho cống Phú Định và xử lý lún nền đê Cần Giuộc thay vì nâng cao trình đỉnh đê tràn lan).
  • Khuyến nghị chính sách: Kiến nghị Bộ Nông nghiệp và PTNT, Bộ Xây dựng sửa đổi các quy chuẩn thiết kế công trình thủy lợi (TCVN 9901, TCVN 9902), cho phép áp dụng mức bảo đảm an toàn phân kỳ theo giai đoạn phát triển kinh tế thay vì ấn định cố định theo cấp hành chính.

Limitations và Future Research

  1. Hạn chế về dữ liệu lún nền động học: Nghiên cứu chưa tích hợp mô hình biến dạng lún nền kiến tạo và lún do khai thác nước ngầm quá mức (tốc độ $10 \div 30\text{ mm/năm}$) vào biến ngẫu nhiên của hàm trạng thái giới hạn theo thời gian.
  2. Giả định tính độc lập thống kê giữa một số cơ chế sự cố: Trong phân tích cây sự cố, một số biến cố thành phần được giả thiết độc lập tương đối, chưa xét triệt để tương quan chéo (cross-correlation) giữa trường thấm và trường ứng suất biến dạng phức tạp trong đất dính.
  3. Giới hạn không gian của hàm thiệt hại vi mô: Dữ liệu điều tra kinh tế - xã hội tập trung chủ yếu tại TP.HCM, cần mở rộng hiệu chỉnh khi áp dụng cho các vùng nông thôn thuần túy thuộc hạ lưu sông Đồng Nai (Đồng Nai, Bình Dương, Long An).

Chương trình nghiên cứu tương lai:

  • Phát triển mô hình độ tin cậy phụ thuộc thời gian (Time-variant Reliability Analysis) có xét đến sự suy giảm độ bền vật liệu, hiện tượng xói mòn lòng dẫn và tốc độ lún chìm mặt đất tích lũy đến năm 2050–2100.
  • Ứng dụng kỹ thuật máy học (Machine Learning Surrogate Models như Kriging, Artificial Neural Networks) để tăng tốc độ hội tụ mô phỏng Monte Carlo đối với các hàm trạng thái giới hạn ngầm định phi tuyến cao.
  • Tích hợp rủi ro phi tiền tệ (tổn thất sinh mạng, gián đoạn giao thông, ô nhiễm môi trường và dịch bệnh) vào hàm tối ưu hóa đa mục tiêu (Multi-objective Risk Optimization).

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

  • Tác động học thuật: Luận án tạo tiền đề cho hàng loạt nghiên cứu chuyên sâu về độ tin cậy công trình thủy công tại Việt Nam; kết quả nghiên cứu đã được công bố trên các tạp chí khoa học uy tín chuyên ngành Thủy lợi và Công trình Thủy.
  • Chuyển đổi ngành công nghiệp xây dựng: Thúc đẩy các tổng công ty tư vấn thiết kế thủy lợi (HEC, WASECO) tiếp cận các công cụ phân tích độ tin cậy hiện đại (VAP, OpenFTA, PLAXIS Reliability) trong lập dự án đầu tư.
  • Ảnh hưởng chính sách công: Cung cấp luận cứ khoa học vững chắc để UBND TP.HCM và các Bộ ngành Trung ương rà soát, điều chỉnh tiến độ và quy mô vốn đầu tư của Quy hoạch 1547, tránh lãng phí hàng nghìn tỷ đồng đầu tư vượt quá mức cần thiết ở giai đoạn ngắn hạn.
  • Lợi ích xã hội: Giảm thiểu thiệt hại ngập úng trực tiếp cho hàng trăm nghìn hộ dân vùng trũng quận 8, Bình Chánh, bảo vệ hạ tầng giao thông huyết mạch phía Tây Nam thành phố.

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

  • Nghiên cứu sinh và học viên sau đại học: Cung cấp khung lý thuyết, thuật toán lập trình FORM/MCS và quy trình phân tích cây sự cố chuẩn mực để mở rộng nghiên cứu cho các công trình hạ tầng khác (đập thủy điện, cảng biển, đê chắn sóng).
  • Các nhà khoa học và giảng viên ngành Kỹ thuật Công trình Thủy: Sử dụng luận án như tài liệu tham khảo giảng dạy cao học và chuyên khảo chuyên sâu về quản lý rủi ro thiên tai.
  • Kỹ sư tư vấn thiết kế và quản lý dự án: Nắm vững công cụ định lượng điểm thiết kế MPP và hệ số nhạy cảm $\alpha_i^2$ để tối ưu hóa kết cấu cọc, bản đáy cống và mái dốc đê kè.
  • Cơ quan hoạch định chính sách đô thị: Sở Kế hoạch & Đầu tư, Sở Xây dựng, Trung tâm Chống ngập có cơ sở khoa học để phân bổ ngân sách chống ngập theo thứ tự ưu tiên rủi ro.

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 và tích hợp thành công lý thuyết tối ưu hóa rủi ro kinh tế của J.K. Vrijling (1998) với lý thuyết độ tin cậy kết cấu Hasofer-Lind (1974) vào một khung tính toán thống nhất cho hệ thống kiểm soát ngập đa thành phần vùng hạ du sông triều nhiệt đới. Luận án đã giải quyết trọn vẹn mối liên kết hàm số giữa xác suất phá hoại kết cấu phi tuyến đa biến ($P_{f,sys}$) và mô hình thiệt hại kinh tế không gian vi mô ($PV(P_f \cdot D)$), xác lập cơ sở khoa học định lượng thay thế cho phương pháp chọn cấp công trình cảm tính truyền thống.

2. Sự đổi mới về mặt phương pháp luận so với các nghiên cứu trước đây?

So với các nghiên cứu trước đây tại Việt Nam (thường chỉ đánh giá độ tin cậy cho từng mặt cắt đơn lẻ bằng công thức giải tích đơn giản hoặc khảo sát thiệt hại bằng thống kê hành chính gộp), luận án đổi mới vượt bậc qua 3 điểm:

  1. Sử dụng Cây sự cố (Fault Tree Analysis) kết hợp giải thuật FORM Cấp độ II và MCS Cấp độ III để mô hình hóa toàn diện 17 cơ chế hư hỏng liên kết;
  2. Kết nối chuỗi mô hình động lực học thủy lực 1D/2D (MIKE) với hệ thông tin địa lý (ArcGIS) để mô phỏng chiều sâu ngập lụt chi tiết;
  3. Xây dựng đường cong thiệt hại phi tuyến dựa trên bộ dữ liệu khảo sát thực địa vi mô $1.139$ phiếu điều tra có trọng số kinh tế rõ ràng.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ và phản trực giác nhất là: Nguy cơ gây mất an toàn hàng đầu cho hệ thống KSNL TP.HCM không phải là hiện tượng nước triều dâng tràn qua đỉnh công trình (vốn là mối lo ngại trực quan phổ biến), mà là các cơ chế phá hoại ngầm về địa kỹ thuật. Cụ thể, tại cống Phú Định, xác suất sự cố do mất khả năng chịu tải trọng ngang của cọc cao gấp $66$ lần so với xác suất tràn đỉnh; tại đê Cần Giuộc, nguy cơ lún thân đê và xói ngầm đáy móng chiếm trên $60%$ tổng rủi ro phá hủy của tuyến đê.

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?

Hoàn toàn có. Luận án cung cấp đầy đủ và chi tiết toàn bộ các phương trình trạng thái giới hạn toán học, bảng tổng hợp đặc trưng thống kê các biến ngẫu nhiên (kỳ vọng, độ lệch chuẩn, dạng hàm phân phối), sơ đồ cây sự cố logic, các bước phối hợp giữa các phần mềm chuyên dụng (BestFit $\to$ VAP $\to$ OpenFTA $\to$ PLAXIS $\to$ MIKE $\to$ ArcGIS) và quy trình tối ưu hóa kinh tế. Bất kỳ nhóm nghiên cứu nào có bộ dữ liệu địa phương tương ứng đều có thể tái lập hoàn toàn quy trình này.

5. Định hướng chương trình nghiên cứu 10 năm tiếp theo được phác thảo như thế nào?

Chương trình nghiên cứu dài hạn tập trung vào 3 trụ cột:

  • Trụ cột 1: Xây dựng mô hình độ tin cậy phụ thuộc thời gian tích hợp biến dạng lún nền 3D và kịch bản nước biển dâng cực đoan RCP8.5;
  • Trụ cột 2: Tự động hóa giải thuật tối ưu hóa rủi ro bằng trí tuệ nhân tạo (AI-driven Reliability-based Design Optimization);
  • Trụ cột 3: Mở rộng mô hình kinh tế sang đánh giá rủi ro xã hội toàn diện (Social Vulnerability Index) và khả năng chống chịu đô thị (Urban Resilience).

Kết luận

  1. Xác lập cơ sở khoa học tiên phong: Luận án là công trình khoa học đầu tiên tại Việt Nam ứng dụng toàn diện Lý thuyết độ tin cậy kết cấu và Phương pháp phân tích rủi ro kinh tế để giải quyết bài toán xác định Mức bảo đảm an toàn cho hệ thống KSNL đa thành phần vùng hạ du sông triều chịu tác động kết hợp của lũ và triều cường.
  2. Lượng hóa chính xác độ tin cậy hệ thống: Nghiên cứu đã chứng minh rằng phương pháp thiết kế truyền thống theo hệ số an toàn cục bộ đã đánh giá quá cao mức độ an toàn thực tế; khi xét theo hệ thống nối tiếp đa cơ chế, xác suất sự cố của hệ thống IA1-3 lên tới $P_{f,sys} = 6,85 \times 10^{-2}$ ($\beta = 1,49$).
  3. Phát hiện các điểm nghẽn địa kỹ thuật then chốt: Xác định rõ tải trọng ngang tác dụng lên cọc cống ngăn triều và hiện tượng lún cố kết - xói ngầm đê bao là những cơ chế phá hoại nguy hiểm nhất cần ưu tiên xử lý trong thiết kế và duy tu.
  4. Xác định MBĐAT kinh tế tối ưu có phân kỳ: Đã chứng minh bằng thực nghiệm mô phỏng rằng MBĐAT tối ưu kinh tế cho vùng IA1-3 là $1/50$ năm cho giai đoạn hiện trạng 2010, nâng lên $1/100$ năm vào giai đoạn quy hoạch 2025, và đạt $1/500$ năm vào năm 2100.
  5. Mở ra các nhánh nghiên cứu mới: Đặt nền móng cho các hướng nghiên cứu về an toàn công trình thích ứng với biến đổi khí hậu, thiết kế công trình thủy lợi theo độ tin cậy mục tiêu, và quy hoạch hạ tầng đô thị dựa trên quản trị rủi ro thông minh.
  6. Giá trị chuyển giao và di sản thực tiễn: Cung cấp bộ công cụ phương pháp luận chuẩn xác hỗ trợ các nhà quản lý đô thị TP.HCM và các tỉnh đồng bằng ven biển đưa ra quyết định đầu tư công chính xác, an toàn, hiệu quả và bền vững trước các thách thức thiên tai trong tương lai.