Tổng quan luận án

Tính cấp thiết của đề tài

Đồng bằng sông Hồng (ĐBSH) là vùng lãnh thổ hạ lưu sông Hồng thuộc miền Bắc Việt Nam, bao gồm 10 tỉnh và thành phố: Hà Nội, Bắc Ninh, Hải Dương, Hưng Yên, Hải Phòng, Hà Nam, Nam Định, Ninh Bình, Thái Bình và Vĩnh Phúc. Đây là trung tâm chính trị, kinh tế, văn hóa và xã hội đặc biệt quan trọng của cả nước, trong đó có thủ đô Hà Nội. Lũ lụt tại khu vực này luôn là loại hình thiên tai gây thiệt hại nghiêm trọng về người và tài sản, điển hình như trận lũ lịch sử năm 1971 do vỡ đê sông Hồng.

Hệ thống công trình phòng chống lũ bảo vệ vùng châu thổ sông Hồng gồm các hồ chứa cắt lũ ở thượng nguồn và hệ thống đê sông Hồng, sông Thái Bình với tổng chiều dài khoảng 3.000 km. Do đặc điểm lịch sử xây dựng qua nhiều thời kỳ, thân đê được tôn cao và áp trúc nhiều lần để chống lũ tràn, trong khi nền đê hầu như chưa được xử lý triệt để. Khi xuất hiện lũ lớn kéo dài, các đoạn đê có nền địa chất yếu thường phát sinh nhiều sự cố biến dạng thấm nguy hiểm như đùn sủi, xói ngầm.

Trong bối cảnh biến đổi khí hậu (BĐKH) và nước biển dâng (NBD), các hình thái thời tiết cực đoan gia tăng, mưa lớn bất thường ở thượng lưu làm giảm hiệu quả cắt lũ của các hồ chứa, đồng thời mực nước biển dâng cản trở khả năng thoát lũ hạ du. Áp lực mực nước lũ duy trì cao và kéo dài làm gia tăng nguy cơ mất an toàn cho hệ thống đê điều, đòi hỏi phải có cơ sở khoa học để đánh giá định lượng độ an toàn của tuyến đê Hữu Hồng nhằm đề xuất các giải pháp công trình và phi công trình phù hợp.

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

Luận án thực hiện nhằm đánh giá mức độ an toàn của hệ thống đê Hữu sông Hồng đoạn từ Sơn Tây về Phú Xuyên (thành phố Hà Nội) trong điều kiện biến đổi khí hậu; trọng tâm nghiên cứu tập trung vào xác định điều kiện biên thủy lực, nghiên cứu quy luật quá trình xói ngầm dưới nền đê, và đánh giá an toàn hệ thống đê kết hợp cống dưới đê theo lý thuyết độ tin cậy và phân tích rủi ro nhằm đề xuất giải pháp giảm thiểu ngập lụt tối ưu.

Câu hỏi nghiên cứu

  1. Ảnh hưởng của biến đổi khí hậu theo các kịch bản khí hậu tác động như thế nào đến chế độ thủy văn trên sông Hồng trong phạm vi nghiên cứu?
  2. Cơ chế sự cố nào có nguy cơ gây mất ổn định hệ thống đê sông Hồng?
  3. Quá trình xói ngầm dưới nền đê sông Hồng diễn biến và tác động như thế nào thông qua mô hình vật lý trong phòng và phân tích thống kê hiện trường?
  4. Mức độ an toàn của hệ thống đê sông Hồng hiện tại trong điều kiện biến đổi khí hậu và nước biển dâng được định lượng như thế nào khi áp dụng lý thuyết độ tin cậy và phân tích rủi ro?

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

  • Đối tượng nghiên cứu: Hệ thống đê Hữu sông Hồng trên địa bàn thành phố Hà Nội, xét trong điều kiện địa chất công trình, địa chất thủy văn và điều kiện biên thủy lực có ảnh hưởng của biến đổi khí hậu.
  • Phạm vi không gian: Tuyến đê Hữu Hồng thuộc địa phận Hà Nội từ Km30+550 (khu vực Sơn Tây) đến Km117+900 (khu vực Phú Xuyên) với tổng chiều dài khoảng 100 km (bao gồm 37,7 km đê cấp đặc biệt và các đoạn đê cấp I), 17 đoạn đê và 10 cống dưới đê tiêu biểu. Khu vực phân tích rủi ro ngập lụt chuyên sâu là đoạn đê Sen Chiểu thuộc huyện Phúc Thọ (Km31+500 ÷ Km33+500).
  • Phạm vi nội dung: Đánh giá an toàn hệ thống đê qua các cơ chế mất ổn định nền và tràn đỉnh (không xét sự cố do phá hỏng kết cấu thân đê phía trên) dưới tác động của các kịch bản BĐKH và NBD cập nhật năm 2016.

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

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

Các nghiên cứu trên thế giới về an toàn đê điều và công trình phòng lũ tập trung vào hai hướng tiếp cận chính: hoàn thiện phương pháp luận đánh giá an toàn dựa trên rủi ro và nghiên cứu cơ chế phá hoại do thấm và xói ngầm.

  • Về phương pháp luận và lý thuyết độ tin cậy: Tại Hà Lan, sau trận lũ năm 1953 (chương trình Delta), các phương pháp phân tích độ tin cậy và rủi ro được phát triển mạnh mẽ. Năm 2005, chương trình VNK05 đã đánh giá lại toàn bộ an toàn hệ thống phòng lũ trong bối cảnh BĐKH, định kỳ 5 năm tái đánh giá. Năm 2008, Ủy ban nghiên cứu về đồng bằng Hà Lan đề xuất kế hoạch ứng phó tới năm 2100 với kịch bản nước biển dâng từ 0,55 m đến 1,30 m. Các tác giả Vrijling và Van Gelder đã xây dựng các mô hình xác định cận xác suất phá hỏng của hệ thống nối tiếp. Lopez De La Cruz nghiên cứu mối quan hệ giữa chỉ số độ tin cậy yêu cầu ($\beta$) và chiều dài hệ thống. Tại Mỹ, Hiệp hội Kỹ sư Quân đội Mỹ (USACE) phát triển các tiêu chuẩn tính toán thiết kế đê, đập dựa trên lý thuyết rủi ro sau thảm họa bão Katrina năm 2005 tại New Orleans.
  • Về cơ chế xói ngầm và biến dạng thấm:
    • Bligh (1910) và Lane (1935) đề xuất các phương pháp kinh nghiệm kiểm toán xói ngầm dựa trên chiều dài đường thấm trọng số và cột nước tới hạn.
    • Tsugaev (1998) xác định gradient áp lực thấm cho phép dựa trên việc ngăn chặn dòng thấm tập trung trong các khuyết tật của thân và nền đất.
    • Sellmeijer (1988) thiết lập mô hình xói ngầm dựa trên điều kiện cân bằng giới hạn của các hạt cát trong ống xói, xét đến ảnh hưởng của đường kính hạt, bề dày tầng chứa nước và hệ số thấm. Vera van (2011) tiếp tục thử nghiệm mô hình Sellmeijer trên mô hình vật lý.
    • Schmertmann (2000) tại Đại học Florida (Mỹ) tiến hành thí nghiệm xói ngầm xác định ảnh hưởng của hệ số không đều hạt ($C_u$), độ chặt và gradient thủy lực. Các nghiên cứu tương tự tại Đại học Okayama (Nhật Bản, 1987) và tại Đức cũng làm rõ mối quan hệ giữa thành phần hạt và gradient tới hạn.
    • Wan và Fell (Australia) nghiên cứu tốc độ xói ngầm trong đất dính. V. Kozlova (1934) và Ixtomina (1952, 1957) xác định gradient cho phép dựa trên tỷ lệ đường kính hạt và hệ số không đều hạt $C_u$.

Các nghiên cứu trong nước

  • Nguyễn Trấn (1980s) và Nghiêm Hữu Hạnh nghiên cứu áp lực thủy động của dòng thấm có áp một chiều dưới đê sông Hồng bằng mô hình tương tự điện thủy động và hàm truyền áp $R(\lambda)$.
  • Nguyễn Hữu An và các cộng sự ứng dụng phương pháp giải tích, sai phân hữu hạn và phần mềm Geo-Studio (Seep/W) để tính thấm ổn định phẳng cho đê nền cát thông với sông.
  • Phạm Văn Quốc nghiên cứu mô hình thấm không ổn định qua đê bằng mô hình khe hẹp; chứng minh đường bão hòa thấm và áp lực thấm lớn nhất tại chân đê phía đồng xuất hiện chậm sau đỉnh lũ từ 1,5 đến 2 ngày; thiết lập biểu thức tính hệ số an toàn chống đẩy bục tầng phủ có xét lực dính.
  • Tô Xuân Vu phân loại các hình thức biến dạng thấm (đùn đất, cát chảy, xói ngầm) theo tỷ số $\xi = m_p/d$, xác định gradient giới hạn và kiến nghị các giải pháp tường chống thấm, sân phủ, giếng giảm áp.
  • Bùi Văn Trường phân chia cấu trúc nền đê sông Hồng (Thái Bình) thành 5 kiểu, 11 phụ kiểu, 23 dạng cấu trúc; tiến hành thí nghiệm hiện trường xác định 3 giai đoạn phát triển biến dạng thấm và đề xuất cừ bản nhựa.
  • Trịnh Minh Thụ và các cộng sự sử dụng mô hình Seep 3D mô phỏng xử lý đùn sủi tại Mai Động - Đức Hiệp (Hưng Yên), chứng minh giếng giảm áp có hiệu quả xử lý áp lực thấm cao hơn tường hào bentonite.
  • Nguyễn Quốc Đạt lập bản đồ phân loại nền đê Hữu Hồng, tả Đáy tại Hà Nam; thử nghiệm công nghệ cọc vây xi măng đất và khoan phụt vữa hóa chất kết hợp xi măng xử lý thấm khẩn cấp.
  • Phạm Quang Tú ứng dụng lý thuyết độ tin cậy đánh giá đê Hữu Hồng theo biên thủy lực của Hà Văn Khối và mô hình xói ngầm của Bùi Văn Trường, tuy nhiên chưa tiến hành thực nghiệm phát triển ống xói theo thời gian và chưa đưa cống dưới đê vào bài toán hệ thống.

Khoảng trống nghiên cứu

Các nghiên cứu trong nước trước đây chủ yếu tiếp cận hiện tượng xói ngầm đơn lẻ hoặc tính toán thấm theo phương pháp tiền định truyền thống, chưa xét đến sự biến đổi chiều dài ống xói theo thời gian thực tế trong một trận lũ. Đồng thời, các nghiên cứu chưa tích hợp đầy đủ điều kiện biên thủy lực cực đoan do BĐKH, chưa đưa toàn bộ các cống dưới đê vào đánh giá độ tin cậy hệ thống và chưa áp dụng nguyên lý phân tích rủi ro kinh tế để lựa chọn phương án gia cố đê tối ưu.

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

Khung lý thuyết và khái niệm sử dụng

  1. Lý thuyết độ tin cậy (LTĐTC): Đánh giá an toàn công trình thông qua hàm trạng thái giới hạn $Z = R(x_i) - L(y_i)$, trong đó $R(x_i)$ là hàm sức kháng của công trình và $L(y_i)$ là hàm tải trọng tác động. Xác suất sự cố được định nghĩa là $P_f = P(Z < 0) = \Phi(-\beta)$, với $\beta$ là chỉ số độ tin cậy ($\beta = \mu_Z / \sigma_Z$). Luận án sử dụng phương pháp độ tin cậy Cấp độ II (phương pháp xấp xỉ bậc nhất FORM/SORM) và Cấp độ III (mô phỏng ngẫu nhiên Monte Carlo), đồng thời phân tích hệ thống nối tiếp xét đến hiệu ứng độ dài của tuyến đê ($P_{fLi}$).
  2. Các cơ chế sự cố được mô hình hóa:
    • Sự cố tràn đỉnh đê và tràn đỉnh cống ($Z = Z_{đđ} - (Z_{TL} + \Delta h)$).
    • Sự cố đẩy bục tầng phủ hạ lưu đê và hạ lưu cống ($Z = H_c - H$).
    • Sự cố xói ngầm dưới nền đê và dưới nền cống ($Z = [J] - J$ hoặc mô hình Sellmeijer hiệu chỉnh).
    • Sự cố trượt mái dốc hạ lưu đê ($Z = FS - 1$) và sự cố nền cống mất khả năng chịu tải.
  3. Mô hình xói ngầm phát triển theo thời gian: Xây dựng phương trình mô tả sự phát triển của chiều dài ống xói $L_{x,t}$ theo thời gian $T$ và chênh cao cột nước $\Delta H$, từ đó tính toán sự suy giảm chiều dài đường thấm hiệu dụng $L_t = L_0 - L_{x,t}$.
  4. Nguyên lý rủi ro và tối ưu hóa: Rủi ro ($Risk$) là tích số giữa xác suất sự cố ($P_f$) và tổng thiệt hại kinh tế do sự cố gây ra ($D$). Phương án thiết kế tối ưu được xác định tại điểm cực tiểu của hàm tổng chi phí bao gồm chi phí đầu tư xây dựng và chi phí rủi ro tổn thất.

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

  • Phương pháp thực nghiệm mô hình vật lý trong phòng: Chế tạo mô hình thí nghiệm thấm ngang và thấm đứng tại Phòng Thí nghiệm Địa kỹ thuật (Trường Đại học Thủy Lợi) để quan sát diễn biến xói ngầm, đo đạc gradient áp lực thấm tới hạn trên các mẫu cát lấy từ hiện trường.
  • Phương pháp khảo sát địa vật lý hiện trường: Sử dụng phương pháp đo sâu điện địa vật lý dọc các tuyến đo tại đoạn đê Sen Chiểu để xác định các dị thường cấu trúc và vùng cát xốp dưới nền đê.
  • Phương pháp phân tích ngược và thống kê hiện trường: Thu thập dữ liệu 65 điểm sự cố đùn sủi trên đê Hữu Hồng (1978–2010) và dữ liệu vỡ đê Vân Cốc năm 1986 để xác định hệ số mô hình $\alpha$ hiệu chỉnh giữa kết quả thực nghiệm trong phòng và thực tế hiện trường.
  • Phương pháp mô hình số: Sử dụng phần mềm Geo-Studio 2016 (mô-đun Seep/W và Slope/W) để phân tích trường thấm và hệ số ổn định trượt mái đê trong các điều kiện thủy lực biến đổi.
  • Phương pháp phân tích hệ thống và mô phỏng xác suất: Xây dựng cây sự cố (Fault Tree) cho hệ thống đê Hữu Hồng gồm 17 đoạn đê và 10 cống dưới đê.

Dữ liệu phục vụ nghiên cứu

  • Dữ liệu địa chất công trình, địa tầng mặt cắt dọc và ngang của tuyến đê Hữu Hồng từ Sơn Tây đến Phú Xuyên.
  • Dữ liệu thủy văn thực đo tại trạm Sơn Tây và trạm Hà Nội; quy trình vận hành liên hồ chứa trên lưu vực sông Hồng (Quyết định 1622/QĐ-TTg).
  • Kịch bản biến đổi khí hậu và nước biển dâng cho Việt Nam do Bộ Tài nguyên và Môi trường công bố năm 2016 (trọng tâm là kịch bản RCP4.5).
  • Bản đồ hiện trạng sử dụng đất năm 2016 của huyện Phúc Thọ và mô hình số độ cao (DEM) vùng nghiên cứu Sen Chiểu.

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

Chương 1: Tổng quan về nghiên cứu an toàn đê trong bối cảnh biến đổi khí hậu

Chương này tổng kết lịch sử phát triển của hệ thống đê điều hạ lưu sông Hồng từ thời phong kiến (năm 1077, 1108, 1248), qua thời kỳ Pháp thuộc (1873–1945), giai đoạn 1945–1954, 1954–1965, trận lũ lịch sử 1971, giai đoạn 1975–2000, đến giai đoạn từ năm 2000 đến nay.

Về hệ thống hồ chứa thượng lưu, chương tổng hợp thông số 05 hồ chứa lớn đa mục tiêu trên lưu vực sông Hồng:

  • Hồ Thác Bà (sông Chảy): MNDBT +58,0 m; MNDGC +61,0 m; dung tích phòng lũ hạ du 0,45 tỷ $\text{m}^3$.
  • Hồ Hòa Bình (sông Đà): MNDBT +117,0 m; MNDGC +122,0 m; dung tích phòng lũ hạ du 2,90 tỷ $\text{m}^3$.
  • Hồ Tuyên Quang (sông Gâm): MNDBT +120,0 m; MNDGC +122,55 m; dung tích phòng lũ hạ du 1,00 tỷ $\text{m}^3$.
  • Hồ Sơn La (sông Đà): MNDBT +215,0 m; MNDGC +217,83 m; dung tích phòng lũ hạ du (kết hợp với Hòa Bình) 7,00 tỷ $\text{m}^3$.
  • Hồ Lai Châu (sông Đà): MNDBT +295,0 m; mực nước lớn nhất thiết kế +297,0 m.

Theo Quyết định số 1622/QĐ-TTg (Quy trình vận hành liên hồ chứa sông Hồng), việc vận hành hệ thống hồ chứa phải bảo đảm khống chế mực nước sông Hồng tại trạm thủy văn Hà Nội không vượt quá +13,10 m đối với lũ chu kỳ 300 năm và không vượt quá +13,40 m đối với lũ chu kỳ 500 năm tại Sơn Tây.

Mực nước khống chế tại Hà Nội Dung tích cắt lũ hồ Sơn La ($\text{tỷ m}^3$) Dung tích cắt lũ hồ Hòa Bình ($\text{tỷ m}^3$) Tổng dung tích cắt lũ ($\text{tỷ m}^3$)
$Z_{HN} \le +11,5\text{ m}$ 1,00 (mực nước hồ 197,3 ÷ 203 m) 1,13 (mực nước hồ 101 ÷ 107 m) 2,13
$Z_{HN} > (+13,1 \div +13,4)\text{ m}$ 3,00 (mực nước hồ 203 ÷ 217,2 m) 1,87 (mực nước hồ 107 ÷ 117 m) 4,87

Khảo sát hiện trạng tuyến đê Hữu Hồng (từ Km30+550 đến Km117+900) cho thấy cao trình đỉnh đê hiện tại cao hơn mực nước lũ thiết kế chu kỳ 500 năm từ 0,95 m đến 2,40 m (ngoại trừ đoạn đê Vân Cốc). Khảo sát 10 cống dưới đê tiêu biểu (2 cống lộ thiên, 8 cống ngầm) ghi nhận tình trạng suy thoái vật liệu và thiết bị đóng mở xuống cấp. Phân tích 65 trường hợp sự cố đùn sủi ghi nhận từ năm 1978 đến 2010 cho thấy gradient áp lực thấm tính toán trung bình $J = H/L$ xuất hiện đùn sủi chỉ dao động từ $0,006$ đến $0,124$ (không đùn cát: $J = 0,015 \div 0,073$; có đùn cát: $J = 0,006 \div 0,124$). Giá trị $J$ thực tế rất nhỏ này cho thấy có sự suy giảm đáng kể chiều dài đường thấm hiệu dụng dưới nền đê do sự phát triển của ống xói ngầm.

Chương 2: Cơ sở khoa học nghiên cứu, đánh giá an toàn đê và giải pháp tăng cường ổn định đê

Chương này thiết lập cơ sở lý thuyết và công cụ phân tích an toàn đê:

  • Điều kiện biên thủy lực: Phân tích quan hệ lưu lượng – mực nước $Q = f(H)$ và biến đổi lòng dẫn tại các trạm Sơn Tây và Hà Nội. Thiết lập sơ đồ tính toán dòng chảy lũ hệ thống sông Hồng có xét đến kịch bản biến đổi khí hậu (RCP4.5) và quy trình vận hành liên hồ chứa.
  • Cấu trúc nền đê: Phân chia địa tầng nền đê Hữu Hồng từ Sơn Tây đến Phú Xuyên thành các mô hình cấu trúc đại diện dựa trên mức độ nguy hiểm về thấm và sự tồn tại của các tầng cát thấm nước mạnh bên dưới lớp phủ dính.
  • Đánh giá an toàn theo phương pháp truyền thống: Thiết lập quy trình đánh giá dựa trên tiêu chuẩn thiết kế hiện hành, kiểm tra hệ số an toàn ổn định mái $[K]$ và gradient thấm cho phép $[J]$.
  • Đánh giá an toàn theo lý thuyết độ tin cậy:
    • Định nghĩa các bất định trong địa kỹ thuật (bất định số liệu đất nền, bất định mô hình tính toán, biến thiên tải trọng thủy văn).
    • Phân tích độ tin cậy cấp độ II (FORM) để xác định chỉ số độ tin cậy $\beta$ và hệ số ảnh hưởng $\alpha_i$ của các biến ngẫu nhiên.
    • Phân tích độ tin cậy cấp độ III (Monte Carlo) để kiểm chứng và giải quyết các bài toán phi tuyến.
    • Xây dựng mô hình xác suất sự cố của hệ thống nối tiếp có xét đến hiệu ứng độ dài đê theo công thức: $$P_{f,Li} = 1 - (1 - P_{f,sec_i})^{L_i / l_i}$$
    • Phân tích sự khác biệt cơ bản giữa phương pháp hệ số an toàn truyền thống và phương pháp xác suất theo lý thuyết độ tin cậy.
  • Cơ sở mô hình vật lý xói ngầm và nguyên lý rủi ro: Đề xuất nguyên lý mô phỏng thấm trong phòng và quy trình phân tích rủi ro kinh tế phục vụ thiết kế tối ưu công trình bảo vệ.

Chương 3: Nghiên cứu xói ngầm dưới nền đê bằng mô hình vật lý trong phòng

Chương này trình bày chi tiết quá trình thí nghiệm mô hình vật lý trong phòng để nghiên cứu quy luật phát triển của ống xói ngầm:

  • Thiết kế mô hình thí nghiệm: Xây dựng hai mô hình vật lý gồm mô hình thấm ngang và mô hình thấm đứng. Mẫu thí nghiệm sử dụng cát lấy tại hai khu vực trọng điểm đùn sủi trên tuyến đê Hữu Hồng là Sen Chiểu (Phúc Thọ) và Hoàng Mai. Cát được thí nghiệm ở ba trạng thái độ chặt: xốp ($D_c$ thấp), chặt vừa, và chặt.
  • Đường dâng nước thí nghiệm: Thiết lập đường dâng nước mô phỏng theo đường lũ lên chu kỳ 500 năm tương ứng với kịch bản biến đổi khí hậu RCP4.5.
  • Kết quả thí nghiệm thấm ngang:
    • Ghi nhận diễn biến phát triển của ống xói qua các thời điểm: $t_s$ (thời điểm quan sát thấy ống xói bắt đầu hình thành ở hạ lưu), $t_1$ (thời điểm cát bắt đầu dịch chuyển mạnh), $t_2$ (thời điểm ống xói đạt chiều dài giới hạn $L_{x\text{-}gh}$), và $t_f$ (thời điểm ống xói nối thông hoàn toàn từ hạ lưu lên thượng lưu gây phá hoại tổng thể).
    • Xây dựng các đường quan hệ giữa chênh áp thủy lực và thời gian ($\Delta H - T$), tương quan chiều dài ống xói phát triển theo thời gian ($L_x - T$) cho cát Sen Chiểu và cát Hoàng Mai ở các cấp độ chặt khác nhau.
  • Kết quả thí nghiệm thấm đứng: Khảo sát hiện tượng cát sôi và bục tầng phủ thẳng đứng; xác định mối quan hệ giữa gradient áp lực thấm $J$, chênh áp $\Delta H$ và thời gian duy trì cột nước $T$ đối với từng độ chặt của cát.
  • Kết luận thực nghiệm: Tốc độ phát triển của ống xói tỷ lệ thuận với độ chênh cột nước $\Delta H$ và tỷ lệ nghịch với độ chặt tương đối $D_c$ của đất cát. Khi ống xói phát triển từ hạ lưu ngược về thượng lưu, chiều dài đường thấm thực tế $L_t$ bị rút ngắn dần, dẫn đến việc gradient thủy lực cục bộ tại đầu ống xói tăng cao đột ngột, thúc đẩy quá trình phá hoại diễn ra nhanh chóng.
Loại mẫu cát Trạng thái độ chặt Hiện tượng quan sát Diễn biến phát triển ống xói ($L_x - T$)
Cát Sen Chiểu Xốp Xuất hiện $t_s$ sớm, cát dịch chuyển nhanh Ống xói phát triển nhanh, đạt $L_{x\text{-}gh}$ trong thời gian ngắn
Cát Sen Chiểu Chặt vừa Cần gradient $\Delta H$ cao hơn để khởi phát Ống xói phát triển ổn định theo từng cấp tăng áp lực
Cát Sen Chiểu Chặt Thời điểm $t_s$ xuất hiện muộn nhất Ống xói chỉ mở rộng khi $\Delta H$ đạt giá trị tới hạn cao
Cát Hoàng Mai Xốp / Chặt vừa / Chặt Hiện tượng xói ngầm tương tự cát Sen Chiểu Tốc độ phát triển phụ thuộc vào thành phần hạt và độ chặt $D_c$

Chương 4: Độ tin cậy của hệ thống đê sông Hồng trong bối cảnh biến đổi khí hậu

Chương 4 tích hợp toàn bộ các kết quả lý thuyết, số liệu thực nghiệm và khảo sát hiện trường vào tính toán an toàn cho hệ thống đê Hữu Hồng:

  • Đánh giá theo phương pháp truyền thống: Tính toán kiểm tra gradient áp lực thấm $J$ so với $[J]$ và hệ số ổn định mái $K_{min}$ so với $[K]$ cho 17 đoạn đê và 10 cống dưới đê Hữu Hồng.
  • Phân tích độ tin cậy của hệ thống đê:
    • Xây dựng cây sự cố chi tiết cho toàn hệ thống đê Hữu Hồng, từng đoạn đê và từng cống dưới đê.
    • Tính toán xác suất sự cố ($P_f$) và chỉ số độ tin cậy ($\beta$) cho 4 cơ chế phá hoại chính: tràn đỉnh, đẩy bục tầng phủ, xói ngầm, và trượt mái.
    • Xác định tỷ lệ ảnh hưởng ($\alpha_i$) của từng biến ngẫu nhiên (mực nước lũ $H$, hệ số thấm $K$, khối lượng riêng $\rho$, góc ma sát trong $\varphi$, lực dính $c$, chiều dày tầng phủ $d$).
    • Đánh giá xác suất sự cố toàn hệ thống khi xét đến hiệu ứng độ dài đê.
  • Phân tích an toàn dưới trận lũ thiết kế và mô hình ống xói:
    • Khảo sát địa vật lý tại đoạn đê Sen Chiểu (Km32+322 ÷ Km32+512) qua 3 tuyến đo T1, T2, T3; dựng hình ảnh 3D mô phỏng vùng cát xốp cục bộ.
    • Phân tích ngược sự cố vỡ đê Vân Cốc năm 1986 để xác định phân bố xác suất của hệ số mô hình $\alpha$ dự báo xói ngầm.
    • Thiết lập phương trình dự báo chiều dài ống xói hiện trường ($L_{x_ht}$) theo thời gian ứng với đường quá trình lũ tần suất $P = 0,2%$ (chu kỳ 500 năm) tại Sơn Tây.
    • Xác định độ trễ thời gian $T_{tr\tilde{e}}$ giữa thời điểm mực nước lũ đạt đỉnh và thời điểm xác suất sự cố xói ngầm đạt giá trị cực đại.
  • Đề xuất giải pháp tăng cường ổn định theo nguyên lý rủi ro:
    • Phân tích vùng nghiên cứu ngập lụt Sen Chiểu (huyện Phúc Thọ) bằng mô hình số độ cao DEM, lập bản đồ ngập lụt ứng với các độ sâu 1m, 3m, 5m.
    • Xây dựng đường cong hàm thiệt hại cho từng loại đất sử dụng theo số liệu hiện trạng năm 2016.
    • So sánh các phương án kỹ thuật và lựa chọn giải pháp tối ưu: bố trí hào thu nước kết hợp tầng lọc hạ lưu đê, giúp giảm thiểu tổng chi phí đầu tư và tổn thất rủi ro.

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

Đóng góp mới về mặt khoa học và lý luận

  1. Thiết lập mô hình xói ngầm theo thời gian có cơ sở thực nghiệm: Xây dựng phương trình toán học mô tả sự phát triển chiều dài ống xói ngầm $L_{x,t}$ dưới nền đê theo thời gian thực của một trận lũ, tích hợp hệ số hiệu chỉnh mô hình $\alpha$ từ phân tích ngược sự cố lịch sử (vỡ đê Vân Cốc 1986).
  2. Xác định quy luật trễ của áp lực thấm và xói ngầm: Làm rõ hiện tượng suy giảm chiều dài đường thấm hiệu dụng và độ trễ pha giữa đỉnh lũ ngoài sông với thời điểm áp lực thấm và xác suất sự cố xói ngầm đạt cực đại trong thân và nền đê.
  3. Phát triển phương pháp đánh giá độ tin cậy hệ thống đê sông hoàn chỉnh: Tích hợp đồng thời các đoạn đê và hệ thống cống dưới đê vào mô hình cây sự cố; kết hợp phương pháp tính toán độ tin cậy Cấp II, Cấp III và hiệu ứng độ dài đê trong điều kiện biên thủy lực chịu tác động của biến đổi khí hậu (kịch bản RCP4.5).
  4. Áp dụng thành công nguyên lý rủi ro trong thiết kế tối ưu công trình phòng lũ: Kết hợp bản đồ ngập lụt số, đường cong thiệt hại kinh tế của các loại hình sử dụng đất và chi phí xây dựng công trình để tìm ra giải pháp kỹ thuật có tổng chi phí nhỏ nhất.

Kiến nghị và giải pháp đề xuất

  • Giải pháp công trình: Đề xuất bố trí hệ thống hào thu nước giảm áp kết hợp tầng lọc ngược ở hạ lưu đoạn đê Sen Chiểu (Phúc Thọ) nhằm tiêu thoát nước thấm an toàn, khống chế quá trình phát triển ống xói ngầm mà không gây xáo trộn môi trường xung quanh.
  • Giải pháp quản lý và vận hành: Kiến nghị các cơ quan quản lý đê điều cần đặc biệt chú trọng công tác kiểm tra, tuần tra canh gác đê trong giai đoạn lũ rút (sau đỉnh lũ từ 1,5 đến 2 ngày) do nguy cơ phát sinh đùn sủi và trượt mái dốc trong thời gian này là cao nhất.

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

  • Hạn chế: Luận án tập trung đánh giá an toàn qua các cơ chế sự cố dưới nền đê (biến dạng thấm, xói ngầm, đẩy bục) và cơ chế tràn đỉnh, chưa xét đến các hư hỏng kết cấu bề mặt thân đê do sóng gió va đập, nứt nẻ do co ngót nhiệt ẩm cục bộ hoặc tác động của động đất.
  • Hướng nghiên cứu tiếp: Mở rộng nghiên cứu áp dụng phương pháp lý thuyết độ tin cậy và phân tích rủi ro cho toàn bộ các tuyến đê còn lại trên lưu vực sông Hồng - Thái Bình; nghiên cứu mô hình hóa 3D không gian tương tác giữa cống ngầm và thân đê dưới tác động của dòng thấm biến đổi phi ổn định.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu cho các nhóm đối tượng:

  • Cơ quan quản lý nhà nước và đơn vị khai thác: Cục Quản lý Đê điều và Phòng chống Thiên tai (Bộ Nông nghiệp và PTNT), Chi cục Đê điều và Phòng chống lụt bão Hà Nội tham khảo phương pháp phân đoạn rủi ro và xác định vị trí xung yếu để lập kế hoạch tu bổ, hộ đê hàng năm.
  • Đơn vị tư vấn thiết kế công trình thủy lợi: Ứng dụng quy trình tính toán độ tin cậy (FORM, Monte Carlo) và mô hình xói ngầm theo thời gian để tối ưu hóa giải pháp xử lý thấm nền đê (hào thu nước, giếng giảm áp, cừ chống thấm).
  • Giảng viên, nghiên cứu sinh và học viên cao học: Các cơ sở đào tạo chuyên ngành Địa kỹ thuật công trình, Xây dựng Thủy lợi - Thủy điện (Trường Đại học Thủy Lợi, Đại học Mỏ - Địa chất, Đại học Bách khoa) sử dụng làm tài liệu nghiên cứu về cơ chế biến dạng thấm và phương pháp phân tích độ tin cậy công trình đất.

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

1. Mục tiêu và phạm vi nghiên cứu chính của luận án là gì?

Luận án nghiên cứu đánh giá mức độ an toàn của hệ thống đê Hữu sông Hồng đoạn qua Hà Nội (từ Sơn Tây đến Phú Xuyên, Km30+550 ÷ Km117+900) trong điều kiện biến đổi khí hậu (kịch bản RCP4.5). Trọng tâm nghiên cứu là xác định điều kiện biên thủy lực, quy luật phát triển xói ngầm dưới nền đê qua mô hình vật lý trong phòng và tính toán xác suất sự cố hệ thống (gồm 17 đoạn đê và 10 cống dưới đê) theo lý thuyết độ tin cậy và phân tích rủi ro.

2. Mô hình thí nghiệm vật lý trong phòng ở Chương 3 đã khảo sát những nội dung gì?

Chương 3 thiết kế hai mô hình vật lý trong phòng (mô hình thấm ngang và mô hình thấm đứng) sử dụng mẫu cát lấy từ hai khu vực đùn sủi thực tế là Sen Chiểu (Phúc Thọ) và Hoàng Mai ở ba trạng thái độ chặt (xốp, chặt vừa, chặt). Thí nghiệm cho phép quan sát các giai đoạn phát sinh ống xói ($t_s, t_1, t_2, t_f$), đo đạc gradient áp lực thấm giới hạn và thiết lập quy luật phát triển chiều dài ống xói theo thời gian ($L_x - T$) dưới tác động của đường dâng nước lũ chu kỳ 500 năm.

3. Đánh giá an toàn đê theo lý thuyết độ tin cậy có điểm gì khác biệt so với phương pháp truyền thống?

Phương pháp truyền thống sử dụng các đại lượng đơn lẻ mang tính tiền định để tính hệ số an toàn $K$ hoặc gradient $[J]$, không phản ánh được mức độ biến động ngẫu nhiên của các yếu tố đầu vào. Phương pháp lý thuyết độ tin cậy coi các thông số địa kỹ thuật (hệ số thấm $K$, khối lượng riêng $\rho$, góc ma sát trong $\varphi$, lực dính $c$) và mực nước lũ $H$ là các biến ngẫu nhiên; tính toán hàm trạng thái giới hạn $Z = R - L$ để định lượng xác suất sự cố ($P_f$) và chỉ số độ tin cậy ($\beta$), đồng thời cho phép tính toán an toàn cho toàn bộ hệ thống nối tiếp xét đến hiệu ứng độ dài đê.

4. Luận án đã đề xuất giải pháp cụ thể nào cho đoạn đê trọng điểm Sen Chiểu (Phúc Thọ)?

Thông qua khảo sát địa vật lý xác định vùng cát xốp, mô phỏng số độ cao ngập lụt và lập đường cong hàm thiệt hại cho từng loại đất sử dụng năm 2016 tại huyện Phúc Thọ, luận án áp dụng nguyên lý rủi ro để so sánh chi phí và đề xuất giải pháp tối ưu là bố trí hệ thống hào thu nước giảm áp kết hợp tầng lọc ngược tại chân đê hạ lưu đoạn Sen Chiểu (Km31+500 ÷ Km33+500).

Kết luận

Luận án của nghiên cứu sinh Đặng Quốc Tuấn đã giải quyết trọn vẹn bài toán đánh giá an toàn hệ thống đê Hữu Hồng đoạn qua Hà Nội trong điều kiện biến đổi khí hậu bằng việc kết hợp giữa thực nghiệm mô hình vật lý, phân tích thống kê hiện trường và lý thuyết độ tin cậy. Nghiên cứu đã làm sáng tỏ cơ chế phát triển của ống xói ngầm theo thời gian thực của trận lũ và chứng minh sự suy giảm của chiều dài đường thấm là nguyên nhân gây ra hiện tượng đùn sủi ở gradient thủy lực thấp. Các kết quả định lượng xác suất sự cố cho 17 đoạn đê và 10 cống dưới đê cùng phương án bố trí hào thu nước tối ưu tại đoạn Sen Chiểu cung cấp luận cứ khoa học quan trọng phục vụ công tác quy hoạch, gia cố đê điều và ứng phó với thiên tai ngập lụt tại thủ đô Hà Nội.