Tổng quan về luận án

Hệ thống đê biển khu vực Bắc Bộ chiếm hơn 65% trong tổng số 2.600 km đê biển của Việt Nam, đóng vai trò sống còn trong việc bảo vệ dân cư và hạ tầng trước các cơn bão nhiệt đới cực đoan. Tuy nhiên, đặc thù đê biển Bắc Bộ có cao trình đỉnh phổ biến từ 5,0 đến 5,5 m, mái nghiêng ngoài dốc $m = 3 \div 4$, khiến hiện tượng nước dềnh bão kết hợp triều cường thường xuyên đẩy sóng lớn áp sát chân công trình, làm phát sinh đồng thời sóng tràn qua đê và tương tác sóng phản xạ cục bộ phức tạp. Các quan trắc thực địa xác nhận nguy cơ xói lở chân đê diễn ra rất nhanh và nghiêm trọng; điển hình như ghi nhận: "Đo đạc thực tế mặt cắt ngang trước và sau bão Damrey 2005 tại Nam Định của Viện Cơ học và hình ảnh bãi trước chân đê bị hạ thấp, độ sâu xói chân đê sau cơn bão lên tới 1m, gây nên mất ổn định cho toàn bộ công trình."

Khoảng trống học thuật cốt lõi (research gap) nằm ở chỗ các mô hình thủy động lực học ven bờ hiện hành (như UNIBEST-TC, DUROSTA, SBEACH) chủ yếu giả định biên công trình thẳng đứng hoặc biên bãi tự nhiên mà bỏ qua quy luật biến đổi không gian của hệ số phản xạ $K_r(x)$ từ kết cấu đê mái nghiêng có tường đỉnh, cũng như chưa định lượng được ảnh hưởng kết hợp của sóng phản xạ và sóng tràn đến cấu trúc dòng phản hồi (undertow) và nồng độ bùn cát lơ lửng.

Nghiên cứu giải quyết 4 câu hỏi khoa học trọng tâm:

  1. Quy luật biến đổi không gian của hệ số phản xạ sóng $K_r$ trên mặt cắt ngang bãi trước chân đê mái nghiêng diễn ra như thế nào?
  2. Sóng phản xạ từ mái đê và tường đỉnh biến đổi cấu trúc vận tốc dòng phản hồi $u(z)$ và vận tốc trung bình độ sâu $u_0$ theo cơ chế nào?
  3. Sự tương tác giữa sóng phản xạ và chuyển động rối chi phối nồng độ bùn cát lơ lửng $C(z)$ và gradient vận chuyển bùn cát ngang bờ $S(x)$ ra sao?
  4. Khả năng ứng dụng và độ chính xác của mô hình toán cập nhật khi mô phỏng hố xói chân đê trong điều kiện bão thực tế tại Nam Định đạt mức độ nào?

Khung lý thuyết của công trình tích hợp lý thuyết bức xạ sóng của Svendsen (1984), mô hình dòng phản hồi logarit - tuyến tính của Steetzel (1993), phương pháp phân tách phổ sóng mảng thẳng đứng của Sheremet (2002) và hàm tương quan không gian của Klopman & Van der Meer (1999). Đóng góp đột phá của công trình là thiết lập thành công mô hình toán Wadibe-TC cải tiến, tích hợp hàm suy giảm phản xạ cục bộ $F_x(x/L)$, cho phép mô phỏng chính xác trắc diện xói chân đê với sai số định lượng dưới 8% so với thực nghiệm và tái hiện trung thực hố xói sâu 1,0 m tại đê biển Thịnh Long (Nam Định).

Literature Review và Positioning

Các nghiên cứu kinh điển về thủy động lực học vùng sóng vỡ xác lập rằng dòng phản hồi là động lực chính chi phối vận chuyển bùn cát ngang bờ ra khơi trong bão. Svendsen (1984) và Stive & Wind (1986) chứng minh dòng phản hồi hình thành do sự mất cân bằng giữa thông lượng động lượng biến đổi theo chiều sâu (gradient ứng suất bức xạ $S_{xx}$) và gradient áp suất thủy tĩnh đồng nhất do nước dềnh tạo bởi sóng. Steetzel (1993) phát triển nghiệm giải tích trắc diện vận tốc $u(z)$ dạng kết hợp tuyến tính và logarit, trong khi Battjes & Janssen (1978, 2008) hoàn thiện mô hình tiêu tán năng lượng sóng vỡ ngẫu nhiên $D_w$.

Tuy nhiên, tồn tại một cuộc tranh luận học thuật sâu sắc giữa hai trường phái:

  1. Trường phái thứ nhất (Neshaei, 2009; Goda & Suzuki, 1976) lập luận rằng sự gia tăng hệ số phản xạ $K_r$ làm suy giảm thông lượng năng lượng sóng tới hướng bờ, từ đó làm giảm độ lớn của dòng phản hồi tổng cộng theo phương trình bảo toàn khối lượng, dẫn đến giảm xói lở đáy.
  2. Trường phái thứ hai (Fredsøe & Deigaard, 2012; Zanuttigh, 2008) chỉ ra rằng sóng phản xạ giao thoa với sóng tới làm khuấy động mạnh mẽ lớp biên đáy, gia tăng ứng suất tiếp rối $\tau(z)$ và hệ số xáo trộn bùn cát $\varepsilon_s(z)$, khiến nồng độ bùn cát lơ lửng $C(z)$ tăng vọt theo phương đứng.

Công trình này định vị chính xác ở điểm giao thoa của hai luồng quan điểm trên. Khi so sánh với nghiên cứu của Neshaei (2009) vốn chỉ thử nghiệm sóng đều trong máng nhỏ với 3 kịch bản bãi tự nhiên, nghiên cứu này mở rộng ra phổ sóng ngẫu nhiên JONSWAP tương tác với đê mái nghiêng $m = 4$ có tường đỉnh trong điều kiện có sóng tràn. Đồng thời, so với tổng hợp 4.000 thí nghiệm của Zanuttigh (2008) về hệ số phản xạ tổng thể, nghiên cứu đã tiến thêm một bước đột phá khi lượng hóa phân bố cục bộ $K_{r,x}(x/L)$ và hàm tương quan không gian $F_x(x/L)$, giải quyết triệt để sự thiếu hụt dữ liệu dòng chảy trắc diện mà các nghiên cứu trong nước trước đây (Vũ Minh Cát, 2008; Thiều Quang Tuấn, 2008; Phạm Văn Lập, 2019) chưa thực hiện được.

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

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

Công trình mở rộng trực tiếp lý thuyết động lượng sóng ven bờ của Svendsen (1984) và mô hình dòng phản hồi của Steetzel (1993) bằng cách bổ sung thành phần sóng phản xạ cục bộ vào hệ phương trình chủ đạo. Luận án trích dẫn rõ bản chất vật lý: "Dòng phản hồi là dòng chảy trung bình hướng ra biển, ở phía dưới mực nước chân sóng trong vùng sóng vỡ... được hình thành do sự mất cân bằng theo phương thẳng đứng giữa thông lượng động lượng biến đổi theo độ sâu và gradient áp suất đồng nhất theo độ sâu do sự dâng nước bởi sóng."

Mô hình lý thuyết đề xuất xác lập mối liên hệ giữa thông lượng khối lượng sóng tới $(m_w + m_r)_i$ và thông lượng sóng phản xạ $(m_w)_r$:

$$m_t = (m_w + m_r)_i - (m_w)_r$$

Đi kèm với phương trình phân bố vận tốc dòng chảy ngang bờ $u(z)$:

$$u(z) = u_0 + K_{lin} z + K_{log} \ln \left[ 1 + \frac{az}{\varepsilon_0} \right]$$

Mệnh đề lý thuyết mới khẳng định: Sóng phản xạ tuy làm giảm vận tốc dòng phản hồi trung bình $u_0$ từ 12% đến 18% ngay sát chân đê, nhưng lại làm tăng hệ số xáo trộn bùn cát $\varepsilon_s(z)$ lên 25 - 40%, dẫn đến tích số $u(z) \cdot C(z)$ tăng mạnh và tạo ra gradient vận chuyển bùn cát $\partial S / \partial x > 0$ cục bộ, thúc đẩy quá trình đào sâu hố xói chân đê vượt trội so với trường hợp không xét phản xạ.

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

Khung phân tích tích hợp đồng thời 4 trụ cột lý thuyết:

  • Lý thuyết bảo toàn động lượng và ứng suất bức xạ sóng (Longuet-Higgins & Stewart, 1964; Svendsen, 1984).
  • Cơ chế tiêu tán năng lượng cuộn sóng mặt (Wave Roller Theory của Fredsøe & Deigaard, 1992).
  • Lý thuyết khuếch tán bùn cát lơ lửng 1DV với nồng độ tham chiếu đáy $C_0$ (Steetzel, 1993; Van Rijn, 1993).
  • Phương pháp phân tách năng lượng sóng vào - ra mảng thẳng đứng (Vertical Array Analysis của Sheremet, 2002).

Phạm vi biên (boundary conditions) được xác định tường minh: áp dụng cho vùng sóng vỡ ven bờ có độ dốc bãi $1/100 \le \tan \beta \le 1/40$, kết cấu đê mái nghiêng $m = 3 \div 4$, chiều cao tường đỉnh $h_w \le 1,0$ m, thông số lưu không tương đối $-1,0 \le R_c/H_s \le 0,5$ và vùng ảnh hưởng phản xạ giới hạn trong khoảng cách không gian $x/L \le 0,4$.

+-------------------------------------------------------------------------+
|                  KHUNG PHAN TICH THUY DONG LUC WADIBE-TC                |
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
| 1. MODULE SONG (Battjes & Janssen 1978; Klopman & Van der Meer 1999)     |
|    - Cân bang nang luong song: d(c_g E)/dx + D_w + D_f = 0              |
|    - Phan tach song mang thang dung: K_r^2 = F^- / F^+                  |
|    - Phan bo song cuc bo: H_m0(x)/H_m0i(x) = [1 + K_r^2 + 2K_r F_x]^0.5  |
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
| 2. MODULE DONG CHAY (Svendsen 1984; Steetzel 1993)                      |
|    - Can bang dong luong: dS_xx/dx + rho*g*(d + eta_bar)*d(eta_bar)/dx = 0 |
|    - Thong luong khoi luong tinh: m_t = (m_w + m_r)_i - (m_w)_r          |
|    - Phan bo van toc undertow u(z) = u_0 + K_lin*z + K_log*ln[1+az/eps_0]|
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
| 3. MODULE VAN CHUYEN BUN CAT (Steetzel 1993; Van Rijn 1993)             |
|    - Khuech tan bun cat: dC/dt = omega_s*(dC/dz) + d/dz[eps_s*(dC/dz)]  |
|    - Nong do day C_0 va he so xao tron eps_s cap nhat song phan xa      |
|    - Suat chuyen cat ngang bo: S = int [u(z) * C(z)] dz                 |
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
| 4. MODULE BIEN DOI HINH THAI (Phuong trinh Exner)                       |
|    - Bao toan khoi luong cat day: (1 - n)*dz_b/dt = -dS/dx              |
|    - Du tinh vi tri va chieu sau ho xoi lon nhat z_max truoc chan de     |
+-------------------------------------------------------------------------+

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ triết lý thực chứng (positivism) kết hợp kiểm chứng thực nghiệm đa tầng (triangulation between physical experiments and numerical modeling). Thiết kế nghiên cứu đa cấp độ gồm:

  • Thí nghiệm mô hình vật lý thu nhỏ 2D trong máng sóng chuyên dụng tại Phòng Thí nghiệm Thủy lực Tổng hợp - Đại học Thủy lợi.
  • Mô hình toán số trị giải phương trình vi phân từng bước theo thời gian 1DV (Time-dependent Cross-shore).
  • Kiểm chứng mô hình hiện trường dựa trên số liệu thực địa bão Damrey (9/2005) tại đê biển Thịnh Long, Nam Định.

Mô hình vật lý tuân theo tiêu chuẩn tương tự Froude với tỷ lệ hình học $N_L$. Cát thí nghiệm sử dụng cát thạch anh tự nhiên có đường kính trung vị $d_{50} = 0,18$ mm, độ rỗng $n = 0,4$, khối lượng riêng hạt $\rho_s = 2.650\text{ kg/m}^3$.

Quy trình nghiên cứu rigorous

Hệ thống thí nghiệm bao gồm hai chuỗi mô hình:

  1. Mô hình lòng cứng (Fixed-bed): Đo đạc đồng thời trắc diện sóng bằng hệ thống đầu đo sóng điện dung (Wave Gauges - WG) và phân bố vận tốc tức thời $u(z,t)$ bằng máy đo vận tốc âm học Doppler đa điểm (Acoustic Doppler Velocimetry - ADV/CG).
  2. Mô hình lòng động (Mobile-bed): Đánh giá biến dạng đáy bãi và kích thước hố xói trước đê với hai cấu hình độ dốc thềm bãi $\tan \beta = 1/100$ và $1/40$, mái đê $m = 4$.

Các cấu hình kết cấu đê khảo sát:

  • Đê cao không tràn ($CW_0$).
  • Đê thấp không tường đỉnh ($TW_0$).
  • Đê thấp có tường đỉnh cao 10 cm trong mô hình ($TW_{10}$).

Các kịch bản thủy lực bão mô phỏng:

  • $D_{65}H_{15}T_{19}$ (độ sâu $D = 65$ cm, chiều cao sóng $H = 15$ cm, chu kỳ đỉnh phổ $T_p = 1,9$ s).
  • $D_{65}H_{17}T_{16}$ ($D = 65$ cm, $H = 17$ cm, $T_p = 1,6$ s).
  • $D_{70}H_{15}T_{19}$ ($D = 70$ cm, $H = 15$ cm, $T_p = 1,9$ s).
  • $D_{70}H_{19}T_{165}$ ($D = 70$ cm, $H = 19$ cm, $T_p = 1,65$ s).

Phương pháp tách sóng phản xạ áp dụng kỹ thuật mảng thẳng đứng của Sheremet (2002) thông qua tích phân phổ mật độ phương sai sóng $S_{\eta\eta}(f,x)$, phổ mật độ phương sai dòng $S_{uu}(f,x)$ và phổ hiệp phương sai $S_{\eta u}(f,x)$:

$$E^{\pm}(f,x) = \frac{1}{4} \left[ S_{\eta\eta}(f,x) + \frac{d}{g} S_{uu}(f,x) \pm 2\sqrt{\frac{d}{g}} S_{\eta u}(f,x) \right]$$

$$F^{\pm}(x) = \int_{f_L}^{f_U} E^{\pm}(f,x) \sqrt{gd} , df \quad \Longrightarrow \quad K_r^2(x) = \frac{F^-(x)}{F^+(x)}$$

Data và phân tích

Mô hình số Wadibe-TC (Wave, Dike and Beach – Time dependent Cross-shore) được nâng cấp toàn diện từ nền tảng UNIBEST-TC của Viện Thủy lực Delft (Delft Hydraulics). Mô hình tích hợp số liệu tái phân tích trường sóng biển Đông từ mô hình Wavewatch III kết hợp dữ liệu khí tượng bão ven bờ Bắc Bộ.

Các kiểm định độ nhạy và tính vững (robustness checks) được thực hiện với các bước lưới không gian $\Delta x = 0,05 \div 0,2$ m và bước thời gian $\Delta t = 0,01$ s, đảm bảo điều kiện ổn định Courant-Friedrichs-Lewy ($CFL < 0,5$). Chỉ số Nash-Sutcliffe Efficiency ($NSE$) giữa kết quả mô phỏng số và dữ liệu thực nghiệm đạt $NSE = 0,91 \div 0,95$ đối với chiều cao sóng $H_{m0}$ và $NSE = 0,87 \div 0,92$ đối với vận tốc dòng phản hồi $u(z)$.

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

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

Thứ nhất, hệ số phản xạ sóng $K_r$ không phân bố đồng đều mà biến đổi phi tuyến mạnh mẽ trong vùng nước sát đê. Nghiên cứu xác lập hàm suy giảm thực nghiệm $F_x(x/L)$ phụ thuộc vào khoảng cách tương đối $x/L$:

$$F_x(x/L) = \exp \left[ -a \left( \frac{x}{L} \right)^b \right]$$

với $a = 6,25$ và $b = 1,18$ ($R^2 = 0,94$). Vùng ảnh hưởng của sóng phản xạ từ đê mái nghiêng bị triệt tiêu hoàn toàn khi khoảng cách $x/L > 0,4$.

Thứ hai, sự xuất hiện của tường đỉnh làm gia tăng hệ số phản xạ tại chân công trình $K_{r,0}$ từ $0,18 \div 0,22$ (đối với đê thấp $TW_0$) lên $0,35 \div 0,42$ (đối với đê có tường đỉnh $TW_{10}$), tương ứng mức tăng từ 65% đến 90%.

Thứ ba, phát hiện hiện tượng đối nghịch thủy động lực học: Sóng phản xạ làm giảm vận tốc dòng phản hồi trung bình $u_0$ khoảng 15%, nhưng làm gia tăng nồng độ bùn cát đáy $C_0$ lên 1,4 - 1,8 lần do sóng phản xạ giao thoa tạo các xoáy rối mạnh sát đáy. Kết quả là thể tích hố xói lớn nhất tăng 22% đến 35% so với kịch bản không tính đến phản xạ sóng.

Thứ tư, ảnh hưởng của sóng tràn qua đê ($q$): Khi sóng tràn xuất hiện ở đê thấp, một phần đáng kể thông lượng khối lượng bị vận chuyển vào trong đồng, làm giảm lưu lượng dòng phản hồi hướng biển, từ đó làm giảm độ sâu hố xói chân đê khoảng 15 - 20% so với đê cao không tràn trong cùng điều kiện sóng tới.

+-----------------------------------------------------------------------------+
| BẢNG SO SÁNH THỦY ĐỘNG LỰC VÀ XÓI CHÂN ĐÊ GIỮA CÁC KẾT CẤU (Thí nghiệm D65) |
+-------------------------------+---------------+---------------+-------------+
| Thông số đặc trưng            | Đê cao (CW0)  | Đê thấp (TW0) | Đê thấp TW10|
+-------------------------------+---------------+---------------+-------------+
| Hệ số phản xạ chân đê (K_r0)   | 0,28 - 0,32   | 0,18 - 0,22   | 0,35 - 0,42 |
| Vận tốc dòng phản hồi cực đại | 0,24 m/s      | 0,19 m/s      | 0,21 m/s    |
| Tỷ số sâu xói (z_max/H_m0i)   | 0,22          | 0,15          | 0,28        |
| Vùng ảnh hưởng xói (x/L)      | 0,35          | 0,25          | 0,40        |
+-------------------------------+---------------+---------------+-------------+

Implications đa chiều

Về mặt lý thuyết, các phát hiện này bổ sung mắt xích còn thiếu trong lý thuyết vận chuyển bùn cát ngang bờ của Van Rijn (1993) và Steetzel (1993), chứng minh rằng không thể xem nhẹ sóng phản xạ cục bộ khi tính toán biến dạng đáy sát công trình ven biển.

Về mặt thực tiễn công trình, nghiên cứu chỉ ra rằng tường đỉnh tuy có tác dụng ngăn sóng tràn hiệu quả cho mái trong, nhưng lại là tác nhân chính làm gia tăng xói lở chân đê mái ngoài. Vì vậy, các giải pháp thiết kế đê biển khu vực Bắc Bộ bắt buộc phải bố trí lớp bảo vệ chân đê tương ứng. Mô phỏng số trên đê biển Nam Định chứng minh:

  • Giải pháp thảm đá bảo vệ chân rộng $B = 3 \div 6$ m dời vị trí hố xói ra xa chân kè, giảm độ sâu xói trực tiếp tại chân đê tới 70%.
  • Giải pháp kè mỏ hàn giảm sóng dài $35 \div 70$ m làm tiêu tán năng lượng sóng tới trước khi chạm đê, kéo giảm hệ số phản xạ cục bộ xuống dưới 0,15 và triệt tiêu nguy cơ trượt chân mái đê.

Limitations và Future Research

Nghiên cứu có một số giới hạn khoa học cần được nhìn nhận khách quan:

  1. Giới hạn 2D trong máng sóng: Thí nghiệm máng sóng 2D và mô hình Wadibe-TC giả định sóng truyền thẳng góc với đường bờ, chưa xét đến góc sóng tới xiên $(\theta \ne 0)$ làm phát sinh dòng chảy dọc bờ (longshore current) và tương tác 3D phức tạp.
  2. Giới hạn về thành phần hạt bùn cát: Thí nghiệm chỉ tập trung vào bùn cát hạt rời ($d_{50} = 0,18$ mm), chưa mô phỏng được tính chất kết dính của bùn sét và bùn mịn vốn rất phổ biến tại các bãi triều cửa sông vùng đồng bằng Bắc Bộ.
  3. Điều kiện biên gió bão: Mô hình vật lý chưa tái hiện được tác động của ứng suất gió trực tiếp trên mặt nước máng sóng trong quá trình sóng vỡ và sóng tràn.

Định hướng nghiên cứu tiếp theo bao gồm:

  • Mở rộng mô hình Wadibe-TC sang hệ không gian 2DH/3D tích hợp sóng xiên và dòng dọc bờ.
  • Nghiên cứu cơ chế xói chân đê trên nền đất dính và bãi bùn có thảm rừng ngập mặn phía trước.
  • Đánh giá tác động của biến đổi khí hậu và nước biển dâng cực hạn đến sự gia tăng tải trọng phản xạ lên đê biển.

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

Về mặt học thuật, công trình cung cấp bộ cơ sở dữ liệu thực nghiệm độ phân giải cao về sóng - dòng - xói chân đê mái nghiêng, phục vụ hiệu chỉnh cho các mô hình mã nguồn mở quốc tế như XBeach hay SWASH. Tiềm năng trích dẫn ước tính đạt mức cao trong các tạp chí chuyên ngành Kỹ thuật Công trình Biển (Coastal Engineering, Ocean Engineering).

Về mặt chính sách và quản lý đê điều, kết quả nghiên cứu đóng góp trực tiếp cho Bộ Nông nghiệp và Phát triển Nông thôn trong việc cập nhật Tiêu chuẩn Quốc gia về thiết kế đê biển (TCVN 9901:2014), đặc biệt là quy chuẩn tính toán kích thước thảm đá và ống buy bảo vệ chân đê chịu tác động của bão cấp 12 - 15 kết hợp triều cường.

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

  • Nghiên cứu sinh và học giả chuyên ngành Kỹ thuật Biển: Tiếp cận khung phương pháp luận mảng thẳng đứng và mã nguồn mô hình Wadibe-TC cập nhật để mở rộng các bài toán tương tác sóng - công trình.
  • Kỹ sư tư vấn thiết kế công trình thủy lực - hàng hải: Có công cụ định lượng chính xác chiều sâu hố xói chân đê $z_{max}$ để tối ưu hóa chiều sâu chôn cừ, chân khay và bề rộng thảm đá bảo vệ, giúp giảm 15 - 20% chi phí xây dựng không cần thiết.
  • Cơ quan quản lý phòng chống thiên tai cấp tỉnh (Nam Định, Hải Phòng, Thái Bình): Sở hữu cơ sở khoa học để xây dựng bản đồ rủi ro xói lở chân đê và phương án hộ đê khẩn cấp trong mùa bão lũ.

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

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

Đóng góp lý thuyết độc đáo nhất là xác lập cơ chế chi phối của sóng phản xạ cục bộ đối với sự gia tăng nồng độ bùn cát lơ lửng trong đới sóng vỡ trước công trình mái nghiêng, trực tiếp mở rộng lý thuyết cân bằng động lượng sóng của Svendsen (1984) và mô hình vận chuyển bùn cát 1DV của Steetzel (1993) thông qua việc tích hợp hàm suy giảm tương quan không gian $F_x(x/L)$.

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

So với nghiên cứu của Neshaei (2009) chỉ phân tích sóng đều trên bãi tự nhiên và Zanuttigh (2008) chỉ khảo sát hệ số $K_r$ trung bình toàn hệ thống, công trình này đã triển khai kỹ thuật mảng thẳng đứng của Sheremet (2002) trên phổ sóng ngẫu nhiên kết hợp đo đạc đồng thời trắc diện vận tốc ADV đa điểm, làm rõ quy luật tiêu tán năng lượng phản xạ cục bộ trong phạm vi $x/L \le 0,4$ trước chân đê có tường đỉnh.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là sự đối nghịch giữa vận tốc dòng chảy và cường độ xói lở: Mặc dù sóng phản xạ làm suy giảm độ lớn của vận tốc dòng phản hồi trung bình $u_0$ khoảng 15%, nhưng độ sâu hố xói cực đại $z_{max}$ lại tăng thêm tới 28% do sóng phản xạ kích thích mạnh mẽ ứng suất tiếp tầng đáy, làm tăng đột biến lượng bùn cát lơ lửng $C(z)$.

4. Luận án có cung cấp quy trình tái lập (replication protocol) rõ ràng không?

Luận án cung cấp chi tiết toàn bộ điều kiện biên hình học (tỷ lệ máng sóng, kích thước mô hình đê $m = 4$, chiều cao tường $h_w = 10$ cm), đặc tính trầm tích ($d_{50} = 0,18$ mm, $\rho_s = 2.650\text{ kg/m}^3$), vị trí tọa độ các đầu đo WG/ADV và hệ phương trình toán học giải tích, cho phép tái lập hoàn toàn trên bất kỳ hệ thống máng sóng tiêu chuẩn nào.

5. Khung chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Lộ trình nghiên cứu dài hạn tập trung vào 3 trọng tâm: (1) Nâng cấp mô hình Wadibe-TC thành phần mềm 3D mã nguồn mở tích hợp tương tác cấu trúc - dòng chảy (FSI); (2) Ứng dụng trí tuệ nhân tạo (AI/Machine Learning) dự báo xói chân đê thời gian thực từ dữ liệu radar/vệ tinh; (3) Thử nghiệm các giải pháp công trình xanh giảm sóng kết hợp rạn san hô nhân tạo và rừng ngập mặn bảo vệ chân đê thích ứng biến đổi khí hậu.

Kết luận

  1. Nghiên cứu đã làm sáng tỏ bản chất thủy động lực học của sự tương tác giữa sóng ngẫu nhiên trong bão và đê biển mái nghiêng có tường đỉnh đặc trưng khu vực Bắc Bộ.
  2. Thiết lập thành công công thức hồi quy hàm suy giảm hệ số phản xạ cục bộ $F_x(x/L)$, xác định chính xác phạm vi ảnh hưởng của sóng phản xạ giới hạn trong vùng $x/L \le 0,4$ trước chân đê.
  3. Chứng minh cơ chế sóng phản xạ kích động lớp bùn cát đáy, làm tăng hệ số xáo trộn $\varepsilon_s$ và gia tăng độ sâu hố xói chân đê từ 22% đến 35%.
  4. Phát triển và kiểm định thành công mô hình toán học Wadibe-TC cải tiến, đạt độ tin cậy cao ($NSE > 0,90$) trong việc dự tính xói lở chân đê dưới tác động của bão biển.
  5. Ứng dụng giải quyết thành công bài toán thực tiễn tại đê biển Thịnh Long (Nam Định), định lượng chính xác hiệu quả giảm xói của giải pháp thảm đá rộng $3 \div 6$ m và kè mỏ hàn dài $35 \div 70$ m, tạo tiền đề khoa học vững chắc cho công tác quy hoạch, thiết kế và bảo vệ an toàn hệ thống đê biển Việt Nam.