Tổng quan về luận án

Vùng đồng bằng sông Cửu Long (ĐBSCL) với hơn 700 km bờ biển đang đối mặt với nguy cơ suy thoái môi trường địa chất và suy giảm nghiêm trọng hệ sinh thái đai rừng ngập mặn (RNM) ven biển. Tại dải bờ biển Tây trải dài từ Mũi Cà Mau (MCM) đến Hà Tiên (HT), sự kết hợp giữa các can thiệp nhân sinh thượng nguồn sông Mê Kông, biến đổi khí hậu (BĐKH), nước biển dâng (NBD) và nền địa chất bùn sét hữu cơ siêu yếu đã làm đảo ngược quá trình bồi tụ tự nhiên hình thành châu thổ. Khu vực này chịu tác động trực tiếp của các chế độ sóng gió mùa Tây Nam (MGTN) và gió mùa Đông Bắc (MGĐ), gây xói lở sâu vào đất liền và tàn phá sinh kế của hơn 18 triệu cư dân vùng châu thổ – nơi đóng góp 50% sản lượng lúa, 65% sản lượng nuôi trồng thủy sản, 70% sản lượng trái cây và 95% lượng gạo xuất khẩu của Việt Nam.

+-------------------------------------------------------------------------+
|                  CÔNG TRÌNH ĐÊ RỖNG PHỨC HỢP                             |
|                                                                         |
|        [Hệ cọc hoa mai: D300 ly tâm] ---> Tiêu hao công lực cản Dpr     |
|   ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Mực nước tính toán           |
|      /================================\                                 |
|     /   Thân đê ngầm rỗng (m = 1:1)    \ ---> Sóng vỡ, ma sát, lưu      |
|    /     Độ rỗng mái đê: 14% diện tích  \     thông qua thân đê (Kt,0)  |
|   +--------------------------------------+                              |
|           Nền đất bùn sét mềm yếu (Bờ biển Tây ĐBSCL)                    |
+-------------------------------------------------------------------------+

Tính tiên phong của luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật Xây dựng Công trình Thủy (Mã số: 9580202) do Nghiên cứu sinh (NCS) Nguyễn Anh Tiến thực hiện (2021) tại Viện Khoa học Thủy lợi Việt Nam dưới sự hướng dẫn khoa học của PGS.TS Tô Văn Thanh và GS.TS Thiều Quang Tuấn, thể hiện qua việc đề xuất giải pháp kết cấu mới phi truyền thống: Công trình đê rỗng phức hợp bảo vệ bờ biển và thiết lập phương pháp giải tích - bán thực nghiệm định lượng chính xác hệ số truyền sóng ($K_t$).

Khoảng trống nghiên cứu (Research Gap) cốt lõi xuất phát từ thực tiễn: Các công trình đê giảm sóng (ĐGS) truyền thống bằng đá đổ hoặc kè cứng trọng lực áp dụng tại các nước phát triển (Mỹ, Nhật Bản, Châu Âu) khi chuyển giao vào dải bờ biển bùn cát mịn và RNM Tây Nam Bộ thường xuyên bị lún chìm, bồi lấp cửa xả, ngăn cản trao đổi thủy động lực nội - ngoại đê làm chết thảm thực vật RNM, hoặc công trình cọc ly tâm chèn đá hộc ($17\text{ km}$ tại Cà Mau) bộc lộ sự thiếu bền vững và chi phí đắt đỏ. Các nghiên cứu kinh điển thế giới (d'Angremond et al., 1996; Van der Meer et al., 2005; Briganti et al., 2003) chỉ giải quyết bài toán đê ngầm đá đổ đồng nhất hoặc tường cọc đơn lẻ (Hayashi et al., 1966; Wiegel et al., 1960), hoàn toàn khuyết thiếu mô hình toán - lý thủy lực tích hợp đồng thời cơ chế tiêu hao năng lượng của thân đê rỗng ngập nước kết hợp hệ cọc thanh mảnh bên trên.

Luận án tập trung giải quyết 3 câu hỏi nghiên cứu và hệ giả thuyết khoa học tương ứng:

  • RQ1: Cơ chế thủy động lực và mức độ tiêu hao năng lượng sóng ngẫu nhiên tương tác với tổ hợp khối thân đê ngầm rỗng và hệ cọc bên trên diễn ra như thế nào?
    • Hypothesis $H_1$: Tổng năng lượng sóng truyền qua công trình phức hợp là sự tích hợp phi tuyến giữa quá trình sóng vỡ, ma sát dòng chảy qua thân rỗng và công của lực cản thủy động ($q_D$) do hệ cọc hình trụ gây ra.
  • RQ2: Các tham số hình học vô thứ nguyên nào chi phối trực tiếp đến hiệu quả giảm sóng ($K_t$) và có thể chuẩn hóa toán học qua định lý Buckingham hay không?
    • Hypothesis $H_2$: Hệ số truyền sóng $K_t$ là hàm số phụ thuộc vào độ ngập tương đối ($R_c/H_{m0,i}$), bề rộng tương đối ($B/L_m, B/H_{m0}$), chỉ số sóng vỡ Iribarren ($\xi_{0p}$) và diện tích cản thủy lực của hệ cọc ($X_b/L_p$).
  • RQ3: Giải pháp kết cấu đê rỗng phức hợp có đáp ứng triệt để bài toán ổn định công trình trên nền bùn yếu và yêu cầu sinh thái bồi tụ phục hồi RNM theo tinh thần Nghị quyết 120/NQ-CP về phát triển bền vững "thuận thiên"?
    • Hypothesis $H_3$: Khối cấu kiện bê tông rỗng đúc sẵn định hình liên kết cọc ly tâm dự ứng lực vừa giảm áp lực truyền xuống nền móng bùn sét, vừa duy trì độ lưu thông nước đạt 14% diện tích mái đê, đảm bảo giảm chiều cao sóng $H_s$ xuống dưới ngưỡng xói lở cho phép RNM tái sinh.

Khung lý thuyết (Theoretical Framework) của luận án tích hợp:

  1. Lý thuyết phân tích thứ nguyên ($\Pi$-theorem của Buckingham).
  2. Định luật tương tự thủy lực Froude cho mô hình vật lý (MHVL).
  3. Phương trình bảo toàn và cân bằng năng lượng sóng ngẫu nhiên (Spectral Wave Energy Balance).
  4. Lý thuyết sóng phi tuyến và công thức lực cản Morison trên hệ cọc trụ tròn.

Ý nghĩa và quy mô nghiên cứu: Dữ liệu thực nghiệm được xây dựng từ $260$ kịch bản thử nghiệm trên máng sóng thủy lực HR Wallingford (Anh) tại 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ích hợp số liệu khí tượng hải văn $10$ năm (2008–2017) từ $3$ trạm NOAA nước sâu (NOAA1, NOAA2, NOAA3), tính toán lan truyền sóng bằng SWANONE và ứng dụng thực chứng tại công trình thử nghiệm Đề tài KHCN cấp Quốc gia mã số ĐTĐL.CN-09/17 tại cống Kênh Mới, xã Khánh Hải, huyện Trần Văn Thời, tỉnh Cà Mau.

Literature Review và Positioning

Các công trình nghiên cứu về ĐGS bảo vệ bờ đã trải qua hơn 7 thập kỷ phát triển với các cột mốc lý thuyết quan trọng:

+-----------------------------------------------------------------------------+
| TIẾN TRÌNH NGHIÊN CỨU ĐÊ GIẢM SÓNG (LITERATURE EVOLUTION)                   |
|                                                                             |
| 1950s-1960s : Johnson (1951), Wiegel (1960), Hayashi (1966), Goda (1969)   |
|               -> Sóng đều, đê đỉnh thấp/cọc đơn lẻ, ma sát & co hẹp dòng    |
| 1980s-1990s : Ahrens (1987), Van der Meer (1990), d'Angremond et al. (1996)|
|               -> Sóng ngẫu nhiên, phổ JONSWAP, chỉ số Iribarren, đê đá đổ   |
| 2000s-2010s : DELOS Project (2002), Calabrese (2002), Briganti et al. (2003)|
|               -> Hệ cơ sở dữ liệu 3.327 thí nghiệm, đê đỉnh rộng B/Hs > 10 |
| 2021 (Luận án): Nguyễn Anh Tiến (2021)                                      |
|               -> ĐỘT PHÁ: Đê rỗng phức hợp (Đê ngầm rỗng m=1 + Cọc ly tâm)   |
+-----------------------------------------------------------------------------+
  1. Giai đoạn nền tảng sóng đều và đê ngập nước (1951–1980): Johnson et al. (1951) lần đầu tiên ứng dụng lý thuyết dòng năng lượng sóng để đề xuất công thức tính $K_t$ qua dải ngầm ngập nước phụ thuộc độ sâu ngập $R_c$ và bước sóng. Goda et al. (1969) phát triển phương trình thực nghiệm cho thùng chìm thẳng đứng dựa trên tỷ số $R_c/H_s$, khẳng định $K_t$ nghịch biến với chiều rộng đỉnh $B$. Đối với công trình cọc, Morison et al. (1950) thiết lập phương trình lực cản thủy động $q_D = \frac{1}{2} C_D \rho \varnothing u |u|$, làm nền tảng cho Wiegel (1960, 1961) và Hayashi et al. (1966, 1968) nghiên cứu truyền sóng qua hàng cọc dựa trên hiệu ứng co hẹp tia nước qua khe hở $b_i/\varnothing$.
  2. Giai đoạn chuẩn hóa sóng ngẫu nhiên và phổ năng lượng (1980–2000): Allsop (1983), Ahrens (1987) bổ sung thông số hình học diện tích mặt cắt đê tương đối ($A/(h \cdot L_p)$) và đường kính hạt $D_{50}$. Bước đột phá thuộc về d'Angremond et al. (1996) khi tổng hợp dữ liệu thực nghiệm, đưa tham số độ dốc sóng và số Iribarren $\xi_{0p} = \tan\alpha / \sqrt{s_{0p}}$ vào phương trình: $$K_t = -0,40 \frac{R_c}{H_{s,i}} + 0,64 \left(\frac{B}{H_{s,i}}\right)^{-0,31} \left(1 - e^{-0,5 \xi_{0p}}\right)$$ Seabrook & Hall (1998) với $633$ thí nghiệm máng sóng đã chứng minh độ thấm hạt lõi đá $D_{n50}$ ảnh hưởng không đáng kể so với chiều rộng đỉnh $B/H_{s,i}$.
  3. Giai đoạn dự án quy mô lớn DELOS và cấu kiện phức hợp (2000–nay): Dự án DELOS (Châu Âu, 1998–2002) với $18$ viện nghiên cứu đối tác (TU Delft, UCA Cantabria, UPC, Aalborg) đã chuẩn hóa cơ sở dữ liệu $3.327$ thí nghiệm máng sóng (Goda & Ahrens, 2008). Briganti et al. (2003) và Van der Meer et al. (2005) xác lập ranh giới áp dụng cho đê đỉnh hẹp ($B/H_s < 8$) và đê đỉnh rất rộng ($B/H_s > 10$).

Tranh luận khoa học cốt lõi:

  • Trường phái thứ nhất (Wiegel, 1960; Hayashi, 1968) cho rằng khoảng hở tương đối giữa các cọc ($b_i/\varnothing$) là tham số thống trị việc tiêu tán sóng qua hệ cọc.
  • Trường phái thứ hai (Truitt & Herbich, 1987; Herbich & Douglas, 1988; Sonia et al., 2015) phản biện rằng độ sâu ngập tương đối ($d/H_s$) và chu kỳ sóng $T_p$ mới là yếu tố chi phối năng lượng sóng truyền qua hàng cọc kép.

Định vị học thuật của luận án: Luận án của NCS. Nguyễn Anh Tiến đã giải quyết điểm xung đột lý thuyết này bằng cách tích hợp mô hình truyền sóng 2 giai đoạn: Giai đoạn 1 – Sóng tương tác với thân đê ngầm rỗng có mái dốc nghiêng ($m=1$) với độ rỗng mái $14%$; Giai đoạn 2 – Dòng sóng sau đê tiếp tục bị tiêu hao qua hệ cọc thanh mảnh bố trí so le hình hoa mai ($b_i = l_i = \varnothing = 0,3\text{ m}$).

So sánh quốc tế:

  • So với các công trình đê ngầm bãi rạn Aquareef tại Nhật Bản (Hirose et al., 2002 với $1.063$ thí nghiệm) vốn đòi hỏi lượng bê tông khối đặc khổng lồ, giải pháp đê rỗng phức hợp của Việt Nam giảm tới $40%$ trọng lượng công trình, thích ứng hoàn hảo cho địa tầng bùn lỏng Tây Nam Bộ.
  • So với các giải pháp cọc ngập nước nghiên cứu tại Đại học Cantabria (Tây Ban Nha - UCA) trong dự án DELOS, cấu trúc đê rỗng phức hợp có khả năng tiêu tán đồng thời sóng ngắn chu kỳ cao và kiểm soát sóng dài chu kỳ thấp nhờ cấu trúc tán xạ 3 chiề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 năng lượng sóng ngẫu nhiên của d'Angremond et al. (1996) và Van der Meer et al. (2005), đồng thời phát triển lý thuyết lực cản thủy động trên mảng cọc của Morison (1950) và Hayashi (1968) trong môi trường sóng chuyển tiếp nước nông.

+-----------------------------------------------------------------------------+
| MÔ HÌNH CÂN BẰNG NĂNG LƯỢNG SÓNG TỔNG QUÁT (THEORETICAL MODEL)              |
|                                                                             |
|  Năng lượng sóng tới (Ei)                                                   |
|                                                                             |
|  PHƯƠNG TRÌNH TỔNG QUÁT:  Kt = m * sqrt( (Kt,0)^2 - Dpr )                  |
|  (Trong đó: m = 0,94; Kt,0: Đê ngầm rỗng; Dpr: Năng lượng tiêu hao do cọc)  |
+-----------------------------------------------------------------------------+

Mô hình lý thuyết được lượng hóa thông qua phương trình cân bằng năng lượng tổng quát: $$E_i = E_r + E_t + E_d$$ Chia cả hai vế cho năng lượng sóng tới $E_i$, ta có biểu thức hệ số truyền năng lượng: $$K_E = K_t^2 = 1 - C_r^2 - D_p$$ Trong đó:

  • $C_r = H_{m0,r} / H_{m0,i}$ là hệ số sóng phản xạ.
  • $D_p$ là tổng năng lượng sóng tiêu hao qua công trình.

Luận án xác lập bước chuyển dịch mô hình (Paradigm Shift): Tách biệt hàm số tiêu hao năng lượng $D_p$ thành hai thành phần độc lập tuyến tính tương tác: $$K_t = m \cdot \sqrt{K_{t,0}^2 - D_{pr}}$$

  • Mệnh đề 1 (Proposition 1): Hệ số truyền sóng qua thân đê ngầm rỗng không cọc ($K_{t,0}$) chịu sự kiểm soát bởi chỉ số sóng vỡ tương đối $\xi_{0m}$, bề rộng tương đối $B/L_m$ và độ sâu ngập $R_c/H_{m0,i}$.
  • Mệnh đề 2 (Proposition 2): Năng lượng sóng tiêu hao do hệ cọc ($D_{pr}$) tỷ lệ thuận với mật độ cản thủy lực, bề rộng hệ cọc tương đối $X_b/L_p$, chỉ số sóng vỡ cục bộ $H_{m0,i}/R_c$ và nghịch biến với độ dốc sóng địa phương $s_p$.
  • Hệ số điều chỉnh mô hình thực nghiệm: Được tối ưu hóa bằng thuật toán hồi quy phi tuyến đa biến trong môi trường MATLAB: $$m = 0,94 \quad (R^2 = 0,935)$$

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

Khung phân tích của luận án tích hợp liên ngành 3 hệ thống lý thuyết:

  1. Lý thuyết phân tích thứ nguyên Buckingham ($\Pi$-Theorem): Thiết lập ma trận thứ nguyên gồm $11$ biến vật lý ($H_{m0,i}, T_p, L_p, h, R_c, B, X_b, \varnothing, b_i, \rho, g$), rút gọn thành 6 nhóm tham số phi thứ nguyên không thứ nguyên chi phối: $\left(\frac{R_c}{H_{m0,i}}, \frac{B}{L_m}, \frac{B}{H_{m0,i}}, \xi_{0p}, \frac{X_b}{L_p}, \frac{H_{m0,i}}{h}\right)$.
  2. Cơ học thủy khí phi tuyến: Mô tả tương tác xoáy cuộn và tiêu tán ma sát khi dòng tia chất lỏng chui qua $14%$ diện tích lỗ rỗng mái đê kết hợp sự co hẹp dòng chảy giữa các hàng cọc so le.
  3. Thủy lực bờ biển ứng dụng sinh thái (Nature-based Coastal Engineering): Lấy nguyên lý tiêu tán năng lượng nhiều tầng của rễ cây mắm, cây đước ven biển làm cảm hứng kỹ thuật để tối ưu hóa khoảng cách cọc $b_i = l_i = 30\text{ cm}$.

Điều kiện biên (Boundary Conditions):

  • Cao trình đỉnh đê ngầm rỗng luôn ngập dưới mực nước biển tính toán ($R_c < 0$).
  • Đỉnh của hệ cọc luôn nhô cao hơn mực nước đỉnh triều thiết kế ($R_{c,\text{pile}} > 0$).
  • Tỷ số độ ngập tương đối nằm trong giới hạn $-3,0 \le R_c/H_{m0,i} \le 0$.
  • Tỷ số bề rộng đỉnh đê tương đối $0,5 \le B/L_p \le 4,0$.

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 học thực chứng (Positivism Paradigm) kết hợp chủ nghĩa hiện thực phê phán (Critical Realism), lấy dữ liệu quan trắc thực nghiệm vật lý làm bằng chứng tối thượng để kiểm chứng và hiệu chỉnh các mô hình toán học giải tích.

+-------------------------------------------------------------------------+
| QUY TRÌNH THỰC NGHIỆM MHVL (HR WALLINGFORD WAVE FLUME)                  |
|                                                                         |
|  [Hệ tạo sóng Piston]                                                   |
|        (~) Phổ JONSWAP (γ = 3.3)                                        |
|        [WG1] -> Đo sóng tới ban đầu (Hm0,i)                             |
|        [WG2 - WG3 - WG4] -> Tách sóng phản xạ Mansard & Funke (Cr)      |
|     |=========|  [ĐÊ RỖNG PHỨC HỢP THÍ NGHIỆM]                          |
|        [WG5 - WG6] -> Đo phổ sóng truyền sau công trình (Hm0,t)         |
|  [Hệ thống tiêu sóng bãi thoải chống phản xạ cuối máng]                 |
+-------------------------------------------------------------------------+

Thiết kế nghiên cứu đa tầng (Multi-level Experimental Design):

  • Cấp độ 1 (Nguyên hình - Prototype): Nghiên cứu chế độ khí tượng - hải văn, địa hình đáy biển, bùn cát và diễn biến xói lở dải bờ biển từ Mũi Cà Mau đến Hà Tiên giai đoạn 1997–2017.
  • Cấp độ 2 (Mô hình vật lý thu nhỏ - Physical Scale Model): Mô hình máng sóng 2D được thiết kế theo định luật tương tự Froude với tỷ lệ hình học $1:N_L$.
    • Tỷ lệ chiều dài: $N_L = 10 \div 15$.
    • Tỷ lệ thời gian: $N_T = \sqrt{N_L} = 3,16 \div 3,87$.
    • Tỷ lệ vận tốc: $N_V = \sqrt{N_L}$.
    • Tỷ lệ lực: $N_F = N_L^3$.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu thực nghiệm:

  1. Thiết bị máng sóng: Máng sóng thủy lực chuyên dụng HR Wallingford (Anh Quốc) có chiều dài $40\text{ m}$, rộng $1,2\text{ m}$, sâu $1,5\text{ m}$, trang bị hệ thống piston tạo sóng ngẫu nhiên điều khiển số hóa, có khả năng mô phỏng chính xác phổ sóng JONSWAP chuẩn ($\gamma = 3,3$) và Bretschneider.
  2. Hệ thống cảm biến đo sóng: Cụm 6 đầu đo sóng điện dung cao tốc (Wave Gauges - WG) độ phân giải $100\text{ Hz}$.
    • Cụm $WG_1, WG_2, WG_3$: Đặt trước công trình để phân tách thành phần sóng tới ($H_{m0,i}$) và sóng phản xạ ($H_{m0,r}$) theo phương pháp tách sóng 3 điểm của Mansard & Funke (1980).
    • Cụm $WG_4, WG_5, WG_6$: Đặt phía sau công trình tại các khoảng cách $1L, 2L, 3L$ để đo sóng truyền ($H_{m0,t}$), triệt tiêu nhiễu động dòng cục bộ.
  3. Quy trình kiểm định độ tin cậy (Validity & Triangulation):
    • Construct Validity: Kiểm tra độ lặp lại ($3$ lần đo cho mỗi kịch bản sóng), hệ số biến thiên đo lường $CV < 2,5%$.
    • Method Triangulation: Đối chiếu kết quả MHVL với mô hình số hóa lan truyền sóng 1D SWANONE và mô hình thủy động lực máng sóng số họ Navier-Stokes IH2-VOF.

Data và phân tích

Ma trận thực nghiệm gồm $260$ thí nghiệm MHVL độc lập được phân bổ chi tiết:

  • Kịch bản sóng không có công trình (đo kiểm chuẩn đường cơ sở bãi tự nhiên): $20$ kịch bản.
  • Kịch bản sóng truyền qua thân đê ngầm rỗng không cọc: $120$ kịch bản (biến thiên $R_c = -0,02\text{ m} \div -0,10\text{ m}$; $H_{m0,i} = 0,04\text{ m} \div 0,16\text{ m}$; $T_p = 1,0\text{ s} \div 2,5\text{ s}$; $B = 0,2\text{ m} \div 0,8\text{ m}$).
  • Kịch bản sóng truyền qua đê rỗng phức hợp (thân đê + cọc): $120$ kịch bản (biến thiên số hàng cọc $n_c = 1, 2, 3$ hàng; cọc D300 ly tâm bố trí hoa mai).
Thông số thủy động lực & hình học Ký hiệu vô thứ nguyên Khoảng giá trị thực nghiệm trong MHVL
Chiều cao sóng tới tương đối $H_{m0,i}/h$ $0,15 \div 0,55$
Độ sâu ngập đỉnh đê tương đối $R_c/H_{m0,i}$ $-2,85 \div -0,15$
Bề rộng đỉnh đê tương đối $B/L_p$ $0,05 \div 0,65$
Bề rộng hệ cọc tương đối $X_b/L_p$ $0,02 \div 0,35$
Chỉ số sóng vỡ Iribarren $\xi_{0p}$ $0,25 \div 1,85$
Độ dốc sóng đỉnh phổ $s_p = H_{m0,i}/L_p$ $0,012 \div 0,065$

Phương pháp phân tích dữ liệu:

  • Phân tích phổ năng lượng sóng Fast Fourier Transform (FFT) bằng ngôn ngữ lập trình MATLAB để trích xuất diện tích phổ mô-men bậc không $m_0 = \int_0^\infty S(f)df$, tính $H_{m0} = 4\sqrt{m_0}$.
  • Phân tích hồi quy phi tuyến tính (Non-linear Least Squares Regression) với khoảng tin cậy $95%$ ($CI = 95%$), kiểm định $F$-test và hệ số xác định $R^2$.

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

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

+-----------------------------------------------------------------------------+
| TỔNG HỢP HIỆU QUẢ GIẢM SÓNG (FINDINGS SUMMARY)                              |
|                                                                             |
|  1. Đê ngầm rỗng đơn lẻ (Kt,0)        : Giảm 30% - 45% chiều cao sóng       |
|  2. Bổ sung hệ cọc hoa mai (nc = 3)   : Tăng thêm 25% - 35% hiệu quả        |
|  3. Đê rỗng phức hợp hoàn chỉnh (Kt)  : Triệt tiêu tới 65% - 78% sóng tới   |
|     (Hệ số Kt đạt 0,22 - 0,35; Duy trì trao đổi nước 14% qua lỗ rỗng)       |
+-----------------------------------------------------------------------------+
  1. Hiệu ứng cộng hưởng giảm sóng vượt trội của hệ cọc hoa mai trên đỉnh đê rỗng: Thực nghiệm chứng minh khi bố trí thêm $3$ hàng cọc bê tông ly tâm so le hình hoa mai trên đỉnh đê ngầm rỗng, năng lượng sóng tiêu hao tăng thêm $25% \div 35%$ so với trường hợp chỉ có thân đê ngầm rỗng. Hệ số truyền sóng toàn phần $K_t$ giảm từ mức $0,55 \div 0,70$ (đê ngầm đơn lẻ) xuống chỉ còn $0,22 \div 0,35$ đối với đê rỗng phức hợp (triệt tiêu tới $65% \div 78%$ chiều cao sóng tới).
  2. Quy luật tiêu tán năng lượng phi tuyến $D_{pr}$ của hệ cọc: Năng lượng sóng tiêu hao qua hệ cọc $D_{pr}$ không tuân theo hàm bậc nhất tuyến tính mà tương quan chặt chẽ với lũy thừa của bề rộng hệ cọc tương đối $X_b/L_p$ và nghịch biến với độ dốc sóng địa phương $s_m$: $$D_{pr} = f\left(\frac{X_b}{L_p}, \frac{H_{m0,i}}{R_c}, s_p\right)$$ Khi số hàng cọc tăng từ $n_c = 1$ lên $n_c = 2$, mức giảm sóng tăng rõ rệt ($18%$), nhưng khi tăng từ $n_c = 2$ lên $n_c = 3$, hiệu quả tiêu tán gia tăng thêm $12%$. Việc bố trí vượt quá 3 hàng cọc mang lại hiệu quả biên giảm dần trong khi làm tăng chi phí móng.
  3. Phát hiện phản trực giác về độ ngập sâu $R_c/H_{m0,i}$: Trái ngược với đê đá đổ đặc truyền thống (nơi mà khi triều dâng cao làm $|R_c|$ tăng thì $K_t$ tăng vọt làm mất tác dụng bảo vệ), đê rỗng phức hợp vẫn duy trì hệ số giảm sóng ổn định ($K_t < 0,45$) ngay cả ở đỉnh triều cao nhất. Nguyên nhân: Khi mức nước dâng cao, thân đê ngầm giảm hiệu ứng sóng vỡ nhưng toàn bộ cột nước ngập lại kích hoạt tối đa công của lực cản ma sát và lực cản hình dạng $q_D$ của hệ cọc ly tâm bên trên.
  4. Cơ chế lưu thông nước và bồi tụ bùn cát mịn: Cấu trúc lỗ rỗng $14%$ trên mái đê cho phép trao đổi lưu lượng nước và phù sa mịn liên tục, ngăn chặn hiện tượng tạo hố xói cục bộ chân khay do sóng phản xạ ($C_r < 0,20$) – điều mà các đê đá đổ khối đặc hoàn toàn thất bại tại bờ biển Cà Mau.

Implications đa chiều

  • Về mặt lý thuyết: Mở rộng lý thuyết thủy lực công trình thủy, xác lập phương trình giải tích - bán thực nghiệm chuẩn xác tính $K_t$ cho loại hình đê rỗng phức hợp – một đóng góp chưa từng xuất hiện trong các tài liệu chuẩn quốc tế như Rock Manual (CIRIA, CUR, CETMEF, 2007) hay Coastal Engineering Manual (USACE, 2006).
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn thí nghiệm máng sóng kết hợp giải tích thứ nguyên $\Pi$-Buckingham và hiệu chỉnh mô hình số học trên MATLAB, có thể chuyển giao áp dụng cho các nghiên cứu đê phức hợp tại các vùng biển nhiệt đới khác.
  • Về mặt thực tiễn kỹ thuật: Cung cấp bộ công cụ tính toán thủy lực cho các kỹ sư thiết kế:
    1. Xác định nhanh chiều rộng đỉnh đê $B$.
    2. Lựa chọn số hàng cọc $n_c$.
    3. Tính toán chính xác cao trình đặt đỉnh đê ngầm để đạt mục tiêu $K_t \le 0,35$ (ngưỡng động lực học cho phép cây ngập mặn tiên phong như mắm trắng, đước đôi bén rễ và tái sinh tự nhiên).
  • Về mặt chính sách: Đóng góp luận cứ khoa học trực tiếp cho Chương trình mục tiêu ứng phó BĐKH và Nghị quyết 120/NQ-CP của Chính phủ, giúp tiết kiệm $30% \div 40%$ vốn đầu tư công so với phương án kè cứng đá đổ truyền thống.

Limitations và Future Research

  1. Giới hạn mô hình 2D máng sóng: Thí nghiệm MHVL được thực hiện trong máng sóng 2 chiều với sóng truyền vuông góc tuyệt đối với tuyến đê ($90^\circ$). Hiện tượng sóng truyền xiên góc, sóng khúc xạ, nhiễu xạ không gian 3D tại hai đầu phân đoạn đê chưa được khảo sát trực tiếp trong máng sóng.
  2. Điều kiện dòng chảy bùn cát liên hợp: Thí nghiệm tập trung vào tương tác thủy động lực sóng - công trình, chưa mô phỏng trực tiếp tương tác 3 pha (nước - sóng - bùn cát lơ lửng) và quá trình bồi lắng trầm tích hạt mịn kết dính ($d < 0,005\text{ mm}$) trong phòng lab.
  3. Giới hạn kết cấu cọc: Nghiên cứu giới hạn ở cọc bê tông ly tâm tiết diện tròn $\varnothing 300\text{ mm}$, chưa mở rộng cho các hình thái cọc tiết diện vuông, cọc rỗng đục lỗ tiêu năng hoặc cọc composite.

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

  • Mở rộng thí nghiệm trên bể sóng 3 chiều (Wave Basin) để đánh giá hệ số truyền sóng dưới tác động của phổ sóng đa hướng và góc truyền xiên $\theta = 0^\circ \div 60^\circ$.
  • Ứng dụng mô hình số trị Navier-Stokes 3D đa pha (OpenFOAM / IHFOAM) mô phỏng chi tiết trường ứng suất cắt đáy biển và cơ chế bồi lắng phù sa sau đê.
  • Nghiên cứu suy giảm năng lượng sóng khi kết hợp đê rỗng phức hợp với dải RNM phục hồi theo thời gian (tích hợp động lực học sinh thái - Bio-hydrodynamics).

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

+-----------------------------------------------------------------------------+
| TÁC ĐỘNG VÀ LAN TỎA (IMPACT & INFLUENCE)                                    |
|                                                                             |
|  [Học thuật]     : Công bố Scopus/Springer, chuẩn hóa lý thuyết ĐGS rỗng    |
|  [Kinh tế]       : Giảm 30% - 40% chi phí xây dựng so với kè đá đổ          |
|  [Môi trường]    : Khôi phục hàng ngàn hecta đai rừng ngập mặn Cà Mau       |
|  [Chính sách]    : Cung cấp luận cứ kỹ thuật cho Nghị quyết 120/NQ-CP       |
+-----------------------------------------------------------------------------+
  • Tác động học thuật (Academic Impact): Luận án công bố các bài báo khoa học trên các tạp chí chuyên ngành uy tín trong danh mục Scopus (Springer Nature) và các tạp chí đầu ngành trong nước (Tạp chí Khoa học và Công nghệ Thủy lợi, Tạp chí Nông nghiệp & PTNT, Tạp chí Khoa học và Công nghệ Biển). Kết quả là tài liệu tham khảo nền tảng cho các viện nghiên cứu và trường đại học đào tạo kỹ sư công trình thủy.
  • Tác động công nghiệp và ngành xây dựng (Industry Transformation): Thay đổi phương thức thi công công trình biển Tây Nam Bộ từ đắp đá thủ công sang công nghiệp hóa đúc sẵn (Precast modular concrete). Khối đê ngầm rỗng và cọc ly tâm được chế tạo hàng loạt tại nhà máy, kiểm soát tuyệt đối chất lượng bê tông bền sunfat, lắp ghép nhanh trên biển bằng hệ sà lan cẩu, rút ngắn $50%$ thời gian thi công ngoài hiện trường.
  • Tác động chính sách và xã hội (Policy & Societal Impact): Giải pháp đê rỗng phức hợp bảo vệ an toàn cho hệ thống đê biển Tây, ngăn chặn sạt lở bờ biển tỉnh Cà Mau và Kiên Giang, bảo vệ trực tiếp tính mạng, tài sản và sinh kế nuôi tôm - cua sinh thái dưới tán rừng của hàng trăm ngàn hộ dân ven biển.

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

  • Nghiên cứu sinh & Nhà nghiên cứu trẻ: Tiếp cận phương pháp phân tích thứ nguyên Buckingham kết hợp mô hình vật lý máng sóng chuẩn quốc tế và các mã nguồn tính toán thủy lực trên MATLAB.
  • Giảng viên & Nhà khoa học cao cấp: Bộ dữ liệu thực nghiệm $260$ kịch bản sóng ngẫu nhiên chuẩn hóa làm cơ sở kiểm định các mô hình số thủy động lực mới.
  • Kỹ sư tư vấn thiết kế & Ban quản lý dự án hàng hải/thủy lợi: Sở Nông nghiệp và PTNT các tỉnh ĐBSCL (Cà Mau, Kiên Giang, Bạc Liêu, Sóc Trăng) có cẩm nang và công thức giải tích tính toán nhanh thông số kích thước hình học đê giảm sóng.
  • Nhà hoạch định chính sách cấp Bộ & Chính phủ: Bộ Nông nghiệp và PTNT, Bộ Khoa học và Công nghệ có cơ sở khoa học phê duyệt các dự án đầu tư công trình bảo vệ bờ biển bền vững, chi phí thấp, thân thiện môi trườ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ì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập công thức giải tích - bán thực nghiệm dạng tích hợp bảo toàn năng lượng: $$K_t = m \cdot \sqrt{K_{t,0}^2 - D_{pr}}$$ Công trình đã mở rộng trực tiếp lý thuyết truyền sóng ngẫu nhiên qua đê ngầm của d'Angremond et al. (1996) và Van der Meer et al. (2005) bằng cách lượng hóa tường minh thành phần tiêu hao năng lượng do lực cản thủy động của hệ cọc trụ tròn ($D_{pr}$) phía trên đỉnh đê ngầm rỗng dựa trên nguyên lý lực cản Morison (1950), giải quyết trọn vẹn bài toán tương tác sóng với kết cấu rỗng phức hợp đa tầng mà các nghiên cứu trước đây chưa thực hiện được.

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

So với nghiên cứu của Briganti et al. (2003) (chỉ phân tích thống kê thuần túy trên cơ sở dữ liệu đê đá đổ DELOS mà không tách biệt các cơ chế tiêu tán thành phần) và nghiên cứu của Hayashi et al. (1966, 1968) (chỉ khảo sát sóng đều qua hàng cọc đơn/đôi trên đáy phẳng không có thân đê ngầm):

  • Luận án đã kết hợp phương pháp phân tích thứ nguyên $\Pi$-Buckingham với quy trình thí nghiệm máng sóng sóng ngẫu nhiên phổ JONSWAP phân tách sóng 3 điểm (Mansard & Funke, 1980).
  • Phương pháp luận của luận án cho phép tách bạch độc lập năng lượng sóng tiêu tán do sóng vỡ/lưu thông qua thân đê ngầm rỗng ($K_{t,0}$) và năng lượng tiêu tán do lực ma sát cản dòng của mảng cọc hoa mai ($D_{pr}$), sau đó tích hợp phi tuyến có hiệu chỉnh mô hình ($m=0,94$).
+-----------------------------------------------------------------------------+
| SO SÁNH PHƯƠNG PHÁP NGHIÊN CỨU VỚI QUỐC TẾ                                  |
|                                                                             |
| Nghiên cứu quốc tế:                                                         |
| - Hayashi (1966, 1968) : Sóng đều, đáy phẳng, chỉ xét khe co hẹp giữa cọc   |
| - Briganti et al.(2003): Thống kê thực nghiệm thuần túy đê đá đổ DELOS      |
|                                                                             |
| Luận án Nguyễn Anh Tiến (2021):                                             |
| -> Sóng ngẫu nhiên phổ JONSWAP trên máng sóng HR Wallingford                |
| -> Tách sóng tới/phản xạ 3 điểm (Mansard & Funke, 1980)                     |
| -> Tách biệt cơ chế: Thân đê rỗng (Kt,0) + Mảng cọc Morison (Dpr)           |
| -> Tích hợp phi tuyến chuẩn hóa toán học m = 0,94 (R2 = 0,935)              |
+-----------------------------------------------------------------------------+

3. Phát hiện bất ngờ nhất (counter-intuitive finding) có dữ liệu minh chứng là gì?

Phát hiện bất ngờ nhất là khả năng duy trì hiệu quả giảm sóng ($K_t \le 0,45$) của đê rỗng phức hợp ngay cả khi mực nước triều dâng cao nhất ($R_c/H_{m0,i} \to -2,85$). Trong khi ở đê đá đổ ngầm thông thường, khi độ ngập $R_c$ tăng cao, sóng truyền thẳng qua đỉnh đê làm $K_t \to 0,85 \div 0,95$ (mất hoàn toàn khả năng bảo vệ bờ), thì ở đê rỗng phức hợp, khi mức nước dâng ngập sâu thân đê, vận tốc hạt nước tương tác trực tiếp với toàn bộ chiều cao của hệ cọc ly tâm D300 nhô cao bên trên, kích hoạt công suất tiêu tán năng lượng của lực cản Morison ($q_D \sim u^2$), bù đắp hoàn hảo cho phần suy giảm hiệu ứng sóng vỡ của thân đê ngầm rỗng.

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

Có. Luận án cung cấp đầy đủ và chi tiết toàn bộ quy trình tái lập nghiên cứu:

  • Bản vẽ cấu tạo hình học khối đê ngầm rỗng (mái $m=1$, tỷ lệ lỗ rỗng $14%$) và bố trí hệ cọc hoa mai ($b_i = l_i = \varnothing$).
  • Ma trận $260$ kịch bản thí nghiệm với thông số sóng chi tiết ($H_{m0}, T_p, h, R_c, B, X_b$).
  • Mã nguồn lập trình MATLAB xử lý tín hiệu sóng, thuật toán phân tách phổ FFT và quy trình hồi quy phi tuyến.
  • Hồ sơ thiết kế kỹ thuật công trình thử nghiệm thực tế tại cống Kênh Mới, xã Khánh Hải, huyện Trần Văn Thời, tỉnh Cà Mau.

5. Chương trình nghiên cứu 10 năm được phác thảo như thế nào?

Chương trình 10 năm tập trung vào:

  1. Nâng cấp mô hình từ máng sóng 2D sang bể sóng không gian 3D chịu tác động của triều, sóng xiên và dòng chảy ven bờ.
  2. Ứng dụng trí tuệ nhân tạo (Machine Learning/Deep Neural Networks) huấn luyện trên cơ sở dữ liệu lớn để dự báo tức thời $K_t$ theo thời gian thực.
  3. Giám sát độ bền vật liệu bê tông rỗng polyme/tro bay bền sunfat và độ ổn định xói lở móng công trình sau 10 năm vận hành thực tế tại bờ biển Cà Mau.

Kết luận

  1. Đề xuất giải pháp công trình đột phá: Luận án đã nghiên cứu đề xuất thành công kết cấu Đê rỗng phức hợp dạng lắp ghép mô-đun bê tông đúc sẵn định hình kết hợp hệ cọc ly tâm dự ứng lực D300 bố trí hoa mai, khắc phục triệt để các nhược điểm lún chìm và cản trở môi trường của đê đá đổ truyền thống trên nền đất yếu bờ biển Tây ĐBSCL.
  2. Xác lập công thức tính toán tổng quát có độ chính xác cao: Xây dựng thành công phương pháp tính toán và phương trình bán thực nghiệm tổng quát xác định hệ số truyền sóng: $$K_t = m \cdot \sqrt{K_{t,0}^2 - D_{pr}} \quad (\text{với } m = 0,94)$$ được kiểm chứng qua $260$ thí nghiệm máng sóng thủy lực chuẩn quốc tế với hệ số tương quan cao ($R^2 = 0,935$).
  3. Làm rõ bản chất cơ chế tiêu tán năng lượng thủy động lực: Chứng minh cơ chế giảm sóng tích hợp đa tầng: thân đê ngầm rỗng chịu trách nhiệm gây sóng vỡ và tiêu hao dòng thấm ($14%$ diện tích), trong khi hệ cọc hoa mai bên trên triệt tiêu động năng sóng thông qua công của lực cản hình dạng và ma sát xoáy cuộn.
  4. Mở ra 3 hướng nghiên cứu mới:
    • Thủy lực công trình rỗng phức hợp 3 chiều dưới tác động của sóng xiên và dòng triều liên hợp.
    • Tương tác động lực học giữa công trình đê giảm sóng thấm nước với quá trình bồi tụ trầm tích bùn cát mịn kết dính.
    • Mô hình tối ưu hóa đa mục tiêu kết hợp công trình bảo vệ bờ với phục hồi hệ sinh thái rừng ngập mặn thích ứng BĐKH.
  5. Khả năng ứng dụng thực chứng và di sản khoa học: Kết quả nghiên cứu đã được ứng dụng thiết kế, thi công thực tế tại công trình thử nghiệm cấp Quốc gia ĐTĐL.CN-09/17 tại bờ biển Cà Mau, chứng minh hiệu quả giảm sóng vượt trội, ổn định nền móng tuyệt đối và tạo bãi bồi phục hồi thành công thảm rừng ngập mặn bản địa, đặt nền móng kỹ thuật vững chắc cho chiến lược bảo vệ bờ biển "thuận thiên" bền vững tại đồng bằng sông Cửu Long.