Tổng quan về luận án

Bảo vệ và phát triển bền vững vùng bờ các đảo ngoài khơi là một bài toán kỹ thuật ven biển đặc biệt phức tạp, giao thoa giữa động lực học chất lỏng phi tuyến, biến động hình thái địa mạo bãi và kỹ thuật công trình thủy. Đảo Phú Quý (tỉnh Bình Thuận) giữ vị trí tiền tiêu chiến lược về an ninh quốc phòng và kinh tế biển tại vùng biển Nam Trung Bộ. Tuy nhiên, hình thái bờ đảo đang chịu tác động dữ dội từ chế độ sóng gió mùa và bão biển cực đoan, dẫn đến sạt lở nghiêm trọng làm suy giảm quỹ đất tự nhiên hạn hẹp. Luận án tiến sĩ kỹ thuật "Nghiên cứu tương tác sóng - đảo Phú Quý và đề xuất giải pháp giảm sóng quanh đảo" do nghiên cứu sinh Mai Quang Khoát thực hiện tại Trường Đại học Thủy lợi dưới sự hướng dẫn khoa học của PGS.TS. Mai Văn Công và GS.TS. Phạm Ngọc Quý (Chuyên ngành: Kỹ thuật Xây dựng Công trình thủy, Mã số: 9580202, năm 2021) đã giải quyết trọn vẹn bài toán động lực học tương tác sóng - bờ đảo và quy hoạch công trình giảm sóng phân đoạn bảo vệ bờ biển ngoài khơi.

                                  KHUNG PHÂN TÍCH TỔNG THỂ
                                             │
    ┌────────────────────────────────────────┼────────────────────────────────────────┐
    ▼                                        ▼                                        ▼
[Động lực học biển]               [Hình thái - Địa mạo]                   [Kỹ thuật Công trình]
- Sóng ngoài khơi (MIKE 21-SW)   - Giải đoán ảnh Landsat (2000-2020)     - Đê ngầm/nổi phân đoạn
- Tương tác gần bờ (MIKE 21-BW)  - Trắc lượng 24 mặt cắt ngang (MC1-MC24)- Khối phủ Haro & Tetrapod
- Khúc xạ, nhiễu xạ, phản xạ     - Mô hình thống kê DSAS (LRR, EPR)      - Thùng chìm Caisson BTCT

Nghiên cứu giải quyết khoảng trống học thuật then chốt: các công trình bảo vệ bờ biển truyền thống tại Việt Nam chủ yếu tiếp cận theo dạng bảo vệ bị động (kè ốp mái cứng dọc tuyến bờ) hoặc áp dụng cho vùng bờ lục địa có bãi thoải. Những nghiên cứu chuyên sâu về trường sóng biển tương tác 3D (khúc xạ, nhiễu xạ, sóng phản xạ, sóng đứng và sóng cuộn vỡ) bao quanh một đảo cô lập ngoài khơi và cơ chế suy giảm năng lượng sóng qua hệ thống đê phá sóng tách bờ xa bờ còn thiếu vắng cơ sở định lượng.

Nghiên cứu thiết lập hệ thống 3 câu hỏi nghiên cứu ($RQ$) và 3 giả thuyết khoa học ($H$):

  • RQ1: Quy luật phân bố không gian và biến đổi năng lượng của trường sóng quanh đảo Phú Quý dưới tác động của các hướng gió mùa thịnh hành (Đông Bắc - Tây Nam) và bão thiết kế diễn ra như thế nào?
  • RQ2: Động lực tương tác sóng - dòng chảy chi phối tốc độ và không gian bồi xói đường bờ đảo Phú Quý trong chu kỳ 20 năm (2000–2020) theo các cấu trúc địa chất bãi khác nhau ra sao?
  • RQ3: Cấu hình công trình giảm sóng nào (đê nhô, đê ngầm, bố trí liên tục hay phân đoạn gián đoạn) đạt hiệu quả tiêu tán năng lượng sóng tối ưu, vừa ổn định bờ bãi vừa duy trì lưu thông thủy vực?
  • H1: Sự bất đối xứng trong trường sóng quanh đảo chịu sự chi phối áp đảo của hiện tượng hội tụ khúc xạ và giao thoa nhiễu xạ tại các mũi nhô địa hình phía Đông Bắc và Tây Nam.
  • H2: Quá trình xói lở cục bộ đường bờ đảo tương quan nghịch với chỉ số tần suất xuất hiện nước ($WFI$) và phụ thuộc trực tiếp vào mức độ can thiệp của các công trình nhân tạo không đồng bộ.
  • H3: Giải pháp đê giảm sóng phân đoạn ngầm kết hợp nổi có khả năng triệt tiêu trên 50% chiều cao sóng có ý nghĩa ($H_s$) tác động vào bờ, vượt trội so với đê liên tục về mặt trao đổi lưu chất và bảo vệ cảnh quan bãi tắm.

Khung lý thuyết của luận án tích hợp: Lý thuyết sóng tuyến tính Airy (1845), phương trình độ dốc thoải Mild-Slope Equation (Berkhoff, 1972), hệ phương trình Boussinesq mở rộng phi tuyến (Madsen & Sørensen, 1992), và lý thuyết động lực trầm tích vùng sóng vỡ CERC (1984). Phạm vi nghiên cứu bao quát toàn bộ 24 mặt cắt địa hình đại diện quanh chu vi bờ đảo Phú Quý, chuỗi số liệu khí tượng hải văn 40 năm (1979–2019), cùng chiến dịch đo đạc sóng thực địa liên tục từ 18/12/2012 đến 03/01/2013 tại trạm ngoài khơi P1.

Literature Review và Positioning

Động lực học tương tác giữa sóng biển và các cấu trúc đảo cô lập đã được đặt nền móng từ các nghiên cứu kinh điển của Arthur (1946, 1951) khi khảo sát sự khúc xạ sóng quanh các đảo tròn và xác định sự thâm nhập năng lượng sóng vào vùng khuất sóng thông qua lý thuyết phân tích trực giao năng lượng dựa trên nguyên lý Fermat và định luật Snell. Như Harrison (1848) và William Morris Davis từng mô tả trực quan: "vận tốc của phần đầu tiên chạm tới vùng cạn sẽ giảm đi, toàn bộ bánh xe sóng quay tròn, và gãy gần như vuông góc trên bãi biển". Tiếp sau đó, Penney & Price (1944) cùng Putnam & Arthur (1948) đã vận dụng nghiệm nhiễu xạ Sommerfeld trong quang học để giải bài toán sóng nước va đập vào chướng ngại vật bán vô hạn.

Trong trường phái mô hình hóa số trị hiện đại, trường sóng quanh các đảo san hô và đảo núi lửa đã được kiểm chứng thông qua các nghiên cứu quốc tế nổi bật. Chao và cộng sự (2006) khi mô phỏng tương tác sóng nội và sóng mặt tại đảo san hô Dongsha (phía Bắc Biển Đông) bằng mô hình số trị 3 chiều đã chứng minh rằng lực Coriolis và độ dốc sườn đảo tạo ra sự bất đối xứng sâu sắc giữa hai nhánh sóng khúc xạ phía Bắc và phía Nam sau khi hợp nhất ở hạ lưu đảo. Nghiên cứu của Jia và cộng sự (2018) sử dụng ảnh vệ tinh radar khẩu độ tổng hợp tiên tiến (Envisat ASAR và TerraSAR-X) đã xác thực cơ chế tái kết nối của các nhánh sóng nhiễu xạ phía sau đảo san hô, chỉ ra rằng mật độ phân tầng đại dương và dòng nền chi phối mạnh mẽ vận tốc pha và vị trí nút sóng giao thoa.

                           TIẾN TRÌNH TIẾP CẬN LÝ THUYẾT SÓNG - BỜ ĐẢO
┌─────────────────────────┐     ┌─────────────────────────┐     ┌─────────────────────────┐
│ Lý thuyết Cổ điển       │     │ Thống kê & Phổ sóng     │     │ Số trị Phi tuyến        │
│ - Airy (1845)           │ ──► │ - Pierson (1950s)       │ ──► │ - Berkhoff (1972)       │
│ - Arthur (1946, 1951)   │     │ - CERC/SPM (1984)       │     │ - Madsen/MIKE (1990s)   │
│ - Penney & Price (1944) │     │ - Zenkovich (1950)      │     │ - Chao/Jia (2006, 2018) │
└─────────────────────────┘     └─────────────────────────┘     └─────────────────────────┘

Tại Việt Nam, các nghiên cứu sóng ven bờ của Lương Phương Hậu (1992), Nguyễn Văn Cư (1995), Nguyễn Thọ Sáo (1995) và Phạm Văn Ninh (2000) tập trung vào các cửa sông ven bờ lục địa (Nam Triệu, Cửa Đại, Soài Rạp). Đối với hải đảo, các đề tài của Lê Đức Tố (2000, 2005), Trịnh Việt An (1995), Nguyễn Hoa Thịnh và Hoàng Xuân Lượng (2005) bước đầu khảo sát hải văn và xói lở quần đảo Trường Sa hay cụm đảo Vịnh Bắc Bộ nhưng dừng lại ở mức độ thống kê mô tả hoặc giải pháp công trình cục bộ.

Luận án định vị chính xác khoảng trống khoa học: kết hợp chuỗi dữ liệu ảnh viễn thám đa thời gian 20 năm với mô hình phổ sóng phổ biến MIKE 21-SW và mô hình sóng phân giải pha Boussinesq MIKE 21-BW độ phân giải cao để làm sáng tỏ cơ chế thủy động lực tương tác 3D tại một hòn đảo có nền địa chất bazan xen kẹp cát san hô đặc thù như đảo Phú Quý.

Đó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 sóng nước nông bằng việc chuẩn hóa cơ chế suy giảm năng lượng sóng khi truyền qua bãi rạn san hô nông bao quanh đảo núi lửa. Luận án kiểm chứng giới hạn vỡ sóng theo tiêu chí Miche và McCowan ($\gamma = H/h = 0{,}78$), chứng minh trong điều kiện địa hình bãi biển Phú Quý có độ dốc thay đổi đột ngột từ đới nước sâu vào thềm san hô, chỉ số sóng vỡ $\gamma$ biến thiên phi tuyến trong khoảng từ $0{,}4$ đến $1{,}2$.

Nghiên cứu phát triển mô hình truyền sóng phi tuyến tích hợp hiệu ứng tiêu tán năng lượng do sóng vỡ bề mặt (surface roller concept) và ma sát đáy thềm san hô. Kết quả làm rõ hiện tượng sóng đứng (clapotis) hình thành trước các tuyến kè dốc đứng hiện trạng gây xói sâu chân khay, giải thích bằng cơ chế phản xạ sóng toàn phần với vận tốc hạt nước phương đứng đạt cực đại tại bụng sóng.

                           KHUNG PHÂN TÍCH THỦY ĐỘNG LỰC HỌC TƯƠNG TÁC
                     ┌──────────────────────────────────────────────────┐
                     │            Trường sóng ngoài khơi                │
                     │          (Phổ JONSWAP / MIKE 21-SW)             │
                     └────────────────────────┬─────────────────────────┘
                                              │
                    ┌─────────────────────────┴─────────────────────────┐
                    ▼                                                   ▼
      ┌───────────────────────────┐                       ┌───────────────────────────┐
      │   Khúc xạ & Vỡ sóng       │                       │   Nhiễu xạ & Phản xạ      │
      │   - Snell Law / Biến sâu  │                       │   - Vật cản tự nhiên      │
      │   - Tiêu tán ma sát đáy   │                       │   - Tương tác mép đảo     │
      └─────────────┬─────────────┘                       └─────────────┬─────────────┘
                    │                                                   │
                    └─────────────────────────┬─────────────────────────┘
                                              │
                                              ▼
                     ┌──────────────────────────────────────────────────┐
                     │     Hệ thống đê giảm sóng phân đoạn (HT02)       │
                     │     - Tiêu tán năng lượng sóng qua đỉnh đê       │
                     │     - Giảm sóng truyền qua: Kt ≤ 0.45            │
                     │     - Duy trì tuần hoàn dòng ven bờ              │
                     └──────────────────────────────────────────────────┘

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

Khung phân tích của luận án kết hợp đồng bộ 3 trụ cột lý thuyết:

  1. Lý thuyết lan truyền và biến đổi sóng nước sâu vào nước nông (Wave Transformation Theory).
  2. Lý thuyết cân bằng hình thái đường bờ và biến động bãi ven biển dưới tác động trường sóng (Morphodynamic Shoreline Equilibrium Theory).
  3. Lý thuyết tương tác sóng - công trình thủy ber (Wave-Structure Interaction Theory).

Cách tiếp cận này cho phép định lượng hóa hệ số truyền sóng $K_t = H_t / H_i$ qua đê giảm sóng ngầm và nhô đỉnh, từ đó xác lập mối quan hệ hàm số giữa cao trình đỉnh đê $Z_d$, bề rộng đỉnh $B$, độ dốc lưu chất và khả năng ổn định động của đường bờ phía sau công trình. Điều kiện biên của khung phân tích được xác lập chặt chẽ cho vùng biển hở, chịu chế độ nhật triều không đều với biên độ triều trung bình $1{,}2 - 1{,}6\text{ m}$.

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

Thiết kế nghiên cứu

Nghiên cứu vận dụng thế giới quan thực chứng (positivism) kết hợp phương pháp mô hình hóa toán học số trị tiên tiến (deterministic numerical modeling) và quan trắc thực địa. Thiết kế nghiên cứu đa tầng bao gồm: quan trắc viễn thám không gian lớn (toàn đảo), mô phỏng trường sóng khu vực biển rộng (MIKE 21-SW, miền tính $45 \times 40\text{ km}$), mô phỏng phân giải pha chi tiết vùng ven bờ (MIKE 21-BW, lưới so le $5 \times 5\text{ m}$), và trắc đạc hình thái bãi tại 24 mặt cắt ngang thực tế.

QUY TRÌNH MÔ HÌNH HÓA ĐA CẤP (MULTI-SCALE MODELING WORKFLOW)
─────────────────────────────────────────────────────────────────────────────────────────
[Dữ liệu Viễn thám Landsat / GEE] ──► [Phân tích DSAS (2000-2020)] ──► [Bản đồ Bồi - Xói]
                                                                               │
[Số liệu Gió/Sóng 40 năm (1979-2019)] ──► [Mô hình Toàn thể MIKE 21-SW]        │
                                                   │                           │
                                          (Biên sóng chi tiết)                 │
                                                   │                           │
                                                   ▼                           ▼
[Hiệu chỉnh Thực đo Trạm P1 (2012-2013)] ─► [Mô hình Thu hẹp MIKE 21-BW] ──► [Tối ưu Đê HT02]
─────────────────────────────────────────────────────────────────────────────────────────

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu tuân thủ các chuẩn mực quốc tế:

  • Xử lý viễn thám: Sử dụng nền tảng Google Earth Engine (GEE) truy xuất toàn bộ ảnh vệ tinh Landsat 5, 7, 8 trong 20 năm (2000–2020). Áp dụng chỉ số nước biến tính $mNDWI$ và chỉ số tần suất nước ($WFI$) để tự động trích xuất đường mép nước thực tế theo mực nước triều chuẩn.
  • Phân tích diễn biến đường bờ: Sử dụng công cụ DSAS (Digital Shoreline Analysis System) của Cục Khảo sát Địa chất Hoa Kỳ (USGS), thiết lập hàng trăm đường trực giao (transects) cách nhau $20\text{ m}$ vuông góc với đường bờ gốc (baseline) để tính toán tốc độ dịch chuyển đường bờ theo hai phương pháp thống kê: Tốc độ điểm cuối ($EPR$) và Tốc độ hồi quy tuyến tính ($LRR$) với khoảng tin cậy 90% ($LCI_{90}$).
  • Đo đạc kiểm định: Sử dụng trạm đo sóng tự ghi P1 đặt tại tọa độ ngoài khơi đảo Phú Quý đo đạc liên tục từ 18/12/2012 đến 03/01/2013 thu thập phổ sóng thực, chiều cao sóng có ý nghĩa $H_s$, chu kỳ đỉnh phổ $T_p$, hướng sóng tới $\theta$.
                              THÔNG SỐ MÔ HÌNH VÀ DỮ LIỆU ĐẦU VÀO
┌───────────────────────────────┬────────────────────────────────────────────────────────┐
│ Thành phần                    │ Thông số kỹ thuật / Giá trị chi tiết                   │
├───────────────────────────────┼────────────────────────────────────────────────────────┤
│ Chuỗi số liệu sóng thống kê   │ 40 năm (1979–2019) tại trạm Phú Quý                    │
│ Đo đạc thực địa kiểm định P1  │ 18/12/2012 – 03/01/2013 (Đo phổ sóng thực tế)          │
│ Miền tính sóng lớn (SW)       │ MIKE 21-SW, lưới phi cấu trúc tam giác linh hoạt       │
│ Miền tính sóng nhỏ (BW)       │ MIKE 21-BW, lưới so le trực giao Δx = Δy = 5.0 m       │
│ Thuật toán giải số trị        │ ADI (Alternating Direction Implicit) ẩn luân hướng     │
│ Phân tích viễn thám không gian│ Landsat 5, 7, 8 (2000–2020), chỉ số mNDWI & WFI        │
│ Mô hình thống kê đường bờ     │ DSAS v5.0, trắc lượng khoảng cách 20m, LRR, EPR, LCI90 │
│ Trắc diện bãi ven bờ          │ 24 mặt cắt ngang đặc trưng (MC1 đến MC24)              │
└───────────────────────────────┴────────────────────────────────────────────────────────┘

Data và phân tích

Mô hình lan truyền sóng MIKE 21-SW được thiết lập trên lưới tam giác phi cấu trúc bao phủ toàn bộ vùng biển quanh đảo, giải phương trình cân bằng tác động sóng phổ:

$$\frac{\partial N}{\partial t} + \nabla \cdot (\vec{v} N) = \frac{S}{\sigma}$$

Trong đó $N(\vec{x}, \sigma, \theta, t)$ là mật độ tác động sóng, $S$ là các số hạng nguồn (sinh sóng do gió, tương tác phi tuyến 4 sóng quadruplet, tiêu tán do sóng vỡ nước sâu và ma sát đáy).

Mô hình thu hẹp chi tiết vùng quanh đảo sử dụng MIKE 21-BW giải hệ phương trình Boussinesq 2D cải tiến bằng thuật toán ẩn luân hướng ADI (Alternating Direction Implicit):

$$\frac{\partial \zeta}{\partial t} + \frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} = 0$$

$$\frac{\partial P}{\partial t} + \frac{\partial}{\partial x}\left(\frac{P^2}{d}\right) + \frac{\partial}{\partial y}\left(\frac{PQ}{d}\right) + g d \frac{\partial \zeta}{\partial x} + \psi_x = 0$$

$$\frac{\partial Q}{\partial t} + \frac{\partial}{\partial x}\left(\frac{PQ}{d}\right) + \frac{\partial}{\partial y}\left(\frac{Q^2}{d}\right) + g d \frac{\partial \zeta}{\partial y} + \psi_y = 0$$

Trong đó $P, Q$ là thông lượng dòng theo phương $x, y$; $\zeta$ là cao độ mặt nước; $d = h + \zeta$ là tổng độ sâu nước; $\psi_x, \psi_y$ là các số hạng tán sắc bậc cao Boussinesq.

Kết quả kiểm định mô hình MIKE 21-SW với số liệu đo thực tế trạm P1 cho thấy hệ số tương quan $R^2 > 0{,}88$, sai số tuyệt đối trung bình của chiều cao sóng có ý nghĩa $MAE < 0{,}18\text{ m}$, khẳng định độ tin cậy khoa học cao của bộ công cụ mô phỏng.

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

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

Thứ nhất, dữ liệu viễn thám và quan trắc địa hình 24 mặt cắt ngang bóc tách thực trạng xói lở nghiêm trọng trên toàn đảo. Luận án trích dẫn số liệu lịch sử: "Tổng chiều dài sạt lở trên toàn đảo Phú Quý tới 5688m chiếm 1/4 chu vi của đảo", tập trung nặng nề nhất tại xã Tam Thanh (Thôn 4, Thôn 5) và Ngũ Phụng (Bãi Lăng), nơi sóng biển từng xóa sổ hơn 100 ngôi nhà và phá hủy hoàn toàn $2\text{ km}$ đường ven biển.

Thứ hai, mô hình dự báo xâm thực dài hạn cảnh báo mức độ mất đất thảm khốc nếu không có công trình can thiệp quy mô tổng thể: "Sau 100 năm, tổng diện tích đất bị mất ở đảo Phú Quý là 0,81km², chiếm 4,82% tổng diện tích toàn đảo. Nước dâng xâm thực của biển lấn sâu vào đất của đảo Phú Quý tới 47,5m".

                         DỰ BÁO XÂM THỰC BỜ ĐẢO PHÚ QUÝ THEO THỜI GIAN
┌────────────┬─────────────────────────┬─────────────────────────┬────────────────────────┐
│ Thời gian  │ Diện tích đất mất (km²) │ Tỷ lệ diện tích đảo (%) │ Xâm thực lấn sâu bờ (m)│
├────────────┼─────────────────────────┼─────────────────────────┼────────────────────────┤
│ 20 năm     │ 0.067                   │ 0.40%                   │ 5.8 m                  │
│ 30 năm     │ 0.745                   │ 4.43%                   │ 8.5 m                  │
│ 50 năm     │ 0.754                   │ 4.49%                   │ 15.8 m                 │
│ 100 năm    │ 0.810                   │ 4.82%                   │ 47.5 m                 │
└────────────┴─────────────────────────┴─────────────────────────┴────────────────────────┘

Thứ ba, mô phỏng trường sóng bão thiết kế hướng Đông Bắc ($H_s = 2{,}5\text{ m}, T_p = 8\text{ s}$) và Tây Nam (chu kỳ lặp 2% cực đoan) làm sáng tỏ hiện tượng bất đối xứng trường sóng. Do địa hình mũi đá nhô và bãi cạn phía Bắc và Đông Bắc, năng lượng sóng bị khúc xạ hội tụ mạnh mẽ, tạo chiều cao sóng ven bờ tăng cục bộ lên tới $3{,}2 - 3{,}8\text{ m}$ tại mặt cắt MC13 và MC14 trước khi vỡ dạng nhào ($plunging$), gây xói lở chân bãi dữ dội.

    TRƯỜNG SÓNG HIỆN TRẠNG (HT00)             TRƯỜNG SÓNG CÓ ĐÊ PHÂN ĐOẠN (HT02)
   Sóng tới bão ĐB (Hs = 2.5m)                   Sóng tới bão ĐB (Hs = 2.5m)
        │     │     │                                 │     │     │
        ▼     ▼     ▼                                 ▼     ▼     ▼
    ═════════════════════                         ═══════  ═══════  ═══════ (11 đê phân đoạn)
    [Bãi cạn thềm san hô]                         ░░░░░░░  ░░░░░░░  ░░░░░░░ (Vùng giảm sóng)
        │     │     │                                 │        │        │
        ▼     ▼     ▼                                 ▼        ▼        ▼
    Hội tụ sóng vỡ nhào                           Sóng truyền qua: Hs ≤ 1.1m (Giảm > 55%)
    Hs_bờ = 3.2 - 3.8 m                           Duy trì tuần hoàn dòng ven - Ổn định bãi cát
    [BỜ ĐẢO BỊ XÓI SẠT NẶNG]                      [BỜ ĐẢO ĐƯỢC BẢO VỆ TOÀN DIỆN]

Thứ tư, kết quả so sánh giữa hai kịch bản bảo vệ: Hệ thống đê giảm sóng liên tục (HT01) và Hệ thống 11 tuyến đê giảm sóng phân đoạn gián đoạn (HT02). Kết quả mô phỏng số trị tại MC13, MC14 chỉ ra phương án đê phân đoạn HT02 làm giảm chiều cao sóng sau đê từ $2{,}5\text{ m}$ xuống dưới $1{,}1\text{ m}$ (hệ số truyền sóng $K_t \le 0{,}44$), triệt tiêu năng lượng phá hủy của sóng bão mà không gây ứ đọng môi trường nước vùng sau đê như phương án đê liên tục HT01.

Implications đa chiều

  • Lý thuyết: Bổ sung cơ sở khoa học định lượng về quá trình sóng nhiễu xạ qua các khe hở giữa các đoạn đê ngầm và tái cấu trúc phổ sóng vùng nước sau đê quanh đảo.
  • Phương pháp: Thiết lập quy trình chuẩn kết hợp viễn thám GEE/DSAS đa thời gian với chuỗi mô hình số trị MIKE 21 (SW và BW) cho các bài toán công trình thủy hải đảo.
  • Thực tiễn & Chính sách: Đóng góp luận cứ khoa học để UBND tỉnh Bình Thuận và Bộ Nông nghiệp & PTNT chuyển dịch tư duy công trình từ "bảo vệ bị động tại mép bờ" (kè cứng) sang "bảo vệ chủ động từ xa" (đê phá sóng ngầm/nhô phân đoạn ngoài khơi kết hợp nuôi bãi nhân tạo).
                 SO SÁNH CÁC GIẢI PHÁP BẢO VỆ BỜ BIỂN TẠI ĐẢO PHÚ QUÝ
┌──────────────────────┬────────────────────────┬──────────────────────┬────────────────────────┐
│ Tiêu chí đánh giá    │ Kè ốp mái cứng (GĐ1/2) │ Đê liên tục (HT01)   │ Đê phân đoạn (HT02)    │
├──────────────────────┼────────────────────────┼──────────────────────┼────────────────────────┤
│ Cơ chế tác động      │ Bị động tại mép bờ     │ Chủ động từ xa       │ Chủ động tiêu tán      │
│ Hiệu quả giảm sóng bờ│ Kém (gây sóng phản xạ) │ Rất cao (Kt ≤ 0.2)   │ Tối ưu (Kt = 0.35–0.45)│
│ Duy trì bãi tắm      │ Phá hủy hoàn toàn bãi  │ Bồi lấp tĩnh lặng    │ Nuôi bãi, phục hồi bãi │
│ Lưu thông thủy vực   │ Không ảnh hưởng        │ Kém (ô nhiễm bồn cảng)│ Tốt (qua các cửa đê)   │
│ Khả năng tàu bè neo  │ Không                  │ Có (hạn chế luồng)   │ Rất tốt đa hướng       │
│ Chi phí & Cảnh quan  │ Thấp / Xấu cảnh quan   │ Rất cao / Chắn tầm   │ Hợp lý / Giữ cảnh quan │
└──────────────────────┴────────────────────────┴──────────────────────┴────────────────────────┘

Limitations và Future Research

  1. Giới hạn trường sóng 3D và dòng rối: Mô hình MIKE 21-BW giải phương trình Boussinesq 2D tích phân theo chiều sâu, chưa mô phỏng chi tiết cấu trúc dòng rối 3 chiều ($3D\text{ Navier-Stokes}$/LES) và dòng xoáy thẳng đứng cục bộ tại chân khay đê giảm sóng.
  2. Giới hạn số liệu động lực bùn cát san hô: Do trầm tích quanh đảo Phú Quý là hỗn hợp cát thạch anh, mảnh vụn san hô và vỏ sinh vật có dung trọng và hình thái hạt phức tạp, mô hình vận chuyển trầm tích chủ yếu kế thừa công thức vùng bờ cát tiêu chuẩn CERC.
  3. Độ phân giải mô hình địa hình đáy: Độ sâu đáy biển ở một số vùng rạn nông phía Bắc chủ yếu nội suy từ hải đồ và viễn thám, cần được đo sâu quét hồi âm đa tia (Multibeam Echo Sounder) chi tiết hơn trong tương lai.
  4. Hướng nghiên cứu tương lai:
    • Ứng dụng mô hình thủy động lực 3D phân giải cao (OpenFOAM / DualSPHysics) mô phỏng cấu trúc lỗ rỗng của khối phủ Haro và Tetrapod.
    • Nghiên cứu mô hình kết hợp sinh thái rạn san hô nhân tạo (eco-engineering breakwaters) trên thân đê ngầm để gia tăng đa dạng sinh học biển.

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

Luận án mang lại giá trị học thuật và thực tiễn to lớn trong chiến lược biển quốc gia:

  • Tác động học thuật: Thiết lập chuẩn mực phương pháp luận nghiên cứu thủy động lực học sóng quanh hải đảo cô lập cho các viện nghiên cứu chuyên ngành tại Việt Nam.
  • Tác động kinh tế - xã hội: Bảo vệ an toàn tính mạng và tài sản cho hơn 30.000 dân cư huyện đảo Phú Quý; bảo tồn diện tích $17\text{ km}^2$ đất đảo tự nhiên trước nguy cơ mất đất dự báo $4{,}82%$ sau 100 năm.
  • Bảo tồn cảnh quan du lịch: Hệ thống đê giảm sóng phân đoạn giúp khôi phục các bãi cát tự nhiên (Bãi Lăng, Triều Dương), tạo động lực phát triển du lịch sinh thái biển đảo bền vững.
  • An ninh - Quốc phòng: Đảm bảo tĩnh lặng cho âu thuyền tránh trú bão tàu cá công suất lớn đến 1000CV và khu vực cảng quân sự bảo vệ chủ quyền biển đảo Nam Trung Bộ.

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

  • Nghiên cứu sinh & Nhà khoa học chuyên ngành Thủy lực - Công trình biển: Kế thừa bộ cơ sở dữ liệu mô phỏng sóng 40 năm, phương pháp luận tích hợp MIKE 21 và DSAS.
  • Các Ban Quản lý Dự án & Kỹ sư Thiết kế Công trình Thủy: Ứng dụng quy trình tính toán kết cấu đê giảm sóng, kích thước khối phủ Haro, Tetrapod ($3\text{T} - 10\text{T}$) và công nghệ thùng chìm caisson bê tông cốt thép ($20 \times 8 \times 8\text{ m}$).
  • Chính quyền địa phương (UBND tỉnh Bình Thuận, huyện Phú Quý): Sử dụng bản đồ phân vùng rủi ro xói lở bờ biển để quy hoạch hạ tầng đô thị ven biển và dân sinh.
  • Ngư dân và doanh nghiệp khai thác thủy sản: Hưởng lợi trực tiếp từ vùng nước tĩnh lặng phục vụ neo đậu tàu thuyền và phát triển lồng bè nuôi trồng hải sản an toàn.

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 độc đáo nhất là chuẩn hóa cơ chế truyền sóng phi tuyến và chuyển dịch năng lượng sóng bão qua hệ thống bãi rạn san hô bao quanh đảo cô lập ngoài khơi. Luận án mở rộng Lý thuyết tương tác sóng của Arthur (1951) và Penney & Price (1944) bằng việc tích hợp phương trình tán sắc Boussinesq phi tuyến bậc cao với sự tiêu tán do vỡ sóng nhào cục bộ, giải quyết bài toán hội tụ năng lượng sóng tại các mũi nhô địa hình đảo tự nhiên.

2. Điểm mới về phương pháp nghiên cứu so với các công trình quốc tế tiền nhiệm?

So với nghiên cứu thuần túy viễn thám của Jia và cộng sự (2018) hay nghiên cứu mô hình số thủy tĩnh đơn thuần của Chao và cộng sự (2006), luận án đã kết nối hoàn chỉnh chuỗi phương pháp: phân tích lịch sử bồi xói 20 năm bằng viễn thám Landsat/DSAS $\rightarrow$ mô hình sóng phổ SW ngoài khơi $\rightarrow$ mô hình Boussinesq BW phân giải pha tương tác công trình $\rightarrow$ kiểm định bằng số liệu đo đạc trạm P1 thực tế.

                           SO SÁNH PHƯƠNG PHÁP VỚI CÁC CÔNG TRÌNH QUỐC TẾ
┌──────────────────────┬────────────────────────┬──────────────────────┬────────────────────────┐
│ Tiêu chí             │ Chao et al. (2006)     │ Jia et al. (2018)    │ Luận án Mai Quang Khoát│
├──────────────────────┼────────────────────────┼──────────────────────┼────────────────────────┤
│ Đối tượng            │ Đảo san hô Dongsha     │ Đảo Dongsha          │ Đảo Phú Quý            │
│ Phương pháp chính    │ Mô hình số 3D thủy tĩnh│ Viễn thám SAR (ASAR) │ Đa mô hình tích hợp    │
│ Tương tác công trình │ Không có               │ Không có             │ Đê ngầm/nổi phân đoạn  │
│ Dữ liệu kiểm định    │ Ảnh vệ tinh định tính  │ Vệ tinh TerraSAR-X   │ Trạm đo phổ sóng P1    │
│ Ứng dụng thực tiễn   │ Lý thuyết thuần túy    │ Động lực học sóng nội│ Quy hoạch công trình bờ│
└──────────────────────┴────────────────────────┴──────────────────────┴────────────────────────┘

3. Phát hiện nào gây bất ngờ nhất và có bằng chứng định lượng ra sao?

Phát hiện bất ngờ nhất là các công trình kè ốp mái nghiêng bảo vệ trực tiếp bờ giai đoạn 1 (tại Thôn 5 Tam Thanh và Bãi Lăng) tuy bảo vệ được khu dân cư ngay sau kè nhưng lại gia tăng xói lở cục bộ ở hai đầu tuyến kè và làm hạ thấp cao trình bãi cát trước kè tới $0{,}8 - 1{,}2\text{ m}$ do hiện tượng sóng đứng phản xạ. Bằng chứng định lượng từ mô hình chỉ ra rằng chỉ có giải pháp đê giảm sóng tách bờ phân đoạn HT02 mới duy trì được hệ số truyền sóng $K_t \approx 0{,}40$, giảm chiều cao sóng từ $2{,}5\text{ m}$ xuống $1{,}0\text{ m}$ tại mặt cắt MC13 và khôi phục bồi tụ bãi biển tự nhiên.

4. Luận án có cung cấp quy trình tái lập (replication protocol) đầy đủ không?

Luận án cung cấp tường minh toàn bộ quy trình thiết lập mô hình toán: thông số bước lưới $\Delta x = \Delta y = 5\text{ m}$, thuật toán ADI, hệ số ma sát Manning đáy biển, điều kiện biên hấp thụ xốp (sponge layer), tọa độ 24 mặt cắt trích xuất sóng và bảng tính kết cấu khối phủ Haro/Tetrapod chi tiết tại Phụ lục 1, 2, 3.

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 nghiên cứu 10 năm tập trung vào 3 hướng trọng tâm: (1) Số hóa mô hình tương tác sóng - đảo theo thời gian thực kết nối hệ thống cảnh báo sớm thiên tai biển thông minh; (2) Phát triển các kết cấu đê tiêu sóng sinh thái (ecologically engineered breakwaters) kích thích san hô bám dính; (3) Khảo sát tương tác sóng - dòng triều - nước biển dâng theo các kịch bản biến đổi khí hậu cực đoan RCP4.5 và RCP8.5.

Kết luận

Luận án "Nghiên cứu tương tác sóng - đảo Phú Quý và đề xuất giải pháp giảm sóng quanh đảo" của NCS. Mai Quang Khoát đã tạo lập một dấu ấn khoa học xuất sắc với 5 đóng góp cốt lõi:

  1. Xác lập quy luật trường sóng đảo cô lập: Làm sáng tỏ bản chất các quá trình khúc xạ, nhiễu xạ, phản xạ sóng và hội tụ năng lượng quanh đảo Phú Quý dưới các hướng gió mùa Đông Bắc và Tây Nam.
  2. Định lượng hóa biến động hình thái bờ đảo 20 năm: Ứng dụng thành công công nghệ viễn thám Landsat/GEE và mô hình DSAS chứng minh quy mô sạt lở toàn đảo đạt $5688\text{ m}$ ($1/4$ chu vi đảo) và cảnh báo nguy cơ mất $0{,}81\text{ km}^2$ đất ($4{,}82%$) sau 100 năm.
  3. Xây dựng chuỗi mô hình số trị đa quy mô: Làm chủ quy trình mô phỏng tương tác sóng từ ngoài khơi (MIKE 21-SW) vào vùng ven bờ phân giải pha (MIKE 21-BW), được hiệu chỉnh tin cậy bằng chuỗi số liệu thực đo trạm P1.
  4. Đề xuất giải pháp công trình đột phá: Thiết kế hoàn chỉnh hệ thống 11 tuyến đê giảm sóng phân đoạn (HT02) ngầm kết hợp nổi, đạt hệ số truyền sóng tối ưu $K_t \le 0{,}44$, vượt trội về khả năng trao đổi lưu chất và bảo vệ bãi biển so với đê liên tục truyền thống.
  5. Chuẩn hóa thiết kế kết cấu công trình biển đảo: Hoàn thiện các thông số mặt cắt ngang đê, cao trình đỉnh $Z_d$, lớp phủ khối Haro, Tetrapod và kết cấu thùng chìm caisson bê tông cốt thép phù hợp với điều kiện thi công ngoài khơi.

Công trình không chỉ mở ra hướng tiếp cận hiện đại trong kỹ thuật xây dựng công trình thủy ven đảo tại Việt Nam mà còn là tài liệu tham khảo mẫu mực cho các chương trình phát triển kinh tế biển gắn liền với bảo vệ an ninh quốc phòng vùng hải đảo Tổ quốc.