Tổng quan về luận án

Nghiên cứu động lực học bờ biển và tương tác giữa sóng bão với các công trình bảo vệ bờ là một trong những bài toán cốt lõi của kỹ thuật công trình thủy và hải văn cận bờ. Tại Việt Nam, hệ thống đê biển trải dài hơn 2.600 km, trong đó khu vực Bắc Bộ chiếm trên 65% tổng chiều dài toàn quốc. Hệ thống này giữ vai trò huyết mạch trong phòng chống thiên tai, bảo vệ dân sinh và cơ sở hạ tầng trước tác động ngày càng gia tăng của bão nhiệt đới và nước biển dâng. Tuy nhiên, các công trình đê biển hiện nay phổ biến có cao trình đỉnh từ +5,0 m đến +5,5 m, bề rộng đỉnh 3 – 6 m, mái ngoài dốc m = 3 – 4 và cơ đê rộng 3 – 5 m. Khi bão lớn đổ bộ kết hợp triều cường, sóng lớn tiến sát bờ gây sóng tràn với lưu lượng hàng trăm lít/giây trên một mét dài, đồng thời tạo ra hiện tượng xói lở chân đê nghiêm trọng. Thực tế ghi nhận trong bão Damrey (9/2005) tại bờ biển Nam Định, độ sâu hố xói trước chân đê hạ thấp tới 1,0 m sau 3 cơn bão liên tiếp, gây sạt lở và sụp đổ kết cấu mái kè.

Khoảng trống học thuật (research gap) trọng tâm nằm ở chỗ: các mô hình thủy động lực học và biến đổi hình thái ven bờ truyền thống (như LITPROF, UNIBEST-TC, DUROSTA hay XBeach) phần lớn giả định biên công trình là biên phản xạ cứng toàn phần hoặc bãi biển tự nhiên đơn thuần, chưa mô hình hóa đầy đủ cơ chế tương tác phức hợp giữa sóng phản xạ từ đê mái nghiêng (có hoặc không có tường đỉnh) tới cấu trúc dòng phản hồi (undertow) và gradient vận chuyển bùn cát ngang bờ (cross-shore sediment transport) trong điều kiện sóng tràn cục bộ.

Nhằm giải quyết triệt để vấn đề này, luận án tập trung giải quyết 4 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 4 giả thuyết khoa học (Hypotheses - H):

  • RQ1 & H1: Hệ số phản xạ sóng ($K_r$) biến đổi phi tuyến theo không gian trước chân đê và suy giảm theo hàm tương quan khoảng cách tương đối $x/L$, chịu sự chi phối mạnh của chỉ số Iribarren ($\xi_0$) và độ cao lưu không tương đối ($R_c/H_s$).
  • RQ2 & H2: Sự hiện diện của sóng phản xạ làm suy giảm thông lượng động lượng sóng tới nhưng lại khuấy động trường ứng suất bức xạ ($S_{xx}$), làm biến dạng trắc diện vận tốc dòng phản hồi $u(z)$ theo phương thẳng đứng.
  • RQ3 & H3: Sóng phản xạ làm gia tăng hệ số xáo trộn bùn cát ($\varepsilon_s$) và nồng độ bùn cát lơ lửng sát đáy ($C_0$), dẫn đến gia tăng độ sâu và quy mô hố xói chân đê cục bộ dù vận tốc dòng phản hồi trung bình có xu hướng giảm.
  • RQ4 & H4: Việc tích hợp module sóng phản xạ không gian vào mô hình toán Wadibe-TC cho phép mô phỏng chính xác động thái hình thành hố xói chân đê trong bão và tối ưu hóa giải pháp công trình bảo vệ chân kè tại thực địa Nam Định.

Nghiên cứu ứng dụng khung lý thuyết cân bằng động lượng và thông lượng khối lượng sóng vỡ (Svendsen, 1984; Longuet-Higgins & Stewart, 1964) kết hợp lý thuyết khuếch tán bùn cát lơ lửng (Steetzel, 1993; Van Rijn, 1993). Luận án được thực hiện với quy mô thí nghiệm vật lý máng sóng toàn diện trên 2 mô hình lòng cứng (đo sóng và vận tốc) và lòng động (đo địa hình xói bồi với độ dốc bãi 1/100 và 1/40), kiểm định qua các trận bão thực tế (Damrey 2005, Doksuri 2017), mang lại giá trị định lượng đột phá cho công tác thiết kế đê biển.

Literature Review và Positioning

Nghiên cứu về dòng phản hồi và vận chuyển bùn cát vùng sóng vỡ đã trải qua hơn một thế kỷ phát triển, bắt đầu từ những phát hiện định tính ban đầu của Johnson (1919) và các thí nghiệm máng sóng của Bagnold (1940). Nền tảng lý thuyết định lượng về dòng phản hồi chỉ thực sự hình thành khi Nielsen và Sorensen (1970) cùng Svendsen (1984) chứng minh rằng dòng phản hồi là hệ quả tất yếu của sự mất cân bằng theo phương thẳng đứng giữa thông lượng động lượng sóng (gradient ứng suất bức xạ $\partial S_{xx}/\partial x$) và gradient áp suất thủy tĩnh do nước dềnh sóng ($\partial \bar{\eta}/\partial x$). Dòng chảy hướng ra biển dưới chân sóng đóng vai trò đền bù cho thông lượng khối lượng ($m_t$) do chuyển động sóng và cuộn sóng (wave rollers) vận chuyển vào bờ.

Trong trường phái mô hình hóa vận chuyển bùn cát ngang bờ, Steetzel (1993) phát triển công thức trắc diện vận tốc dòng phản hồi dạng logarit-tuyến tính: $$u(z) = u_0 + K_{lin} z + K_{log} \ln \left[ 1 + \frac{az}{\varepsilon_0} \right]$$ đồng thời tích hợp vào điều kiện biên bảo toàn khối lượng cục bộ $\int_{0}^{d_t} u(z) dz = -m_t/\rho$. Song song đó, Fredsøe & Deigaard (1992) và Van Rijn (1993) hoàn thiện lý thuyết lớp biên sóng - dòng chảy và phương trình khuếch tán nồng độ bùn cát $C(z) = C_0 f_c(z)$.

Tuy nhiên, tồn tại một tranh luận học thuật lớn giữa hai nhóm quan điểm:

  • Nhóm thứ nhất (Neshaei et al., 2009; Goda & Suzuki, 1976): Khẳng định khi hệ số phản xạ sóng $K_r$ tăng, thông lượng sóng tiến vào bờ giảm ($m_t = (m_w + m_r)_i - (m_w)_r$), do đó độ lớn dòng phản hồi bắt buộc phải suy giảm, kéo theo sự suy giảm cường độ vận chuyển cát ra xa bờ.
  • Nhóm thứ hai (Vũ Minh Cát, 2008; Zanuttigh, 2008; Van der Meer, 2005): Các quan sát thực nghiệm chỉ ra rằng sóng phản xạ từ kết cấu đê mái nghiêng tạo ra giao thoa sóng mạnh, gây chuyển động rối quy mô lớn sát đáy, làm tăng đột biến hệ số xáo trộn bùn cát $\varepsilon_s$ và dẫn đến xói lở chân đê sâu hơn dù dòng phản hồi trung bình suy giảm.

So sánh với các công trình quốc tế, nghiên cứu của Neshaei et al. (2009) chỉ dừng lại ở các kịch bản sóng đều trong máng sóng nhỏ và chưa tính đến tường đỉnh hay hiện tượng sóng tràn; trong khi nghiên cứu của Klopman & van der Meer (1999) xác lập hàm suy giảm sóng phản xạ $F_x(x/L)$ nhưng chỉ áp dụng cho đê đá đổ nước sâu. Luận án này đã định vị chính xác khoảng trống khoa học bằng cách kết hợp cơ chế phản xạ sóng của đê mái nghiêng có tường đỉnh trong điều kiện nước nông vùng sóng vỡ, giải quyết mâu thuẫn giữa độ lớn dòng chảy và nồng độ bùn cát, tạo bước tiến đột phá cho mô hình hóa hình thái ven bờ tại Việt Nam và khu vực châu Á.

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

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

Luận án mở rộng trực tiếp lý thuyết cân bằng động lượng của Svendsen (1984) và lý thuyết vận chuyển bùn cát của Steetzel (1993) sang miền tương tác sóng - công trình phức hợp. Luận án đóng góp 3 định đề lý thuyết cốt lõi:

  1. Định đề 1 (Quy luật suy giảm không gian của $K_r$): Hệ số phản xạ không phải là một hằng số tại chân công trình mà là một trường biến thiên không gian $K_{r}(x)$, suy giảm theo hàm mũ - đa thức liên hệ chặt chẽ với khoảng cách tương đối $x/L$ và chỉ số Iribarren cục bộ $\xi_0$.
  2. Định đề 2 (Biến dạng trường động lượng và dòng phản hồi): Sự xuất hiện của sóng phản xạ bổ sung thành phần ứng suất bức xạ phản xạ $(S_{xx})_r$, làm thay đổi độ dốc mặt nước trung bình $\partial \bar{\eta}/\partial x$ và phân bố lại hệ số nhớt rối $\nu_t(z)$, dẫn đến sự trễ pha và dịch chuyển điểm cực đại của trắc diện $u(z)$ lên cao hơn so với điều kiện bãi tự nhiên.
  3. Định đề 3 (Nghịch lý bùn cát - xói lở chân đê): Cơ chế xói chân đê trong bão không do vận tốc dòng phản hồi đơn thuần quyết định, mà bị chi phối bởi gradient vận chuyển bùn cát lơ lửng $\partial S_x/\partial x$, vốn gia tăng đột biến do tương tác sóng phản xạ - sóng tới kích thích nồng độ bùn cát tham chiếu $C_0$ sát chân đê.
+-------------------------------------------------------------+
|   SÓNG BÃO TIẾN VÀO BỜ (Hs, Tp, Nước dềnh triều cường)      |
+-------------------------------------------------------------+
                              |
                              v
+-------------------------------------------------------------+
| TƯƠNG TÁC ĐÊ MÁI NGHIÊNG CÓ TƯỜNG ĐỈNH (m=4, Rc/Hs, f_r)     |
+-------------------------------------------------------------+
         |                                           |
         v (Phản xạ sóng Kr(x))                      v (Sóng tràn q)
+---------------------------------+         +------------------------+
| Ứng suất bức xạ (Sxx = Sxxi+Sxxr)|         | Tổn thất khối lượng mt |
+---------------------------------+         +------------------------+
         |                                           |
         +--------------------+----------------------+
                              |
                              v
+-------------------------------------------------------------+
| CẤU TRÚC DÒNG PHẢN HỒI u(z) & HỆ SỐ NHỚT RỐI \nu_t(z)        |
+-------------------------------------------------------------+
                              |
                              v
+-------------------------------------------------------------+
| XÁO TRỘN BÙN CÁT: Nồng độ C(z) = C0.fc(z) & \epsilon_s(z)   |
+-------------------------------------------------------------+
                              |
                              v
+-------------------------------------------------------------+
| GRADIENT VẬN CHUYỂN BÙN CÁT LƠ LỬNG \partial Sx/\partial x   |
+-------------------------------------------------------------+
                              |
                              v
+-------------------------------------------------------------+
| HỐ XÓI CHÂN ĐÊ BIỂN: \partial zb/\partial t = -(1/(1-n))dSx/dx|
+-------------------------------------------------------------+

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

Khung phân tích của luận án tích hợp liền giác 3 lý thuyết nền tảng:

  • Lý thuyết sóng ngẫu nhiên phi tuyến trong vùng sóng đổ (Battjes & Janssen, 1978; Sheremet, 2002).
  • Lý thuyết cấu trúc dòng chảy 3 lớp (dòng trôi mặt, dòng phản hồi giữa cột nước, dòng chảy biên đáy).
  • Định luật bảo toàn thể tích bùn cát đáy xốp của Exner-Van Rijn: $\frac{\partial z_b}{\partial t} = -\frac{1}{1-n} \frac{\partial S_x}{\partial x}$.

Điều kiện biên xác định (boundary conditions) được thiết lập tường minh: tại chân đê ($x=0$), hệ số phản xạ đạt cực đại $K_{r0}$, năng lượng sóng phản xạ suy giảm dần ra phía ngoài khơi và triệt tiêu hoàn toàn khi $x/L \ge 0,4$; lưu lượng sóng tràn qua đỉnh đê làm giảm tương ứng thông lượng bù dòng hướng bờ.

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

Thiết kế nghiên cứu

Luận án áp dụng thế giới quan thực chứng (positivism) với phương pháp tiếp cận đa mức độ (multi-level design), kết hợp chặt chẽ giữa nghiên cứu thực nghiệm vật lý trong máng sóng thuỷ lực và mô hình hóa toán học số trị nâng cao.

Thiết kế thí nghiệm bao gồm:

  • Tiêu chuẩn tương tự Froude: Tỷ lệ hình học không méo $N_L = 1/15$ đến $1/25$, bảo đảm đồng dạng về lực trọng trường và lực quán tính trong dòng chảy sóng vỡ.
  • Kịch bản kết cấu đê: Gồm 3 cấu hình điển hình Bắc Bộ: Đê cao không tràn (CW0), Đê thấp không tường đỉnh (TW0), và Đê thấp có tường đỉnh cao tương ứng 10 cm trong mô hình (TW10); mái ngoài $m=4$, mái trong $m=3$.
  • Kịch bản địa hình bãi: Hai độ dốc bãi trước đê gồm bãi thoải $1/100$ và bãi dốc $1/40$.
BẢNG THIẾT KẾ CÁC KỊCH BẢN THỦY ĐỘNG LỰC TRỌNG TÂM:
+---------------+---------------+---------------+---------------+
| Tên kịch bản  | Mực nước D(cm)| Chiều cao H(cm)| Chu kỳ T(s)  |
+---------------+---------------+---------------+---------------+
| D65H15T19     | 65            | 15            | 1,9           |
| D65H17T16     | 65            | 17            | 1,6           |
| D70H15T19     | 70            | 15            | 1,9           |
| D70H19T165    | 70            | 19            | 1,65          |
+---------------+---------------+---------------+---------------+

Quy trình nghiên cứu rigorous

Quy trình đo đạc và phân tích tuân thủ các chuẩn mực quốc tế khắt khe:

  • Tách sóng tới và sóng phản xạ: Áp dụng phương pháp Mảng thẳng đứng (Vertical Array method của Sheremet, 2002), đồng bộ hóa tín hiệu cao độ mặt nước $\eta(x,t)$ từ đầu đo sóng điện dung (Wave Gauge - WG) và tín hiệu vận tốc hạt nước $u(x,t)$ từ đầu đo vận tốc âm học ADV (Current Gauge - CG) bố trí trên cùng trục đứng.
  • Tính toán phổ năng lượng: Phổ mật độ phương sai của sóng $S_{\eta\eta}(f)$, dòng chảy $S_{uu}(f)$ và mật độ phổ hiệp phương sai $S_{\eta u}(f)$ được tích phân qua miền tần số $[f_L, f_U]$ để tính toán thông lượng năng lượng vào $F^+$ và ra $F^-$, từ đó xác định $K_r^2 = F^- / F^+$.
  • Đo trắc diện xói lòng động: Sử dụng thiết bị đo quét địa hình đáy tự động với độ chính xác $\pm 0,5$ mm, quét trắc diện dọc trước và sau khi tác dụng sóng liên tục trong thời gian tương đương 4 giờ bão thực tế.
                   SƠ ĐỒ BỐ TRÍ THÍ NGHIỆM MÁNG SÓNG VẬT LÝ
                                                        [Tường đỉnh]
[Máy tạo sóng] ===> ~~~~~ Sóng tới ~~~~~>                 _||_
                     <--- Sóng phản xạ <---             / | Đê
~~~~~ Bãi cát trước đê (Độ dốc 1/100 hoặc 1/40) ~~~~~ /   | m=4
____________________________________________________/     |
  [WG1]      [WG2/CG2]     [WG3/CG3]      [WG4/CG4] [Chân đê]
  (x/L > 1)   (x/L = 0,4)   (x/L = 0,2)   (x/L = 0,05) (x=0)

Data và phân tích

Mô hình toán học Wadibe-TC (Wave, Dike and Beach – Time dependent Cross-shore) được phát triển cập nhật từ nền tảng Unibest-TC thông qua 4 module tương tác:

  1. Module sóng: Giải phương trình cân bằng năng lượng sóng có xét tiêu tán do sóng vỡ Battjes-Janssen và sóng phản xạ: $\frac{d(c_g E)}{dx} + D_w + D_f = 0$.
  2. Module dòng chảy: Tính toán phân bố dòng phản hồi $u(z)$ giải tích - số trị có tích hợp hàm tương quan phản xạ $F_x(x/L)$ và hệ số suy giảm sóng tràn $f_r = 0,2 R_c/H_s + 0,9$.
  3. Module bùn cát: Giải phương trình khuếch tán với nồng độ đáy $C_0$ phụ thuộc ứng suất tiếp tới hạn và tham số Shields hiệu chỉnh theo năng lượng rối phản xạ.
  4. Module hình thái: Cập nhật biến dạng đáy theo sơ đồ sai phân hữu hạn hiện/ẩn đan lưới vững chắc.

Mô hình được hiệu chỉnh tham số độ nhám đáy Nikuradse $k_s$, hệ số nhớt rối $\nu_t$, và hệ số sóng vỡ $\gamma$, đạt chỉ số Nash-Sutcliffe $NSE > 0,88$ và sai số quân phương $RMSE < 0,04$ m đối với cả trắc diện sóng và cao độ xói đáy.

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 khoa học đột phá có giá trị thực nghiệm sâu sắc:

  1. Quy luật suy giảm không gian của sóng phản xạ: Luận án chứng minh hệ số phản xạ $K_r$ suy giảm phi tuyến từ chân đê ra ngoài khơi theo khoảng cách tương đối $x/L$. Tỷ số chiều cao sóng tổng cộng trên sóng tới tuân theo hàm tương quan không gian: $$\frac{H_{m0,x}}{H_{m0i,x}} = \sqrt{1 + K_{r,0}^2 + 2 K_{r,0} F_x(x/L)}$$ với $F_x(0) = 1$ tại chân đê và $F_x(x/L) \approx 0$ khi $x/L \ge 0,4$. Điều này xác nhận vùng ảnh hưởng thủy động lực trực tiếp của sóng phản xạ tập trung trong phạm vi 0,4 lần bước sóng trước chân đê.
  2. Ảnh hưởng kép của tường đỉnh: Sự hiện diện của tường đỉnh (kịch bản TW10) làm tăng hệ số phản xạ sóng $K_{r0}$ thêm 15 – 25% so với đê không tường đỉnh (TW0). Sóng dội từ tường đỉnh kết hợp với gió bão làm tăng đột biến lượng sóng bắn và sóng tràn qua đê.
  3. Nghịch lý dòng phản hồi và xói lở: Mặc dù sóng phản xạ làm giảm vận tốc dòng phản hồi trung bình ngang bờ $u_0$ từ 8 – 14% do suy giảm thông lượng động lượng sóng tới, nhưng năng lượng rối giao thoa làm tăng nồng độ bùn cát lơ lửng sát đáy $C_0$ lên 1,4 – 1,8 lần. Do đó, tổng lượng vận chuyển cát ra xa bờ $S_x$ tăng mạnh, làm hố xói chân đê sâu hơn từ 20 – 35% so với trường hợp không xét phản xạ.
  4. Quy luật độ sâu hố xói cực đại: Trong điều kiện sóng bão ngẫu nhiên kéo dài 4 giờ, chiều sâu hố xói lớn nhất trước chân đê đạt giá trị xấp xỉ bằng $1/5$ chiều cao sóng trước hố xói ($z_{scour} \approx 0,2 H_{m0}$), và có tương quan thuận tuyến tính với chỉ số sóng vỡ Iribarren $\xi_0$.
  5. Ảnh hưởng của độ dốc bãi: Bãi dốc $1/40$ tạo ra hố xói hẹp nhưng sâu và áp sát chân kè hơn so với bãi thoải $1/100$, do đới sóng vỡ và vùng tập trung năng lượng phản xạ bị co hẹp sát chân công trình.
BIẾN THIÊN HỆ SỐ PHẢN XẠ VÀ CHIỀU SÂU HỐ XÓI THEO KẾT CẤU ĐÊ:
+--------------------------+--------------------+------------------------+
| Loại kết cấu đê          | Hệ số phản xạ Kr0  | Chiều sâu hố xói z_max |
+--------------------------+--------------------+------------------------+
| Đê cao không tràn (CW0)  | 0,35 - 0,42        | 1,25 * z_base          |
| Đê thấp không tường (TW0)| 0,20 - 0,28        | 1,00 * z_base (chuẩn)  |
| Đê thấp có tường (TW10)  | 0,28 - 0,36        | 1,22 * z_base          |
+--------------------------+--------------------+------------------------+

Implications đa chiều

  • Về mặt lý thuyết: Luận án giải quyết triệt để tranh luận khoa học kéo dài về vai trò của sóng phản xạ trong động lực học bãi trước công trình, thiết lập cơ sở liên kết giữa lý thuyết sóng phi tuyến và mô hình xáo trộn bùn cát 3 chiều rút gọn.
  • Về phương pháp luận: Xây dựng quy trình thực nghiệm đồng bộ mảng thẳng đứng ADV-WG kết hợp mô hình số Wadibe-TC, có khả năng chuyển giao áp dụng cho các dạng công trình biển phức tạp khác như đê chắn sóng xốp, kè giảm sóng ngầm.
  • Về ứng dụng thực tiễn: Mô phỏng kiểm định tại tuyến đê biển Thịnh Long (Nam Định) trong bão Damrey 2005 cho thấy mô hình tái hiện chính xác vị trí và độ sâu hố xói 1,0 m thực tế. Luận án đề xuất giải pháp kỹ thuật có tính ứng dụng cao:
    • Bố trí thảm đá gia cố chân đê rộng từ 3 – 6 m (tương ứng với chiều rộng hố xói $L_{scour}$ dự báo).
    • Kết hợp hệ thống kè mỏ hàn dài 35 – 70 m hướng dòng, giảm vận tốc dòng phản hồi và bồi lắng bùn cát bảo vệ chân kè bền vững.

Limitations và Future Research

Luận án thẳng thắn chỉ ra 4 giới hạn nghiên cứu xuất phát từ điều kiện kỹ thuật và giả thiết mô hình:

  1. Giới hạn trường sóng 2D: Thí nghiệm máng sóng và mô hình Wadibe-TC tập trung mô phỏng bài toán biến đổi ngang bờ 1D/2DV, chưa xét đến góc sóng tới xiên bãi lớn ($>30^\circ$) vốn tạo ra dòng dọc bờ (longshore currents) kết hợp phức tạp.
  2. Đặc tính hạt bùn cát đồng nhất: Mô hình vật lý sử dụng cát thạch anh có cấp phối hạt $d_{50}$ chuẩn, chưa mô phỏng hết tính chất dính bám của bùn sét và trầm tích hạt mịn đặc trưng của một số vùng bãi bồi ven biển Bắc Bộ.
  3. Giả thiết điều kiện biên máng sóng: Mặc dù đã tính đến lưu lượng sóng tràn, nhưng mô hình vật lý chưa tái hiện hoàn chỉnh chu trình hồi lưu nước sau đê ra biển qua các cống tiêu thoát lớn trong bão.

Định hướng nghiên cứu tương lai:

  • Nâng cấp mô hình Wadibe-TC lên mô hình hình thái 3D tích hợp sóng - dòng tổng hợp (Wave-Current 3D coupling).
  • Nghiên cứu động lực học xói lở chân đê sinh thái kết hợp rừng ngập mặn phía trước bãi và mái đê phủ thảm thực vật giảm sóng.
  • Ứng dụng kỹ thuật đo đạc ảnh số PIV (Particle Image Velocimetry) độ phân giải cao để quan sát chi tiết cấu trúc vi xoáy rối sát lớp biên đáy chân đê.

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

  • Tác động học thuật: Cung cấp bộ dữ liệu thực nghiệm chuẩn hóa (benchmark datasets) về sóng phản xạ và vận tốc dòng phản hồi cho cộng đồng nghiên cứu kỹ thuật bờ biển trong và ngoài nước; dự kiến đóng góp vào các công bố uy tín thuộc danh mục ISI/Scopus (Q1/Q2) trong lĩnh vực Coastal Engineering và Ocean Engineering.
  • Tác động ngành và công nghệ: Chuyển giao phần mềm mô hình Wadibe-TC cho các viện nghiên cứu, trường đại học chuyên ngành thủy lợi - hải văn và các tập đoàn tư vấn xây dựng công trình thủy, nâng cao độ chính xác trong tính toán thiết kế đê kè biển.
  • Tác động chính sách và an ninh ven biển: Cung cấp luận cứ khoa học vững chắc cho Bộ Nông nghiệp và Phát triển nông thôn cùng UBND các tỉnh ven biển Bắc Bộ (Nam Định, Hải Phòng, Thái Bình, Ninh Bình) trong việc rà soát, điều chỉnh Tiêu chuẩn Quốc gia (TCVN) về thiết kế đê biển thích ứng với biến đổi khí hậu.
  • Lợi ích kinh tế - xã hội: Giảm thiểu nguy cơ vỡ đê do xói chân, tiết kiệm hàng trăm tỷ đồng kinh phí sửa chữa, duy tu công trình hàng năm và bảo vệ an toàn tính mạng, tài sản của hàng triệu người dân vùng duyên 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 tách sóng mảng thẳng đứng tiên tiến, quy trình thí nghiệm lòng động và mã nguồn mô hình toán Wadibe-TC để phát triển các đề tài chuyên sâu.
  • Giảng viên và nhà khoa học chuyên ngành công trình biển: Kế thừa khung lý thuyết hoàn chỉnh về tương tác sóng - dòng - bùn cát trước công trình để giảng dạy và mở rộng nghiên cứu.
  • Kỹ sư thiết kế và chuyên gia tư vấn: Sở hữu công cụ định lượng chính xác kích thước hố xói chân đê ($z_{scour}, L_{scour}$) để tính toán quy mô thảm đá gia cố chân và cọc cừ bảo vệ tối ưu, tránh lãng phí vật liệu hoặc thiếu an toàn.
  • Cơ quan quản lý đê điều và phòng chống thiên tai: Có cơ sở khoa học để lập bản đồ phân vùng rủi ro xói lở chân đê trong bão và xây dựng phương án hộ đê chủ động trong mùa mưa bão.

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 làm sáng tỏ cơ chế nghịch lý giữa suy giảm vận tốc dòng phản hồi và gia tăng xói lở chân đê do sóng phản xạ. Luận án đã mở rộng trực tiếp lý thuyết cân bằng động lượng của Svendsen (1984) và mô hình trắc diện dòng phản hồi của Steetzel (1993) bằng cách tích hợp trường sóng phản xạ biến thiên không gian $K_r(x)$ và hàm tương quan $F_x(x/L)$ vào phương trình phân bố ứng suất bức xạ và phương trình khuếch tán bùn cát lơ lửng.

2. Điểm đổi mới về phương pháp nghiên cứu so với các công trình trước đây là gì?

So với nghiên cứu của Neshaei et al. (2009) (chỉ dùng sóng đều cho đê đứng) và Vũ Minh Cát (2008) (chỉ đo địa hình xói, không đo trường dòng chảy), luận án đã thiết lập quy trình thực nghiệm đồng bộ: đo đồng thời trường sóng và trường vận tốc 3D bằng mảng thẳng đứng ADV-WG trên mô hình đê mái nghiêng có tường đỉnh; đồng thời cập nhật toàn diện module sóng - dòng - bùn cát trong mô hình toán Wadibe-TC có kiểm định chéo giữa phòng thí nghiệm và hiện trường.

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

Phát hiện phản trực giác nhất là: Mặc dù đê có tường đỉnh làm giảm vận tốc dòng phản hồi trung bình $u_0$ (do sóng phản xạ triệt tiêu một phần thông lượng sóng tới), nhưng độ sâu hố xói chân đê lại tăng thêm 22% so với đê không tường đỉnh. Nguyên nhân do sóng dội từ tường đỉnh tạo ra năng lượng rối giao thoa cực mạnh tại chân đê, đẩy nồng độ cát lơ lửng $C_0$ lên gấp 1,5 lần, khiến suất vận chuyển cát ra ngoài khơi tăng vọt.

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

Có. Luận án cung cấp chi tiết toàn bộ thông số hình học mô hình (tỷ lệ $1/25$, mái $m=4$, tường đỉnh 10 cm), đặc tính hạt cát ($d_{50}$, trọng lượng riêng), tọa độ vị trí bố trí đầu đo (WG, CG), các kịch bản sóng (JONSWAP, $H_s, T_p$) và các phương trình vi phân sai phân trong mô hình Wadibe-TC, bảo đảm khả năng tái lập độc lập 100% trong các phòng thí nghiệm thủy lực chuẩn quốc tế.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Chương trình 10 năm định hướng vào 3 trụ cột: (1) Phát triển mô hình Wadibe-3D tích hợp trí tuệ nhân tạo (AI/Machine Learning) dự báo xói chân đê thời gian thực theo dữ liệu bão vệ tinh; (2) Nghiên cứu giải pháp công trình xanh giảm sóng phản xạ bằng bê tông xốp sinh thái và rạn san hô nhân tạo; (3) Quan trắc tự động liên tục xói lở chân đê bằng cảm biến thông minh IoT tại các điểm nóng xói lở Bắc Bộ.

Kết luận

  1. Xác lập quy luật phân bố không gian của sóng phản xạ: Chứng minh hệ số $K_r(x)$ suy giảm phi tuyến từ chân đê ra biển và vùng ảnh hưởng thủy động lực tập trung trong phạm vi $x/L < 0,4$.
  2. Làm rõ cơ chế động lực học dòng phản hồi: Định lượng sự biến dạng trắc diện vận tốc $u(z)$ và phân bố ứng suất bức xạ dưới tác động đồng thời của sóng phản xạ và sóng tràn qua đê mái nghiêng.
  3. Khám phá bản chất tương tác xói lở chân đê: Khẳng định xói chân đê bị chi phối bởi sự gia tăng nồng độ bùn cát lơ lửng do năng lượng rối phản xạ, với độ sâu xói cực đại đạt xấp xỉ $1/5$ chiều cao sóng trước chân đê ($z_{scour} \approx 0,2 H_{m0}$).
  4. Nâng cấp và làm chủ công cụ mô hình số Wadibe-TC: Tích hợp thành công các module sóng phản xạ, dòng phản hồi và biến đổi địa hình đáy, đạt độ tin cậy cao ($NSE > 0,88$).
  5. Đề xuất giải pháp bảo vệ đê biển Nam Định hiệu quả: Khẳng định tính khả thi của giải pháp thảm đá gia cố chân rộng 3 – 6 m kết hợp hệ thống kè mỏ hàn 35 – 70 m để triệt tiêu hố xói trong bão lớn.
  6. Mở ra hướng tiếp cận mới trong tiêu chuẩn thiết kế đê biển: Cung cấp cơ sở khoa học then chốt để chuyển đổi từ phương pháp thiết kế tĩnh truyền thống sang phương pháp thiết kế động lực học hình thái hiện đại tại Việt Nam và quốc tế.