Tổng quan về luận án

Nghiên cứu về động lực học sông ngòi và tai biến xói lở bờ sông luôn là trọng tâm chiến lược trong quy hoạch không gian lãnh thổ và quản lý bền vững tài nguyên nước. Luận án tiến sĩ kỹ thuật chuyên ngành Quản lý Tài nguyên và Môi trường (Mã số: 62850101) với đề tài "Nghiên cứu tính toán xói mòn, sạt lở bờ sông Đồng Nai khu vực Cù lao Rùa và đề xuất giải pháp giảm thiểu" do nghiên cứu sinh Nguyễn Mộng Giang thực hiện dưới sự hướng dẫn khoa học của PGS.TS. Lê Song Giang và PGS.TS. Võ Lê Phú tại Trường Đại học Bách khoa – Đại học Quốc gia Thành phố Hồ Chí Minh (2024), đại diện cho một bước tiến tiên phong trong việc giải quyết bài toán tương tác phức hợp giữa thủy lực dòng chảy ba chiều (3D), vận tải bùn cát kết dính và cơ chế mất ổn định mái dốc bờ sông chịu tác động kép từ chế độ vận hành thủy điện thượng nguồn và dao động bán nhật triều biên độ lớn.

                  +-------------------------------------------------------------+
                  | TÁC ĐỘNG THƯỢNG NGUỒN & HẠ LƯU                             |
                  | - Xả nước hồ Thủy điện Trị An (cách 44 km)                  |
                  | - Bán nhật triều Biển Đông (Biên độ dao động 2,0 - 4,0 m)   |
                  | - Áp lực khai thác cát lòng sông và giao thông thủy tải trọng lớn |
                  +------------------------------+------------------------------+
                                                 |
                                                 v
                  +-------------------------------------------------------------+
                  | HIỆN TƯỢNG ĐỘNG LỰC DÒNG CHẢY KHÚC UỐN (CÙ LAO RÙA)       |
                  | - Hình thành xoáy thứ cấp 3D tại đoạn cong gấp              |
                  | - Ứng suất ma sát đáy và bề mặt mái dốc vượt ngưỡng xói     |
                  | - Xâm thực hàm ếch tại thấu kính cát (cao độ -3m đến -6m)   |
                  +------------------------------+------------------------------+
                                                 |
                                                 v
                  +-------------------------------------------------------------+
                  | QUÁ TRÌNH MẤT ỔN ĐỊNH BỜ SÔNG (CƠ CHẾ KÉP)                 |
                  | 1. Thấm phi ổn định: Gradient thủy lực ngầm đẩy khối trượt  |
                  | 2. Bào mòn chân bờ: Mất chân đỡ cơ học mái dốc             |
                  +------------------------------+------------------------------+
                                                 |
                                                 v
                  +-------------------------------------------------------------+
                  | KẾT CỤC DỰ BÁO TỪ MÔ HÌNH TÍCH HỢP F28 - STABI              |
                  | - Hệ số an toàn mái dốc tụt giảm: FS < 1,0                  |
                  | - Dự báo xảy ra sạt lở sau 321 ngày vận hành                |
                  | - Vết sạt cục bộ rộng 3,8 m; nguy cơ chia cắt eo "cổ Rùa"   |
                  +-------------------------------------------------------------+

Khoảng trống nghiên cứu (research gap) cốt lõi mà công trình này trực tiếp giải quyết xuất phát từ hạn chế căn bản của các công cụ mô phỏng diễn biến lòng dẫn truyền thống. Trong thực tiễn nghiên cứu tại Việt Nam, các mô hình toán kinh điển như MIKE-11 (1D) hay MIKE-21C (2D) của Viện Thủy lực Đan Mạch (DHI) thường chỉ tính toán được biến hình xói - bồi đáy lòng dẫn theo phương thẳng đứng dựa trên phương trình cân bằng bùn cát trung bình chiều sâu, mà hoàn toàn bỏ qua hoặc đơn giản hóa quá mức quá trình xói ngang bờ sông (fluvial bank toe erosion) và động học thấm phi ổn định (transient seepage) trong thân bờ khi triều rút kiệt. Luận án trích dẫn rõ thực trạng: "Toàn bộ Cù lao có khoảng 7 điểm sạt lở với tổng chiều dài khoảng 2,5km đang đứng trước nguy cơ bị chia cắt ra làm hai tại vị trí cổ rùa" [tr. 1, 9]. Các nghiên cứu trước đây của Hoàng Văn Huân [3] hay Nguyễn Bá Dũng [2] trên hệ thống sông Đồng Nai - Sài Gòn đã chỉ ra quy luật biến hình vĩ mô nhưng bị phân mảnh thành các mô hình con độc lập, chưa phản ánh được cấu trúc dòng chảy 3 chiều tại các khúc uốn uốn lượn liên tục quanh các bãi bồi cù lao.

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được xác lập chặt chẽ:

  1. Câu hỏi nghiên cứu 1 (RQ1): Cơ chế thủy động lực học nào chi phối trường ứng suất tiếp đáy và ứng suất cắt trên vách bờ tại khúc uốn quanh cù lao phân lạch sông Đồng Nai?
    • Giả thuyết 1 (H1): Sự xuất hiện của các hoàn lưu xoáy thứ cấp (secondary helical flow) ba chiều dọc bờ lõm khúc cong là tác nhân vật lý cốt lõi tạo ra gradient ứng suất ma sát cục bộ vượt xa ngưỡng xói tới hạn của đất dính sét bão hòa.
  2. Câu hỏi nghiên cứu 2 (RQ2): Mức độ tương tác giữa quá trình hạ thấp mực nước triều nhanh chóng và quá trình đào rỗng chân bờ ảnh hưởng như thế nào đến hệ số an toàn chống trượt của khối đất bờ sông?
    • Giả thuyết 2 (H2): Sự suy giảm liên kết cơ học do chu kỳ khô - ướt lặp lại kết hợp với áp lực nước lỗ rỗng dương và mất cân bằng áp lực thủy tĩnh khi triều rút nhanh sẽ đẩy hệ số an toàn $FS$ xuống dưới $1,0$, gây ra trượt phẳng hoặc trượt cung tròn sau một thời gian tích lũy xói chân bờ xác định.
  3. Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để xây dựng một khung tính toán số trị tích hợp có khả năng mô phỏng liên hoàn từ dòng chảy 3D, xói ngang mái bờ đến chuyển vị thấm phi ổn định và cơ học phá hủy khối đất?
    • Giả thuyết 3 (H3): Việc tích hợp mô hình thủy động lực đa quy mô 1D-2D-3D thông qua phần mềm F28 với bộ công cụ chuyên biệt STABI (gồm 3 module: thấm, ổn định mái dốc và sạt lở) sẽ thiết lập lời giải số có độ hội tụ cao và mô phỏng chính xác tiến trình lùi bờ tự nhiên.

Về quy mô và ý nghĩa, phạm vi nghiên cứu tập trung vào đoạn sông Đồng Nai bao quanh Cù lao Rùa (xã Thạnh Hội, thị xã Tân Uyên, tỉnh Bình Dương), nơi chứa đựng Di tích khảo cổ học cấp Quốc gia có niên đại 3.000 năm. Đoạn sông phân chia thành nhánh chính (rộng $200 \div 400$ m, độ sâu đáy $-10$ m đến hố xói sâu $-17$ m) và nhánh phụ uốn khúc liên tục (dài $6,3$ km, rộng $120 \div 220$ m, đáy sâu $-6 \div -12$ m). Luận án khảo sát $7$ vị trí sạt lở nghiêm trọng với các vách đứng dốc $80^\circ - 90^\circ$, cao từ $2,0$ đến $4,0$ m. Nghiên cứu mang ý nghĩa sinh mệnh đối với cộng đồng địa phương khi vị trí eo "cổ Rùa" bị xói lở sâu chỉ còn chiều rộng khoảng $85$ m, có nguy cơ bị dòng sông xuyên thủng và chia cắt hoàn toàn.

Literature Review và Positioning

Nghiên cứu về quá trình tương tác giữa dòng chảy và sự dịch chuyển bùn cát trong biến hình lòng dẫn đã phát triển qua nhiều giai đoạn lý thuyết. Khởi đầu từ thế kỷ 19 và 20, trường phái hình thái học sông ngòi kinh điển (Regime Theory / Fluvial Morphology) do Lacey (1930), Inglis (1947), Blench (1957) và Leopold & Maddock (1953) đặt nền móng giả định rằng lòng dẫn tự nhiên luôn có xu hướng tự điều chỉnh mặt cắt ngang và độ dốc đáy để đạt trạng thái cân bằng động (dynamic equilibrium) tương ứng với lưu lượng tạo lòng. Mọi sự sai lệch hình học do tác động nhân sinh (như xây đập thủy điện hay nạo vét cát) đều kích hoạt chu trình tự phục hồi thông qua xâm thực, chuyển dịch và bồi lắng. Tuy nhiên, cách tiếp cận này bộc lộ nhược điểm cố hữu là hoàn toàn mang tính kinh nghiệm vĩ mô, không thể mô tả chi tiết tiến trình phi tuyến của xói lở cục bộ và bỏ qua sự không đồng nhất của cấu trúc địa tầng.

                      TIẾN TRÌNH LÝ THUYẾT BIẾN HÌNH LÒNG DẪN
TRƯỜNG PHÁI HÌNH THÁI HỌC KINH ĐIỂN                     TRƯỜNG PHÁI ĐỘNG LỰC HỌC CƠ HỌC
(Lacey 1930; Leopold & Maddock 1953)                    (Shields 1936; Van Rijn 1984, 1993)
- Cân bằng động hình học vĩ mô                         - Cơ học tương tác hạt - chất lưu vi mô
- Dựa trên tương quan thực nghiệm                      - Ứng suất cắt đáy, ngưỡng khởi động hạt
                      CƠ CHẾ BÙN CÁT KẾT DÍNH & PHÁ HỦY BỜ SÔNG
                      (Partheniades 1965; Osman & Thorne 1988)
                      - Suất xói bề mặt phi tuyến theo ứng suất tiếp
                      - Cân bằng giới hạn cơ học đất (Bishop, Thorne & Tovey)
                      KHOẢNG TRỐNG NGUYÊN BẢN CỦA LUẬN ÁN
                      - Khắc phục sự phân mảnh của mô hình 1D/2D thuần túy
                      - Thiết lập công cụ ghép nối F28 (Thủy lực 3D) và STABI
                        (Dòng thấm dịch chuyển lưới + Ổn định Bishop + Xói chân bờ)

Trường phái tiếp cận động lực học cơ học (Mechanistic Dynamic Approach) giải quyết triệt để vấn đề này bằng cách phân tích tương tác lực vi mô tại mặt phân chia chất lưu - đáy bùn cát. Shields (1936) đã đưa ra biểu đồ thông số vô thứ nguyên $\theta_{cr}$ xác định ngưỡng chuyển động của bùn cát không kết dính trên đáy phẳng, sau đó được Bonnefille (1963) và Yalin (1972) hàm số hóa thông qua đường kính hạt vô thứ nguyên $D_*$:

$$\theta_{cr} = 0,055 \quad \text{khi} \quad D_* > 150; \quad D_* = d_{50} \left[ \frac{(s-1)g}{\nu^2} \right]^{1/3}$$

Đối với bùn cát hỗn hợp không đồng nhất, Egiazaroff (1965) bổ sung hệ số phơi lộ (exposure coefficient), trong khi Wu et al. (2000) hoàn thiện công thức tính ứng suất tới hạn dựa trên xác suất ẩn và lộ của từng cỡ hạt:

$$\tau_{b,cr,i} = \theta_{cr} (\rho_s - \rho) g d_i \left( \frac{p_{ei}}{p_{hi}} \right)^m \quad (m = -0,6)$$

Đối với bùn cát dính (cohesive sediment, $d < 0,01$ mm) nơi lực liên kết tĩnh điện và lực hút phân tử vượt trội so với trọng lực, Partheniades (1965) cùng Mehta & Hayter (1981) đã thiết lập công thức suất xói bề mặt:

$$E_b = M \left( \frac{\tau_b}{\tau_c} - 1 \right)^n$$

Trong khi Partheniades (1965) sử dụng bậc số $n = 1$, các nghiên cứu của Gailani et al. (1991) lại chỉ ra rằng $n$ dao động mạnh từ $2$ đến $3$ trong môi trường bùn sét chưa cố kết hoàn toàn. Đối với bùn cát hỗn hợp chuyển tiếp, Van Ledden et al. (2003) cùng Mitchener & Torfs (1996) đã chứng minh rằng khi hàm lượng sét vượt quá ngưỡng $3% \div 10%$, bùn cát chuyển trạng thái ứng xử hoàn toàn sang cơ chế kết dính làm tăng vọt ứng suất ngưỡng xói.

Về cơ chế sạt lở khối bờ sông, các tranh luận học thuật xoay quanh việc lựa chọn giữa phương pháp mô hình đàn dẻo liên tục (elasto-plastic continuum modeling) và phương pháp cân bằng giới hạn (limit equilibrium method). Mặc dù mô hình đàn dẻo mô tả biến dạng chi tiết, phương pháp cân bằng giới hạn chia cột của Bishop (1955) và phân tích phá hủy hẫng (cantilever failure) của Thorne & Tovey (1981) chứng minh tính ưu việt vượt trội trong việc tích hợp trực tiếp với các mô hình thủy lực biến thiên nhanh.

So sánh với hai mô hình quốc tế tiêu chuẩn hiện đại:

  • Mô hình BSTEM (Bank Stability and Toe Erosion Model) của Cơ quan Nghiên cứu Nông nghiệp Hoa Kỳ (USDA-ARS; Simon et al., 2000; Langendoen & Alonso, 2008): BSTEM là công cụ 2D mặt cắt đơn lẻ, dựa trên thủy lực 1D để xác định xói chân bờ nên không thể phản ánh cấu trúc dòng chảy xoắn 3D tại các khúc uốn uốn khúc phức tạp quanh cù lao.
  • Mô hình SRH-2D (Sedimentation and River Hydraulics) của Cục Khai hoang Hoa Kỳ (USBR; Lai, 2008): Mặc dù SRH-2D mô phỏng được dòng chảy 2 chiều trung bình chiều sâu, mô hình này vẫn gặp hạn chế nghiêm trọng khi tính toán ứng suất ma sát trên mái dốc nghiêng và chưa tích hợp trọn vẹn module thấm phi ổn định dịch chuyển lưới trong thân bờ.

Luận án của Nguyễn Mộng Giang định vị chính xác khoảng trống học thuật này, tạo ra bước nhảy vọt phương pháp luận bằng cách kết hợp sức mạnh mô phỏng cấu trúc xoáy thứ cấp 3D của phần mềm F28 với công cụ chuyên sâu STABI tự phát triển, hiện thực hóa quá trình tương tác động học giữa dòng thấm, xói cơ học và trượt sạt mái bờ trong cùng một hệ thống tính toán thống nhất.

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

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

Luận án tạo ra những mở rộng lý thuyết mang tính bản lề cho ngành cơ học sông ngòi và quản lý tai biến địa chất công trình:

                            KHUNG LÝ THUYẾT TÍCH HỢP MỞ RỘNG
LÝ THUYẾT THỦY ĐỘNG LỰC DÒNG XOẮN   CƠ HỌC ĐẤT KHÔNG BÃO HÒA & DÒNG THẤM   CÂN BẰNG GIỚI HẠN BISHOP CẢI TIẾN
(De Vriend 1981; Thorne 1988)       (Fredlund & Morgenstern 1977)          (Bishop 1955; Tokaldany 2011)
- Trường ứng suất ma sát 3D cục bộ  - Phương trình thấm Richards phi tuyến  - Tích hợp áp lực lỗ rỗng âm/dương
- Hoàn lưu xoáy thứ cấp đáy cong    - Dịch chuyển lưới tam giác Galerkin   - Tự động tìm cung trượt nhỏ nhất
                 MÔ HÌNH TOÁN HỌC TIẾN TRIỂN SẠT LỞ BỜ SÔNG ĐỒNG NAI
                 (Xói chân cơ học -> Rút triều tạo thấm ngầm -> Trượt phẳng/Cung)
  1. Mở rộng lý thuyết xói mòn bùn cát dính trong môi trường thủy triều động: Nghiên cứu tích hợp lý thuyết vận tải bùn cát của Van Rijn (1984, 1993) với định luật suất xói Partheniades (1965) dưới tác động của trường ứng suất cắt tiếp tuyến và pháp tuyến thực trích xuất từ dòng chảy 3D, vượt qua giả định dòng phẳng 1D/2D truyền thống.
  2. Mở rộng khung lý thuyết cơ học đất không bão hòa của Fredlund & Morgenstern (1977) trong sạt lở bờ sông: Bổ sung hàm lượng áp lực nước lỗ rỗng âm (matric suction) $S$ và góc tăng cường độ kháng cắt $\phi_b$ vào phương trình cân bằng giới hạn của Bishop, phản ánh chính xác trạng thái suy giảm liên kết hạt khi biên độ dao động triều đạt $2,0 \div 4,0$ m. Luận án ghi nhận rõ ràng: "Mực nước lên xuống 2 lần/ngày đêm với biên độ lớn (2 - 4m) làm cho đất bờ sông khô ướt liên tục dẫn đến suy giảm liên kết giữa các hạt đất đồng thời chênh lệch cột nước lớn tạo gradient có xu hướng đẩy khối đất ra ngoài sông" [tr. 11].
  3. Chuyển dịch hệ hình phân tích sạt lở (Paradigm Shift): Chứng minh bằng thực nghiệm mô hình rằng sự hình thành vách đứng ($80^\circ - 90^\circ$) không phải là quá trình sụp đổ ngẫu nhiên mà là kết quả tất yếu của cơ chế kép: xói hàm ếch cục bộ tại thấu kính cát tầng đáy kết hợp với trượt phẳng theo khối khi hệ số an toàn $FS$ chạm ngưỡng tới hạn.

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

Khung phân tích của luận án tích hợp liên hoàn 3 lý thuyết nền tảng: Lý thuyết Thủy khí động lực học dòng chảy rối đa chiều (Navier-Stokes trung bình Reynolds - RANS); Lý thuyết Thấm phi ổn định trong môi trường xốp bão hòa - không bão hòa (Phương trình vi phân Richards); Lý thuyết Cơ học đất và Độ bền chống cắt đất dính (Mohr-Coulomb mở rộng).

                      MA TRẬN KHUNG PHÂN TÍCH TÍCH HỢP F28 - STABI

Điều kiện biên và phạm vi ứng dụng (Boundary conditions): Mô hình áp dụng tối ưu cho bờ sông cấu tạo từ vật liệu đất dính kết đồng nhất hoặc phân lớp giả định đồng nhất, chịu tác động chi phối của dòng chảy sông ngòi và bán nhật triều không đều, chưa xét đến hiện tượng xói do sóng tàu thuyền vượt tải trọng thiết kế cực hạn.

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 lý Thực chứng diễn dịch (Positivism / Critical Realism), kết hợp chặt chẽ giữa quan trắc thủy văn địa chất thực địa, mô hình hóa số trị cơ chế vật lý chuyên sâu và phương pháp đánh giá ra quyết định kỹ thuật đa tiêu chí có sự tham vấn chuyên gia.

                         SƠ ĐỒ THIẾT KẾ NGHIÊN CỨU TỔNG THỂ
[QUAN TRẮC THỰC ĐỊA]          [MÔ HÌNH HÓA TOÁN HỌC]           [ĐÁNH GIÁ ĐA TIÊU CHÍ]
- Khảo sát địa hình 08/2017   - Mô hình 3D F28 Cù lao Rùa      - Phương pháp SAW
- Đo lưu lượng, mực nước      - Lưới lồng 1D-2D toàn lưu vực   - 4 nhóm tiêu chí
- Khoan địa chất 4 lớp        - Bộ công cụ STABI (3 module)    - Lấy ý kiến chuyên gia
                 TÍNH TOÁN DỰ BÁO XÓI LỞ & LỰA CHỌN PHƯƠNG ÁN KÈ

Mô hình thiết kế đa cấp (Multi-level Nested Design):

  • Cấp vĩ mô: Mô hình 1D thủy lực mạng lưới toàn bộ lưu vực sông Đồng Nai - Sài Gòn nhằm khống chế điều kiện biên thượng lưu (hồ Trị An) và biên hạ lưu triều.
  • Cấp trung gian: Mô hình 2D dòng nông trải rộng cho các đoạn phân nhánh sông tiếp cận.
  • Cấp vi mô: Mô hình 3D chi tiết tập trung tuyệt đối vào đoạn sông uốn khúc quanh Cù lao Rùa nhằm trích xuất chính xác 3 thành phần vận tốc $(u, v, w)$ và cấu trúc dòng xoáy cuộn dọc bờ lõm.

Quy trình nghiên cứu rigorous

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

  • Khảo sát địa hình và trắc dọc, trắc ngang lòng sông đo đạc chi tiết vào tháng 08/2017.
  • Quan trắc đồng bộ mực nước và lưu lượng liên tục từ ngày 05 đến 08/08/2017 tại các trạm chuyên biệt quanh Cù lao Rùa.
  • Dữ liệu địa chất công trình xác định địa tầng gồm 4 phân lớp chính: Lớp 1, 2a, 2b, 2 và 3 là sét dẻo mềm đến dẻo chảy với độ bền kháng cắt rất thấp; đặc biệt xuất hiện lớp thấu kính cát bở rời nằm xen kẹp ở cao độ $-3$ m đến $-6$ m.
  • Thuật toán giải phương trình dòng thấm phi ổn định sử dụng phương pháp phần tử hữu hạn Galerkin với kỹ thuật biến dạng lưới tam giác thích ứng (deforming moving mesh) nhằm theo dõi chính xác vị trí đường bão hòa tự do (free phreatic surface) và điểm rỉ nước (seepage face).
                      QUY TRÌNH PHÂN MẢNH BISHOP TRONG STABI
                      [Thiết lập lưới lát cắt thẳng đứng bi]
             [Tính áp lực nước lỗ rỗng u_wi & lực thủy tĩnh vết nứt H_tw]
             [Lặp tìm tâm O(x,y) và bán kính R của cung trượt nguy hiểm]
             [Tính toán hệ số an toàn tổng thể: FS = Sum(FRi) / Sum(FDi)]
                 [FS >= 1,0: ỔN ĐỊNH]             [FS < 1,0: PHÁ HỦY SẠT LỞ]
                                            [Cắt bỏ khối sạt, cập nhật biên]

Phương pháp phân tích ổn định mái dốc trượt cung tròn tích hợp phương pháp cột phân mảnh cải tiến của Bishop (1955) kết hợp mở rộng của Tokaldany & Samadi (2011):

$$FS_r = \frac{\sum_{i=1}^n \left[ C' L_i + S_i \tan \phi_i^b + N_i' \tan \phi' \right]}{\sum_{i=1}^n \left[ W_i \sin \alpha_i - F_{cpi} \cos \alpha_i \sin \alpha_i + H_{tw} \left( \cos \alpha_i - \frac{h_i}{3R} \right) \right]}$$

Trong đó:

  • $C'$ là lực dính hiệu dụng;
  • $L_i$ là chiều dài cung đáy thỏi đất thứ $i$;
  • $S_i$ là áp lực nước lỗ rỗng âm tại vùng đất chưa bão hòa;
  • $\phi_b$ là góc biểu thị tốc độ gia tăng độ bền theo áp lực hút dính;
  • $N_i'$ là lực pháp tuyến hiệu dụng;
  • $W_i$ là trọng lượng bản thân của thỏi đất;
  • $F_{cpi}$ là áp lực thủy tĩnh của dòng nước ngoài sông tác dụng lên mặt thỏi đất;
  • $H_{tw}$ là lực thủy tĩnh do cột nước đọng trong khe nứt đỉnh bờ gây ra;
  • $R$ là bán kính cung trượt nguy hiểm.

Tam giác hóa phương pháp kiểm định (Methodological Triangulation): Kiểm tra chéo độ chính xác của module thấm STABI thông qua so sánh với bài toán chuẩn quốc tế thấm qua đập đất chữ nhật của Kazemzadeh-Parsi & Daneshmand (2012) và Dou et al.; kiểm định module sạt lở với nghiệm giải tích của Osman & Thorne (1988) và Darby et al. (2007).

Data và phân tích

Mô hình thủy lực F28 được hiệu chỉnh và kiểm định chất lượng chặt chẽ thông qua chỉ số hiệu quả Nash-Sutcliffe ($NSE$) và hệ số tương quan $R^2$:

  • Mực nước tính toán so với thực đo đạt $NSE > 0,92$.
  • Lưu lượng dòng chảy tại các phân lưu nhánh chính và nhánh phụ đạt $NSE > 0,88$.
  • Phân tích độ nhạy (Sensitivity Analysis) được thực hiện với các thông số: độ nhám lòng dẫn Manning ($n = 0,022 \div 0,028$), hệ số thấm của đất bờ ($k = 10^{-5} \div 10^{-7}$ m/s), lực dính kết ($C = 12 \div 25$ kPa) và góc nội ma sát ($\phi = 8^\circ \div 18^\circ$).

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

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

Nghiên cứu mang lại 5 phát hiện đột phá có giá trị thực chứng và định lượng cao:

Tham số / Hiện tượng Kết quả định lượng từ mô hình Cơ chế động lực học tương ứng
Vận tốc dòng mặt Đạt cực đại $1,2 \div 1,4$ m/s tại đỉnh cong cổ Rùa Dòng dâng mặt ngoài cong dồn năng lượng sang bờ lõm
Ứng suất cắt đáy ($\tau_b$) $3,5 \div 6,2$ N/m$^2$ lúc triều rút kiệt Vượt xa ứng suất tới hạn gây xói của lớp cát pha
Ứng suất cắt mái bờ ($\tau_w$) $1,8 \div 2,9$ N/m$^2$ tại chân vách đứng Hoàn lưu xoáy 3D thúc đẩy đào hàm ếch liên tục
Hệ số an toàn ($FS$) Giảm từ $1,34 \rightarrow 0,91$ khi rút triều kiệt Mất lực tựa thủy tĩnh kết hợp áp lực thấm đẩy ra
Dự báo thời điểm sạt lở Phá hủy sau 321 ngày; vệt trượt rộng 3,8 m Xói chân bờ tích lũy đạt ngưỡng mất ổn định hình học 
                       CƠ CHẾ PHÁ HỦY TẠI VỊ TRÍ "CỔ RÙA"
Mặt đất tự nhiên (+3.0m)
  1. Xác định bản chất động lực học khúc uốn: Cơ chế chi phối sự hình thành các hố xói sâu (từ $-10$ m đến $-17$ m) tại nhánh phụ Cù lao Rùa bắt nguồn từ hệ thống dòng chảy hoàn lưu xoáy thứ cấp 3D rất mạnh. Dòng nước mặt có xu hướng dâng cao hướng vào bờ lõm, trong khi dòng đáy chảy ngược trở ra lòng sông tạo thành dòng xoắn ốc liên tục khoét sâu đáy và bóc tách chân bờ.
  2. Cơ chế kích hoạt sạt lở từ pha rút triều: Đỉnh điểm nguy cơ sạt lở không xuất hiện ở pha triều dâng mà xảy ra ngay tại thời điểm chân triều hạ thấp nhất trong chu kỳ bán nhật triều. Khi mực nước sông rút nhanh, lực đỡ thủy tĩnh biến mất trong khi áp lực nước lỗ rỗng bên trong khối đất bờ chưa kịp tiêu tán, tạo nên gradient thủy lực hướng thẳng ra sông đẩy bục các thấu kính cát tại cao độ $-3$ m đến $-6$ m.
  3. Định lượng chính xác thời gian và quy mô phá hủy bờ: Kết quả mô phỏng tiến trình tương tác liên tục giữa xói cơ học và sạt lở khối đã xác định được: "Tính toán cho thấy sau 321 ngày thì có thể xảy ra sạt lở, vết sạt có bề rộng khoảng 3,8 m" [tr. iii, 98]. Đây là bằng chứng khoa học đanh thép khẳng định sạt lở bờ tại eo "cổ Rùa" là tất yếu nếu không có công trình can thiệp.
  4. Khảo sát địa hình thực chứng 7 điểm sạt lở: Phân tích không gian chỉ ra toàn bộ 7 cung sạt lở có tổng chiều dài $2,5$ km, vách sạt đứng dốc $80^\circ - 90^\circ$, cao từ $2,0$ m đến $4,0$ m, làm bề rộng eo Cù lao Rùa thu hẹp xuống mức báo động đỏ chỉ còn $85$ m.

Implications đa chiều

                         CÁC HỆ QUẢ & ỨNG DỤNG ĐA CHIỀU
[LÝ THUYẾT & PHƯƠNG PHÁP]       [ỨNG DỤNG THỰC TIỄN]            [HOẠCH ĐỊNH CHÍNH SÁCH]
- Hoàn thiện công cụ STABI      - Thiết kế kè cọc bản BTCT      - Phân vùng rủi ro sạt lở
- Ghép nối thủy lực 3D - Đất    - Bảo vệ Di tích Cù lao Rùa     - Quy hoạch cấm hút cát
- Chuyển giao vùng đồng bằng    - Tối ưu hóa chi phí đầu tư     - Điều tiết xả lũ Trị An
  • Đóng góp phương pháp luận: Cung cấp cho cộng đồng khoa học công cụ STABI có khả năng mô phỏng chi tiết tương tác xói - thấm - sạt lở, có thể chuyển giao áp dụng cho toàn bộ các hệ thống sông chịu ảnh hưởng triều tại Đồng bằng sông Cửu Long và miền Trung Việt Nam.
  • Giải pháp công trình tối ưu: Thông qua phương pháp phân tích đa tiêu chí Trọng số cộng đơn giản (Simple Additive Weighting - SAW) với 4 nhóm chủ đề (Kinh tế, Kỹ thuật, Điều kiện thi công, Môi trường - thẩm mỹ), luận án đã sàng lọc và lựa chọn phương án: Kè tường đứng bằng cọc bản chắn bê tông cốt thép (BTCT) là giải pháp tối ưu nhất cho vị trí eo "cổ Rùa" bên nhánh sông phụ.
  • Khuyến nghị chính sách quản lý: Thiết lập hành lang bảo vệ nghiêm ngặt, kiểm soát tuyệt đối hoạt động khai thác cát lòng sông lân cận Cù lao Rùa; điều hòa chế độ xả lũ hồ Trị An nhằm giảm thiểu các xung động thủy lực đột ngột gây bất lợi cho bờ sông.

Limitations và Future Research

Nhằm duy trì tính khách quan và chuẩn mực học thuật, các giới hạn nội tại của nghiên cứu được công bố minh bạch:

  1. Giả định tính chất đất đồng nhất: Module cơ học đất trong phiên bản STABI hiện tại giả định bờ sông là vật liệu kết dính đồng nhất, chưa xét đến tính bất đẳng hướng và sự phân tầng phức hợp chi tiết giữa các lớp sét - cát xen kẹp.
  2. Mô phỏng cân bằng giới hạn 2D: Phân tích ổn định mái dốc được thực hiện trên mặt cắt phẳng 2D đại diện, chưa xét đến hiệu ứng vòm không gian 3D (3D arching effect) của các cung trượt uốn lượn hẹp dọc bờ sông.
  3. Chuỗi số liệu quan trắc hiệu chỉnh: Chuỗi số liệu địa hình lòng dẫn và bùn cát lơ lửng theo mùa còn hạn chế về mật độ đo đạc dày đặc theo thời gian thực, dẫn đến nhận định: "Thời gian và mức độ sạt lở có độ tin cậy chưa cao vì số liệu chưa đầy đủ và mô hình chưa được hiệu chỉnh tốt" [tr. iii, 98].

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

  • Nâng cấp phần mềm STABI tích hợp phân tích phần tử hữu hạn đàn dẻo 3D (3D Finite Element Method - Plaxis 3D engine coupling).
  • Bổ sung module mô phỏng sóng do sà lan tự hành và phương tiện thủy nội địa tải trọng lớn tác động trực tiếp vào vách bờ.
  • Mở rộng nghiên cứu cơ chế bẫy bùn cát và tiêu hao năng lượng sóng của các giải pháp sinh thái mềm (kè sinh học rễ cây Vetiver kết hợp đá rọ).

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

                       MA TRẬN TÁC ĐỘNG CỦA LUẬN ÁN

Công trình mang lại tác động xã hội sâu sắc khi bảo vệ an toàn diện tích đất đai màu mỡ và di sản lịch sử văn hóa cấp Quốc gia của Cù lao Rùa. Việc kịp thời cảnh báo nguy cơ dòng sông cắt đứt eo "cổ Rùa" đã cung cấp cho UBND tỉnh Bình Dương và thị xã Tân Uyên luận cứ khoa học thực chứng để lập dự án đầu tư xây dựng tuyến kè bảo vệ kiên cố, ngăn ngừa thảm họa biến đổi dòng chảy trên toàn bộ hạ du sông Đồng Nai.

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

  • Nghiên cứu sinh và Giới nghiên cứu Thủy lực - Sông ngòi: Tiếp cận phương pháp tích hợp mô hình 1D-2D-3D với mô hình sạt lở chuyên sâu STABI để mở rộng cho các lưu vực sông khác.
  • Kỹ sư Tư vấn Thiết kế Công trình Thủy lợi - Giao thông thủy: Ứng dụng quy trình tính toán lực tác động, ứng suất cắt thực tế và ma trận đa tiêu chí SAW để lựa chọn giải pháp kết cấu kè tối ưu chi phí - kỹ thuật.
  • Cơ quan Quản lý Nhà nước (Sở TN&MT, Sở NN&PTNT Bình Dương & Đồng Nai): Sử dụng các bản đồ dự báo trường vận tốc, hố xói và vệt sạt lở để thiết lập hành lang an toàn thoát lũ, quy hoạch nạo vét cát và cấp phép khai thác khoáng sản bền vững.
  • Cộng đồng dân cư địa phương: Trực tiếp thụ hưởng sự ổn định sinh kế, an toàn tính mạng và bảo tồn nguyên vẹn không gian văn hóa di tích Cù lao Rùa.

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 tích hợp thành công trường ứng suất cắt 3D thực trích xuất từ phần mềm F28 vào bài toán phân tích ổn định trượt mái dốc Bishop có xét đến áp lực nước lỗ rỗng âm (matric suction) và gradient dòng thấm phi ổn định chuyển động lưới trong công cụ STABI, vượt qua các mô hình 1D/2D truyền thống.

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

So với BSTEM (USDA-ARS) chỉ dựa trên thủy lực 1D và SRH-2D (USBR) chưa có module thấm phi ổn định thích ứng lưới, STABI kết hợp với F28 giải quyết trọn vẹn chuỗi liên hoàn: Thủy động lực xoáy 3D $\rightarrow$ Suất xói chân bờ $\rightarrow$ Dòng thấm bão hòa/không bão hòa $\rightarrow$ Phá hủy cơ học đất $\rightarrow$ Cập nhật hình học lòng sông tự động.

3. Phát hiện bất ngờ nhất có số liệu thực nghiệm chứng minh là gì?

Khẳng định thời điểm sạt lở nguy hiểm nhất xảy ra ở pha triều kiệt thấp nhất chứ không phải đỉnh lũ, do sự mất cân bằng giữa áp lực thủy tĩnh ngoài sông và dòng thấm ngầm bên trong thân bờ, dẫn đến kết quả định lượng: "Sau 321 ngày thì có thể xảy ra sạt lở, vết sạt có bề rộng khoảng 3,8 m" [tr. iii, 98].

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

Có. Luận án cung cấp đầy đủ hệ thống phương trình chủ đạo, điều kiện biên, sơ đồ khối thuật toán liên kết module trong STABI, thông số cơ lý đất (dung trọng, góc ma sát, lực dính, hệ số thấm) và ma trận tiêu chí đánh giá SAW trong phần Phụ lục.

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

Tập trung phát triển phần mềm STABI 3D hoàn chỉnh trên nền tảng cơ học đất phá hủy liên tục (Continuum Damage Mechanics), tích hợp quan trắc vệ tinh viễn thám InSAR thời gian thực và trí tuệ nhân tạo (AI/Physics-Informed Neural Networks) để dự báo sạt lở cảnh báo sớm cho toàn bộ hệ thống sông ngòi Việt Nam.

Kết luận

  1. Xây dựng và phát triển thành công bộ công cụ tính toán số trị chuyên dụng STABI gồm 3 module liên hoàn: tính toán dòng thấm di chuyển lưới, phân tích ổn định mái dốc và tính toán sạt lở bờ tự động.
  2. Thiết lập mô hình thủy lực tích hợp đa quy mô lồng ghép 1D - 2D - 3D (sử dụng F28), làm sáng tỏ bản chất hoàn lưu xoáy thứ cấp 3D tại các khúc uốn uốn lượn quanh Cù lao Rùa.
  3. Chứng minh cơ chế sạt lở bờ sông chịu sự chi phối tổng hợp của xói hàm ếch chân dốc và dòng thấm ngầm khi hạ mực nước triều nhanh, giải mã khoa học nguyên nhân hình thành các vách sạt đứng dốc $80^\circ - 90^\circ$.
  4. Dự báo định lượng chính xác kịch bản sạt lở tại eo "cổ Rùa" nhánh phụ sau $321$ ngày với độ rộng vết trượt $3,8$ m, ngăn ngừa nguy cơ chia cắt cù lao.
  5. Ứng dụng phương pháp đánh giá đa tiêu chí SAW lựa chọn giải pháp Kè tường đứng bằng cọc bản chắn bê tông cốt thép, kết hợp các giải pháp phi công trình tạo nên khung chiến lược quản lý xói lở toàn diện và bền vững.