Tổng quan về luận án

Việt Nam sở hữu hơn 3.260 km bờ biển, nằm trọn trong vùng khí hậu nhiệt đới gió mùa và chịu tác động trực tiếp của từ 6 đến 10 cơn bão mỗi năm. Các hình thái thiên tai cực đoan với cường độ gió lớn, sóng biển dữ dội cùng nước dâng do bão thường xuyên uy hiếp an toàn của hệ thống đê kè ven biển Bắc Bộ và Trung Bộ. Thực tiễn vận hành công trình tại tuyến đê biển Cát Hải (thành phố Hải Phòng) và nhiều địa phương cho thấy một hiện tượng phá hủy nghiêm trọng: đá hộc rải thảm bảo vệ chân kè bị sóng đánh bật khỏi vị trí, liên tục trượt lên xuống theo dòng leo và dòng rút trên mái kè, tạo ra ma sát chà xát cơ học làm thủng lớp bê tông bảo vệ, dẫn tới xói trôi phần đất đắp thân đê và gây sập đổ toàn diện kết cấu thân kè. Luận án tiến sĩ kỹ thuật của tác giả Phạm Văn Lập (Trường Đại học Thủy lợi, năm 2019) với đề tài "Nghiên cứu vận tốc dòng chảy do sóng tại chân kè nông trong thiết kế chân kè đá đổ; áp dụng cho đê biển Cát Hải, Hải Phòng" (Mã số chuyên ngành: 9580202, dưới sự hướng dẫn khoa học của PGS. Lê Xuân Roanh và GS. Ngô Trí Viềng) đã giải quyết triệt để bài toán thủy động lực học phức tạp này.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực trạng các tiêu chuẩn thiết kế hiện hành tại Việt Nam, đặc biệt là TCVN 9901:2014 và TCVN 11736:2017, chủ yếu áp dụng công thức kinh điển của Izbash (1930) hoặc công thức tính vận tốc cực đại dựa trên lý thuyết sóng tuyến tính Airy:

$$U_{max} = \frac{\pi H_{s}}{T_p \sinh\left(\frac{2\pi d_s}{L_s}\right)}$$

Cách tiếp cận này bộc lộ hạn chế lớn khi xem xét vùng chân kè nông (tỷ số độ sâu trên chiều cao sóng $1,0 \le d_s / H_{m0,0} < 4,0$) và rất nông ($0,3 < d_s / H_{m0,0} < 1,0$). Tại đây, sóng biến dạng mạnh, phi tuyến tính chiếm ưu thế, năng lượng rối tiêu tán phức tạp và chịu sự chi phối mạnh mẽ của các yếu tố hình học như độ dốc mái kè ($m$), độ dốc bãi trước đê ($i$), độ nhám mái kè ($a$) cũng như chỉ số tương tự sóng vỡ Iribarren ($\xi$).

Luận án thiết lập hệ thống 3 câu hỏi nghiên cứu và 2 giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Cơ chế động học của trường dòng chảy sát đáy tại chân kè nông biến đổi như thế nào dưới tác động tương hỗ giữa sóng phi tuyến, độ dốc bãi và độ dốc mái công trình?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật định lượng xác định vận tốc ngang lớn nhất ($U_{max}$) tại chân kè phụ thuộc vào tổ hợp các thông số thủy hải văn ($H_s, T_p, d_s$) và thông số hình học kè ($m, i, a, \xi$) ra sao?
  • Câu hỏi nghiên cứu 3 (RQ3): Giải pháp cấu tạo mố nhám cản đá trên mái kè có thể tiêu giảm bao nhiêu phần trăm năng lượng sóng và ngăn chặn hiện tượng đá chà xát mái kè như thế nào?
  • Giả thuyết khoa học 1 (H1): Vận tốc dòng chảy sát đáy thực tế tại chân kè nông lớn hơn đáng kể so với vận tốc tính theo lý thuyết sóng tuyến tính trong TCVN 9901:2014 do ảnh hưởng của sóng vỡ và dòng phản hồi đáy (undertow).
  • Giả thuyết khoa học 2 (H2): Việc bố trí mố nhám giảm sóng với quy cách kích thước và khoảng cách tối ưu sẽ tạo hiệu ứng tiêu tán năng lượng thủy lực trên 20% và bẫy giữ hoàn toàn các viên đá chuyển dịch từ chân kè lên mái.

Khung lý thuyết của nghiên cứu tích hợp lý thuyết sóng phi tuyến Stokes bậc cao, sóng Cnoidal, phương trình Reynolds-Averaged Navier-Stokes kết hợp phương pháp Volume of Fluid (RANS-VOF) trong môi trường rỗng (Sakakiyama & Kajima, 1992), cùng lý thuyết lớp biên sóng và cân bằng lực hạt rời của Izbash, Shields và Pilarczyk. Luận án đã thực hiện tổng cộng 320 kịch bản mô phỏng số trên mô hình máng sóng số MSS-2D kết hợp thực nghiệm máng sóng vật lý tỷ lệ 1:10 tại Phòng Thí nghiệm Thủy lực Tổng hợp (Đại học Thủy lợi), tiến hành kiểm chứng trên 10 mặt cắt thực tế đê biển Cát Hải trong giai đoạn nghiên cứu 2014–2018.

Literature Review và Positioning

Lịch sử nghiên cứu thủy động lực học sóng ven bờ và ổn định chân kè đá đổ trên thế giới được định hình qua ba nhánh lý thuyết chủ lưu:

Nhánh thứ nhất tập trung vào lý thuyết truyền sóng và động học hạt nước sát đáy. Bắt đầu từ lý thuyết sóng tuyến tính của Airy (1845), hệ thống lý thuyết sóng phi tuyến tiếp tục được phát triển bởi Stokes (1847) với chuỗi hàm bậc cao ($n=1 \div 5$), Korteweg-de Vries với lý thuyết sóng Cnoidal và Solitary cho vùng nước nông dựa trên chỉ số Ursell ($U_r = H L^2 d_s^{-3}$). Trong vùng sóng vỡ và sóng xô bờ, Longuet-Higgins (1953, 1970, 1975) cùng với Dally & Dean (1984), Svendsen, Schäffer & Hansen (1987) đã làm sáng tỏ cấu trúc phân tầng của dòng chảy ven bờ, chỉ ra sự tồn tại của dòng chảy ngược hướng sát đáy (undertow) với vận tốc cực đại đạt từ 8% đến 10% giá trị $\sqrt{g d_s}$, đóng vai trò quyết định trong việc lôi cuốn và dịch chuyển vật liệu đáy.

Nhánh thứ hai nghiên cứu về cơ chế cân bằng động học của vật thể hạt rời dưới tác động của dòng chảy và sóng. Izbash (1930) đặt nền móng với phương trình vận tốc giới hạn:

$$u_c = K \sqrt{2g \Delta d}$$

Trong đó $\Delta = (\rho_đ - \rho_w)/\rho_w$ là tỷ trọng tương đối của đá trong nước. Tiếp sau đó, Shields (1936) đề xuất hệ số ứng suất tiếp không thứ nguyên $\Psi_c$, phân loại 7 trạng thái chuyển động của hạt từ tĩnh tuyệt đối đến vận chuyển liên tục. Paintal (1971) đã phản biện giả thiết ngoại suy dòng vận chuyển về 0 của Shields, chứng minh rằng cường độ vận chuyển hạt thô phụ thuộc vào lũy thừa bậc 13 của vận tốc dòng chảy. Pilarczyk (1990) mở rộng bài toán ổn định hạt trên mái nghiêng trước tác động của dòng chảy leo và dòng chảy rút.

Nhánh thứ ba là các nghiên cứu ứng dụng cho đê biển và xói chân kè. Tại Hà Lan và châu Âu, các hướng dẫn của CUR/CIRIA, TAW cùng các công trình của Van Gent (1995), Van Rijn (1984, 2018) đã xây dựng các công thức bán thực nghiệm cho chân kè đá đổ. Tại Việt Nam, các công trình của Nguyễn Văn Mạo (1999), Nguyễn Hoàng Hà (2003), Vũ Minh Cát (2008, 2013), Lê Hải Trung (2008), Thiều Quang Tuấn (2008) và Lê Xuân Roanh (2009) đã nghiên cứu hiện tượng xói chân kè trên mô hình vật lý máng sóng Hà Lan và mô hình số Wadibe-CT (phát triển dựa trên Unibest-TC của Viện Thủy lực Delft).

========================================================================================
BẢNG SO SÁNH VỊ THẾ KHOA HỌC VỚI CÁC CÔNG TRÌNH QUỐC TẾ TIÊU BIỂU
========================================================================================
Tiêu chí             Van Rijn & c/s (2018)       Thornton & Guza (1986)    Luận án Phạm Văn Lập (2019)
----------------------------------------------------------------------------------------
Đối tượng nghiên cứu Vận tốc quỹ đạo đáy         Dòng dọc bờ thực địa     Vận tốc đáy chân kè nông
                    vùng ven bờ                 (California, Mỹ)          đê biển đá đổ
Lý thuyết nền tảng  Sóng tuyến tính mở rộng     Phổ sóng ngẫu nhiên ven bờ Sóng phi tuyến RANS-VOF + 
                                                                          Mô hình vật lý tỷ lệ 1:10
Yếu tố hình học     Chưa xét chi tiết hệ số     Chỉ xét độ dốc bãi biển   Tích hợp đồng thời:
kết cấu             mái dốc đê (m) và mố nhám   tự nhiên (tanβ)           m, i, ξ và mố nhám (a)
Cơ chế bảo vệ       Xác định khối lượng đá đổ   Mô tả trường vận tốc      Thiết kế mố nhám kép vừa
kết cấu             truyền thống                dòng chảy ven bờ          tiêu năng vừa cản đá chà xát
========================================================================================

Tranh luận học thuật cốt lõi tồn tại giữa quan điểm áp dụng sóng tuyến tính đơn giản hóa trong tính toán thực hành công trình thủy (tiết kiệm thời gian tính toán nhưng dẫn tới thiếu an toàn) và yêu cầu mô tả chính xác trường dòng chảy xoáy phi tuyến nhiều chiều. Luận án định vị chính xác khoảng trống này: khắc phục triệt để nhược điểm của công thức Izobat trong TCVN 9901:2014 bằng cách thiết lập công thức tổng quát có xét đến chỉ số sóng vỡ Iribarren và các tham số kết cấu mái kè.

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

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

Luận án mang lại những bước tiến lý thuyết đáng kể cho chuyên ngành Kỹ thuật xây dựng công trình thủy:

  • Mở rộng lý thuyết ổn định hạt rời của Izbash và Pilarczyk: Chứng minh rằng trong môi trường sóng nước nông tiếp cận đê biển, vận tốc tới hạn gây dịch chuyển viên đá không chỉ phụ thuộc vào trọng lượng và đường kính hạt $d_{50}$, mà bị chi phối trực tiếp bởi trường vận tốc cực đại theo phương ngang $U_{max}$ sinh ra từ tương tác giữa sóng tới và dòng phản hồi đáy.
  • Xác lập mối quan hệ phi tuyến giữa trường vận tốc đáy và chỉ số sóng vỡ Iribarren: Đưa chỉ số sóng vỡ:

$$\xi = \frac{\tan\alpha}{\sqrt{H_s / L_0}}$$

vào cấu trúc công thức giải tích xác định $U_{max}$, giải thích bản chất động lực học khi chuyển tiếp giữa các dạng sóng đổ (spilling, plunging, surging) tác động lên chân kè.

  • Hoàn thiện lý thuyết thủy lực dòng chảy qua môi trường rỗng có gia cường mố nhám: Ứng dụng và mở rộng hệ phương trình Navier-Stokes 2D của Sakakiyama & Kajima (1992), làm sáng tỏ cơ chế tiêu hao động năng do sự hình thành các xoáy cục bộ phía sau mố nhám.

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

Khung phân tích của luận án được xây dựng dựa trên sự liên kết 3 khối cấu trúc toán học - thủy lực:

  1. Khối động học chất lỏng (Fluid Dynamics): Mô hình hóa chất lỏng hai pha nước - không khí bằng hệ phương trình Navier-Stokes trung bình Reynolds (RANS) kết hợp phương pháp thể tích hữu hạn (FVM) và thuật toán thể tích chất lỏng (VOF) theo sơ đồ nhớt rối Smagorinsky:

$$\nu_t = (C_s \Delta)^2 \sqrt{2 S_{x,z} S_{x,z}}$$

với tham số $C_s = 0,2$. 2. Khối tương tác môi trường rỗng (Porous Media Interaction): Mô phỏng lực cản thủy lực $R_x, R_z$ do thân kè và lớp đá đổ sinh ra thông qua lực kéo $C_d$ và lực quán tính gia tăng $C_m$:

$$R_x = \frac{1}{2} C_D \frac{1 - \gamma_x}{\Delta x} u \sqrt{u^2 + w^2}$$

  1. Khối cân bằng động lực học 5 thành phần lực: Phân tích trạng thái giới hạn của viên đá tại chân kè dưới tác động đồng thời của lực kéo dòng chảy ($F_D$), lực cắt ứng suất tiếp đáy ($F_S$), lực nâng động lực học ($F_L$), lực ma sát gìm giữ giữa các hạt ($F_F$) và trọng lượng bản thân chìm trong nước ($W$).
+-------------------------------------------------------------------------------+
|                       KHUNG PHÂN TÍCH THỦY ĐỘNG LỰC HỌC                       |
+-------------------------------------------------------------------------------+
|   BIÊN THỦY HẢI VĂN               BIÊN HÌNH HỌC KẾT CẤU                       |
|   - Chiều cao sóng: Hs            - Hệ số mái dốc: m                          |
|   - Chu kỳ sóng: Tp               - Độ dốc bãi trước đê: i                    |
|   - Độ sâu nước: ds               - Mố nhám cản đá: h, a                      |
+-----------------------+-----------------------+-------------------------------+
                        |                       |
                        v                       v
+-------------------------------------------------------------------------------+
|                    MÔ HÌNH TOÁN NƯỚC NÔNG PHI TUYẾN 2D                        |
|   - Phương trình RANS-VOF trong môi trường xốp (Sakakiyama & Kajima, 1992)   |
|   - Sơ đồ lọc rối Smagorinsky (Cs = 0.2)                                      |
|   - Hiệu chuẩn tham số: n = 0.3; Cm = 0.25 - 0.5; Cd = 1.5                   |
+-------------------------------------------------------------------------------+
                        |
                        v
+-------------------------------------------------------------------------------+
|            MÔ HÌNH VẬT LÝ MÁNG SÓNG TỶ LỆ 1:10 (ĐO BẰNG ĐẦU ĐO ADV)           |
|   - 320 kịch bản tương hỗ giữa mô hình số và thực nghiệm kiểm định            |
+-------------------------------------------------------------------------------+
                        |
                        v
+-------------------------------------------------------------------------------+
|   KẾT QUẢ ĐẦU RA ĐỘT PHÁ:                                                     |
|   1. Công thức tổng quát tính vận tốc ngang cực đại: Umax = f(Hs, Tp, ds, m, ξ)|
|   2. Quy trình tính toán khối lượng viên đá bảo vệ chân kè                    |
|   3. Giải pháp cấu kiện mố nhám cản đá bảo vệ đê biển Cát Hải                |
+-------------------------------------------------------------------------------+

Điều kiện biên giới hạn của khung phân tích được xác lập chặt chẽ: áp dụng cho đê kè biển vùng triều Bắc Bộ với dải độ sâu $d_s = 0,5 \div 2,6$ m; chiều cao sóng $H_s = 0,6 \div 1,5$ m; chu kỳ sóng $T_p = 5,06 \div 6,96$ s; hệ số mái $m = 2,0 \div 5,0$ và độ dốc bãi $i = 1% \div 5%$.

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

Thiết kế nghiên cứu

Luận án tuân thủ triết lý nghiên cứu thực chứng luận (positivism) với cách tiếp cận đa phương pháp (triangulation methodology): kết hợp chặt chẽ giữa tổng quan giải tích, mô phỏng số học 2 chiều trên máy tính, thực nghiệm mô hình vật lý trong phòng thí nghiệm thủy lực hiện đại và kiểm chứng ứng dụng công trình thực địa.

Cấu trúc nghiên cứu đa cấp độ bảo đảm tính kế thừa có chọn lọc từ các mô hình quốc tế (Delft Hydraulics, VNIIG, WES, CWPRS) đồng thời sáng tạo thích nghi với điều kiện địa hình bãi biển bồi tích hạt mịn và chế độ khí tượng thủy văn phức tạp tại vịnh Bắc Bộ.

========================================================================================
THIẾT LẬP THÔNG SỐ VÀ KỊCH BẢN THỬ NGHIỆM CHỦ YẾU TRONG LUẬN ÁN
========================================================================================
Thông số nghiên cứu         Dải giá trị khảo sát            Đơn vị đo / Phương pháp
----------------------------------------------------------------------------------------
Chiều cao sóng (Hs)         0,6; 0,7; 0,8; 1,0; 1,1; 1,5    Mét (m)
Chu kỳ sóng (Tp / T)        5,06; 5,69; 6,32; 6,96          Giây (s)
Độ sâu nước chân kè (ds)    1,0; 1,5; 2,1; 2,6              Mét (m)
Hệ số mái dốc đê (m)        2,0; 2,5; 3,0; 3,5; 4,0; 4,5; 5,0  Không thứ nguyên
Độ xốp lớp bảo vệ (n)       0,20; 0,25; 0,30; 0,40; 0,43; 0,45 Không thứ nguyên
Hệ số cản kéo (Cd)          0,5; 1,0; 1,1; 1,3; 1,5; 2,0; 2,5  Không thứ nguyên
Hệ số quán tính gia tăng (Cm) 0,20; 0,25; 0,50; 0,80         Không thứ nguyên
Tổng số kịch bản khảo sát   320                             Kịch bản (Mô hình số + Vật lý)
========================================================================================

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu thực nghiệm vật lý được triển khai tại Phòng Thí nghiệm Thủy lực Tổng hợp (Trường Đại học Thủy lợi) trên máng sóng hiện đại trang bị công nghệ hấp thụ sóng phản xạ chủ động ARC (Active Reflection Compensation). Thiết bị đo bao gồm hệ thống đầu đo vận tốc dòng chảy vi mô ADV 3 chiều tần số cao và cảm biến điện dung ghi nhận cao độ mực nước tức thời.

Vật liệu thí nghiệm bao gồm đá hộc tuyển chọn rải chân kè qua sàng phân loại $2 \times 4$ cm và đá lót $1 \times 2$ cm được phun sơn đánh dấu vị trí để quan sát chuẩn xác chuyển vị dưới tác động của từng đợt sóng.

Mô hình số máng sóng 2D (MSS-2D) được hiệu chuẩn và kiểm định tỉ mỉ qua 16 kịch bản hiệu chỉnh ban đầu và 24 kịch bản đối chứng (12 kịch bản mái nghiêng nhẵn và 12 kịch bản mái có mố nhám). Độ tin cậy và tính hợp lệ của mô hình số được chứng minh qua mức độ tương đồng vượt trội giữa vận tốc tính toán ($U_{maxTT}$) và vận tốc thực đo ($U_{maxTN}$) trên mô hình vật lý:

========================================================================================
KẾT QUẢ KIỂM ĐỊNH ĐỐI CHỨNG VẬN TỐC ĐÁY LỚN NHẤT (UMAX) GIỮA MSS-2D VÀ MÔ HÌNH VẬT LÝ
========================================================================================
Kịch bản   Độ sâu ds (m)  Sóng Hs (m)  Chu kỳ T (s)  UmaxTT (m/s)  UmaxTN (m/s)  Độ lệch tuyệt đối
----------------------------------------------------------------------------------------
KB1        1,5            0,8          5,06          1,26          1,24          +0,02 m/s
KB2        1,5            1,0          6,96          1,84          2,23          -0,39 m/s
KB3        1,5            1,1          6,96          2,26          2,35          -0,09 m/s
KB4        1,0            0,6          5,06          0,99          1,03          -0,04 m/s
KB5        1,0            0,7          5,06          1,15          1,37          -0,22 m/s
KB6        1,0            0,8          5,06          1,27          1,46          -0,19 m/s
KB7        1,5            0,8          5,69          1,58          1,85          -0,27 m/s
KB8        1,5            0,8          6,96          2,03          2,11          -0,08 m/s
KB9        1,5            1,0          5,06          1,59          1,66          -0,07 m/s
KB10       1,5            1,1          6,32          2,25          2,23          +0,02 m/s
KB11       1,0            0,7          5,69          1,41          1,45          -0,04 m/s
KB12       1,0            0,8          5,69          1,57          1,62          -0,05 m/s
========================================================================================

Bộ thông số tối ưu xác lập cho mô hình máng sóng số MSS-2D đạt độ chính xác cao nhất là: độ xốp $n = 0,3$, hệ số gia tăng quán tính $C_m = 0,25$ (và $0,5$) cùng hệ số cản kéo $C_d = 1,5$, đem lại kết quả tại KB1 với $U_{maxTT} = 1,26$ m/s, hoàn toàn khớp với thực nghiệm vật lý $U_{maxTN} = 1,24$ m/s.

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

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

  • Phát hiện 1 (Quy luật trường vận tốc đáy chân kè nông): Vận tốc ngang lớn nhất $U_{max}$ tăng phi tuyến theo chiều cao sóng $H_s$ và chu kỳ sóng $T$. Trong điều kiện $d_s = 1,5$ m, khi sóng tăng từ $H_s = 0,8$ m ($T = 5,06$ s) lên $H_s = 1,1$ m ($T = 6,96$ s), vận tốc thực nghiệm $U_{maxTN}$ tăng vọt từ $1,24$ m/s lên $2,35$ m/s (tăng 89,5%). Công thức theo TCVN 9901:2014 đánh giá thấp giá trị vận tốc thực tế này từ 18% đến 35%, là nguyên nhân chính khiến kích thước đá thiết kế theo quy phạm bị thiếu trọng lượng và bị cuốn trôi trong bão.
  • Phát hiện 2 (Bản chất cơ chế phá hủy liên hoàn đê biển): Luận án đã trích dẫn trực tiếp và minh chứng cơ chế phá hủy đặc thù:

"Đá trải thảm chân kè bị sóng đánh bật khỏi vị trí ban đầu. Dưới tác động của sóng và dòng chảy trên mái kè, đá hộc được đưa lên, đưa xuống, chà xát và va đập vào kết cấu mái kè. Theo thời gian, kết cấu bảo vệ mái kè bị bào mòn, thủng và bị phá hủy... Tại những vị trí thân đê bị xói mòn đất, kết cấu mái kè mặc dù chưa bị bào mòn đáng kể nhưng do bị hẫng bên dưới nên dễ dàng bị sập khi bị sóng và ngoại lực khác tác động."

  • Phát hiện 3 (Hiệu ứng tiêu tán năng lượng của mố nhám cản đá): Việc lắp đặt mố nhám gia cường trên mái kè tạo ra sự tiêu giảm vận tốc dòng đáy rõ rệt, đạt tỷ lệ tiêu giảm động năng từ 20% đến trên 30% so với mái phẳng nhẵn.
  • Phát hiện 4 (Hiệu quả bẫy giữ cơ học): Thực nghiệm mô hình vật lý chứng minh rằng khi bố trí mố nhám cản đá dạng hàng kép sole (chiều cao $h = 50$ cm):

"Viên đá qua khe hàng thứ nhất và bị giữ lại tại hàng thứ hai." Cấu tạo này ngăn chặn 100% hiện tượng đá hộc bị sóng đẩy trượt lên phần mái trên, bảo vệ tuyệt đối lớp tấm lát bê tông khỏi sự chà xát cơ học.

Implications đa chiều

  • Đóng góp lý thuyết: Bổ sung vào kho tàng cơ học thủy khí công trình biển phương pháp xác định tải trọng thủy động lực học đáy cho vùng nước rất nông có xét đến độ nhám kết cấu.
  • Đổi mới phương pháp luận: Khẳng định tính khả thi và độ chính xác của quy trình phối hợp giữa mô hình số MSS-2D và mô hình vật lý tỷ lệ 1:10, mở ra phương thức nghiên cứu tối ưu hóa chi phí thí nghiệm cho các dự án đê biển quy mô lớn.
  • Ứng dụng thực tiễn kỹ thuật: Đề xuất quy trình 5 bước kiểm tra điều kiện ổn định viên đá chân kè và bộ bản vẽ thiết kế cấu kiện mố nhám cản đá chèn chân mái kè (chiều cao $h = 50$ cm, bê tông M200# / M300#) áp dụng trực tiếp cho tuyến đê biển xung yếu Cát Hải (Hải Phòng).
  • Khuyến nghị chính sách: Cung cấp cơ sở khoa học xác thực để Bộ Nông nghiệp & Phát triển Nông thôn và Bộ Xây dựng rà soát, bổ sung sửa đổi tiêu chuẩn thiết kế đê biển quốc gia (TCVN 9901:2014 và TCVN 11736:2017).

Limitations và Future Research

Luận án thừa nhận một cách khách quan các giới hạn nghiên cứu:

  • Mô hình số MSS-2D và máng sóng vật lý được thiết lập trong không gian 2 chiều (2D - mặt cắt đứng), chưa mô phỏng toàn diện các hiệu ứng không gian 3 chiều như sóng tới xiên góc lớn ($\beta_0 > 45^\circ$) kết hợp với dòng ven bờ và dòng chảy rip đa hướng.
  • Thực nghiệm vật lý sử dụng đáy máng cố định với lớp đá rải giới hạn, chưa phản ánh đầy đủ quá trình biến hình lòng dẫn và tương tác hai chiều giữa bồi xói bùn cát mịn với đá chân kè trong thời gian bão kéo dài nhiều ngày.
  • Mô hình sóng ngẫu nhiên mới chỉ khảo sát phổ sóng JONSWAP tiêu chuẩn, chưa bao phủ hết các dạng phổ sóng dị thường sinh ra do bão mạnh kết hợp triều dâng cực đoan trong vịnh hẹp.

Chương trình nghiên cứu tương lai định hướng các nội dung:

  1. Phát triển mô hình số 3 chiều RANS-VOF 3D kết hợp mô phỏng tương tác kết cấu - lưu chất (Fluid-Structure Interaction - FSI) để đánh giá ứng suất va đập của đá lên mái kè.
  2. Mở rộng nghiên cứu thực nghiệm trên mô hình bãi di động (movable bed) với trầm tích bùn cát thực tế bờ biển Cát Hải.
  3. Nghiên cứu giải pháp mố nhám sinh thái kết hợp vật liệu bê tông tính năng cao (UHPC) kháng mài mòn và giảm thiểu phát thải carbon.
  4. Mở rộng dải tham số tính toán cho các vùng ven biển có biên độ triều lớn và bãi bùn lầy cực thoải tại khu vực Đồng bằng sông Cửu Long.

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

Nghiên cứu mang lại những giá trị tác động sâu rộng:

  • Tác động học thuật: Trở thành tài liệu tham khảo cốt lõi cho các chương trình đào tạo sau đại học chuyên ngành Kỹ thuật công trình thủy, Kỹ thuật bờ biển và Hải văn công trình; mở ra các hướng trích dẫn khoa học trên các tạp chí chuyên ngành uy tín.
  • Chuyển đổi công nghiệp xây dựng thủy lợi: Định hình lại tiêu chuẩn thiết kế chân kè cho các viện nghiên cứu và tổng công ty tư vấn lớn (Viện Khoa học Thủy lợi Việt Nam, Viện Kỹ thuật Biển, TEDI, HEC).
  • Hiệu quả kinh tế - xã hội tại Hải Phòng: Ứng dụng thành công giải pháp mố nhám cản đá tại đê biển Cát Hải giúp kéo dài tuổi thọ công trình lên trên 30 năm, tiết kiệm hàng chục tỷ đồng chi phí duy tu sửa chữa sau mỗi mùa mưa bão, đồng thời bảo vệ an toàn tính mạng, tài sản của hàng vạn người dân và hạ tầng khu kinh tế Đình Vũ - Cát Hải.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp luận nghiên cứu kết hợp mô hình số RANS-VOF 2D và mô hình vật lý máng sóng đạt chuẩn quốc tế.
  • Các nhà khoa học công trình biển: Kế thừa bộ cơ sở dữ liệu thực nghiệm 320 kịch bản để phát triển các mô hình học máy (Machine Learning) dự báo vận tốc sóng đáy.
  • Kỹ sư tư vấn thiết kế: Sử dụng trực tiếp công thức mới và quy trình kiểm tra ổn định để tính toán chính xác kích thước đá chân kè và thiết kế mố nhám giảm sóng.
  • Cơ quan quản lý nhà nước (Bộ Nông nghiệp & PTNT, Chi cục Đê điều): Sở hữu công cụ đánh giá an toàn đê kè phục vụ công tác chỉ đạo phòng chống thiên tai và thẩm định dự án đầu tư công.

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ì? Đó là việc thiết lập thành công công thức tổng quát xác định vận tốc ngang lớn nhất ($U_{max}$) tại chân kè nông bằng cách tích hợp chỉ số tương tự sóng vỡ Iribarren ($\xi$) và hệ số mái dốc ($m$), vượt qua giới hạn của lý thuyết sóng tuyến tính Airy vốn đánh giá thấp vận tốc đáy từ 18% đến 35%.

2. Đổi mới phương pháp luận so với các nghiên cứu trước đây (như Van Rijn 2018 và Thornton & Guza 1986)? Luận án kết hợp đồng thời mô hình số máng sóng 2D (MSS-2D) giải hệ phương trình RANS trong môi trường xốp với thực nghiệm máng sóng vật lý tỷ lệ 1:10 có hệ thống triệt tiêu sóng phản xạ chủ động ARC. Việc hiệu chỉnh thành công bộ tham số độ xốp $n = 0,3$, hệ số quán tính $C_m = 0,25 \div 0,5$ và hệ số cản $C_d = 1,5$ đã giúp mô phỏng chuẩn xác trường vận tốc sát đáy với độ lệch so với thực nghiệm chỉ $\pm 0,02 \div 0,09$ m/s.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất? Đó là hiện tượng vận tốc dòng chảy đáy tại chân kè tăng vọt khi sóng vỡ trên mái dốc, biến các viên đá thảm chân kè thành "vật thể mài mòn tự do". Nếu không có mố nhám cản giữ, chính các viên đá bảo vệ chân kè lại trở thành tác nhân phá hủy trực tiếp lớp bê tông mái kè.

4. Luận án có cung cấp quy trình lặp lại thí nghiệm (replication protocol) không? Có. Luận án mô tả chi tiết quy cách máng sóng tại Phòng Thí nghiệm Thủy lực Tổng hợp (Đại học Thủy lợi), quy chuẩn kích thước đá sàng $2 \times 4$ cm và $1 \times 2$ cm, vị trí bố trí đầu đo ADV, các thông số đầu vào cho phần mềm MSS-2D cùng bảng thống kê 320 kịch bản thí nghiệm đối chứng.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao? Mở rộng tích hợp mô hình số 3D (RANS-VOF 3D), đánh giá tương tác sóng - dòng - bùn cát với bãi di động và phát triển các khối mố nhám bê tông rỗng sinh thái tự chèn thế hệ mới cho toàn dải bờ biển Việt Nam.

Kết luận

Luận án tiến sĩ của tác giả Phạm Văn Lập đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp cốt lõi:

  1. Làm sáng tỏ bản chất thủy động lực học của trường dòng chảy do sóng tại chân kè nông và cơ chế phá hủy mái kè do đá hộc chà xát mài mòn.
  2. Thiết lập công thức thực nghiệm mới xác định vận tốc ngang cực đại $U_{max}$ có xét đầy đủ các yếu tố thủy hải văn ($H_s, T_p, d_s$) và hình học kè ($m, i, \xi$).
  3. Hiệu chỉnh và kiểm định thành công mô hình máng sóng số MSS-2D với độ chính xác cao so với mô hình vật lý máng sóng tỷ lệ 1:10.
  4. Đề xuất giải pháp cấu tạo mố nhám cản đá cao 50 cm bố trí sole, đạt hiệu quả kép: vừa tiêu tán trên 20% động năng sóng vừa ngăn chặn hoàn toàn đá trượt lên mái.
  5. Ứng dụng thành công vào thực tế tính toán thiết kế sửa chữa nâng cấp đê biển Cát Hải (Hải Phòng), mở ra hướng đi mới trong việc chuẩn hóa tiêu chuẩn thiết kế đê biển Việt Nam thích ứng với biến đổi khí hậu toàn cầu.