Tổng quan về luận án

Xói lở bờ biển và suy thoái hệ sinh thái rừng ngập mặn tại các vùng đồng bằng châu thổ là thách thức môi trường cấp bách mang tính toàn cầu. Tại hạ lưu lưu vực sông Mekong thuộc vùng Đồng bằng sông Cửu Long (ĐBSCL), hiện tượng xói lở ven bờ đang diễn ra với quy mô và cường độ nghiêm trọng do sự kết hợp của biến đổi khí hậu, nước biển dâng và tình trạng suy giảm trầm tích nghiêm trọng từ thượng nguồn. Như số liệu từ văn bản luận án đã chỉ rõ: "dưới tác động của các đập thượng nguồn, bùn cát bị giảm sau khi đắp đập, là 60%, 69%, 85%, 94% và 98% tương ứng đối với đồng bằng sông Danube, Mississippi, Rhône, Indus và Nile." Tại vùng ven biển ĐBSCL với tổng chiều dài đường bờ 744 km, có tới 268 km đang trong tình trạng sạt lở mạnh; trong đó: "Tỷ lệ phần trăm xói lở bờ biển so với chiều dài bờ biển toàn vùng khoảng 35,9% và tỷ lệ phần trăm đường bờ xói lở đặc biệt nguy hiểm so với toàn bộ chiều dài đường bờ biển vùng ĐBSCL khoảng 11,34%". Các giải pháp công trình truyền thống như kè lát mái nghiêng (dạng bị động) hay đê chắn sóng trọng lực không thấm nước thường bộc lộ nhiều hạn chế về mặt chi phí, gây lún sụt trên nền đất yếu và ngăn chặn dòng bùn cát tự nhiên nuôi bãi.

                    +-------------------------------------------------+
                    |              BỐI CẢNH KHOA HỌC                  |
                    | 744 km bờ biển ĐBSCL: 35.9% xói lở, 11.34% nguy |
                    | hiểm; thiếu hụt trầm tích thượng nguồn (60-98%) |
                    +------------------------+------------------------+
                                             |
                                             v
                    +-------------------------------------------------+
                    |               RESEARCH GAP CHÍNH                |
                    | Chưa có mô hình lý thuyết & công thức thực      |
                    | nghiệm cho đê cọc ly tâm đá đổ (thay đổi B & n) |
                    +------------------------+------------------------+
                                             |
                     +-----------------------+-----------------------+
                     |                                               |
                     v                                               v
+------------------------------------------+   +------------------------------------------+
|          MÔ HÌNH VẬT LÝ (2D FLUME)       |   |               MÔ HÌNH SỐ                 |
| 234 kịch bản, 7 kim đo sóng WG1-WG7      |   | Phương trình Boussinesq mở rộng 2 lớp    |
| B = 24, 38, 52 cm; D50 = 4-7 cm (n=0.45) |   | rỗng (FDM), kiểm chứng giải tích & MHVL  |
+--------------------+---------------------+   +---------------------+--------------------+
                     |                                               |
                     +-----------------------+-----------------------+
                                             |
                                             v
                    +-------------------------------------------------+
                    |               KẾT QUẢ ĐỘT PHÁ                   |
                    | • Quy chuẩn tối ưu: B=2.7m, H=2.7m, D50=30-50cm |
                    | • Công thức thực nghiệm tính Kt đa thông số      |
                    | • Cơ chế tiêu tán sóng kết hợp bảo tồn bùn cát   |
                    +-------------------------------------------------+

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ việc các công thức tính toán hệ số truyền sóng kinh điển (d’Angremond et al., 1996; Van der Meer et al., 2005) chủ yếu áp dụng cho đê đá đổ mái nghiêng dạng hình thang thông thường, chưa phản ánh chính xác bản chất thủy lực của đê giảm sóng kết cấu rỗng dạng hai hàng cọc bê tông ly tâm xen kẽ lõi đá hộc (gọi tắt là đê cọc ly tâm). Cấu trúc tường đứng rỗng kép này sở hữu tính chất tương tác thủy động lực học dị hướng rất đặc thù, trong đó sự thay đổi đồng thời của bề rộng đê ($B$) và độ rỗng lõi đá ($n$) dưới tác động của phổ sóng ngẫu nhiên chưa từng được mô hình hóa toàn diện.

Nghiên cứu tập trung giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 3 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Tương quan định lượng giữa các tham số hình học đê ($B, R_c$), đặc trưng đá ($D_{50}, n$) và phổ sóng tới ngẫu nhiên ($H_{m0,i}, T_p, S_{0p}$) tác động như thế nào đến hệ số truyền sóng ($K_t$), phản xạ ($K_r$) và tiêu tán ($K_d$)?
    • H1: Hiệu quả giảm sóng $K_t$ của đê cọc ly tâm tỉ lệ nghịch phi tuyến với bề rộng tương đối $B/H_{m0,i}$ và chịu sự chi phối quyết định bởi mức độ lưu không tương đối $R_c/H_{m0,i}$ kết hợp sức cản ma sát rối trong các lỗ rỗng của lõi đá.
  • RQ2: Làm thế nào để thiết lập một công thức thực nghiệm dự báo chính xác hệ số truyền sóng $K_t$ đặc thù cho đê cọc ly tâm trên nền đất yếu ven biển?
    • H2: Tồn tại một hàm hồi quy phi tuyến đa biến tích hợp đầy đủ các nhóm thông số phi thứ nguyên ($R_c/H_{m0,i}, B/H_{m0,i}, S_{0p}$) cho phép tính toán nhanh $K_t$ với hệ số tương quan $R^2 > 0.85$.
  • RQ3: Phương trình Boussinesq mở rộng cho môi trường hai lớp rỗng có thể thay thế hiệu quả các mô hình vật lý tốn kém để mô phỏng tương tác sóng - đê cọc ly tâm hay không?
    • H3: Mô hình số giải bằng phương pháp sai phân hữu hạn (FDM) dựa trên hệ phương trình Boussinesq mở rộng hai lớp rỗng có khả năng tái hiện chính xác trường sóng biến dạng với sai số bình phương trung bình (RMSE) nằm trong giới hạn kiểm định thủy động lực học.

Khung lý thuyết của nghiên cứu được xây dựng trên cơ sở thuyết sóng phi tuyến kết hợp cơ học môi trường rỗng (Porous Media Hydrodynamics). Đóng góp đột phá của luận án là xác định tổ hợp thông số hình học và vật liệu tối ưu: "đề xuất đê có bề rộng đê B=270cm, chiều cao đê H=270cm và đường kính viên đá đổ trong đê từ 30cm÷50cm (độ rỗng đê khoảng 45%) phù hợp cho dạng đê cọc ly tâm ở cả Biển Đông và Biển Tây khu vực ĐBSCL". Phạm vi nghiên cứu bao quát 234 kịch bản thí nghiệm mô hình vật lý máng sóng 2D và chuỗi mô phỏng số mở rộng, mang ý nghĩa nền tảng cho chiến lược công trình bảo vệ bờ biển bền vững tại Việt Nam.

Literature Review và Positioning

Các nghiên cứu quốc tế về hệ số truyền sóng ($K_t$) qua đê giảm sóng bắt đầu từ các công trình nghiên cứu sóng đều qua chướng ngại vật ngầm của Johnson et al. (1951) và Goda et al. (1969). Đến thập niên 1980–1990, Seelig (1980), Allsop (1983) và Ahrens (1987) đã mở rộng sang phổ sóng ngẫu nhiên trên đê đá đổ mái nghiêng. Đặc biệt, chương trình nghiên cứu DELOS (1998–2002) tại châu Âu với hơn 2.400 thí nghiệm mô hình vật lý (Van der Meer & Daemen, 1994; d’Angremond et al., 1996; Briganti et al., 2003; Van der Meer et al., 2005) đã thiết lập các công thức chuẩn mực cho đê đỉnh thấp (Low-Crested Structures - LCS). Cụ thể, d’Angremond et al. (1996) đề xuất phương trình xác định $K_t$ dựa trên chiều cao lưu không tương đối ($R_c/H_s$), bề rộng tương đối ($B/H_s$) và chỉ số Iribarren ($\xi_{0p}$):

$$K_t = -0.4 \frac{R_c}{H_s} + 0.64 \left(\frac{B}{H_s}\right)^{-0.31} \left(1 - e^{-0.5 \xi_{0p}}\right)$$

Tuy nhiên, trong y văn tồn tại hai luồng quan điểm đối lập sâu sắc:

  1. Trường phái Van der Meer (1990): Cho rằng với các đê đỉnh hẹp ($B/H_s < 8$), đại lượng lưu không tương đối $R_c/H_s$ là thông số chi phối gần như tuyệt đối, tương quan $K_t \sim R_c/H_s$ mang tính tuyến tính bậc nhất và có thể bỏ qua ảnh hưởng của đường kính hạt đá $D_{n50}$ hay độ rỗng lõi đê.
  2. Trường phái Seabrook & Hall (1998) và Calabrese et al. (2002): Khẳng định khi bề rộng đê tăng lên ($B/H_s > 10$) hoặc trong cấu trúc rỗng dị tính, sự tiêu hao năng lượng qua ma sát trong lỗ rỗng đá ($D_{n50}$) và chiều rộng đỉnh đê ($B$) đóng vai trò quyết định, lấn át cơ chế sóng vỡ trên đỉnh dốc.
Tiêu chí so sánh Nghiên cứu DELOS (Châu Âu) Nghiên cứu Presque Isle (Mỹ) Nghiên cứu Luận án (Nguyễn Hải Dương, 2024)
Loại kết cấu Đê đá đổ mái nghiêng ngầm/bán ngập Đê chắn sóng đá rời xa bờ Đê cọc ly tâm đôi chèn đá hộc (tường đứng rỗng)
Cơ chế giảm sóng chính Sóng vỡ trên mái nghiêng dốc Sóng vỡ và nhiễu xạ qua khe hở Ma sát rối trong lõi đá + vỡ sóng cưỡng bức qua hàng cọc
Môi trường nền móng Nền cát, đá cứng ven biển châu Âu Nền cát hồ Erie, Pennsylvania Nền bùn sét cực yếu ven biển châu thổ ĐBSCL
Khả năng gây bồi Hạn chế, chủ yếu chống xói cục bộ Gây bồi tạo mũi cát dạng Tombolo Cho phép bùn cát lơ lửng đi qua để khôi phục rừng ngập mặn
Mô hình tính toán Công thức d'Angremond / Van der Meer Phân tích cân bằng trầm tích ven bờ Boussinesq mở rộng 2 lớp rỗng (FDM) + Công thức thực nghiệm riêng

So với các công trình quốc tế như hệ thống đê rời tại Presque Isle (Mỹ) hay các đê ngầm tại Ý, Tây Ban Nha trong dự án DELOS, công trình của luận án định vị một hướng tiếp cận hoàn toàn mới: phát triển hệ kết cấu bán cứng chịu lực trên nền đất bùn lầy ven biển ĐBSCL. Tại Việt Nam, các nghiên cứu trước đây như Thiều Quang Tuấn et al. (2018, 2020) chỉ tập trung vào hàng rào tre hoặc đê rỗng phi truyền thống lắp ghép; Lê Xuân Tú et al. (2020) bước đầu khảo nghiệm đê cọc ly tâm nhưng chưa xét đến biến thiên đồng thời của bề rộng đê $B$ và đường kính đá $D_{50}$. Luận án đã vượt lên các giới hạn này, lấp đầy khoảng trống học thuật bằng việc phân tích đa biến 234 kịch bản thực nghiệm kết hợp mô hình hóa 2 lớp rỗng.

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

                   +-------------------------------------------------------+
                   |                 TÍCH HỢP 3 LÝ THUYẾT                  |
                   |  1. Thuyết sóng phi tuyến Boussinesq mở rộng          |
                   |  2. Định luật cản Forchheimer trong môi trường rỗng   |
                   |  3. Định luật bảo toàn năng lượng sóng (Kt²+Kr²+Kd²=1)|
                   +---------------------------+---------------------------+
                                               |
                                               v
                   +-------------------------------------------------------+
                   |                 MÔ HÌNH HÓA 2 LỚP RỖNG                |
                   | • Lớp trên: Nước tự do hoặc đá rỗng (n1)             |
                   | • Lớp dưới: Khối đá hộc chèn giữa 2 cọc (n2)         |
                   +---------------------------+---------------------------+
                                               |
                                               v
                   +-------------------------------------------------------+
                   |                 CÔNG THỨC THỰC NGHIỆM                 |
                   |       Kt = f(Rc/Hm0,i, B/Hm0,i, S0p, D50, n)          |
                   |         Phạm vi: -1.49 <= Rc/Hm0,i <= 1.48            |
                   +-------------------------------------------------------+

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

Nghiên cứu tạo ra bước phát triển mở rộng đối với khung lý thuyết của d’Angremond et al. (1996) và Van der Meer et al. (2005) khi chứng minh rằng: đối với đê cọc ly tâm, hệ số truyền sóng $K_t$ không tuân theo quy luật suy giảm tuyến tính đơn thuần theo $R_c/H_s$ mà bị chi phối mạnh bởi hiệu ứng cản thủy lực cục bộ qua các khe cọc và cấu trúc ma trận lỗ rỗng của đá hộc.

Khung mệnh đề lý thuyết được xác lập thông qua nguyên lý bảo toàn năng lượng sóng:

$$K_t^2 + K_r^2 + K_d^2 = 1$$

Trong đó, năng lượng sóng tiêu tán ($E_d$ hay $K_d$) chiếm ưu thế áp đảo so với năng lượng sóng phản xạ ($K_r$) nhờ cơ chế xoáy cuộn rối bên trong thân đê. Nghiên cứu xác lập sự chuyển dịch mô thức (paradigm shift) từ tư duy "chắn sóng tuyệt đối" (ngăn cản hoàn toàn động lực biển) sang tư duy "tiêu tán năng lượng thích ứng" (cho phép trao đổi thủy triều và phù sa mịn đi qua đê để nuôi dưỡng hệ sinh thái rừng ngập mặn).

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

Khung phân tích của luận án tích hợp đồng thời 3 hệ lý thuyết:

  1. Lý thuyết sóng Boussinesq bậc cao: Cho phép mô tả hiện tượng phi tuyến và tán sắc sóng nước nông.
  2. Lý thuyết dòng chảy qua môi trường rỗng Forchheimer: Mô tả tổn thất áp lực do ma sát tầng và ma sát rối:

$$-\nabla p = a \mathbf{u} + b |\mathbf{u}|\mathbf{u}$$

  1. Lý thuyết phổ năng lượng sóng ngẫu nhiên JONSWAP: Đánh giá sự chuyển dịch phổ tần số bậc cao khi sóng tương tác với công trình.

Điều kiện biên và phạm vi hữu dụng của khung phân tích được xác định tường minh: áp dụng cho vùng nước chuyển tiếp $0.14 \le h/L \le 0.5$, tỷ số lưu không tương đối $-1.49 \le R_c/H_{m0,i} \le 1.48$, độ dốc sóng tới $0.009 \le S_{0p} \le 0.051$, và tỷ số độ sâu $B/D$ từ $0.48$ đến $1.04$.

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

+---------------------------------------------------------------------------------------+
|                               QUY TRÌNH NGHIÊN CỨU                                    |
+---------------------------------------------------------------------------------------+
| [Thực địa ĐBSCL] -> Khảo sát 744km đường bờ -> Phân tích đa tiêu chí -> Chọn đê cọc  |
|                                         |                                             |
|                                         v                                             |
| [MHVL Máng sóng 2D] -> 234 Kịch bản -> 7 Kim đo sóng (WG1-WG7) -> Tách sóng tới/phản  |
|                                         |                                             |
|                                         v                                             |
| [Hồi quy Thực nghiệm] -> Xây dựng công thức Kt đa biến (R² cao, kiểm tra S0p, Rc, B)  |
|                                         |                                             |
|                                         v                                             |
| [Mô hình số Boussinesq] -> Rời rạc hóa FDM -> 2 Lớp rỗng -> Kiểm chứng giải tích/MHVL |
+---------------------------------------------------------------------------------------+

Thiết kế nghiên cứu

Luận án vận dụng thế giới quan thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phê phán (Critical Realism). Nghiên cứu sử dụng phương pháp hỗn hợp quy chuẩn (Rigorous Mixed Methods) gồm hai trụ cột: thực nghiệm thủy động lực học trên mô hình vật lý máng sóng 2D và mô hình hóa toán học số trị (Numerical Modeling).

Thiết kế đa mức (Multi-level design) được triển khai từ quy mô nguyên hình (Prototype) bờ biển Cà Mau - Bạc Liêu xuống quy mô mô hình (Model) theo định luật đồng dạng Froude:

  • Tỷ lệ hình học chiều dài và chiều cao: $N_L = 1/7$
  • Tỷ lệ thời gian: $N_t = \sqrt{N_L} = 1/\sqrt{7} \approx 0.378$
  • Tỷ lệ vận tốc: $N_v = \sqrt{N_L} = 1/\sqrt{7}$
  • Tỷ lệ khối lượng: $N_m = N_L^3 = (1/7)^3 = 1/343$

Ba quy mô bề rộng đê được kiểm chứng: $B_{24} = 24\text{ cm}$ (nguyên hình $1.7\text{ m}$), $B_{38} = 38\text{ cm}$ (nguyên hình $2.7\text{ m}$), và $B_{52} = 52\text{ cm}$ (nguyên hình $3.7\text{ m}$).

Quy trình nghiên cứu rigorous

Thực nghiệm được tiến hành tại Máng sóng của Phòng Thí nghiệm Thủy động lực Sông biển (Viện Khoa học Thủy lợi miền Nam tại Bình Dương). Máng sóng có hệ thống tạo sóng điện động lực học điều khiển bằng máy tính, tích hợp mái tiêu sóng hấp thụ năng lượng ở cuối máng.

                      SƠ ĐỒ BỐ TRÍ THÍ NGHIỆM TRONG MÁNG SÓNG (2D)
                      
   Bàn tạo sóng          WG1  WG2  WG3          Đê cọc ly tâm        WG6  WG7    Mái tiêu sóng
   [~~~~~~~~~~] -------- [|]--[|]--[|] -------- [||:ĐÁ:||] -------- [|]--[|] -- /////////
     (Piston)             <--Trước đê-->           Bề rộng B        <--Sau đê-->  (Absorber)
                          (Sóng tới/px)           (24,38,52cm)     (Sóng truyền)

Hệ thống đo đạc gồm 7 kim đo sóng điện dung (Wave Gauges - WG):

  • Cụm trước đê (WG1, WG2, WG3, WG4, WG5): Đo và tách sóng tới ($H_{m0,i}$) và sóng phản xạ ($H_{m0,r}$) bằng thuật toán Mansard & Funke (1980).
  • Cụm sau đê (WG6, WG7): Đo chiều cao sóng truyền qua đê ($H_{m0,t}$) để tính toán $K_t$.

Cốt liệu đá thí nghiệm sử dụng đá hộc đồng nhất kích thước $D_{50} = 4\text{--}7\text{ cm}$ (tương ứng $30\text{--}50\text{ cm}$ ngoài thực tế), khối lượng thể tích xốp $\rho_x$ được kiểm định chính xác qua cân đo thể tích, xác lập độ rỗng vật liệu đạt chuẩn $n = 45%$.

Data và phân tích

Chương trình thí nghiệm bao gồm tổng cộng 234 kịch bản đo đạc. Phân tích hồi quy phi tuyến đa biến (Nonlinear Multivariable Regression) được áp dụng để xác lập hàm truyền sóng thực nghiệm.

Mô hình số được xây dựng trên nền tảng phương trình Boussinesq mở rộng cho môi trường 2 lớp rỗng (Extended Boussinesq equations for two porous layers). Phương pháp sai phân hữu hạn (Finite Difference Method - FDM) dạng lưới lệch (staggered grid) kết hợp bước thời gian Crank-Nicolson được sử dụng để giải hệ phương trình. Kết quả mô phỏng số được đối soát với lời giải giải tích và 234 kịch bản thực nghiệm, đánh giá mức độ hội tụ qua chỉ số sai số chuẩn tuyệt đối (RMSE) của hệ số phản xạ $K_r$ và hệ số tiêu tán $K_d$.

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

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

  • Quy chuẩn thông số mặt cắt đê cọc ly tâm tối ưu: Thực nghiệm trên 234 kịch bản chứng minh mặt cắt đê có bề rộng $B = 2.7\text{ m}$ ($B_{38}$ trên mô hình), chiều cao $H = 2.7\text{ m}$ kết hợp lõi đá hộc $D_{50} = 30\text{--}50\text{ cm}$ (độ rỗng $n \approx 45%$) đạt hiệu suất thủy lực và kinh tế tối ưu. Cấu trúc này làm suy giảm từ $70%$ đến $85%$ năng lượng sóng tới trong các điều kiện bất lợi nhất của cả vùng biển Đông và biển Tây ĐBSCL.
  • Quy luật phi tuyến của tỷ số lưu không tương đối ($R_c/H_{m0,i}$): Khi đỉnh đê ở trạng thái nhô ($R_c/H_{m0,i} > 0.5$), hệ số truyền sóng giảm mạnh xuống $K_t < 0.3$. Ngược lại, khi đê chìm sâu trong triều cường ($R_c/H_{m0,i} < -1.0$), sóng tràn qua đỉnh khiến $K_t$ tiệm cận $0.75\text{--}0.80$.
  • Cơ chế tiêu tán năng lượng nội tại vượt trội: Khác với đê tường đứng cứng tạo sóng phản xạ lớn ($K_r > 0.7$), đê cọc ly tâm lõi đá triệt tiêu phần lớn năng lượng thông qua dòng thấm rối nội tại ($K_d = 0.60\text{--}0.85$, $K_r < 0.35$), giúp bảo vệ triệt để đáy biển trước chân đê không bị xói lở cục bộ do sóng phản xạ.
  • Xác lập 2 công thức thực nghiệm dự báo $K_t$: Luận án xây dựng thành công phương trình dự báo $K_t$ có độ tin cậy cao ($R^2 > 0.85$):

$$K_t = f\left(\frac{R_c}{H_{m0,i}}, \frac{B}{H_{m0,i}}, S_{0p}\right)$$

Phương trình phản ánh chính xác tác động đồng thời của độ dốc sóng đỉnh phổ $S_{0p}$ và bề rộng tương đối của đê.

       HỆ SỐ TRUYỀN SÓNG (Kt) THEO LƯU KHÔNG TƯƠNG ĐỐI (Rc/Hm0,i)
    Kt
   1.0 |      *                                      (Trạng thái ngập)
   0.8 |        *  *
   0.6 |             *  *                            (Bán ngập)
   0.4 |                  *  *  
   0.2 |                        *  *  *  *  *        (Trạng thái nhô)
   0.0 +------------------------------------------
      -1.5   -1.0   -0.5    0.0    0.5    1.0   1.5   Rc/Hm0,i

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung nhánh tri thức về thủy lực công trình bảo vệ bờ dạng tường rỗng composite (cọc bê tông + lõi đá) vào lý thuyết công trình biển ngập nước.
  • Về mặt phương pháp luận: Chứng minh mô hình Boussinesq mở rộng 2 lớp rỗng giải bằng sai phân hữu hạn (FDM) hoàn toàn có thể thay thế các thí nghiệm máng sóng đắt tiền để dự báo nhanh trường sóng và hiệu quả giảm sóng.
  • Về mặt thực tiễn và chính sách: Cung cấp cơ sở khoa học trực tiếp cho các Ban Quản lý Dự án Nông nghiệp & Giao thông tại Cà Mau, Bạc Liêu, Kiên Giang trong việc chuẩn hóa thiết kế đê giảm sóng xa bờ, tối ưu hóa hơn $30%$ chi phí đầu tư so với đê khối đúc sẵn và nâng cao tuổi thọ công trình lên trên 40 năm.

Limitations và Future Research

  1. Giới hạn không gian thí nghiệm máng sóng 2D: Thực nghiệm 234 kịch bản mới chỉ khảo sát sóng truyền vuông góc ($90^\circ$). Hiện tượng sóng truyền xiên góc, sóng nhiễu xạ 3D qua các đoạn đê ngắt quãng và dòng chảy dọc bờ chưa được mô phỏng toàn diện trong nghiên cứu này.
  2. Giả định bùn cát và dòng bùn lơ lửng: Nghiên cứu tập trung vào khía cạnh thủy động lực học sóng (Wave Hydrodynamics) mà chưa trực tiếp thả bùn cát hạt mịn (cohesive sediment) vào máng sóng để định lượng tốc độ bồi tụ thực tế theo thời gian.
  3. Hiện tượng xói chân cọc cục bộ dưới tác động kết hợp của dòng triều và sóng cực đoan: Cần được kiểm chứng thêm về mặt cơ học đất địa kỹ thuật lâu dài.

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda)

  • Thí nghiệm mô hình vật lý 3D trong bể sóng lớn (Wave Basin) để đánh giá góc sóng xiên và hiệu ứng đầu mũi đê.
  • Tích hợp module chuyển bùn cát hạt mịn vào mô hình số Boussinesq 2 lớp rỗng nhằm dự báo chính xác trường bồi lắng sau đê.
  • Nghiên cứu độ bền ăn mòn sunfat và clo của bê tông ly tâm dự ứng lực trong môi trường nước biển ĐBSCL dưới tác động của biến đổi khí hậu.

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

  • Học thuật: Cung cấp bộ cơ sở dữ liệu mở gồm 234 kịch bản thực nghiệm chất lượng cao và công thức dự báo $K_t$ chuẩn mực, tạo tiền đề trích dẫn quan trọng cho các công bố quốc tế thuộc danh mục ISI/Scopus chuyên ngành Coastal Engineering và Ocean Engineering.
  • Ngành xây dựng công trình thủy - hàng hải: Định hình tiêu chuẩn thiết kế mới cho kết cấu đê cọc ly tâm chèn đá, khắc phục tình trạng thiết kế theo cảm tính hoặc sao chép rập khuôn đê mái nghiêng châu Âu.
  • Chính sách & Xã hội: Đóng góp trực tiếp vào Đề án tổng thể chống sạt lở bờ sông, bờ biển vùng ĐBSCL theo Quyết định của Thủ tướng Chính phủ, bảo vệ sinh kế cho hơn 17 triệu dân vùng đồng bằng trước hiểm họa nước biển dâng.

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

  • Nghiên cứu sinh & Nhà khoa học trẻ: Kế thừa bộ công thức thực nghiệm, hệ phương trình Boussinesq mở rộng 2 lớp rỗng và phương pháp luận thiết lập máng sóng 2D.
  • Giảng viên & Chuyên gia công trình biển: Sử dụng số liệu kiểm chứng và kết quả phân tích đa tiêu chí làm giáo trình, tài liệu chuyên khảo cao học và tiến sĩ.
  • Kỹ sư tư vấn & Doanh nghiệp xây lắp: Sở hữu thông số thiết kế chuẩn ($B=2.7\text{ m}, H=2.7\text{ m}, D_{50}=30\text{--}50\text{ cm}$) để áp dụng trực tiếp vào hồ sơ thiết kế bản vẽ thi công.
  • Nhà hoạch định chính sách (Bộ Nông nghiệp & PTNT, UBND các tỉnh ĐBSCL): Có luận cứ khoa học tin cậy để thẩm định, phân bổ nguồn vốn đầu tư công trình phòng chống thiên tai hiệu quả và bền vữ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 mở rộng lý thuyết tương tác sóng - môi trường rỗng sang dạng kết cấu thành đứng hỗn hợp (cọc bê tông kết hợp ma trận đá rỗng), chứng minh sự kết hợp giữa hệ số lưu không tương đối $R_c/H_{m0,i}$, bề rộng tương đối $B/H_{m0,i}$ và độ dốc sóng $S_{0p}$ xác lập nên quy luật tiêu tán sóng phi tuyến, vượt khỏi giới hạn của mô hình đê mái nghiêng d’Angremond (1996).

  2. Đột phá về phương pháp luận so với các nghiên cứu trước đây? So với nghiên cứu của Sasikumar et al. (2020) (chỉ dùng RANS 3D phức tạp, tốn tài nguyên tính toán) hay Cea et al. (2009) (phương trình NLSW thiếu tính tán sắc), luận án đã phát triển mô hình Boussinesq mở rộng 2 lớp rỗng giải bằng sai phân hữu hạn (FDM), kết hợp đối soát với 234 kịch bản thực nghiệm máng sóng 2D đo đạc bằng 7 kim đo điện dung, đạt độ chính xác cao tương đương nhưng tốc độ tính toán nhanh hơn vượt bậc.

  3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì? Hiện tượng đê cọc ly tâm khi ở trạng thái ngập hoàn toàn ($R_c/H_{m0,i} < -0.5$) vẫn duy trì khả năng tiêu tán năng lượng sóng đạt trên $50%$ mà không phụ thuộc vào cơ chế sóng vỡ trên đỉnh dốc. Tổn thất áp lực do dòng thấm cưỡng bức qua các khe cọc và mạng lỗ rỗng đá hộc đóng vai trò hấp thụ năng lượng chủ đạo.

  4. Luận án có cung cấp giao thức nhân bản (Replication Protocol) không? Toàn bộ giao thức thí nghiệm được mô tả chi tiết: tỷ lệ đồng dạng Froude ($N_L = 1/7$), sơ đồ bố trí 7 kim đo sóng (WG1-WG7), đặc tính đá hộc ($D_{50} = 4\text{--}7\text{ cm}$, khối lượng thể tích $\rho_x$, độ rỗng $n = 45%$), các bước sóng kiểm tra phổ JONSWAP và thuật toán phân tách sóng phản xạ Mansard & Funke.

  5. Chương trình nghiên cứu 10 năm tới được vạch ra như thế nào? Tập trung giải quyết bài toán biến động đáy biển 3D (3D Morphodynamics), tương tác sóng - dòng chảy - bùn cát mịn, và ứng dụng trí tuệ nhân tạo (Machine Learning/PINNs) để tối ưu hóa hình học cọc ly tâm theo thời gian thực.

Kết luận

Luận án của NCS. Nguyễn Hải Dương đã hoàn thành xuất sắc các mục tiêu nghiên cứu với các đóng góp cốt lõi:

  1. Xác lập cơ sở khoa học thủy động lực học toàn diện cho đê cọc ly tâm bảo vệ bờ biển ĐBSCL.
  2. Thực hiện 234 kịch bản thí nghiệm máng sóng 2D nghiêm ngặt, làm sáng tỏ cơ chế truyền, phản xạ và tiêu tán sóng qua kết cấu rỗng thành đứng.
  3. Đề xuất quy chuẩn kích thước hình học tối ưu: bề rộng đê $B = 2.7\text{ m}$, chiều cao đê $H = 2.7\text{ m}$, đường kính đá hộc $D_{50} = 30\text{--}50\text{ cm}$ (độ rỗng $n \approx 45%$).
  4. Thiết lập công thức thực nghiệm dự báo hệ số truyền sóng $K_t$ với hệ số tương quan $R^2 > 0.85$, ứng dụng thuận tiện cho kỹ thuật bờ biển.
  5. Ứng dụng thành công mô hình số Boussinesq mở rộng 2 lớp rỗng (FDM), mở ra giải pháp công nghệ mô phỏng số hiện đại thay thế mô hình vật lý.
  6. Mở ra 3 hướng nghiên cứu mới về công trình thủy ven biển thích ứng biến đổi khí hậu, cung cấp giải pháp kỹ thuật bền vững bảo vệ 744 km đường bờ và phục hồi rừng ngập mặn vùng Đồng bằng sông Cửu Long.