Tổng quan luận án

Luận án tiến sĩ kỹ thuật với đề tài "Nghiên cứu phương pháp xác định chỉ số an toàn và độ tin cậy yêu cầu cho hệ thống đê vùng đồng bằng sông Hồng theo lý thuyết độ tin cậy và phân tích rủi ro" do nghiên cứu sinh Trần Quang Hoài thực hiện tại Trường Đại học Thủy Lợi dưới sự hướng dẫn khoa học của PGS. Mai Văn Công và GS. Trịnh Minh Thụ. Công trình thuộc chuyên ngành Xây dựng Công trình thủy, mã số 62-58-40-01, bảo vệ năm 2018.

1. Tính cấp thiết và khoảng trống nghiên cứu

Thiên tai lũ lụt tại Việt Nam diễn ra với tần suất cao và cường độ phức tạp. Các sự kiện lũ lịch sử ghi nhận thiệt hại nghiêm trọng: trận lụt năm 1971 tại đồng bằng sông Hồng làm hơn 10.000 người thiệt mạng và 250.000 ha đất nông nghiệp bị ngập úng; lũ năm 1996 làm gần 1.000 người chết tại 6 tỉnh miền Bắc; lũ năm 1999 tại miền Trung làm 750 người chết; và lũ năm 2000 tại Đồng bằng sông Cửu Long làm 539 người chết. Thống kê ước tính thiệt hại kinh tế trung bình hằng năm do thiên tai tại Việt Nam chiếm từ 1% đến 1,5% GDP và làm thiệt mạng trung bình trên 300 người/năm.

Hệ thống đê điều giữ vai trò nòng cốt trong chiến lược phòng chống lũ lụt (được quy định tại Luật Đê điều 2006). Tuy nhiên, phần lớn các tuyến đê vùng đồng bằng sông Hồng (ĐBSH) được hình thành qua nhiều giai đoạn lịch sử bằng phương pháp thủ công, vật liệu đắp tại chỗ không đồng nhất, nền đê chưa qua xử lý kỹ thuật triệt để nên tiềm ẩn nhiều ẩn họa trong thân và nền đê.

Khoảng trống nghiên cứu được tác giả chỉ rõ: Các phương pháp tính toán và đánh giá an toàn hệ thống đê theo quy định truyền thống hiện nay chủ yếu dựa trên hệ số an toàn tĩnh $K$ tại một số mặt cắt đại diện, mang tính định mức cố định. Cách tiếp cận này bộc lộ các hạn chế:

  • Chưa xét đến tính ngẫu nhiên của các yếu tố tải trọng thủy lực và sức chịu tải của vật liệu đất;
  • Chưa xác định được xác suất phá hỏng tổng hợp của toàn hệ thống công trình;
  • Bỏ qua ảnh hưởng của quy mô chiều dài đê (xác suất xảy ra sự cố tăng lên khi chiều dài tuyến đê tăng);
  • Chưa tích hợp định lượng yếu tố giá trị kinh tế - xã hội của vùng được bảo vệ và tác động của biến đổi khí hậu (BĐKH), nước biển dâng (NBD) để xác định tiêu chuẩn an toàn (TCAT) và độ tin cậy yêu cầu (ĐTCYC) tối ưu làm cơ sở ra quyết định đầu tư nâng cấp.

2. Mục tiêu nghiên cứu

Mục đích nghiên cứu của luận án là xây dựng phương pháp đánh giá, xác định chỉ số an toàn hiện tại và độ tin cậy yêu cầu cho hệ thống đê vùng đồng bằng sông Hồng trong điều kiện hiện tại và tương lai, có xét đến biến đổi khí hậu và phát triển kinh tế - xã hội.

3. Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Hệ thống đê bảo vệ vùng hạ du đồng bằng bao gồm các đoạn đê sông, đê cửa sông và đê biển. Hai trường hợp đại diện được lựa chọn tính toán ứng dụng cụ thể:
    1. Hệ thống 1 (HT1): Tuyến đê Hữu Hồng bảo vệ khu vực trung tâm thành phố Hà Nội (đại diện cho hệ thống đê sông vùng đồng bằng bảo vệ đô thị tập trung đông dân cư).
    2. Hệ thống 2 (HT2): Tuyến đê bao bảo vệ huyện Giao Thủy, tỉnh Nam Định (đại diện cho hệ thống đê phức hợp gồm đê sông Hữu Hồng, đê cửa sông và đê biển tạo thành vòng bảo vệ khép kín).
  • Phạm vi nghiên cứu:
    • Đánh giá mức độ an toàn, ổn định của hệ thống đê dưới tác động chính của mực nước lũ sông, mực nước biển, các tổ hợp tải trọng và độ bền thường xuyên.
    • Phân tích độ tin cậy có xét tính ngẫu nhiên của tải trọng và độ bền, không xét mức độ suy giảm độ bền theo thời gian.
    • Không tính toán các tải trọng đặc biệt như động đất và mưa cục bộ; độ tin cậy của các cống và công trình qua đê không đi sâu phân tích trong phần ứng dụng.
    • Mốc thời gian dự báo và tính toán nâng cấp mở rộng đến năm 2050.

Tổng quan tài liệu và vị trí của luận án

1. Các hướng nghiên cứu quốc tế

  • Lý thuyết độ tin cậy (LTĐTC): Khởi xướng từ nghiên cứu của M. Mayer (1920) với việc sử dụng giá trị trung bình thống kê và độ lệch chuẩn để tính toán an toàn công trình. Đến thập niên 1960, E. Freudenthal đề xuất lại khái niệm chỉ số độ tin cậy, sau đó được phát triển thành các phương pháp phân tích độ tin cậy kết cấu bởi Cornell, Ang tại Bắc Mỹ và Châu Âu, đưa chỉ số độ tin cậy $\beta$ thành chỉ số tiêu chuẩn đánh giá mức an toàn.
  • Ứng dụng LTĐTC và phân tích rủi ro (PTRR) trong công trình phòng chống lũ:
    • Hà Lan: Sau thảm họa lũ Biển Bắc năm 1953, Ủy ban Đồng bằng Hà Lan đã xây dựng khung khái niệm tiêu chuẩn an toàn gắn với vốn đầu tư và rủi ro kinh tế - xã hội, ứng dụng LTĐTC từ năm 1970 và chính thức đưa vào quy phạm thiết kế đê điều năm 1990 với tần suất thiết kế từ 1/4.000 đến 1/10.000 năm cho các vòng đê xung yếu.
    • Châu Âu: Triển khai chương trình nghiên cứu chung FLOODsite (2006) và hoàn thiện khung quản lý rủi ro lũ tích hợp theo phương pháp PTRR & LTĐTC vào năm 2010.
    • Hoa Kỳ: Sau sự cố vỡ đê do bão Katrina năm 2005 tại New Orleans, quân đội Mỹ đã áp dụng phương pháp PTRR & LTĐTC để quy hoạch lại hệ thống đê với chu kỳ lặp thiết kế 1/1.000 năm.
    • Châu Á: Nhật Bản và Trung Quốc chuyển dần sang phương thức thiết kế đê điều và hồ đập theo độ tin cậy yêu cầu định trước.
  • Mô hình hóa ngập lụt và xác định hàm thiệt hại: Dutta và cộng sự (2003) xây dựng mô hình tổng hợp tính toán thiệt hại lũ lụt cho lưu vực sông Ichinomiya (Nhật Bản), tích hợp GIS, mô phỏng thủy lực và các đường cong quan hệ giữa độ sâu ngập với mức độ thiệt hại nhà ở, thời gian ngập với thiệt hại cây trồng.

2. Các hướng nghiên cứu trong nước

  • Ổn định đê và địa kỹ thuật: Các công trình nghiên cứu cơ học đất thấm nền đê, thấm không ổn định, biến dạng thấm nền đê Hữu Hồng (Hà Nội) và đê tỉnh Thái Bình bằng giải pháp tường xi măng đất.
  • Mô phỏng ngập lụt: Nghiên cứu ứng dụng mô hình MIKE 21 và GIS tính toán cảnh báo ngập lụt sông Hương (Thừa Thiên - Huế, 2008); lập bản đồ ngập lụt cho kế hoạch ứng phó khẩn cấp (EPP) hồ Phú Ninh, Kẻ Gỗ, Đồng Nghệ.
  • Ứng dụng LTĐTC và PTRR:
    • Mai Văn Công và cộng sự (2008) xây dựng quy trình phân tích rủi ro, thiết lập đường cong tần suất thiệt hại $FN$ cho hệ thống đê biển Việt Nam giai đoạn 1970–2007.
    • Lê Xuân Bảo (2016) áp dụng mô phỏng và đánh giá rủi ro ngập lụt hạ du sông Đồng Nai - Sài Gòn.
    • Nguyễn Văn Mạo (2006–2011) chủ trì đề tài cấp nhà nước nghiên cứu độ tin cậy đập chính hồ chứa miền Trung và Tây Nguyên.
    • Các nghiên cứu khác áp dụng phương pháp độ tin cậy Cấp độ I và II (chuẩn hóa biến ngẫu nhiên, tuyến tính hóa hàm trạng thái) cho cấu kiện đơn lẻ.

3. Khoảng trống luận án lựa chọn giải quyết

Các công trình trong nước trước đây chủ yếu giải bài toán độ tin cậy cho từng cấu kiện hoặc mặt cắt đơn lẻ theo phương pháp gần đúng (Cấp độ I, II) và phân tích rủi ro riêng lẻ cho một nguồn lũ. Luận án lựa chọn giải quyết bài toán ở cấp độ hệ thống theo phương pháp tiếp cận ngẫu nhiên Cấp độ III (mô phỏng Monte Carlo), tích hợp bài toán phân tích độ tin cậy hệ thống với bài toán phân tích rủi ro kinh tế và an toàn sinh mạng, có xét đến hiệu ứng chiều dài đê và tác động của biến đổi khí hậu đến năm 2050 cho hệ thống đê vùng ĐBSH.

Cơ sở lý thuyết và phương pháp nghiên cứu

1. Lý thuyết và khung phân tích

  • Lý thuyết độ tin cậy (LTĐTC):
    • Cơ chế sự cố được diễn toán thông qua hàm trạng thái giới hạn (hàm tin cậy) dạng: $$Z = R - S$$ Trong đó $R$ là khả năng chịu tải (độ bền), $S$ là tải trọng tác dụng. Trạng thái phá hoại xảy ra khi $Z < 0$ ($R < S$), trạng thái an toàn khi $Z > 0$.
    • Xác suất xảy ra sự cố ($P_f$) và xác suất an toàn ($P_s$) liên hệ qua biểu thức $P_s = 1 - P_f$.
    • Xác suất sự cố tổng quát: $$P_f = P(Z < 0) = \iint_{Z < 0} f_{R,S}(R,S) , dR , dS$$ Chỉ số độ tin cậy $\beta$ được xác định qua hàm phân phối chuẩn tích lũy: $P_f = \Phi(-\beta)$.
  • Phương pháp Cấp độ III - Mô phỏng ngẫu nhiên Monte Carlo (MCS): Sử dụng các hàm mật độ xác suất thực tế của các biến ngẫu nhiên đầu vào, rút mẫu trực tiếp để tính xác suất sự cố $P_f$ mà không cần tuyến tính hóa hàm $Z$.
  • Lý thuyết phân tích rủi ro (PTRR):
    • Rủi ro được lượng hóa bằng tích số giữa xác suất sự cố ($P_f$) và hậu quả thiệt hại ($D$).
    • Thiệt hại kinh tế của vùng bảo vệ được xác định qua tích phân diện tích ngập trên các ô lưới: $D = \sum D_i(h_i)$, kết hợp hàm thiệt hại (Damage Function - DF) theo độ sâu ngập $h$.
    • Tiêu chuẩn an toàn tối ưu kinh tế được xác định thông qua cực tiểu hóa hàm tổng chi phí vòng đời: $$C_{tot} = C(I_{P_f}) + PV(M_{P_f}) + PV(P_f \cdot D) \to \min$$ Trong đó $C(I_{P_f})$ là chi phí đầu tư nâng cấp đê; $PV(M_{P_f})$ là giá trị hiện tại của chi phí duy tu bảo dưỡng; $PV(P_f \cdot D)$ là giá trị hiện tại của rủi ro thiệt hại kinh tế theo tỷ lệ lãi suất hiệu quả $r$ và tuổi thọ công trình $T$.
  • Mô hình hóa hệ thống và Cây sự cố (Fault Tree Analysis - FTA): Mô hình hóa hệ thống đê theo cấu trúc nối tiếp (hỏng một đoạn dẫn đến hỏng toàn hệ thống), cấu trúc song song (tất cả các thành phần hỏng mới hỏng hệ thống), hoặc cấu trúc phức hợp.
  • Hiệu ứng chiều dài tuyến đê: Xem xét sự suy giảm độ tin cậy khi chiều dài tuyến đê tăng lên, dựa trên ma trận sự cố $m$ đoạn đê và $n$ cơ chế phá hoại, có xét khoảng cách độc lập và hệ số tương quan không gian của các biến tải trọng và biến độ bền.

2. Phương pháp nghiên cứu và công cụ phần mềm

  • BESTFIT: Phân tích thống kê chuỗi số liệu thực đo của các biến tải trọng và độ bền, kiểm định mức độ phù hợp và xác định quy luật phân phối xác suất (kỳ vọng toán, độ lệch chuẩn).
  • PROB2B: Tính toán chỉ số độ tin cậy và xác suất sự cố cho từng cơ chế phá hoại đơn lẻ của các đoạn đê theo phương pháp mô phỏng Monte Carlo.
  • OpenFTA: Thiết lập và giải bài toán cây sự cố hệ thống, xác định xác suất sự cố tổng thể $P_f$ cho toàn tuyến và toàn vùng bảo vệ.
  • Bộ phần mềm MIKE (MIKE 11, MIKE 21, MIKE Flood): Mô phỏng thủy lực dòng tràn, vỡ đê và lan truyền ngập lụt trong vùng bảo vệ để trích xuất bản đồ phân bố độ sâu ngập lụt.

3. Nguồn tư liệu và dữ liệu

  • Tài liệu khảo sát địa hình, địa chất công trình, hồ sơ mặt cắt ngang đê từ Sở Nông nghiệp và PTNT Hà Nội (2017) và Sở Nông nghiệp và PTNT Nam Định (2015, 2017).
  • Dữ liệu thống kê sự cố đê điều, các điểm đê xung yếu giai đoạn trước mùa lũ của Cục Quản lý đê điều và Phòng chống lụt bão (Vụ Quản lý Đê điều).
  • Chuỗi số liệu thủy văn mực nước lũ sông Hồng, triều dâng bão biển Đông và các kịch bản biến đổi khí hậu đến năm 2050.

Nội dung chính theo từng chương

Chương 1: Tổng quan về hệ thống đê phòng chống lũ, nghiên cứu rủi ro lũ lụt và an toàn đê điều

  • Lịch sử đê điều Việt Nam: Trải dài từ thời An Dương Vương (257 TCN), thời Lý đắp đê Cơ Xá (1088), thời Nguyễn đắp 145 km đê công và 400 km đê tư (1802–1829), thời Pháp thuộc đắp 87 triệu $\text{m}^3$ đất (ghi nhận 8 lần vỡ đê sông Hồng), và giai đoạn sau Cách mạng Tháng Tám năm 1945 theo Sắc lệnh số 70-SL (ngày 22/5/1946) đã đắp trên 260 triệu $\text{m}^3$ đất.
  • Quy mô hệ thống đê hiện nay: Khu vực Bắc Bộ và Khu 4 cũ có khoảng 5.000 km đê sông, đê cửa sông và đê biển (trong đó đê sông từ cấp III đến cấp đặc biệt khoảng 2.000 km).
    • Hà Nội: 20 tuyến đê chính dài khoảng 470 km (gồm hơn 37 km đê Hữu Hồng cấp đặc biệt, 211,5 km đê cấp I, 67,4 km đê cấp II, 87,3 km đê cấp III, 65,8 km đê cấp IV) và 27 tuyến đê bối dài 82,5 km.
    • Nam Định: Hệ thống đê sông và cửa sông dài 421,039 km (8,135 km cấp I, 80 km cấp II, 185,058 km cấp III) và 137,19 km đê biển.
  • Đánh giá hiện trạng an toàn đê ĐBSH:
    • Toàn vùng ĐBSH có khoảng 251 km đê hệ thống sông Hồng (119 km vật liệu kém, 132 km mối hoạt động) và 212 km đê hệ thống sông Thái Bình (176 km vật liệu kém, 36 km mối) ở tình trạng thân đê yếu; 226 km đê đi qua nền đất yếu (194 km sông Hồng, 32 km sông Thái Bình); 137 km yếu cả thân và nền; 75 km có ao đầm sát chân đê. Năm 2017 ghi nhận 168 điểm xung yếu và xảy ra 91 sự cố đê điều.
    • Cấu trúc địa chất nền đê gồm 3 tầng: Tầng 1 (tầng phủ đất dính dày 3–5m, hệ số thấm $k = 10^{-5} \div 10^{-6}\text{ cm/s}$); Tầng 2 (gồm đất dính, bùn sét chịu tải $0,4 \div 0,6\text{ kg/cm}^2$, hoặc tầng cát pha hạt bụi thông nước với sông có $k = 10^{-3} \div 10^{-2}\text{ cm/s}$); Tầng 3 (sét chặt phân bố ở độ sâu 10–30m).
TT Vị trí đê Hệ số rỗng $e_1$ Hệ số rỗng $e_2$ Hệ số rỗng $e_3$
1 Hữu Hồng - Sen Chiểu, Phúc Thọ, Hà Tây 0,726 1,377 1,223
2 Hữu Hồng - Thanh Chiểu, Phúc Thọ, Hà Tây 0,726 1,258 1,339
3 Đê tả Hồng - Mai Động, Đức Hợp, Hưng Yên 0,640 1,210 1,080
4 Đê hữu Hồng K13 - K14, Ba Vì, Hà Tây 0,740 1,360 0,830

Bảng 1-1: Hệ số độ rỗng của một số nền cát

TT Vị trí đoạn đê Thân đê $K\text{ (cm/s)}$ Tầng phủ $K\text{ (cm/s)}$ Tầng cát $K\text{ (cm/s)}$
1 Đê hữu Hồng – Phương Độ $1,8 \times 10^{-5}$ $1,8 \times 10^{-5}$ $39,1 \times 10^{-3}$
2 Đê hữu Hồng – Sen Chiểu $3,0 \times 10^{-6}$ $1,0 \times 10^{-5}$ $2,0 \times 10^{-3}$
3 Đê Hữu Hồng – Thanh Chiểu $3,0 \times 10^{-6}$ $1,0 \times 10^{-5}$ $3,0 \times 10^{-3}$
4 Đê tả Hồng – Mai Động $1,42 \times 10^{-6}$ $1,42 \times 10^{-5}$ $4,1 \times 10^{-3}$
5 Đê Hữu Hồng K13 – K14 $1,0 \times 10^{-5}$ $1,0 \times 10^{-4}$ $1,0 \times 10^{-3}$
6 Đê Hữu Hồng – Nhật Tân, Hà Nội $5,2 \times 10^{-6}$ $3,4 \times 10^{-4}$ $8,4 \times 10^{-3}$
7 Đê Hữu Hồng, Nam Định $6,5 \times 10^{-5}$ $9,2 \times 10^{-4}$ $8,7 \times 10^{-3}$

Bảng 1-2: Hệ số thấm các lớp đất chính của một số đoạn đê

  • Phân tích bất cập của tiêu chuẩn hiện hành: Quy định phân cấp đê dựa vào 3 tiêu chí (số dân/diện tích bảo vệ, lưu lượng lũ thiết kế, độ ngập sâu trung bình) chưa phản ánh đúng giá trị tài sản kinh tế của vùng bảo vệ và chưa lượng hóa tác động biến đổi khí hậu làm tăng điều kiện biên khí tượng thủy văn thêm 5% đến 10% đến năm 2050.

Chương 2: Phương pháp luận phân tích an toàn và xác định độ tin cậy yêu cầu cho hệ thống đê

  • Trình bày cơ sở toán học của LTĐTC theo Cấp độ III (Monte Carlo) để giải phương trình tích phân xác suất kết hợp giữa hàm mật độ xác suất tải trọng $f_S(S)$ và hàm mật độ xác suất độ bền $f_R(R)$.
  • Thiết lập quy trình phân tích rủi ro lũ lụt thông qua tích hợp chuỗi: Mô hình thủy văn $\to$ Mô phỏng ngập lụt thủy lực $\to$ Bản đồ phân bố độ sâu ngập $\to$ Đường cong hàm thiệt hại $f_j(h_i) \to$ Lượng hóa thiệt hại kinh tế kỳ vọng $E(D)$ và rủi ro cá nhân/cộng đồng.
  • Xây dựng mô hình cây sự cố cho các hệ thống liên kết nối tiếp, song song và phức hợp; quy định các cổng logic AND/OR và quy tắc tổ hợp xác suất phá hủy của hệ thống đê nhiều thành phần.

Chương 3: Thiết lập bài toán xác định chỉ số an toàn và độ tin cậy yêu cầu cho hệ thống đê

  • Sơ đồ hóa các cơ chế sự cố đê điều:
    • Đê sông: Chảy tràn đỉnh đê; Mất ổn định cấu kiện bảo vệ mái; Xói ngầm và đẩy trồi nền đê (mạch đùn, mạch sủi theo mô hình Bligh $c = c_B$); Trượt mái đê và mất ổn định tổng thể; Mất ổn định thấm thân và nền đê.
    • Đê biển: Chảy tràn do sóng và triều cường; Mất ổn định kết cấu bảo vệ mái phía biển; Xói lở chân đê biển; Mất ổn định thấm do chênh lệch mực nước triều.
  • Mô hình hóa hiệu ứng chiều dài: Phân chia tuyến đê thành $m$ đoạn đê độc lập tương đối và $n$ cơ chế phá hoại; xác định ma trận sự cố, khoảng cách độc lập không gian của các biến ngẫu nhiên tải trọng (mực nước, gió, sóng) và biến độ bền (cao trình đỉnh, hệ số thấm, cường độ kháng cắt).
  • Thuật toán xác định độ tin cậy yêu cầu (ĐTCYC):
    • Thiết lập hàm mục tiêu tối ưu kinh tế $C_{tot}(P_f) \to \min$.
    • Tích hợp điều kiện biên rủi ro chấp nhận về an toàn sinh mạng theo chuẩn rủi ro cá nhân ($10^{-5} \div 10^{-6}/\text{năm}$) và rủi ro tập thể (đường cong FN).

Chương 4: Ứng dụng phân tích độ tin cậy cho các hệ thống đê điển hình vùng Đồng bằng sông Hồng

  • Hệ thống 1 (HT1 - Đê Hữu Hồng bảo vệ trung tâm Hà Nội):
    • Phân chia tuyến đê thành 3 đoạn tính toán đặc trưng với các mặt cắt đại diện (Đoạn 1, Đoạn 2, Đoạn 3 theo hồ sơ Sở Nông nghiệp & PTNT Hà Nội 2017).
    • Phân tích độ tin cậy từng cơ chế sự cố bằng phần mềm PROB2B; đánh giá hệ số ảnh hưởng của các biến ngẫu nhiên đến cơ chế chảy tràn, mất ổn định bảo vệ mái, xói chân, xói ngầm/đẩy trồi và ổn định trượt mái.
    • Thiết lập cây sự cố HT1 bằng OpenFTA và xác định đường cong xác suất sự cố hệ thống khi có và không xét hiệu ứng chiều dài cho các kịch bản KB1, KB2.
    • Mô phỏng ngập lụt bằng bộ mô hình MIKE cho 3 trường hợp vỡ đê, thiết lập bản đồ ngập lụt và đường cong thiệt hại kinh tế khu vực trung tâm Hà Nội.
    • Tính toán hàm tổng chi phí $C_{tot}$ bao gồm chi phí đầu tư nâng cấp $I_H$, chi phí quản lý vận hành $PV(M)$ và chi phí rủi ro tiềm tàng $PV(P_f \cdot D)$, từ đó xác định tần suất đảm bảo phòng lũ và ĐTCYC tối ưu kinh tế cho hệ thống đê Hữu Hồng đến năm 2050.
  • Hệ thống 2 (HT2 - Hệ thống đê huyện Giao Thủy, tỉnh Nam Định):
    • Phân tích hệ thống phức hợp vòng khép kín gồm đê sông Hữu Hồng, đê cửa sông và đê biển Giao Thủy.
    • Tính toán chỉ số độ tin cậy riêng biệt cho từng tuyến đê (Bảng 4-23, 4-24, 4-25) và tổng hợp xác suất sự cố hệ thống nối tiếp (Bảng 4-26).
    • Thiết lập đường cong thiệt hại lịch sử, tính toán chi phí nâng cấp $I_H$, chi phí bảo dưỡng $PV(M)$ và chi phí rủi ro $R_{Pf}$ để tìm điểm tối ưu $C_{tot}$, đưa ra tần suất thiết kế và ĐTCYC phù hợp cho khu vực ven biển Giao Thủy.
  • Đề xuất giải pháp kỹ thuật: Kiến nghị xây dựng tuyến đê dự phòng tạo thành hệ thống song song 2 lớp đê và thiết lập các tuyến đê phụ phân chia khu vực ngập lụt thành các ô nhỏ để hạn chế lan truyền rủi ro.

Kết quả và những đóng góp mới

1. Những đóng góp mới về mặt lý luận

  • Thiết lập hoàn chỉnh khung phương pháp luận và quy trình tính toán tích hợp giữa Lý thuyết độ tin cậy Cấp độ III (phương pháp mô phỏng ngẫu nhiên Monte Carlo) và Phân tích rủi ro (PTRR) trong đánh giá mức độ an toàn cho hệ thống đê vùng đồng bằng.
  • Mở rộng bài toán phân tích từ phạm vi công trình/mặt cắt đơn lẻ sang phạm vi hệ thống hoàn chỉnh (kết hợp tương tác giữa tuyến đê và vùng được bảo vệ), tích hợp mô hình cây sự cố (FTA) và mô hình hóa định lượng hiệu ứng chiều dài tuyến đê có xét tương quan không gian của các biến ngẫu nhiên.
  • Xây dựng phương pháp xác định độ tin cậy yêu cầu (ĐTCYC) của hệ thống đê nâng cấp dựa trên bài toán tối ưu hóa kinh tế tổng chi phí vòng đời ($C_{tot} \to \min$) kết hợp với chuẩn rủi ro chấp nhận về sinh mạng cộng đồng trong bối cảnh biến đổi khí hậu và phát triển kinh tế - xã hội đến năm 2050.

2. Những đóng góp mới về mặt thực tiễn

  • Lượng hóa chính xác chỉ số độ tin cậy hiện tại và xác định độ tin cậy yêu cầu đến năm 2050 cho hai hệ thống đê điển hình vùng ĐBSH: hệ thống đê Hữu Hồng (bảo vệ trung tâm Hà Nội) và hệ thống đê biển - đê sông Giao Thủy (Nam Định).
  • Đưa ra cơ sở khoa học định lượng giúp các cơ quan quản lý nhà nước xác định quy mô đầu tư, phân kỳ nâng cấp đê hợp lý, tránh lãng phí vốn đầu tư công hoặc thiếu hụt an toàn phòng lũ.
  • Đề xuất các giải pháp kỹ thuật công trình nhằm giảm thiểu rủi ro: xây dựng tuyến đê dự phòng tạo hệ thống song song 2 lớp đê và thiết lập các tuyến đê phụ chia nhỏ lưu vực ngập lụt.

Hạn chế và hướng nghiên cứu tiếp

1. Hạn chế của nghiên cứu

  • Chưa xem xét quá trình suy giảm độ bền cơ lý của vật liệu đắp thân đê và đất nền theo thời gian (quá trình lão hóa, suy thoái địa kỹ thuật).
  • Bỏ qua các tổ hợp tải trọng đặc biệt bất thường như động đất, dư chấn địa chất và tác động xói mòn do mưa lớn cục bộ kéo dài trực tiếp trên thân đê.
  • Chưa đi sâu mô hình hóa xác suất sự cố chi tiết của các công trình ngầm xuyên đê (như cống tiêu dưới đê, trạm bơm, cửa khẩu qua đê) trong bài toán ứng dụng thực tế.

2. Hướng phát triển tiếp theo

  • Mở rộng bài toán phân tích độ tin cậy theo thời gian (Time-variant reliability) có kể đến quy luật suy thoái cơ lý của vật liệu đắp và xói lở nền tích lũy qua các chu kỳ lũ.
  • Tích hợp chi tiết các cơ chế phá hỏng và xác suất sự cố của các cụm công trình qua đê (cống, van điều tiết, khớp nối cống) vào cây sự cố tổng thể của hệ thống.
  • Mở rộng phạm vi ứng dụng phương pháp cho toàn bộ các tuyến đê thuộc hệ thống sông Thái Bình, hệ thống đê sông và đê biển các tỉnh Bắc Trung Bộ và Đồng bằng sông Cửu Long.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị ứng dụng và tham khảo đối với nhiều đối tượng chuyên môn:

  • Cơ quan quản lý nhà nước và hoạch định chính sách: Bộ Nông nghiệp và Phát triển Nông thôn (Vụ Quản lý Đê điều, Tổng cục Phòng chống thiên tai), Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn các tỉnh ĐBSH có thể sử dụng phương pháp luận của luận án để rà soát, điều chỉnh quy định phân cấp đê và lập quy hoạch đầu tư công trình phòng chống lũ thích ứng với BĐKH.
  • Kỹ sư thiết kế, tư vấn công trình thủy lợi - thủy điện: Ứng dụng quy trình tính toán độ tin cậy Cấp độ III (sử dụng phối hợp BESTFIT, PROB2B, OpenFTA và MIKE) để tính toán kiểm tra an toàn và xác định kích thước mặt cắt đê tối ưu thay cho phương pháp hệ số an toàn truyền thống.
  • Nghiên cứu sinh, giảng viên và học viên cao học: Tham khảo khung phân tích rủi ro thiên tai, lý thuyết độ tin cậy hệ thống và mô hình hóa cây sự cố công trình thủy trong giảng dạy và nghiên cứu chuyên sâu ngành Xây dựng Công trình thủy, Kỹ thuật Tài nguyên nước và Địa kỹ thuật công trình.

Câu hỏi thường gặp

1. Sự khác biệt cốt lõi giữa phương pháp đánh giá an toàn truyền thống và phương pháp lý thuyết độ tin cậy trong luận án là gì?

Phương pháp truyền thống đánh giá an toàn dựa trên hệ số an toàn tĩnh $K$ tại một mặt cắt đại diện với các giá trị tải trọng và độ bền cố định ứng với tần suất thiết kế định mức. Phương pháp lý thuyết độ tin cậy (Cấp độ III) của luận án mô hình hóa tải trọng và độ bền dưới dạng các biến ngẫu nhiên theo phân phối xác suất thực tế, giải tích phân xác suất để xác định xác suất xảy ra sự cố tổng hợp $P_f$ và chỉ số độ tin cậy $\beta$ cho toàn bộ hệ thống công trình.

2. Hiệu ứng chiều dài tuyến đê được mô hình hóa và ảnh hưởng thế nào đến an toàn hệ thống đê?

Trong thực tế, khi chiều dài tuyến đê tăng lên, xác suất xuất hiện các điểm xung yếu hoặc xảy ra sự cố phá hoại tại một vị trí bất kỳ trên tuyến đê sẽ gia tăng. Luận án mô hình hóa hiệu ứng chiều dài bằng cách phân chia tuyến đê thành $m$ đoạn đê độc lập tương đối với $n$ cơ chế sự cố, xem xét khoảng cách độc lập không gian và hệ số tương quan giữa các biến tải trọng (thủy văn, sóng, gió) và biến độ bền (địa chất, hình học đê) để xác lập chính xác đường cong phá hủy của toàn tuyến đê.

3. Bộ công cụ và phần mềm mô hình toán được luận án sử dụng gồm những gì?

Luận án sử dụng 4 công cụ mô hình toán chuyên sâu:

  1. BESTFIT: Phân tích thống kê và xác định quy luật phân phối xác suất của các chuỗi số liệu biên tải trọng và độ bền.
  2. PROB2B: Phân tích độ tin cậy cho từng cơ chế sự cố thành phần của mặt cắt đê theo mô phỏng Monte Carlo.
  3. OpenFTA: Phân tích cây sự cố logic để tính toán xác suất sự cố tổng thể của hệ thống.
  4. Bộ mô hình MIKE (MIKE 11, MIKE 21, MIKE Flood): Mô phỏng diễn biến thủy lực ngập lụt khi xảy ra tràn hoặc vỡ đê và xây dựng bản đồ độ sâu ngập.

4. Tiêu chí xác định độ tin cậy yêu cầu (ĐTCYC) tối ưu cho hệ thống đê được thiết lập như thế nào?

ĐTCYC được xác định thông qua bài toán tối ưu hóa kinh tế: tìm giá trị xác suất sự cố $P_f$ sao cho tổng chi phí vòng đời $C_{tot} = C(I_{P_f}) + PV(M_{P_f}) + PV(P_f \cdot D)$ đạt cực tiểu (cân bằng giữa vốn đầu tư xây lắp nâng cấp với giá trị rủi ro thiệt hại kinh tế được giảm thiểu), đồng thời thỏa mãn ràng buộc về mức rủi ro chấp nhận đối với tính mạng con người trong cộng đồng.

Kết luận

Luận án tiến sĩ của tác giả Trần Quang Hoài đã hoàn thành việc xây dựng khung phương pháp luận khoa học và quy trình thuật toán đánh giá an toàn hệ thống đê điều dựa trên sự tích hợp giữa Lý thuyết độ tin cậy Cấp độ III và Phân tích rủi ro. Công trình đã giải quyết bài toán an toàn ở cấp độ hệ thống, tích hợp mô hình cây sự cố, hiệu ứng chiều dài tuyến đê và kịch bản biến đổi khí hậu đến năm 2050. Kết quả tính toán ứng dụng cụ thể cho hệ thống đê Hữu Hồng (Hà Nội) và đê biển Giao Thủy (Nam Định) cung cấp cơ sở dữ liệu định lượng phục vụ công tác quy hoạch, ra quyết định đầu tư nâng cấp và quản lý an toàn đê điều vùng đồng bằng sông Hồng.