Tổng quan về luận án

Xói lở bờ biển và suy thoái các hệ thống đụn cát ven bờ dưới tác động của biến đổi khí hậu cùng các hiện tượng thời tiết cực đoan đang là thách thức nghiêm trọng đối với an toàn đê kè, cơ sở hạ tầng du lịch và sinh kế cộng đồng ven biển miền Trung Việt Nam. Công trình luận án tiến sĩ kỹ thuật của tác giả Nguyễn Ngọc Thể với tiêu đề "Nghiên cứu ảnh hưởng của nước dâng do sóng đến biến động bãi biển và giải pháp bảo vệ bãi biển Cửa Đại, Hội An" (Chuyên ngành: Kỹ thuật xây dựng công trình Biển, Mã số: 9580203; Người hướng dẫn: PGS.TS. Trần Thanh Tùng, GS.TS. Nguyễn Trung Việt; Trường Đại học Thủy Lợi, 2021) là một công trình tiên phong giải quyết toàn diện mối quan hệ tương tác thủy - thạch động lực học giữa hiện tượng nước dâng do sóng (NDDS - Wave setup) và hình thái bãi biển cát.

       ┌────────────────────────────────────────────────────────┐
       │                NĂNG LƯỢNG SÓNG NƯỚC SÂU               │
       │    (Bão, Gió mùa Đông Bắc: Hs > 5m, Chu kỳ 5-7s)      │
       └───────────────────────────┬────────────────────────────┘
                                   │ Lan truyền vào vùng nước nông
                                   ▼
       ┌────────────────────────────────────────────────────────┐
       │            HIỆN TƯỢNG SÓNG ĐỔ & BIẾN ĐỔI               │
       │   Ứng suất bức xạ Gradient (dSxx/dx) tạo NDDS mép bờ   │
       └───────────────────────────┬────────────────────────────┘
                                   │ Nâng cao mực nước hiệu dụng (TWP)
                                   ▼
       ┌────────────────────────────────────────────────────────┐
       │             TÁC ĐỘNG VÙNG SÓNG TRÀN (SWASH)            │
       │   Sóng vươn cao tạt sâu chân đụn cát & bãi cao         │
       └───────────────────────────┬────────────────────────────┘
                                   │ Nứt kéo, trượt phẳng & bào mòn
                                   ▼
       ┌────────────────────────────────────────────────────────┐
       │         VẬN CHUYỂN BÙN CÁT NGANG BỜ RA KHƠI            │
       │    Hạ thấp cao trình bãi, lấn sâu đường bờ >2-3m/bão   │
       └───────────────────────────┬────────────────────────────┘
                                   │ Cơ sở khoa học đề xuất
                                   ▼
       ┌────────────────────────────────────────────────────────┐
       │            GIẢI PHÁP BẢO VỆ TỔ HỢP BỀN VỮNG            │
       │  Đê ngầm giảm sóng xa bờ + Kè mỏ hàn + Nuôi bãi cát   │
       └────────────────────────────────────────────────────────┘

Về mặt khoa học, khoảng trống nghiên cứu (research gap) lớn nhất tại Việt Nam trước công trình này là việc định lượng thành phần NDDS trong các cơn bão hầu như chỉ dựa vào các công thức kinh nghiệm giải tích hoặc điều tra hồi tố vết ngập tĩnh sau bão, dẫn đến độ tin cậy bị suy giảm nghiêm trọng và bỏ qua động thái thực tế của đới sóng tràn (Swash zone). Như Lê Đình Mẫu (2015) đã khẳng định: "từ những năm 1965 trở lại đây vùng ven biển Cửa Đại, Hội An hiện tượng xói lở luôn chiếm ưu thế hơn hiện tượng bồi tụ và cường độ xói lở ngày càng gia tăng, đường bờ ngày càng lấn sâu vào trong đất liền". Đồng thời, nghiên cứu của Nguyễn Kim Đan và Nguyễn Trung Việt (2016-2017) chỉ ra các nhóm nguyên nhân chính gây xói lở là "do thiếu hụt bùn cát, do thay đổi hướng sóng, do tác động từ xây dựng công trình trong khu vực". Tuy nhiên, cơ chế tác động trực tiếp của NDDS làm gia tăng mực nước cực trị, kích hoạt quá trình phá hủy chân đụn cát trong các đợt bão cực đoan vẫn chưa được định lượng rõ nét bằng số liệu thực nghiệm hiện trường liên tục.

Luận án tập trung giải quyết 3 câu hỏi nghiên cứu và giả thuyết khoa học cốt lõi:

  1. Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để quan trắc và trích xuất chuẩn xác dao động mực nước tức thời và giá trị NDDS thực tế trong bão tại dải sóng tràn ven bờ?
    Giả thuyết 1 (H1): Ứng dụng công nghệ xử lý ảnh số Video-Camera kết hợp hệ thống cọc tiêu chia vạch màu RGB cho phép trích xuất liên tục chuỗi dao động mực nước và tách lọc thành phần NDDS với độ chính xác tương đương các thiết bị đo sóng chuyên dụng.
  2. Câu hỏi nghiên cứu 2 (RQ2): Mức độ đóng góp và cơ chế tác động của NDDS đến biến động sạt lở bãi cao và chân đụn cát diễn ra như thế nào dưới các kịch bản triều khác nhau?
    Giả thuyết 2 (H2): NDDS làm nâng cao đáng kể đường mực nước tổng cộng (Total Water Level - TWP), đưa năng lượng sóng tràn vươn sâu tới chân đụn cát gây mất ổn định trượt mái dốc, khiến cường độ xói lở bãi cao tăng vọt so với điều kiện không xét NDDS.
  3. Câu hỏi nghiên cứu 3 (RQ3): Phương án công trình chỉnh trị bờ biển nào có khả năng tiêu tán năng lượng sóng từ xa, kiểm soát NDDS và tái tạo bãi cát bền vững?
    Giả thuyết 3 (H3): Tổ hợp công trình đê ngầm giảm sóng xa bờ kết hợp kè mỏ hàn và nuôi bãi nhân tạo sẽ tạo ra hiệu ứng bồi tụ dạng Tombolo, triệt tiêu dòng xói ngang bờ và bảo vệ toàn diện bờ biển.

Khung lý thuyết của nghiên cứu được xây dựng trên nền tảng Lý thuyết Ứng suất bức xạ sóng (Radiation Stress Theory) của Longuet-Higgins & Stewart (1962, 1964), lý thuyết động lực học sóng tràn của Masselink & Puleo (2006), và mô hình hình thái động lực học cát bờ biển bờ rời. Phạm vi nghiên cứu bao trùm 7.6 km dải bờ biển phía Bắc Cửa Đại từ cửa sông Thu Bồn đến bãi tắm An Bàng, thu thập và xử lý chuỗi số liệu khí tượng bão trong 42 năm (1979–2020), thực hiện đo đạc liên tục qua nhiều chiến dịch hiện trường giai đoạn 2016–2017 với các cơn bão điển hình như Sarika (2016) và Doksuri (2017).

Literature Review và Positioning

Các nghiên cứu thủy động lực học ven bờ trên thế giới về NDDS đã phát triển mạnh mẽ từ các công trình kinh điển của Longuet-Higgins & Stewart (1962, 1964), chứng minh rằng sự biến thiên của tensor ứng suất bức xạ sóng ($S_{xx}, S_{yy}, S_{xy}$) khi sóng tiến vào vùng nước nông và vỡ đổ chính là nguyên nhân cơ bản dẫn đến hiện tượng nước rút do sóng (Wave setdown) bên ngoài đới sóng vỡ và nước dâng do sóng (Wave setup) bên trong đới sóng vỡ tới mép bờ nước mép cát. Bowen et al. (1968) đã kiểm chứng thực nghiệm trong phòng thí nghiệm sự tương thích tuyệt đối giữa lý thuyết ứng suất bức xạ và phân bố mực nước dâng. Trên cơ sở thực nghiệm hiện trường, Hanslow & Nielsen (1993) đã thiết lập công thức tính toán NDDS tại mép bờ: $$\bar{\eta} = 0.048 \sqrt{H_{0} L_{0}}$$ trong khi Stockdon et al. (2006) tổng hợp dữ liệu từ 10 bãi biển tự nhiên quốc tế để đề xuất công thức tính NDDS phụ thuộc vào độ dốc bãi $\beta_f$ và thông số sóng ngoài khơi: $$\langle \eta \rangle = 0.35 \beta_f \sqrt{H_0 L_0}$$

Trong bức tranh học thuật quốc tế tồn tại hai luồng quan điểm và tranh luận sâu sắc:

  • Luồng quan điểm thứ nhất (Funakoshi et al., 2008 khi tích hợp ADCIRC + SWAN) cho rằng NDDS chỉ đóng góp tỷ trọng thứ yếu, dao động từ 10% đến 15% trong tổng mực nước dâng cực trị trong bão.
  • Luồng quan điểm đối nghịch thứ hai (điển hình là Chen et al., 2008 nghiên cứu bão Katrina tại Hoa Kỳ) chứng minh rằng tại các vùng bờ có địa hình đáy biến đổi dốc cục bộ, hiệu ứng sóng ven bờ và NDDS có thể chiếm tới 80% mực nước dâng cực trị, vượt xa ảnh hưởng của áp suất khí quyển và ứng suất gió thuần túy.

Tại Việt Nam, các công trình nghiên cứu trước đây (Nguyễn Xuân Hiển, 2013; Đỗ Đình Chiến, 2014; Vũ Hải Đăng et al., 2017; Nguyễn Ba Thủy, 2017) đã ứng dụng mô hình tích hợp sóng - triều - nước dâng (SuWAT, Delft3D, SWAN) cho thấy NDDS chiếm khoảng 20%–35% tổng lượng nước dâng bão ven biển Bắc Bộ và Trung Bộ. Tuy nhiên, điểm hạn chế căn bản của các nghiên cứu này là thiếu số liệu thực nghiệm đo đạc liên tục đới sóng tràn để hiệu chỉnh mô hình, đồng thời chưa kết nối được hệ quả của NDDS với quá trình biến đổi hình thái mặt cắt ngang bãi biển.

Về cơ chế bồi xói tại bãi biển Cửa Đại, tồn tại sự tranh luận giữa vai trò của vận chuyển bùn cát dọc bờ (Lê Đình Mẫu et al., 2003, sử dụng GENESIS) và vận chuyển bùn cát ngang bờ. Nghiên cứu của Huỳnh Công Hoài et al. (2018, mô hình TOMAWAC) khẳng định: "xói lở nghiêm trọng bãi biển Cửa Đại, Hội An chủ yếu xảy ra về thời kỳ mùa gió Đông Bắc" với chiều cao sóng cực đại $H_s$ lên tới 4.7 m. Luận án của Nguyễn Ngọc Thể đã định vị chính xác vị trí khoa học của mình: kế thừa nhận định của Lê Đình Mẫu (2012) rằng "cơ chế VCBC ngang bờ trong bão chiếm ưu thế hơn VCBC dọc bờ", từ đó đột phá đi sâu chứng minh NDDS là tác nhân thủy động lực trực tiếp điều khiển cơ chế vận chuyển ngang bờ, làm gia tăng biên độ xâm thực bãi cao trong các pha triều cường kết hợp bão.

So sánh với các nghiên cứu quốc tế tương đồng:

  1. Nghiên cứu bãi biển Woolacombe (Vương quốc Anh) của King et al. (1990): Thiết lập tỷ lệ $0.11 < \bar{\eta} / H_{b,rms} < 0.15$ trên bãi phẳng cát mịn dưới điều kiện sóng bão $H_s > 2$ m kéo dài. Luận án của tác giả Nguyễn Ngọc Thể đã mở rộng kiểm chứng trên bãi biển Cửa Đại có độ dốc thay đổi phức tạp và chịu tác động của chế độ bão nhiệt đới dữ dội.
  2. Nghiên cứu tại Grand Popo (Benin) của Raimundo Ibaceta & Rafael Almar (2018, 2020): Sử dụng công nghệ viễn thám Lidar và Camera video để quan trắc quá trình sóng tràn và tiêu tán năng lượng trên bãi cát dốc. Luận án đã phát triển một quy trình quang học Video-Camera kết hợp cọc tiêu vạch màu thích ứng với điều kiện thời tiết khắc nghiệt trong bão tại Việt Nam, mang lại hiệu quả đo đạc tương đương với chi phí vận hành tối ư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 Lý thuyết Ứng suất bức xạ sóng của Longuet-Higgins & Stewart (1964) thông qua việc làm sáng tỏ cơ chế truyền dẫn động lượng sóng trong điều kiện bãi biển có sự hiện diện của các dạng đụn cát nhân tạo và công trình kè bảo vệ cục bộ. Công trình chỉ ra rằng gradient ứng suất bức xạ $\partial S_{xx}/\partial x$ trong đới sóng tràn không chỉ tạo ra độ nâng mực nước dâng tĩnh đơn thuần mà còn sinh ra dao động sóng chu kỳ dài (Infragravity waves), tác động trực tiếp vào vùng chân đụn cát.

Khung lý thuyết bổ sung cơ chế phá hủy đụn cát ven bờ theo mô hình địa kỹ thuật - thủy lực kết hợp:

  1. Mệnh đề 1 ($P_1$): Hiện tượng NDDS làm nâng cao mặt nước trung bình tại mép bờ, dịch chuyển điểm vỡ sóng và vùng sóng tràn sâu hơn vào phía đất liền từ 15 m đến 35 m so với mực nước triều thiên văn thông thường.
  2. Mệnh đề 2 ($P_2$): Dưới tác động ngâm ướt liên tục và lực kéo thủy động của sóng tràn dâng cao, chân đụn cát bị xói ngầm (toe scour), tạo khe nứt kéo trên đỉnh đụn cát do ứng suất kéo vượt quá giới hạn dính kết của đất cát không bão hòa.
  3. Mệnh đề 3 ($P_3$): Sự mất cân bằng cơ học gây sạt lở trượt phẳng (planar shear failure) từng mảng lớn đụn cát; dòng triều rút và sóng rút (backwash) cuốn trôi toàn bộ lượng cát sạt lở ra phía ngoài lạch sâu, hình thành các đới cồn cát ngầm (sandbars) chắn ngoài khơi.
       ┌────────────────────────────────────────────────────────┐
       │   CƠ CHẾ PHÁ HỦY ĐỤN CÁT DO NƯỚC DÂNG DO SÓNG (NDDS)   │
       └───────────────────────────┬────────────────────────────┘
                                   │
                                   ▼
        + Cao trình đỉnh đụn cát
        |       | \
        |       |   \  <= Khe nứt kéo bề mặt do trọng lực & mất chân
        |       |     \
        |       |       \ [Khối trượt sạt lở]
        |       |         \
        +-------+           \ ===> Cát sạt trượt xuống đáy
        | BÃI BIỂN CAO       \
        +~~~~~~~~~~~~~~~~~~~~~\~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
        |  Mực nước dâng NDDS  \ (Xói lở khoét chân đụn cát)
        +-------------------------------------------------------
        |  Mực nước triều bình thường

Sự chuyển dịch hệ hình (paradigm shift) mà luận án mang lại thể hiện ở việc thay đổi hoàn toàn quan niệm truyền thống: từ việc xem xói lở bãi biển trong bão thuần túy là do "sóng to gió lớn đánh trực diện" sang việc xác định NDDS là nhân tố điều kiện cốt lõi (catalyst) nâng nền nước để sóng bão xâm nhập phá hủy các cấu trúc bãi cao vốn an toàn trong điều kiện thủy triều thông thường.

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

Khung phân tích của luận án là sự tích hợp đồng bộ giữa 3 cấu trúc lý thuyết:

  • Lý thuyết Động lực học sóng phi tuyến và lan truyền sóng phổ (Spectral Wave Action Balance).
  • Lý thuyết Thủy lực học vùng sóng tràn (Swash Zone Hydrodynamics & Groundwater Infiltration).
  • Lý thuyết Tiến hóa hình thái động lực học đáy cát (Process-based Coastal Morphodynamics).

Phương pháp tiếp cận giải tích mới kết hợp thuật toán xử lý ảnh số nhận dạng tín hiệu quang học màu RGB trên các cọc tiêu hiện trường với hệ phương trình nước nông phi tuyến giải bằng mô hình 2D XBEACH.

Các khái niệm then chốt được chuẩn hóa và định nghĩa rõ ràng:

  • Nước dâng do sóng (NDDS / Wave setup - $\bar{\eta}$): Độ gia tăng của mực nước biển trung bình trên mực nước tĩnh bên trong đới sóng vỡ do sự chuyển đổi thông lượng động lượng sóng (gradient ứng suất bức xạ) khi sóng vỡ đổ.
  • Mực nước cực trị tổng cộng (Total Water Level - TWP): Mực nước tức thời xác định bởi tổng đại số của mực triều thiên văn ($Z_{tide}$), nước dâng do bão khí áp/gió ($Z_{surge}$) và nước dâng do sóng ($\bar{\eta}$).
  • Vùng sóng tràn (Swash zone): Đới bờ biển xen kẽ ngập nước và lộ bãi giữa giới hạn chân sóng rút thấp nhất và đỉnh sóng vươn cao nhất trên mặt dốc bãi biển.

Điều kiện biên xác định (Boundary Conditions): Khung phân tích áp dụng chuẩn xác cho bờ biển cát thoải đến dốc trung bình ($\tan \beta = 0.02 - 0.08$), chịu chế độ bán nhật triều không đều biên độ nhỏ ($< 1.5$ m), chịu tác động của bão nhiệt đới có chiều cao sóng nước sâu $H_{s0} > 3$ m và chu kỳ sóng đỉnh $T_p = 6 - 12$ s.

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

Thiết kế nghiên cứu

Luận án tuân thủ nghiêm ngặt quan điểm triết học thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism), lấy số liệu thực nghiệm đo đạc tại chỗ làm thước đo chân lý để kiểm chứng và hiệu chỉnh các mô hình số trị lý thuyết.

Thiết kế nghiên cứu đa tầng (Multi-level design) bao gồm 3 không gian tích hợp:

  1. Tầng vĩ mô (Macro-level): Vùng biển mở nước sâu vịnh Đà Nẵng - Quảng Nam (Lưới tính 1, Lưới tính 2 mô hình SWAN), mô phỏng trường sóng truyền từ ngoài khơi vào vùng nước nông.
  2. Tầng trung mô (Meso-level): Dải ven biển Cửa Đại dài 7.6 km từ Cửa Đại đến An Bàng, lưới tính phi cấu trúc mịn $15 \times 15$ m, mô phỏng thủy thạch động lực học tổng thể.
  3. Tầng vi mô (Micro-level): Các mặt cắt ngang bãi biển đặc thù (MCN-01, MCN-02, bãi tắm KS Agribank, Palm Garden, Tân Mỹ; 4 mặt cắt tính toán chi tiết CD01 đến CD04), tập trung phân tích tương tác sóng - đụn cát trong đới sóng tràn.
       ┌────────────────────────────────────────────────────────┐
       │   CẤU TRÚC THIẾT KẾ ĐA TẦNG (MULTI-LEVEL DESIGN)       │
       └───────────────────────────┬────────────────────────────┘
                                   │
        [TẦNG 1: VĨ MÔ]            ▼
        Mô hình SWAN nước sâu (Lưới 15x15m, NCEP Reanalysis 1988-2017)
                                   │ Cung cấp điều kiện biên sóng ngoài khơi
                                   ▼
        [TẦNG 2: TRUNG MÔ]
        Mô hình XBEACH 2D ven bờ (Dải bờ 7.6km Cửa Đại - An Bàng)
                                   │ Truyền tải năng lượng vào đới sóng tràn
                                   ▼
        [TẦNG 3: VI MÔ]
        4 Mặt cắt chi tiết (CD01-CD04), Quan trắc Camera + Cọc C1-C6
        (Phân tích sạt lở bãi cao, trượt đụn cát & kiểm định mô hình)

Dung lượng mẫu và tiêu chí lựa chọn:

  • Chuỗi thống kê 83 cơn bão trong 42 năm (1979–2020) trong bán kính ảnh hưởng 100–200 km; trong đó bão tập trung cao nhất vào tháng 9 (30.12%) và tháng 10 (30.12%), tổng cộng chiếm 60.24%.
  • Hai chiến dịch đo đạc thủy hải văn đồng bộ mùa khô (mùa gió Tây Nam, tháng 10/2016) và mùa mưa bão (mùa gió Đông Bắc, tháng 3/2017).
  • Chiến dịch quan trắc quang học chuyên sâu trong bão số 10 (Doksuri, 12–15/9/2017) tại bãi biển khách sạn Agribank với hệ thống 6 cọc tiêu liên hoàn ($C_1$ đến $C_6$).

Quy trình nghiên cứu rigorous

Quy trình thu thập dữ liệu hiện trường sử dụng các thiết bị đo đạc đạt chuẩn quốc tế:

  • Đo đạc địa hình: Máy toàn đạc điện tử Leica TS02 với độ chính xác đo cạnh $[3\text{ mm} + 2\text{ ppm} \times D]$ và đo góc $\pm 2''$, kết hợp máy đo sâu hồi âm Garmin đo sâu thủy đạc đáy biển ngoài khơi đến độ sâu $h > 10$ m.
  • Đo đạc hải văn sóng và dòng chảy: Sử dụng đồng thời máy đo sóng - dòng chảy âm thanh Doppler chuyên dụng AWAC và AQUAPRO đặt tại 2 trạm SMS01 ($15^\circ 54'14.77''\text{N}, 108^\circ 22'15.69''\text{E}$, sâu $-6.0$ m) và SMS02 ($15^\circ 52'43.79''\text{N}, 108^\circ 25'13.73''\text{E}$, sâu $-6.0$ m); trạm dòng chảy đới sóng vỡ SCRI ($15^\circ 53'30.6''\text{N}, 108^\circ 21'45.37''\text{E}$, sâu $-1.8$ m).
  • Đo mực nước liên tục: Đầu đo áp lực tự ghi Level TROLL 500 (In-Situ, USA) lắp đặt tại cầu cảng Đồn Biên phòng Cửa Đại (Trạm MN01: $15^\circ 52'33.06''\text{N}, 108^\circ 23'20.81''\text{E}$, sâu $-2.3$ m) với tần suất lấy mẫu 10 phút/lần suốt 2 năm liên tục (12/1/2016 – 4/3/2017).
  • Hệ thống quan trắc quang học Camera - Cọc tiêu: Lắp đặt Video-Camera kỹ thuật số độ phân giải cao kết hợp tuyến cọc tiêu chia vạch màu đối phản RGB đóng sâu vào bãi cát tại đới sóng tràn. Mỗi đợt quan sát (obs) kéo dài 500–700 giây (ghi nhận trên 100 chu kỳ sóng liên tiếp), xử lý ma trận ảnh trích xuất cường độ sáng - màu theo thời gian thực để xác định đường dao động mực nước tức thời và NDDS.
       ┌────────────────────────────────────────────────────────┐
       │   QUY TRÌNH XỬ LÝ HÌNH ẢNH TRÍCH XUẤT NDDS QUA CAMERA  │
       └───────────────────────────┬────────────────────────────┘
                                   │
        [BƯỚC 1: THU NHẬN ẢNH VIDEO]
        Camera ghi hình bãi tắm KS Agribank (500-700s/obs, >100 con sóng)
                                   │
                                   ▼
        [BƯỚC 2: PHÂN TÁCH FRAME & ĐỊNH CHUẨN PIXEL]
        Trích xuất chuỗi frame ảnh, xác định ma trận tỷ lệ độ dài/pixel
        trên cọc tiêu C1 - C6 (Kiểm định sai số mm)
                                   │
                                   ▼
        [BƯỚC 3: PHÂN TÍCH TÍN HIỆU MÀU R-G-B]
        Quét cường độ sáng và dải màu RGB dọc thân cọc để nhận dạng
        ranh giới giao thoa không khí - mặt nước tức thời
                                   │
                                   ▼
        [BƯỚC 4: LỌC TÍN HIỆU & TÁCH THÀNH PHẦN THỦY TRIỀU]
        Lọc bỏ sóng thành phần ngắn chu kỳ, tách thành phần triều thiên văn
        thu được từ trạm Level TROLL 500 (MN01) & trạm hải văn Sơn Trà
                                   │
                                   ▼
        [BƯỚC 5: XÁC LẬP GIÁ TRỊ NDDS THỰC ĐO]
        Thu nhận đường dâng mực nước trung bình do sóng (Wave setup)
        chuẩn xác phục vụ kiểm chuẩn mô hình XBEACH

Tam giác đạc dữ liệu (Data Triangulation) được thực hiện chặt chẽ qua 3 nguồn: Số liệu quan trắc thực địa (In-situ sensors/Camera) $\leftrightarrow$ Mô phỏng số trị (SWAN + XBEACH) $\leftrightarrow$ Số liệu tái phân tích khí tượng toàn cầu (NCEP/SEAMOS/JMA).

Data và phân tích

Đặc trưng trầm tích đáy: Mẫu trầm tích lấy tại 14 vị trí phân bố đại diện cho toàn bộ bãi biển Cửa Đại và đáy biển ven bờ; thành phần hạt đồng nhất là cát thạch anh hạt mịn đến hạt trung, đường kính trung vị $d_{50}$ dao động từ 0.20 mm đến 0.35 mm, đưa trực tiếp vào mô-đun vận chuyển bùn cát không bão hòa của XBEACH.

Các kỹ thuật mô hình hóa nâng cao:

  • Mô hình SWAN (Simulating Waves Nearshore): Giải phương trình cân bằng mật độ tác động sóng phổ trên lưới tính phủ kín vùng biển miền Trung. Các tham số hiệu chỉnh tối ưu đạt được: hệ số tiêu tán sóng vỡ sâu, hệ số ma sát đáy Collins ($C_f = 0.015$). Kết quả kiểm định tại 2 trạm SMS01 và SMS02 trong chiến dịch tháng 3/2017 cho thấy mức độ tương đồng vượt trội: sai số tương đối chiều cao sóng $H_s < 8.5%$, sai số chu kỳ sóng đỉnh $T_p < 6.2%$, sai số hướng sóng $< 7^\circ$.
  • Mô hình XBEACH (eXtreme Beach behavior): Mô phỏng quá trình thủy động lực học phi tuyến trong đới sóng vỡ và biến đổi địa hình đáy bãi cát dưới tác động của bão. Tham số hiệu chỉnh mô hình: hệ số nhớt dòng xoáy ngang facsl $= 1.0$, hệ số sạt trượt mái cát khô dưới nước wetslp $= 0.3$, góc nghỉ của cát khô dryslp $= 1.0$, vận tốc lắng hạt cát $w_s$ hiệu chỉnh theo nhiệt độ nước biển. Kết quả kiểm định mô hình với chuỗi số liệu bão Doksuri (tháng 9/2017) tại mặt cắt KS Agribank đạt chỉ số tương quan $R^2 = 0.92$ đối với độ cao sóng, $R^2 = 0.88$ đối với giá trị NDDS và sai số hình thái đáy bãi sau bão $\text{RMSE} < 0.18$ m.

Hệ thống kịch bản kiểm tra độ nhạy và tính vững chắc (Robustness checks) được thiết kế qua 3 kịch bản lớn:

  • Kịch bản 1 (KB1): Bão tác động kết hợp mực triều trung bình ($+0.0$ m), có xét đầy đủ thành phần NDDS.
  • Kịch bản 2 (KB2): Bão tác động kết hợp mực triều cao ($+0.7$ m, triều cường), có xét đầy đủ thành phần NDDS.
  • Kịch bản 3 (KB3): Bão tác động kết hợp mực triều trung bình ($+0.0$ m), loại bỏ thành phần NDDS (giả định truyền thống).

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

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

Công trình mang lại 5 phát hiện khoa học mang tính đột phá với đầy đủ bằng chứng định lượng từ thực nghiệm và mô hình:

       ┌────────────────────────────────────────────────────────┐
       │   5 PHÁT HIỆN ĐỘT PHÁ CỦA CÔNG TRÌNH NGHIÊN CỨU       │
       └───────────────────────────┬────────────────────────────┘
                                   │
        ├── [PHÁT HIỆN 1] Định lượng chính xác NDDS tại mép bờ Cửa Đại
        │   đạt cực đại 0.65m - 1.15m trong bão lớn (chiếm 30-40% TWP).
        │
        ├── [PHÁT HIỆN 2] NDDS là động lực chính phá hủy đụn cát:
        │   KB1 (có NDDS) làm sâu xói tăng 1.8-2.5 lần, lấn sâu bờ 15-28m
        │   so với KB3 (không xét NDDS).
        │
        ├── [PHÁT HIỆN 3] Dòng bùn cát ngang bờ chiếm ưu thế tuyệt đối
        │   trong bão; dòng chảy đới sóng vỡ đổi hướng theo chu kỳ sóng.
        │
        ├── [PHÁT HIỆN 4] Kè cứng bảo vệ bờ trực tiếp thất bại do sóng tràn
        │   vượt đỉnh kết hợp NDDS gây sập mái bê tông & xói sâu chân kè.
        │
        └── [PHÁT HIỆN 5] Tổ hợp Đê ngầm giảm sóng xa bờ + Kè mỏ hàn
            giảm 60-75% sóng bờ, tạo bãi bồi Tombolo phục hồi >25-40m cát.

  1. Định lượng chuẩn xác giá trị và phân bố NDDS tại vùng biển miền Trung:
    Trong các cơn bão có chiều cao sóng nước sâu $H_s > 4.5 - 5.0$ m (như bão Sarika 2016, bão Doksuri 2017), giá trị NDDS lớn nhất tại mép bờ biến thiên từ $0.65$ m đến $1.15$ m tùy thuộc vào độ dốc đáy biển cục bộ. NDDS đóng góp từ 30% đến 42% vào tổng giá trị mực nước dâng cực trị trong bão, chứng minh rằng bỏ qua NDDS sẽ dẫn đến việc đánh giá thấp nghiêm trọng rủi ro ngập lụt và xói lở.

  2. Làm rõ vai trò kích hoạt của NDDS đối với cơ chế sạt lở bãi cao và chân đụn cát:
    So sánh đối chứng giữa kịch bản KB1 (có NDDS) và KB3 (không xét NDDS) trên các mặt cắt tính toán CD01 đến CD04 đã lượng hóa được thiệt hại:

    • Tại mặt cắt CD04 (bãi biển tự nhiên chưa có công trình, khu vực An Bàng - Tân Mỹ), kịch bản KB1 làm gia tăng chiều sâu xói đáy bãi thêm $0.85 - 1.40$ m và làm dịch chuyển đường bờ lấn sâu vào đất liền từ $15$ m đến $28$ m so với KB3.
    • Khi bão đổ bộ trùng thời kỳ triều cường (Kịch bản KB2), mực nước dâng tổng cộng TWP dâng cao vượt qua toàn bộ thềm bãi phẳng, sóng bão tràn trực diện vào chân đụn cát ven bờ gây sạt lở nghiêm trọng với khối lượng cát bị xói cuốn trôi tăng $210% - 260%$ so với điều kiện triều thông thường.

  3. Bằng chứng thực nghiệm khẳng định cơ chế vận chuyển bùn cát ngang bờ chiếm ưu thế:
    Phân tích số liệu trích xuất từ trạm hải văn SCRI trong đới sóng vỡ cho thấy dòng chảy ven bờ trong suốt thời gian bão hoạt động không hình thành dòng chảy dọc bờ ổn định theo một hướng nhất định mà liên tục đảo chiều thuận nghịch vuông góc với đường bờ theo từng chu kỳ sóng tràn và sóng rút. Toàn bộ lượng cát xói lở từ bãi cao bị vận chuyển cưỡng bức ra ngoài khơi lắng đọng tại độ sâu $-3.5$ m đến $-6.0$ m tạo thành các đê cát ngầm chắn sóng tự nhiên.

  4. Lý giải nguyên nhân phá hủy hàng loạt của các công trình kè bảo vệ trực tiếp:
    Nghiên cứu làm sáng tỏ lý do vì sao hàng loạt công trình kè lát mái bê tông nghiêng, kè bao tải cát trọng lực tại khách sạn Fusion Alya, Vinpearl, Palm Garden bị sụp đổ hoàn toàn trong bão năm 2016 và 2020: Các công trình này được thiết kế dựa trên cao trình mực nước dâng bão thuần túy mà không tính đến độ nâng mực nước bổ sung do NDDS. Khi NDDS xuất hiện, sóng bão không bị tiêu hao năng lượng từ xa mà tràn qua đỉnh kè (wave overtopping), đồng thời dòng nước rút tạo hố xói sâu cục bộ ngay tại chân kè ($> 1.8$ m), gây mất ổn định cơ học dẫn đến sập toàn bộ mái kè.

  5. Xác lập giải pháp công trình chỉnh trị tổ hợp tối ưu:
    Mô phỏng 4 phương án bố trí mặt bằng công trình bảo vệ bờ biển chứng minh phương án đê ngầm giảm sóng xa bờ kết hợp hệ thống kè mỏ hàn ngắn tạo hiệu quả vượt trội: Chiều cao sóng sau khi qua tuyến đê ngầm giảm từ $55%$ đến $72%$, triệt tiêu hoàn toàn hiện tượng NDDS tại mép bờ trong, thúc đẩy quá trình bồi tụ cát tạo thành các dải bãi bồi dạng Tombolo mở rộng đường bờ ra phía biển từ $25$ m đến $40$ m sau 1–2 mùa gió.

Implications đa chiều

  • Đóng góp lý thuyết (Theoretical Advances): Hoàn thiện khung lý thuyết động lực học sóng tràn ven bờ tại các vùng biển chịu chế độ bán nhật triều không đều, cung cấp luận cứ khoa học để hiệu chỉnh các công thức kinh nghiệm tính toán NDDS quốc tế phù hợp với đặc thù hình thái địa mạo bờ biển miền Trung Việt Nam.
  • Đổi mới phương pháp luận (Methodological Innovations): Quy trình công nghệ quan trắc dao động mực nước bão bằng Camera xử lý ảnh số RGB kết hợp cọc tiêu mở ra phương pháp tiếp cận kinh tế, chính xác và có thể nhân rộng để thiết lập mạng lưới trạm giám sát xói lở tự động cho toàn bộ dải bờ biển Việt Nam.
  • Ứng dụng thực tiễn (Practical Applications): Cung cấp bộ thông số thiết kế chuẩn xác (cao trình đỉnh kè, độ sâu chân khay chống xói, khoảng cách đặt đê ngầm xa bờ $150 - 250$ m, chiều sâu đỉnh ngầm $-1.0$ m đến $-1.5$ m) cho các dự án chỉnh trị bờ biển Cửa Đại và các khu nghỉ dưỡng cao cấp tại Hội An.
  • Khuyến nghị chính sách (Policy Recommendations): Đề xuất Bộ Nông nghiệp & PTNT cùng UBND tỉnh Quảng Nam chuyển đổi căn bản chiến lược bảo vệ bờ biển: chấm dứt việc xây dựng các tuyến kè cứng bảo vệ bờ trực tiếp mang tính chắp vá; chuyển sang áp dụng giải pháp công trình mềm giảm sóng từ xa (đê ngầm, mỏ hàn mềm, nuôi bãi nhân tạo) kết hợp quản lý tích hợp vùng bờ (ICZM).
  • Điều kiện khái quát hóa (Generalizability): Kết quả và mô hình tính toán có thể áp dụng trực tiếp cho các bãi biển cát miền Trung có cấu trúc địa mạo tương đồng như bãi biển Mỹ Khê (Đà Nẵng), Cảnh Dương (Huế), Tam Quan (Bình Định) và Nha Trang (Khánh Hòa).

Limitations và Future Research

Luận án thẳng thắn chỉ ra 4 giới hạn nghiên cứu:

  1. Độ phân giải hình ảnh Camera trong thời tiết bão cực đoan: Tầm nhìn của hệ thống Camera quang học bị suy giảm cục bộ trong các pha mưa giông bão lớn có mật độ hạt mưa dày và ban đêm, đòi hỏi phải bổ sung thuật toán lọc nhiễu quang học hoặc camera hồng ngoại trong tương lai.
  2. Số chiều mô hình thủy động lực: Mô hình XBEACH 2D dạng averaged-depth tuy mô phỏng xuất sắc quá trình xói lở ngang bờ nhưng chưa phản ánh đầy đủ cấu trúc dòng xoáy cuộn 3 chiều (3D turbulence) xung quanh đầu kè mỏ hàn.
  3. Tính biến thiên theo mùa của cấp phối bùn cát đáy: Dữ liệu thành phần hạt cát $d_{50}$ được đo đạc đại diện theo 2 mùa chính, chưa theo dõi được liên tục sự biến thiên đường phân bố hạt theo từng giờ trong suốt quá trình bão quét qua.
  4. Chi phí nuôi bãi nhân tạo: Luận án tập trung đánh giá hiệu quả kỹ thuật thủy lực của phương án nuôi bãi kết hợp đê ngầm mà chưa lượng hóa đầy đủ bài toán tối ưu hóa kinh tế - tài chính dài hạn trong chu kỳ 20–30 năm.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Phát triển hệ thống quan trắc AI-Vision kết hợp viễn thám Radar mảng pha và thiết bị Lidar cố định ven bờ để tự động hóa hoàn toàn việc trích xuất NDDS và biến đổi đường bờ thời gian thực.
  • Ứng dụng mô hình thủy động lực học 3D hoàn chỉnh (như OpenFOAM hoặc TELEMAC-3D) để nghiên cứu chi tiết trường ứng suất cắt đáy cục bộ quanh các kết cấu đê ngầm giảm sóng có hình dạng phức tạp.
  • Nghiên cứu tác động liên hợp giữa biến đổi khí hậu, nước biển dâng dài hạn và hiện tượng sụt lún địa chất khu vực đồng bằng hạ lưu sông Thu Bồn đến cân bằng bùn cát bãi biển Cửa Đại tầm nhìn 2050–2100.
  • Thiết lập khung đánh giá kinh tế sinh thái lượng hóa dịch vụ bảo vệ bờ của các rạn san hô và bãi cỏ biển nhân tạo tích hợp trên thân đê ngầm giảm sóng.

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

  • Tác động học thuật (Academic Impact): Công trình là tài liệu tham khảo mẫu mực trong đào tạo sau đại học ngành Kỹ thuật công trình biển, Kỹ thuật tài nguyên nước và Hải dương học; mở ra hướng nghiên cứu thực nghiệm hiện trường chuyên sâu về tương tác sóng - công trình tại Việt Nam.
  • Chuyển đổi ngành kỹ thuật công trình thủy (Industry Transformation): Thay đổi hoàn toàn tư duy thiết kế của các viện quy hoạch, công ty tư vấn xây dựng thủy lợi - thủy điện ven biển: bắt buộc phải đưa chỉ tiêu tính toán NDDS vào tiêu chuẩn thiết kế đê biển và công trình bảo vệ bờ.
  • Ảnh hưởng chính sách công (Policy Influence): Cung cấp luận cứ khoa học trực tiếp phục vụ việc phê duyệt và triển khai Dự án đầu tư xây dựng công trình chống xói lở khẩn cấp bờ biển Cửa Đại (gói thầu đê ngầm giảm sóng và tôn tạo bãi cát trị giá hàng trăm tỷ đồng từ nguồn ngân sách Trung ương và vốn ODA).
  • Lợi ích kinh tế - xã hội và du lịch: Bảo vệ trực tiếp dải tài sản nghỉ dưỡng hàng tỷ USD dọc bờ biển Hội An, hồi sinh bãi tắm Cửa Đại - từng được bầu chọn là một trong những bãi biển đẹp nhất châu Á, bảo vệ sinh kế bền vững cho hàng vạn hộ dân làm dịch vụ du lịch tại địa phương.
  • Tầm vóc quốc tế: Đóng góp bộ dữ liệu thực nghiệm bão nhiệt đới quý giá cho cộng đồng nghiên cứu động lực học bờ biển quốc tế tại khu vực biển Đông và Tây Thái Bình Dương.

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

  • Nghiên cứu sinh và học viên cao học ngành Công trình biển: Tiếp cận quy trình nghiên cứu mẫu mực, phương pháp xử lý ảnh Camera trích xuất mực nước và kỹ thuật ghép nối mô hình SWAN - XBEACH.
  • Các nhà khoa học và chuyên gia Hải dương học: Sử dụng bộ dữ liệu thực nghiệm sóng, dòng chảy, NDDS và biến động địa hình mặt cắt ngang bãi biển Cửa Đại để kiểm chứng các lý thuyết động lực học mới.
  • Kỹ sư tư vấn và các doanh nghiệp xây dựng ven biển: Sở hữu thông số thiết kế chuẩn xác và giải pháp kết cấu đê ngầm mỏ hàn tối ưu để áp dụng cho các dự án bảo vệ bờ kè resort, khu đô thị biển.
  • Cơ quan quản lý nhà nước (Bộ NN&PTNT, UBND tỉnh Quảng Nam, Tổng cục Phòng chống Thiên tai): Có căn cứ khoa học vững chắc để xây dựng quy hoạch không gian biển, ban hành tiêu chuẩn kỹ thuật thiết kế đê kè biển thích ứng với thiên tai cực đoan.

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?
Trả lời: Đóng góp độc đáo nhất là việc mở rộng Lý thuyết Ứng suất bức xạ sóng của Longuet-Higgins & Stewart (1964) vào cơ chế phá hủy bãi cát cao và chân đụn cát trong bão nhiệt đới. Luận án đã lượng hóa được rằng gradient ứng suất bức xạ tạo ra NDDS nâng cao nền nước mép bờ thêm $0.65 - 1.15$ m, làm dịch chuyển đới sóng tràn vươn sâu tới chân đụn cát, kích hoạt trạng thái mất cân bằng cơ học nứt kéo và sạt trượt mái cát khô – một cơ chế trước đây thường bị bỏ qua trong các mô hình xói lở bãi biển thuần túy do sóng.

2. Đổi mới phương pháp luận của luận án so với các nghiên cứu trước đây thể hiện ở điểm nào?
Trả lời: So với các nghiên cứu của Nguyễn Xuân Hiển (2013) và Đỗ Đình Chiến (2014) vốn chỉ dựa vào mô hình số trị thuần túy và kiểm chứng bằng vết nước dâng tĩnh sau bão, luận án đã đột phá xây dựng quy trình thực nghiệm đồng bộ: kết hợp trạm đo áp lực liên tục Level TROLL 500 với công nghệ Video-Camera trích xuất dao động mực nước tức thời trên cọc tiêu RGB trong bão lớn (bão Doksuri 2017). Đây là lần đầu tiên tại Việt Nam số liệu NDDS thực đo trực tiếp tại hiện trường được sử dụng để kiểm định mô hình hình thái XBEACH với độ tin cậy vượt bậc ($R^2 > 0.88$).

3. Phát hiện bất ngờ nhất về mặt số liệu thực nghiệm là gì?
Trả lời: Phát hiện bất ngờ nhất là dòng chảy trong đới sóng vỡ tại bãi biển Cửa Đại trong thời gian bão (ghi nhận tại trạm SCRI) hoàn toàn không tuân theo quy luật dòng chảy dọc bờ truyền thống mà đổi hướng liên tục thuận nghịch theo từng chu kỳ sóng tràn/sóng rút, khẳng định cơ chế vận chuyển bùn cát ngang bờ chiếm tỷ trọng chi phối gần như tuyệt đối ($> 85%$) trong quá trình làm mất bãi cát trong bão.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) chi tiết không?
Trả lời: Có. Luận án cung cấp chi tiết toàn bộ quy trình công nghệ gồm: thông số kỹ thuật camera, bảng tra ma trận pixel/độ dài thực tế trên cọc tiêu (Bảng 2.6, 2.7), sơ đồ khối thuật toán phân tách màu RGB, bộ tham số hiệu chỉnh mô hình sóng SWAN (Bảng 2.9, 2.10) và mô hình hình thái XBEACH (Bảng 2.11, 2.12), cho phép bất kỳ nhóm nghiên cứu nào cũng có thể tái lập và ứng dụng chính xác cho các vùng bờ biển khác.

5. Chương trình nghiên cứu 10 năm tới được phác thảo như thế nào?
Trả lời: Chương trình 10 năm tới tập trung vào: (1) Tự động hóa mạng lưới trạm quan trắc Camera-AI đo đạc NDDS và xói lở bờ biển toàn quốc; (2) Tích hợp mô hình biến đổi khí hậu toàn cầu và dự báo bão thời gian thực vào hệ thống cảnh báo sớm xói lở đụn cát; (3) Hoàn thiện tiêu chuẩn quốc gia (TCVN) về tính toán NDDS trong thiết kế công trình bảo vệ bờ biển.

Kết luận

Luận án tiến sĩ kỹ thuật của tác giả Nguyễn Ngọc Thể đã hoàn thành xuất sắc các mục tiêu đề ra với 5 đóng góp cốt lõi:

  1. Xây dựng và chuẩn hóa thành công công nghệ quan trắc, phân tích NDDS từ hình ảnh Video-Camera kết hợp hệ cọc tiêu chia vạch màu RGB tại hiện trường trong bão cực đoan.
  2. Lượng hóa chính xác biên độ NDDS tại bãi biển Cửa Đại, Hội An đạt từ $0.65$ m đến $1.15$ m, chiếm tỷ trọng $30% - 42%$ tổng mực nước dâng cực trị trong bão.
  3. Làm sáng tỏ bản chất thủy - thạch động lực học của quá trình xói lở bãi cao và sạt trượt chân đụn cát do sự kết hợp nguy hiểm giữa NDDS, triều cường và sóng tràn.
  4. Chứng minh một cách thuyết phục cơ chế vận chuyển bùn cát ngang bờ chiếm ưu thế vượt trội so với vận chuyển dọc bờ trong điều kiện bão gió mùa Đông Bắc.
  5. Đề xuất giải pháp công trình chỉnh trị tổng thể: Đê ngầm giảm sóng xa bờ kết hợp kè mỏ hàn ngắn và nuôi bãi nhân tạo, tạo tiền đề khoa học vững chắc để khôi phục và bảo vệ bền vững bãi biển Cửa Đại, Hội An.

Công trình tạo ra bước chuyển dịch hệ hình quan trọng trong kỹ thuật công trình biển tại Việt Nam, mở ra 3 nhánh nghiên cứu học thuật mới (quan trắc quang học đới sóng tràn, mô hình hóa tương tác NDDS - đụn cát, và công nghệ đê ngầm giảm sóng mềm), để lại di sản khoa học và thực tiễn lâu dài phục vụ sự phát triển kinh tế - xã hội bền vững của các đô thị biển Việt Nam.