Tổng quan về luận án

Mất ổn định bờ sông tại các vùng đồng bằng châu thổ là một trong những hiểm họa địa chất công trình và tai biến tự nhiên nghiêm trọng nhất, đe dọa trực tiếp đến hạ tầng dân sinh, phát triển kinh tế - xã hội và an toàn sinh thái. Công trình nghiên cứu tiến sĩ chuyên ngành Kỹ thuật Địa chất (Mã số: 62520501) của tác giả Trần Lê Thế Diễn, được thực hiện tại Trường Đại học Bách khoa – Đại học Quốc gia TP. Hồ Chí Minh dưới sự hướng dẫn khoa học của PGS.TS Tạ Đức Thịnh và TS. Bùi Trọng Vinh, mang tựa đề: "Nghiên cứu dự báo ổn định bờ sông Hậu đoạn chảy qua tỉnh An Giang và đề xuất giải pháp bảo vệ bờ thích hợp". Công trình đặt trọng tâm giải quyết bài toán cấp bách về sạt lở bờ sông vùng đồng bằng hạ lưu sông Mê Kông trên nền trầm tích Đệ Tứ yếu và phức tạp.

                  ┌─────────────────────────────────────────────────────────┐
                  │                 HỆ THỐNG MẤT ỔN ĐỊNH                    │
                  │              BỜ SÔNG HẬU (TỈNH AN GIANG)                │
                  └────────────────────────────┬────────────────────────────┘
                                               │
             ┌─────────────────────────────────┼────────────────────────────────┐
             │                                 │                                │
             ▼                                 ▼                                ▼
┌─────────────────────────┐       ┌─────────────────────────┐      ┌─────────────────────────┐
│ NGUYÊN NHÂN CỐT LÕI (x) │       │ TÁC ĐỘNG TRỰC TIẾP (y)  │      │  YẾU TỐ THÚC ĐẨY (z)    │
│  Cấu trúc nền địa chất  │       │ Thủy động lực & Hìnhthái│      │    Hoạt động nhân sinh  │
├─────────────────────────┤       ├─────────────────────────┤      ├─────────────────────────┤
│• Kiểu I: Đất yếu đồng   │       │• Lưu lượng lũ & vận tốc │      │• Tải trọng công trình   │
│  nhất / phân lớp xen kẹp│       │  cực đại qua Vàm Nao    │      │  ven sông (QL91, đô thị)│
│• Kiểu II: Đất cù lao mới│       │• Hình thái khúc uốn cong│      │• Khai thác cát lòng sông│
│• 5 Phụ kiểu: IA, IB, IC,│       │  (bờ lõm xã Bình Mỹ)    │      │  bất hợp lý             │
│  IIA, IIB (thấu kính cát│       │• Chế độ phân lưu, mở    │      │• Sóng tàu thuyền tải    │
│  mịn - trung xói ngầm)  │       │  rộng lòng dẫn cục bộ   │      │  trọng lớn (H > 0.3m)   │
└────────────┬────────────┘       └────────────┬────────────┘      └────────────┬────────────┘
             │                                 │                                │
             └─────────────────────────────────┼────────────────────────────────┘
                                               │
                                               ▼
                  ┌─────────────────────────────────────────────────────────┐
                  │                   3 CƠ CHẾ PHÁ HỦY                      │
                  ├─────────────────────────────────────────────────────────┤
                  │ 1. Trượt lở sâu (Mái dốc đất yếu + xói thấu kính cát)   │
                  │ 2. Xói lở hàm ếch (Dòng bao & xâm thực chân bờ sông)    │
                  │ 3. Xói mòn bề mặt (Ứng suất sóng tàu thuyền > tau_c)    │
                  └────────────────────────────┬────────────────────────────┘
                                               │
                                               ▼
                  ┌─────────────────────────────────────────────────────────┐
                  │           PHÂN VÙNG DỰ BÁO & GIẢI PHÁP CHỈNH TRỊ        │
                  ├─────────────────────────────────────────────────────────┤
                  │• Vùng Mất ổn định: Châu Đốc, Bình Mỹ, đầu cù lao BT-MHH │
                  │• Vùng Nguy cơ: Khu vực chuyển tiếp, mở rộng lòng dẫn    │
                  │• Vùng Ổn định: An Châu, Bình Long                       │
                  │=> Giải pháp: Cọc xi măng đất, Cừ Larsen, Kè rọ đá, Mỏhàn│
                  └─────────────────────────────────────────────────────────┘

Theo số liệu điều tra thực nghiệm của luận án, "Dọc bờ sông Hậu có 21 đoạn trượt lở, xói lở với tổng chiều dài là 61.720m chiếm gần 66,3% độ dài sông Hậu" chảy qua tỉnh An Giang. Sự xuất hiện lặp đi lặp lại của các đợt biến dạng sập đổ nghiêm trọng, điển hình như biến cố trượt vỡ đứt gãy đê kè Quốc lộ 91 (xã Bình Mỹ, huyện Châu Phú) qua các năm 2010, 2019, 2020 và sạt lở quy mô lớn tại phường Bình Đức, Bình Khánh (TP. Long Xuyên) năm 2012, đặt ra yêu cầu phải làm sáng tỏ cơ chế nội tại thay vì chỉ nhìn nhận từ khía cạnh thủy lực bề mặt đơn thuần.

Khoảng trống học thuật cốt lõi (research gap) mà nghiên cứu định vị là sự thiếu vắng một phương pháp tiếp cận liên ngành địa chất công trình – thủy động lực học định lượng. Các công trình trước đây tại đồng bằng sông Cửu Long chủ yếu tập trung vào quan trắc hình thái sông hoặc mô phỏng thủy văn vĩ mô, xem cấu trúc nền đất là môi trường đồng nhất tĩnh hoặc là nguyên nhân thứ yếu. Luận án đặt ra và giải quyết 3 câu hỏi nghiên cứu tường minh:

  1. Bản chất tương tác cơ học giữa các phân vị địa tầng đất yếu Đệ Tứ chứa thấu kính cát với quá trình mất ổn định bờ sông là gì?
  2. Cơ chế thủy động lực học dòng chảy kết hợp với sóng nội sinh/nhân sinh tác động như thế nào lên từng kiểu cấu trúc nền đặc trưng?
  3. Mô hình dự báo phân vùng nào tích hợp chính xác các trường dữ liệu địa chất, thủy lực và nhân sinh để làm tiền đề cho thiết kế giải pháp công trình thích ứng?

Phạm vi không gian nghiên cứu bao quát tuyến sông Hậu dài gần 100 km qua tỉnh An Giang, với biên khảo sát địa chất công trình mở rộng 200 m theo chiều ngang tính từ mép bờ và chiều sâu nghiên cứu cấu trúc địa tầng đạt 60 m. Đóng góp đột phá của công trình được định lượng qua việc xây dựng thành công bộ cơ sở dữ liệu địa chất tích hợp gồm 200 lỗ khoan thu thập, 56 lỗ khoan khảo sát chuyên đề và 07 lỗ khoan kiểm tra sâu, giải mã nguồn gốc phá hoại bờ sông và thiết lập bản đồ phân vùng dự báo tỷ lệ lớn phục vụ công tác chỉnh trị bền vững.

Literature Review và Positioning

Các nghiên cứu quốc tế về ổn định bờ sông và diễn biến lòng dẫn đã hình thành nhiều trường phái tiếp cận. Trường phái cơ học xói lở bờ tập trung vào tương tác ứng suất cắt đáy dòng chảy và động lực sóng. Điển hình, Sahameddin Mahmoudi Kurdistani et al. (2019) đã phát triển mô hình tính toán thông số bảo vệ bờ chịu tác động đồng thời của dòng chảy và sóng tàu thuyền, chứng minh chiều cao sóng cực đại ($H_{\max}$), chu kỳ sóng và độ xiên là những tham số điều khiển xói mòn hạt mịn. Trước đó, Larissa Laderoute et al. (2013) tại sông Shuswap (Canada) và McConchie et al. (2003) tại sông Waikato (New Zealand) chỉ ra rằng năng lượng sóng do tàu thuyền tạo ra có thể lớn gấp 2 đến hơn 100 lần sóng do gió, trở thành tác nhân xung kích gây phá hủy cấu trúc bãi bồi ven bờ.

Ở khía cạnh tương tác thủy lực khúc uốn và áp lực nước lỗ rỗng, Massimo Rinaldi et al. (2008) khi mô phỏng sông Cecina (Ý) đã phát hiện cơ chế phá hoại bờ xảy ra khốc liệt nhất không phải tại đỉnh lũ mà ở giai đoạn sau đỉnh lũ (nước rút nhanh), khi sự kết hợp giữa áp lực thủy tĩnh giảm đột ngột và áp lực nước lỗ rỗng dư kích hoạt trượt cung tròn. Về mặt thạch học, Kukulak et al. (2016) khi nghiên cứu vùng Podhale (Ba Lan) khẳng định tính dị hướng và sự không đồng nhất của trầm tích lòng dẫn đóng vai trò kiểm soát các mặt trượt tiềm năng.

Tại Việt Nam, các nghiên cứu trên hệ thống sông Cửu Long đã ghi nhận nhiều dấu mốc quan trọng:

  • Nhóm nghiên cứu thủy văn và hình thái sông: Lê Ngọc Bích et al. (1995-1998), Lê Mạnh Hùng et al. (1999-2004), Hoàng Văn Huân et al. (1998-2003) đã thiết lập các mối quan hệ hình thái lòng dẫn, đánh giá biến thiên mặt cắt ngang và dự báo diễn biến bồi xói vĩ mô. Đinh Công Sản et al. (2009) tập trung giải bài toán phân lưu qua sông Vàm Nao từ sông Tiền sang sông Hậu.
  • Nhóm nghiên cứu viễn thám và địa mạo: Hà Quang Hải et al. (2009-2011) và Lâm Đạo Nguyên ứng dụng GIS – viễn thám đo vẽ địa hình đáy sông, làm rõ xu thế dịch chuyển của các bãi bồi, cù lao. Nguyễn Nghĩa Hùng et al. (2013-2016) nghiên cứu quá trình biến động hình thái cù lao và định hướng giải pháp kè bán kiên cố tại Cần Thơ.
  • Nhóm nghiên cứu địa chất công trình hạ lưu: Phạm Văn Tỵ, Đậu Văn Ngọ, Nguyễn Bá Hoằng (1995-2001) đã phân tích môi trường địa chất tác động đến bờ sông Sài Gòn – Đồng Nai. Bùi Trọng Vinh et al. (2009-2020) mở rộng ứng dụng mô hình số trị (MIKE, Flow 3D, GeoSlope/W) tại cửa sông Thị Vải và Tiền Giang.
                                  TIẾN TRÌNH HỌC THUẬT QUỐC TẾ & VIỆT NAM
  ┌───────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │ QUỐC TẾ:                                                                                              │
  │ • Massimo Rinaldi et al. (2008): Cơ chế xói lở sau đỉnh lũ do áp lực nước lỗ rỗng dư (Sông Cecina, Ý) │
  │ • Kurdistani et al. (2019), Laderoute et al. (2013): Thủy động lực học sóng tàu thuyền & xói bờ hạt mịn│
  │ • Kukulak et al. (2016): Tính không đồng nhất thạch học kiểm soát mặt trượt cung tròn (Ba Lan)        │
  └──────────────────────────────────────────────────┬────────────────────────────────────────────────────┘
                                                     │
  ┌──────────────────────────────────────────────────┴────────────────────────────────────────────────────┐
  │ VIỆT NAM (ĐỒNG BẰNG SÔNG CỬU LONG):                                                                   │
  │ • Lê Ngọc Bích, Lê Mạnh Hùng, Hoàng Văn Huân (1995-2004): Thủy lực hình thái & biến động lòng dẫn     │
  │ • Đinh Công Sản (2009): Động lực học phân lưu dòng chảy qua ngã ba sông Vàm Nao                      │
  │ • Hà Quang Hải, Lâm Đạo Nguyên (2009-2011): GIS & Viễn thám biến động đường bờ và bãi bồi             │
  │ • Bùi Trọng Vinh et al. (2009-2020): Mô phỏng số trị thủy lực - ổn định mái dốc ven biển & cửa sông   │
  └──────────────────────────────────────────────────┬────────────────────────────────────────────────────┘
                                                     │
                                                     ▼ [RESEARCH GAP]
  ┌───────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │ ĐỘT PHÁ CỦA LUẬN ÁN TRẦN LÊ THẾ DIỄN (2023):                                                          │
  │ Tích hợp ĐỊA CHẤT CÔNG TRÌNH NỀN ĐẤT YẾU ĐỆ TỨ (Phân lập 2 kiểu, 5 phụ kiểu chứa thấu kính cát)     │
  │ với THỦY ĐỘNG LỰC DÒNG CHẢY 2D (MIKE 21) & CƠ HỌC ĐẤT MÔ PHỎNG SỐ (GeoSlope/W, PLAXIS)                │
  └───────────────────────────────────────────────────────────────────────────────────────────────────────┘

Điểm định vị học thuật vượt trội của luận án Trần Lê Thế Diễn so với các công trình trên là sự kết hợp định lượng cấu trúc nền địa chất công trình trầm tích Đệ Tứ với động lực học dòng chảy và áp lực tải trọng hiện đại. Trong khi các nghiên cứu quốc tế thường giả định tính chất cơ lý nền tương đối đồng nhất hoặc chỉ nghiên cứu xói lở bề mặt đơn lẻ, luận án đã mổ xẻ cấu trúc địa chất dạng phân lớp chứa "thấu kính cát", chứng minh hiện tượng xói ngầm thủy lực (piping/internal erosion) trong các thấu kính này là ngòi nổ kích hoạt sự trượt lở sụp đổ tổng thể của toàn bộ khối trượt.

Đó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 quan trọng vào lý thuyết cơ học đất ứng dụng và địa kỹ thuật công trình sông thông qua việc mở rộng định luật sức kháng cắt Coulomb và lý thuyết cân bằng giới hạn Mohr ($\tau_f = c' + (\sigma - u)\tan\phi'$) trong điều kiện nền địa chất có sự hiện diện của các thấu kính cát không liên tục. Luận án chỉ ra rằng sự suy giảm hệ số ổn định $F_s$ diễn ra qua hai pha: pha suy giảm ứng suất hữu hiệu $\sigma'$ do dao động áp lực nước lỗ rỗng $u$, tiếp nối bởi pha mất liên kết hình thái do xói ngầm thấu kính cát mịn dưới gradien thủy lực dòng ngầm thấm thoát ra mép sông.

Hệ thống luận điểm lý thuyết được chuẩn hóa thông qua 3 mệnh đề cốt lõi:

  1. Mệnh đề 1: Độ ổn định của bờ sông là hàm đa biến đồng thời: $F_s = f(x, y, z)$, trong đó biến số nội sinh $x$ (kiểu cấu trúc nền địa chất) đóng vai trò điều kiện cần và quyết định hình thức phá hủy; biến số động lực $y$ (thủy động lực dòng chảy, năng lượng sóng, hình thái khúc sông cong) và biến số nhân sinh $z$ (gia tải công trình, khai thác cát đáy sông) là điều kiện đủ kích hoạt và đẩy nhanh tốc độ phá hủy.
  2. Mệnh đề 2: Cơ chế mất ổn định bờ sông Hậu là sự chuyển tiếp giữa 3 trạng thái động: trượt lở mái dốc cung tròn sâu (landslide), xói lở xâm thực chân bờ tạo hàm ếch (bank collapse/scour), và xói mòn bóc tách hạt bề mặt (surface erosion).
  3. Mệnh đề 3: Sự tồn tại của thấu kính cát hạt mịn đến hạt trung trong các phân vị bùn sét/sét pha tạo nên đới tập trung ứng suất cắt và đường thoát nước ưu tiên, làm triệt tiêu sức kháng cắt cục bộ và hạ thấp hệ số an toàn tính toán $F_s$ xuống dưới ngưỡng cân bằng giới hạn ($F_s < 1.0$) khi chịu tải trọng động hoặc tải trọng tĩnh bề mặt.
                              KHUNG MÔ HÌNH PHÂN TÍCH TỔNG HỢP: Fs = f(x, y, z)
  ┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │                               BIẾN SỐ NỘI SINH: CẤU TRÚC NỀN ĐỊA CHẤT (x)                              │
  │ • Kiểu I & 5 Phụ kiểu (IA, IB, IC, IIA, IIB)                                                          │
  │ • Sự hiện diện thấu kính cát (TK) hạt mịn-trung trong bùn sét pha (Q2) -> Xói ngầm & tạo rỗng         │
  └──────────────────────────────────────────────────┬─────────────────────────────────────────────────────┘
                                                     │
                                                     ▼
  ┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │                          BIẾN SỐ ĐỘNG LỰC: THỦY LỰC, HÌNH THÁI, ĐỘNG LỰC SÓNG (y)                      │
  │ • Vmax dòng chủ mùa lũ & phân lưu Vàm Nao -> Xâm thực đáy, mở rộng lòng dẫn (hệ số B/h)               │
  │ • Ứng suất cắt sóng tàu thuyền tau_b > tau_c (sức kháng cắt bề mặt từ Jet test) -> Xói mòn bóc vỏ      │
  └──────────────────────────────────────────────────┬─────────────────────────────────────────────────────┘
                                                     │
                                                     ▼
  ┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │                           BIẾN SỐ NHÂN SINH: TẢI TRỌNG & TÁC ĐỘNG NGOẠI TẠI (z)                        │
  │ • Gia tải đô thị, tuyến giao thông sát mép bờ (QL91, P. Bình Đức, P. Bình Khánh) -> Tăng momen trượt   │
  │ • Khai thác cát hạ thấp cao độ đáy sông cục bộ -> Tăng góc nghiêng mái dốc tự nhiên                    │
  └──────────────────────────────────────────────────┬─────────────────────────────────────────────────────┘
                                                     │
                                                     ▼
  ┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │ KẾT QUẢ ĐÁNH GIÁ ĐỊNH LƯỢNG HỆ SỐ AN TOÀN: Fs = M_chống / M_gây (GeoSlope/W, PLAXIS, MIKE 21)          │
  │ • Không thấu kính cát: Fs > 1.25 (Ổn định)                                                             │
  │ • Có thấu kính cát + Xói ngầm + Gia tải: Fs = 0.65 - 0.88 (Mất ổn định trượt sâu tức thời)            │
  └────────────────────────────────────────────────────────────────────────────────────────────────────────┘

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: Lý thuyết địa chất Đệ Tứ và Địa mạo công trình; Lý thuyết cân bằng giới hạn cơ học đất (Limit Equilibrium Method - Bishop, Fellenius); Lý thuyết thủy lực dòng sông hở (phương trình Navier-Stokes trung bình Reynolds dạng 2D shallow water); và Phương pháp xói mòn phản lực tia nước (Jet Index Method - Hanson, 2002).

Công trình đã thiết lập một hệ thống phân loại cấu trúc nền đất yếu bờ sông Hậu gồm 2 kiểu và 5 phụ kiểu:

  • Kiểu I (Cấu trúc nền đất yếu bờ sông tự nhiên và đô thị):
    • Phụ kiểu IA: Nền đất sét dẻo mềm đến dẻo chảy chiều dày lớn, bên dưới là cát hoặc sét pha chứa thấu kính cát ngậm nước (phân bố tại TP. Châu Đốc, TP. Long Xuyên).
    • Phụ kiểu IB: Nền đất gồm lớp bùn sét/sét pha hữu cơ xen kẹp các thấu kính cát hạt mịn đến hạt trung ở độ sâu từ -5 m đến -25 m (phân bố tại xã Bình Mỹ, P. Bình Đức, P. Bình Khánh).
    • Phụ kiểu IC: Nền chuyển tiếp với chiều dày lớp đất yếu mỏng dần, sức chịu tải tăng dần.
  • Kiểu II (Cấu trúc nền đất các cù lao và bãi bồi mới thành tạo):
    • Phụ kiểu IIA: Trầm tích sông (a$Q_2$) và đầm lầy (ab$Q_2$), bề mặt là lớp sét pha bụi có tính trương nở - co ngót mạnh khi thay đổi độ ẩm (đầu cù lao Bình Thạnh).
    • Phụ kiểu IIB: Nền cát pha, sét pha phân lớp mỏng, độ liên kết kém, dễ bị dòng bao xoáy làm xói chân và bóc vỏ bề mặt do sóng (cù lao Mỹ Hòa Hưng).

Biên giới hạn ứng dụng (boundary conditions) của khung lý thuyết này áp dụng chuẩn xác cho các hệ thống sông thuộc châu thổ bồi tích Đệ Tứ phẳng, chịu chế độ bán nhật triều không đều và có sự tương tác thủy lực mạnh từ các trục phân lưu (như sông Vàm Nao).

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

Thiết kế nghiên cứu

Luận án vận dụng triết lý thực chứng khoa học kết hợp quan điểm duy vật biện chứng địa chất công trình, xem đất đá là hệ thống động biến đổi liên tục dưới tác động của các trường lực nội sinh và ngoại sinh. Thiết kế nghiên cứu đa tầng bậc (multi-level research design) kết hợp giữa phương pháp thực địa thực nghiệm, mô phỏng số học (numerical simulation) và phân tích không gian GIS.

                                  QUY TRÌNH THỰC HIỆN NGHIÊN CỨU TÍCH HỢP
  ┌───────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │ 1. DỮ LIỆU ĐỊA CHẤT CÔNG TRÌNH & THỰC NGHIỆM                                                          │
  │ • Thu thập 200 lỗ khoan ĐCCT lịch sử + 56 lỗ khoan chuyên khảo + 07 lỗ khoan kiểm tra sâu 60m         │
  │ • Thí nghiệm trong phòng: Cơ lý đất tiêu chuẩn (TCVN), nén lún, cắt trực tiếp, cắt ba trục (CU, UU)   │
  │ • Thí nghiệm chuyên biệt: Độ trương nở - co ngót thể tích, Thí nghiệm xói Jet Index Test (Hanson 2002)│
  └──────────────────────────────────────────────────┬────────────────────────────────────────────────────┘
                                                     │
                                                     ▼
  ┌───────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │ 2. DỮ LIỆU THỦY VĂN, ĐỘNG LỰC DÒNG CHẢY & SÓNG                                                        │
  │ • Quan trắc thủy lực mùa kiệt (tháng 04/2010) và mùa lũ (tháng 10/2010) tại Châu Đốc, Vàm Nao, B.Khánh │
  │ • Thiết lập hệ thống trạm đo sóng thực tế: Tần suất & năng lượng sóng tàu thuyền (ST), sóng gió (SG)   │
  │ • Giải đoán ảnh viễn thám Landsat 5, 8, 9 (1993, 2011, 2016, 2021) xác định biến động đường bờ        │
  └──────────────────────────────────────────────────┬────────────────────────────────────────────────────┘
                                                     │
                                                     ▼
  ┌───────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │ 3. MÔ PHỎNG SỐ HỌC (NUMERICAL MODELING) & HIỆU CHỈNH                                                  │
  │ • MIKE 21 (HD, ST, MT): Lưới tính phi cấu trúc, mô phỏng trường vận tốc Vmax, hướng dòng & bồi xói đáy│
  │   -> Kiểm định Nash-Sutcliffe (NASH) mực nước và lưu lượng bùn cát đạt độ chính xác cao                │
  │ • GeoSlope/W & PLAXIS 2D/3D: Phân tích ổn định mái dốc cung trượt, biến dạng chuyển vị & áp lực kẽ rỗng│
  │   xét đến kịch bản xói ngầm thấu kính cát và gia tải bề mặt                                           │
  └──────────────────────────────────────────────────┬────────────────────────────────────────────────────┘
                                                     │
                                                     ▼
  ┌───────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │ 4. TỔNG HỢP & DỰ BÁO KHÔNG GIAN (GIS OVERLAY)                                                         │
  │ • Chồng chập các lớp thông tin: Cấu trúc nền, Thủy động lực học, Hình thái sông, Hoạt động nhân sinh  │
  │ => Bản đồ phân vùng dự báo ổn định bờ sông Hậu tỷ lệ lớn (Vùng mất ổn định, Vùng nguy cơ, Vùng ổn định)│
  └───────────────────────────────────────────────────────────────────────────────────────────────────────┘

Mẫu nghiên cứu địa chất bao gồm 200 mặt cắt lỗ khoan thu thập từ các dự án cầu đường, kè sông trong vùng, 56 lỗ khoan khảo sát bổ sung và 07 lỗ khoan đường kính lớn chuyên biệt do chính tác giả thực hiện tại các điểm nóng sạt lở: xã Bình Mỹ (huyện Châu Phú), phường Bình Đức, phường Bình Khánh (TP. Long Xuyên), đầu cù lao Bình Thạnh và cù lao Mỹ Hòa Hưng.

Quy trình nghiên cứu rigorous

Quy trình thu thập và phân tích dữ liệu tuân thủ nghiêm ngặt các tiêu chuẩn kỹ thuật quốc gia (TCVN) và quốc tế (ASTM, IAEG):

  • Khảo sát địa chất và thí nghiệm cơ lý: Mẫu đất nguyên dạng được lấy bằng ống mẫu thành mỏng nòng đôi để hạn chế tối đa phá hủy kết cấu. Thí nghiệm xác định các chỉ tiêu dung trọng ($\gamma$), độ ẩm ($W$), góc ma sát trong ($\phi$), lực dính kết ($c$), hệ số nén lún ($a$), và hệ số thấm ($k$).
  • Thí nghiệm chuyên biệt phản lực tia nước (Jet Test Apparatus): Ứng dụng theo quy trình chuẩn của Hanson (2002) nhằm xác định trực tiếp sức kháng cắt tới hạn bề mặt ($\tau_c$) và hệ số kháng xói mòn ($K_d$) trên các mẫu đất bùn sét bãi bồi cù lao.
  • Thiết lập trạm quan trắc động lực sóng: Hệ thống cảm biến áp suất đáy và cọc đo quang học ghi nhận sóng tàu thuyền (ST) và sóng do gió (SG) tại 4 trạm quan trắc xung quanh cù lao Mỹ Hòa Hưng và Bình Thạnh, đồng bộ hóa với nhật trình lưu thông phương tiện vận tải thủy.
  • Kiểm soát độ tin cậy tam giác hóa (Triangulation): Kết quả khảo sát thạch học từ lỗ khoan được kiểm chứng chéo với tài liệu đo sâu địa chấn nông phân giải cao và số liệu quét địa hình đáy sông bằng máy đo sâu hồi âm đa chùm tia (Multibeam Sonar).
          MÔ HÌNH THÍ NGHIỆM JET TEST THEO HANSON (2002)
   ┌────────────────────────────────────────────────────────────┐
   │ Ống cấp nước áp lực không đổi (H constant)                 │
   └─────────────────────────────┬──────────────────────────────┘
                                 │
                                 ▼
   ┌────────────────────────────────────────────────────────────┐
   │ Vòi phun phản lực nước (Jet Nozzle d = 6.35 mm)           │
   │ Vận tốc dòng tia tác động vuông góc lên bề mặt mẫu        │
   └─────────────────────────────┬──────────────────────────────┘
                                 │
                                 ▼
   ┌────────────────────────────────────────────────────────────┐
   │ Mẫu đất nguyên sinh/chế bị trong buồng ngập nước          │
   │ Quan trắc độ sâu vết xói rãnh theo thời gian t (phút)      │
   ├────────────────────────────────────────────────────────────┤
   │ => Xác định: Ứng suất cắt tới hạn tau_c (Pa)               │
   │              Hệ số tốc độ xói mòn Kd (cm3/N.s)             │
   └────────────────────────────────────────────────────────────┘

Data và phân tích

Mô hình hóa thủy lực 2D dòng chảy sông Hậu được xây dựng trên phần mềm MIKE 21 (DHI) sử dụng lưới tính phi cấu trúc (Flexible Mesh) với hàng chục nghìn phần tử tam giác và tứ giác, cho phép tinh mịn hóa tại các khúc uốn và ngã ba sông:

  • Điều kiện biên: Biên lưu lượng thượng lưu tại Châu Đốc và sông Tiền; biên mực nước hạ lưu tại các trạm kiểm soát triều phụ cận Vàm Cống; biên nhập lưu bổ sung đặc biệt tại ngã ba sông Vàm Nao.
  • Hệ số nhám Manning ($n$): Được hiệu chỉnh theo từng phân vùng đáy từ 0.022 đến 0.035 $\text{s/m}^{1/3}$.
  • Kiểm định mô hình (Model Calibration & Validation): Sử dụng hệ số hiệu quả Nash-Sutcliffe ($E_{\text{NASH}}$) đối chuẩn với chuỗi số liệu thực đo mùa kiệt (tháng 04/2010) và mùa lũ (tháng 10/2010) tại các trạm Châu Đốc, Vàm Nao, Bình Khánh. Kết quả $E_{\text{NASH}}$ đều đạt giá trị từ 0.85 đến 0.94 (đạt cấp độ rất tốt trong tiêu chuẩn thủy lực quốc tế).

Phân tích ổn định mái dốc được tiến hành trên phần mềm GeoSlope/W (module SLOPE/W và SEEP/W) kết hợp mô hình phần tử hữu hạn PLAXIS 2D/3D:

  • Áp dụng các phương pháp phân mảnh lát cắt cổ điển và nâng cao: Bishop đơn giản hóa, Janbu hiệu chỉnh, Morgenstern-Price đáp ứng đồng thời cân bằng lực và cân bằng mômen.
  • Thiết lập bài toán kiểm tra độ nhạy (robustness checks) với các trường hợp: (1) Trạng thái tự nhiên; (2) Mực nước rút nhanh mùa lũ kết hợp áp lực kẽ rỗng dư; (3) Mái dốc bị xói lở chân tạo hàm ếch cục bộ; (4) Khuyết tán rỗng do xói ngầm thấu kính cát; (5) Gia tải công trình xây dựng và đê bao Quốc lộ 91.
       BÀI TOÁN PHÂN TÍCH PHÂN MẢNH LÁT CẮT (MORGENSTERN-PRICE / BISHOP)
                      TRONG GEOSLOPE/W VÀ PLAXIS
                       Tải trọng công trình q (kN/m2)
                       ↓↓↓↓↓↓↓↓↓↓↓↓↓↓↓↓↓↓↓↓↓↓↓↓↓↓
       Mặt đất tự nhiên ──────────────────────┐
                                               \  Mái dốc tự nhiên
                                                \
         Lớp đất phủ sét pha                     \
       ═══════════════════════════════════════════\
         Lớp bùn sét yếu                           \
       ─────────────────────────────────────────────\   Mực nước sông
         [ THẤU KÍNH CÁT (XÓI NGẦM) ]                \ ~~~~~~~~~~~~~~
       ═══════════════════════════════════════════════\──────────────
         Lớp sét dẻo cứng chịu lực                     \ Đáy sông xói
                                                        \__________
         - - - - - - - - - - - - - - - - - - - - - - - /
                 Cung trượt tiềm năng (Mặt trượt tròn sâu)

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

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

Thứ nhất, giải mã nguyên nhân cốt lõi gây sạt lở lặp lại tại Quốc lộ 91 (xã Bình Mỹ) và TP. Long Xuyên. Kết quả khảo sát lỗ khoan và mô phỏng số đã chỉ ra sự tồn tại của phụ kiểu cấu trúc nền IB: các thấu kính cát hạt mịn đến hạt trung nằm xen kẹp trong lớp bùn sét yếu ở cao trình từ -8 m đến -18 m. Dưới tác động của chênh lệch cột nước áp lực và gradien thấm khi lũ rút, hiện tượng xói ngầm thủy lực (piping) xảy ra bên trong thấu kính cát, rửa trôi hạt mịn tạo thành các khoang rỗng ngầm. Khi xuất hiện tải trọng động của phương tiện giao thông và tải trọng tĩnh của công trình ven bờ, hệ số an toàn mái dốc sụt giảm đột ngột từ $F_s = 1.18$ xuống $F_s = 0.68 - 0.74$, kích hoạt trượt lở quy mô lớn dạng cung trượt sâu phá hủy toàn bộ thân đường.

Thứ hai, phát hiện cơ chế xói mòn bóc vỏ bề mặt do sóng tàu thuyền tại các cù lao. Kết quả đo đạc từ thí nghiệm Jet test cho thấy sức kháng cắt tới hạn bề mặt của lớp đất sét pha bụi cù lao Mỹ Hòa Hưng và Bình Thạnh dao động trong khoảng $\tau_c = 0.42 - 0.78\text{ Pa}$. Trong khi đó, các xung sóng tàu thuyền có chiều cao $H > 0.3\text{ m}$ tạo ra ứng suất cắt đáy tức thời $\tau_b$ lên tới $1.25 - 2.10\text{ Pa}$, vượt xa giá trị $\tau_c$. Điều này giải thích tại sao tốc độ xâm thực đường bờ tại mũi đầu cù lao Bình Thạnh đạt mức kỷ lục 22.0 m/năm (giai đoạn 1993 - 2016) và đầu cù lao Mỹ Hòa Hưng đạt 10.0 - 30.0 m/năm.

Thứ ba, định lượng hóa tác động thủy động lực từ sự tiếp nhận lưu lượng sông Vàm Nao. Sông Hậu sau khi tiếp nhận tới 40 - 50% tổng lưu lượng từ sông Tiền chuyển qua sông Vàm Nao đã tạo nên sự gia tăng đột biến về vận tốc dòng chảy ($V_{\max} > 1.8 - 2.2\text{ m/s}$ trong mùa lũ) và tạo ra các xoáy cuộn trục đứng tại các khu vực lòng dẫn hẹp hoặc khúc sông uốn cong (như đoạn qua xã Bình Mỹ và bờ phải phường Bình Đức). Dòng chủ áp sát bờ gây xói chân bờ dốc đứng với góc dốc tự nhiên bị kéo dốc từ $25^\circ$ lên tới $55^\circ - 70^\circ$, phá vỡ thế cân bằng cơ học tự nhiên của mái dốc.

Thứ tư, xác lập bản đồ phân vùng dự báo ổn định bờ sông Hậu bằng kỹ thuật chồng chập đa lớp. Nghiên cứu đã tích hợp 4 trường thông tin địa không gian: (1) Bản đồ cấu trúc nền đất yếu; (2) Bản đồ trường phân bố vận tốc thủy động lực; (3) Bản đồ chỉ số hình thái khúc sông uốn và tỷ lệ mở rộng lòng dẫn $B/h$; (4) Bản đồ mật độ tải trọng kinh tế - xây dựng và khai thác cát. Phân chia chuẩn xác 3 vùng:

  1. Vùng mất ổn định: TP. Châu Đốc, xã Bình Mỹ, mũi đầu cù lao Bình Thạnh, dải cù lao Mỹ Hòa Hưng.
  2. Vùng nguy cơ mất ổn định: Đoạn bờ chuyển tiếp xã Bình Khánh, thị trấn Cái Dầu.
  3. Vùng ổn định: Đoạn bờ xã Bình Long, thị trấn An Châu (nơi cấu trúc nền đồng nhất, không chứa thấu kính cát và lòng sông thoải).
   BẢN ĐỒ PHÂN VÙNG DỰ BÁO ỔN ĐỊNH BỜ SÔNG HẬU (TỈNH AN GIANG) - TỔNG HỢP
  ┌─────────────────────────────────────────────────────────────────────────┐
  │ 1. VÙNG MẤT ỔN ĐỊNH CAO (ĐỎ):                                           │
  │    • TP. Châu Đốc (Phụ kiểu IA + Dòng chảy xiết sát bờ)                 │
  │    • Xã Bình Mỹ - QL91 (Phụ kiểu IB + Thấu kính cát + Khúc uốn + Gia tải)│
  │    • Đầu Cù lao Bình Thạnh & Mỹ Hòa Hưng (Phụ kiểu IIA/IIB + Sóng tàu)  │
  │    • P. Bình Đức, P. Bình Khánh (Phụ kiểu IB + Gia tải đô thị)          │
  ├─────────────────────────────────────────────────────────────────────────┤
  │ 2. VÙNG NGUY CƠ MẤT ỔN ĐỊNH (VÀNG):                                     │
  │    • Đoạn chuyển tiếp từ kênh Cây Dương đến Phà Năng Gù                 │
  │    • Đuôi cù lao Mỹ Hòa Hưng & dải bờ giáp ranh Cần Thơ                 │
  ├─────────────────────────────────────────────────────────────────────────┤
  │ 3. VÙNG ỔN ĐỊNH (XANH):                                                 │
  │    • Bờ sông qua Thị trấn An Châu (Huyện Châu Thành)                    │
  │    • Đoạn bờ sông qua Xã Bình Long (Huyện Châu Phú)                     │
  └─────────────────────────────────────────────────────────────────────────┘

Implications đa chiều

  • Đóng góp học thuật: Xác lập phương pháp luận chuẩn cho địa kỹ thuật công trình sông đồng bằng, chuyển trọng tâm từ phòng chống thụ động sang dự báo dựa trên nhận diện cấu trúc địa chất 3D.
  • Đổi mới phương pháp luận: Cung cấp quy trình tích hợp giữa mô phỏng thủy văn 2D (MIKE 21), động lực học thấm - trượt dốc (GeoSlope/W, PLAXIS) và kiểm tra sức bền vi mô (Jet test).
  • Ứng dụng thực tiễn và giải pháp kỹ thuật thích ứng:
    1. Khu vực xã Bình Mỹ (Phụ kiểu IB chứa thấu kính cát): Kiến nghị giải pháp cọc xi măng đất (Soil-cement mixing piles) đường kính D800 - D1000 cắm sâu qua lớp bùn sét vào tầng cát chịu lực nhằm đông kết các thấu kính cát ngậm nước, kết hợp xây dựng hệ mỏ hàn chỉnh trị hướng dòng chảy xa bờ.
    2. Khu vực đầu cù lao Bình Thạnh và Mỹ Hòa Hưng (Cơ chế xói mòn sóng và dòng bao): Ứng dụng thảm đá/kè rọ đá (gabions) đa tầng kết hợp vải địa kỹ thuật lọc ngược để tiêu tán năng lượng sóng tàu thuyền và dòng xoáy, bảo vệ chân bờ bãi bồi.
    3. Khu vực đô thị TP. Long Xuyên (Phụ kiểu IA, IB chịu tải trọng công trình): Áp dụng tường cừ Larsen thép kết hợp neo mềm trong đất nhằm giữ ổn định thành đứng bờ kè, giảm thiểu biến dạng ngang. Đồng thời điều chỉnh quy hoạch khai thác cát khơi thông lạch phụ cù lao An Thạnh Trung và cù lao Phó Ba để giải tỏa áp lực lưu lượng cho nhánh chính.
  • Khuyến nghị chính sách: Cung cấp luận cứ khoa học để UBND tỉnh An Giang và Bộ Nông nghiệp & PTNT ban hành hành lang an toàn sạt lở, phân loại cấp phép khai thác khoáng sản cát lòng sông dựa trên độ nhạy cảm địa kỹ thuật.
                  GIẢI PHÁP KỸ THUẬT BẢO VỆ BỜ THEO TỪNG CƠ CHẾ
  ┌─────────────────────────────────────────────────────────────────────────┐
  │ CƠ CHẾ 1: TRƯỢT SÂU NỀN ĐẤT YẾU CHỨA THẤU KÍNH CÁT (BÌNH MỸ, LONG XUYÊN)│
  │  => CỌC XI MĂNG ĐẤT (SOIL CEMENT PILES D800-D1000)                      │
  │     • Đông kết thấu kính cát, triệt tiêu xói ngầm                       │
  │     • Tạo tường chắn composite gia cường sức kháng cắt                  │
  ├─────────────────────────────────────────────────────────────────────────┤
  │ CƠ CHẾ 2: DÒNG CHỦ TẤN CÔNG & MỞ RỘNG LÒNG DẪN (BỜ LÕM, KHÚC CONG)     │
  │  => HỆ THỐNG MỎ HÀN CHỈNH TRỊ (GROYNES) + NẠO VÉT KHƠI THÔNG LẠCH PHỤ   │
  │     • Chuyển hướng dòng lưu tốc lớn ra xa mép bờ dốc                    │
  │     • Khơi thông lạch cù lao An Thạnh Trung / Phó Ba                    │
  ├─────────────────────────────────────────────────────────────────────────┤
  │ CƠ CHẾ 3: XÓI MÒN BỀ MẶT DO SÓNG TÀU THUYỀN (CÙ LAO BÌNH THẠNH, MHH)   │
  │  => KÈ RỌ ĐÁ TIÊU NĂNG (GABIONS) + VẢI ĐỊA KỸ THUẬT LỌC NGƯỢC           │
  │     • Triệt tiêu năng lượng sóng va đập (H > 0.3m)                      │
  │     • Giữ ổn định hạt sét pha bụi dưới ứng suất cắt dao động            │
  └─────────────────────────────────────────────────────────────────────────┘

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, công trình vẫn tồn tại một số giới hạn nghiên cứu:

  • Giới hạn dữ liệu thời gian thực: Các đợt khảo sát động lực học thủy văn tập trung điển hình vào tháng 04 (mùa kiệt) và tháng 10 (mùa lũ) năm 2010. Do đó, mô hình chưa phản ánh liên tục chuỗi biến động thủy văn đa năm dưới tác động của các đập thủy điện thượng nguồn sông Mê Kông và biến đổi khí hậu cực đoan trong giai đoạn 2021 - 2026.
  • Giới hạn không gian phân tích 3D: Các mô hình phân tích ổn định mái dốc chủ yếu thực hiện dưới dạng mặt cắt phẳng 2D (2D limit equilibrium/FEM), chưa mô phỏng hoàn chỉnh hiệu ứng vòm không gian 3D tại các vị trí uốn khúc phức tạp.
  • Giới hạn về đặc tính lưu biến của bùn sét: Mô hình vật liệu trong PLAXIS chủ yếu sử dụng Mohr-Coulomb và Soft Soil Model, chưa tích hợp đầy đủ tính lưu biến (creep behavior) và phá hủy mỏi của đất yếu dưới tác động lặp liên tục của tải trọng sóng tàu thuyền.

Định hướng nghiên cứu mở rộng trong tương lai bao gồm:

  1. Phát triển hệ thống giám sát tự động IoT kết hợp cảm biến áp lực nước lỗ rỗng và độ dịch chuyển sâu (inclinometer) truyền dữ liệu thời gian thực tại các cung trượt xung yếu như xã Bình Mỹ.
  2. Ứng dụng mô hình liên kết thủy lực - địa kỹ thuật 3D (Coupled Hydro-Mechanical 3D Modeling) mô phỏng tương tác động giữa dòng chảy rối 3 chiều và sự phá hủy lũy tiến của bờ sông.
  3. Nghiên cứu ảnh hưởng của suy giảm bùn cát mịn từ thượng nguồn Mê Kông đến sự thay đổi độ nhám đáy và động lực biến đổi lòng dẫn hạ du.
  4. Đánh giá khả năng ứng dụng các giải pháp kỹ thuật sinh thái (Bio-engineering), kết hợp thảm thực vật ngập nước bản địa với kè mềm để giảm giá thành và phục hồi cảnh quan ven sông.

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

Nghiên cứu mang lại giá trị tác động học thuật và thực tiễn to lớn:

  • Tác động học thuật: Đặt nền móng cho hướng nghiên cứu địa kỹ thuật sông chuyên sâu tại Việt Nam. Các phân loại kiểu cấu trúc nền (Type I, Type II) và cơ chế xói ngầm thấu kính cát trở thành tài liệu tham khảo chuẩn mực cho các chương trình đào tạo sau đại học và các đề tài cấp Nhà nước về tai biến đồng bằng. Ước tính công trình có tiềm năng trích dẫn cao trong các tạp chí ISI/Scopus chuyên ngành (như Engineering Geology, Geomorphology, Catena).
  • Tác động kinh tế – xã hội: Cung cấp phương án kỹ thuật tối ưu giúp tiết kiệm hàng trăm tỷ đồng cho ngân sách nhà nước trong việc xử lý sự cố sạt lở Quốc lộ 91, thay thế các giải pháp đắp đá đổ tạm bợ kém hiệu quả bằng hệ cọc xi măng đất và mỏ hàn nắn dòng bài bản.
  • Tác động chính sách: Là cơ sở khoa học để Sở Tài nguyên và Môi trường tỉnh An Giang, Cục Quản lý Đê điều và Phòng chống Thiên tai quy hoạch các cụm tuyến dân cư an toàn, di dời hơn 1.500 hộ dân ra khỏi vùng có nguy cơ trượt lở cao tại TP. Long Xuyên và huyện Châu Phú.
  • Ý nghĩa quốc tế: Phương pháp luận nghiên cứu hoàn toàn có thể chuyển giao và áp dụng hiệu quả cho các quốc gia hạ lưu sông xuyên biên giới như Campuchia (sông Mekong), Bangladesh (sông Jamuna - Padma), và Myanmar (sông Ayeyarwady).

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

                                      MA TRẬN ĐỐI TƯỢNG HƯỞNG LỢI
  ┌─────────────────────────┬───────────────────────────────────────────┬───────────────────────────────────────────┐
  │    Nhóm đối tượng       │              Giá trị tiếp nhận            │            Lợi ích định lượng             │
  ├─────────────────────────┼───────────────────────────────────────────┼───────────────────────────────────────────┤
  │ NCS & Giới học thuật    │ Khung lý thuyết tích hợp Fs = f(x,y,z),   │ Tiết kiệm 40% thời gian xây dựng mô hình; │
  │                         │ bộ dữ liệu 200 lỗ khoan chuẩn hóa         │ Khai phá 3 hướng nghiên cứu mới           │
  ├─────────────────────────┼───────────────────────────────────────────┼───────────────────────────────────────────┤
  │ Kỹ sư thiết kế nền móng │ Bộ chỉ tiêu cơ lý chi tiết 5 phụ kiểu;   │ Giảm thiểu rủi ro biến dạng công trình;   │
  │                         │ Giải pháp cọc đất - xi măng, cừ Larsen    │ Tối ưu hóa 25-30% chi phí kết cấu bảo vệ  │
  ├─────────────────────────┼───────────────────────────────────────────┼───────────────────────────────────────────┤
  │ Cơ quan quản lý & Đô thị│ Bản đồ phân vùng dự báo tỷ lệ lớn;       │ Giảm thiểu thiệt hại hạ tầng giao thông;  │
  │ (UBND, Sở TN&MT, GTVT)  │ Vị trí điều chỉnh luồng tàu & mỏ cát      │ Đảm bảo an toàn sinh mạng cho >1.500 hộ   │
  └─────────────────────────┴───────────────────────────────────────────┴───────────────────────────────────────────┘

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc mở rộng định luật sức kháng cắt Coulomb - Mohr trong cơ học đất bão hòa sang môi trường địa chất trầm tích Đệ Tứ chứa "thấu kính cát không liên tục". Luận án đã chứng minh hiện tượng mất ổn định mái dốc không chỉ do vượt ứng suất cắt tiếp tuyến thuần túy mà là quá trình hủy hoại liên hoàn do xói ngầm cơ học (piping) bên trong thấu kính cát dưới gradien thấm, làm sụp đổ rỗng xốp tầng đỡ và chuyển đổi trạng thái ứng suất cục bộ sang biến dạng dẻo trượt phẳng/trượt cung tròn sâu.

2. Điểm mới về phương pháp luận so với các nghiên cứu quốc tế tương tự là gì?

So với các nghiên cứu của Massimo Rinaldi et al. (2008) tại Ý hay Larissa Laderoute et al. (2013) tại Canada, điểm mới của luận án là phương pháp tiếp cận "Địa kỹ thuật công trình liên ngành": tích hợp đồng thời trường dữ liệu địa chất cấu trúc 3 chiều (phân chia 2 kiểu, 5 phụ kiểu nền), trường thủy động lực 2D phi cấu trúc (MIKE 21), thí nghiệm xói mòn phản lực tia nước (Jet Index Test) và phân tích cân bằng giới hạn biến dạng (GeoSlope/W, PLAXIS). Phương pháp này giúp khắc phục hoàn toàn sự phiến diện giữa trường phái thuần túy thủy lực và trường phái thuần túy địa chất.

3. Phát hiện nào trong dữ liệu thực nghiệm gây bất ngờ nhất và được giải thích ra sao?

Phát hiện bất ngờ nhất là tại khu vực đầu cù lao Bình Thạnh và cù lao Mỹ Hòa Hưng, năng lượng phá hủy bề mặt đường bờ chủ yếu đến từ sóng tàu thuyền có tải trọng vận tải lớn ($H > 0.3\text{ m}$) chứ không phải từ lưu tốc dòng chảy tự nhiên mùa lũ. Dữ liệu thực nghiệm Jet test chứng minh sức kháng cắt bề mặt đất bùn sét bãi bồi ($\tau_c = 0.42 - 0.78\text{ Pa}$) thấp hơn nhiều so với ứng suất cắt động do sóng tàu tạo ra ($\tau_b = 1.25 - 2.10\text{ Pa}$), gây bóc tách lớp phủ bề mặt liên tục với tốc độ lùi bờ lên tới 22 - 30 m/năm.

4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) cho các lưu vực sông khác không?

Luận án cung cấp một quy trình nhân bản 4 bước chuẩn hóa cao độ: (1) Thiết lập mạng lưới khoan khảo sát và phân chia kiểu cấu trúc nền theo tiêu chuẩn địa tầng Đệ Tứ; (2) Thu thập dữ liệu thủy lực mùa kiệt/mùa lũ và thiết lập mô hình 2D hiệu chỉnh hệ số $E_{\text{NASH}} > 0.85$; (3) Thí nghiệm cơ lý chuyên sâu kết hợp Jet test xác định ngưỡng xói $\tau_c, K_d$; (4) Tích hợp các lớp thông tin lên nền GIS để thành lập bản đồ phân vùng dự báo ổn định phục vụ thiết kế công trình bảo vệ.

5. Khung chương trình nghiên cứu 10 năm (10-year research agenda) kế thừa luận án được vạch ra như thế nào?

Khung chương trình 10 năm bao gồm:

  • Giai đoạn 1 (Năm 1 - 3): Xây dựng mạng lưới quan trắc biến dạng bờ sông thời gian thực bằng công nghệ viễn thám radar khẩu độ tổng hợp giao thoa (InSAR) kết hợp cảm biến địa kỹ thuật IoT dọc tuyến sông Hậu và sông Tiền.
  • Giai đoạn 2 (Năm 4 - 6): Mở rộng mô hình thủy cơ liên hợp 3D (3D Hydro-Mechanical Coupled Model) trên toàn bộ hệ thống phân lưu sông Cửu Long, đánh giá tác động lũy tích của chuỗi đập thủy điện dòng chính Mê Kông.
  • Giai đoạn 3 (Năm 7 - 10): Phát triển vật liệu sinh thái địa kỹ thuật mới (Geopolymers, vi sinh kết tủa canxi MICP) gia cố nền đất yếu và hoàn thiện khung thể chế quản trị thích ứng sạt lở sông ngòi cho toàn vùng châu thổ Nam Bộ.

Kết luận

Công trình nghiên cứu tiến sĩ của tác giả Trần Lê Thế Diễn đã giải quyết trọn vẹn, khoa học và hệ thống bài toán ổn định bờ sông Hậu qua tỉnh An Giang với các đóng góp cốt lõi:

  1. Xác lập phương pháp luận liên ngành địa chất công trình – thủy động lực học, làm sáng tỏ bản chất nguyên nhân mất ổn định bờ sông là hàm đa biến $F_s = f(x, y, z)$.
  2. Lần đầu tiên phân chia và mô tả chi tiết hệ thống cấu trúc nền đất yếu gồm 2 kiểu và 5 phụ kiểu (IA, IB, IC, IIA, IIB), phát hiện vai trò kích nổ của các thấu kính cát trong cơ chế trượt sâu.
  3. Nhận diện và định lượng 3 cơ chế phá hoại: trượt lở mái dốc, xói lở chân bờ và xói mòn bóc vỏ do sóng tàu thuyền ($H > 0.3\text{ m}$).
  4. Xây dựng bản đồ phân vùng dự báo ổn định bờ sông Hậu tích hợp công nghệ GIS và kiểm định thủy lực đạt hệ số $E_{\text{NASH}} > 0.85 - 0.94$.
  5. Đề xuất hệ giải pháp công trình chỉnh trị thích ứng theo từng phân vùng: cọc xi măng đất đông kết thấu kính cát cho xã Bình Mỹ, kè rọ đá tiêu năng cho các cù lao, cừ Larsen cho đô thị Long Xuyên, và mỏ hàn nắn dòng phân lưu.
  6. Mở ra hướng tiếp cận quản lý rủi ro thiên tai hiện đại, đóng góp thiết thực cho sự nghiệp phát triển bền vững kinh tế - xã hội vùng đồng bằng sông Cửu Long.