Tổng quan nghiên cứu
Bờ biển Việt Nam trải dài khoảng 3.260 km, giữ vị trí chiến lược đặc biệt quan trọng đối với phát triển kinh tế biển và an ninh quốc phòng. Vịnh Nha Trang thuộc tỉnh Khánh Hòa có diện tích tự nhiên khoảng 500 km², được che chắn bởi hệ thống 19 đảo lớn nhỏ, trong đó đảo Hòn Tre có diện tích lớn nhất lên tới 36 km². Vùng biển này thường xuyên chịu tác động phức tạp của các trường sóng ven bờ, gây ra hiện tượng xói lở bãi cát và biến đổi địa mạo đới bờ nghiêm trọng. Khi sóng biển truyền từ khơi vào vùng nước nông có độ sâu từ 5 m đến 10 m, chiều cao sóng đạt ngưỡng giới hạn so với độ sâu đáy dẫn đến hiện tượng sóng đổ vỡ. Vùng sóng tràn (swash zone)—dải đất hẹp nằm giữa mép nước triều kiệt và điểm cao nhất sóng leo lên—là nơi diễn ra quá trình tiêu tán năng lượng rối cực mạnh cùng sự di chuyển bùn cát với vận tốc dòng chảy biến thiên từ 0,01 m/s đến hơn 1,76 m/s.
Luận văn Thạc sĩ Khoa học chuyên ngành Hải dương học của tác giả Trần Văn Mỹ (Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội, năm 2015) tập trung giải quyết bài toán: phân tích cấu trúc trường dòng chảy, ước tính năng lượng rối và ứng dụng mô hình số vỡ đập (Dam-Break) để mô phỏng động lực học sóng sau đới sóng đổ tại bãi biển Nha Trang. Nghiên cứu khai thác dữ liệu từ 2 đợt khảo sát thực địa chuyên sâu vào tháng 5 và tháng 12 năm 2013 kết hợp chuỗi dữ liệu tái phân tích sóng WaveWatch-III giai đoạn 1988–2012 (24 năm), mang lại cơ sở khoa học định lượng chuẩn xác giúp kiểm soát xói lở bờ biển và tối ưu hóa thiết kế công trình bảo vệ bờ.
Cơ sở lý thuyết và phương pháp nghiên cứu
Khung lý thuyết áp dụng
Nghiên cứu được xây dựng trên nền tảng thủy động lực học sóng nước nông và lý thuyết chuyển động rối trong cơ học chất lỏng. Quá trình lan truyền sóng ngoài khơi tuân theo định luật bảo toàn thông lượng năng lượng $P = E \cdot c_g = \text{const}$, trong đó $E$ là năng lượng sóng và $c_g$ là vận tốc nhóm sóng. Khi tiến vào vùng nước nông ven bờ, vận tốc truyền sóng $c_g$ giảm dần làm gia tăng chiều cao sóng cho đến khi đạt tỷ lệ giới hạn và sụp đổ. Động lực học vùng sóng tràn được phân loại theo lý thuyết của Short (1999) gồm 3 dải tần số chính: tần số cao (0,07–0,20 Hz, chu kỳ 5–15 s), tần số trung bình (0,03–0,07 Hz, chu kỳ 15–30 s) và tần số thấp (0,003–0,03 Hz, chu kỳ 30–300 s).
Hình thái bãi biển được phân loại theo mô hình của Wright và Short (1984) gồm: bãi phản xạ (độ dốc lớn, năng lượng sóng tiêu tán thấp), bãi khuếch tán (độ dốc thoai thoải, nhiều bar chắn, năng lượng tiêu tán rất lớn) và bãi trung gian. Quá trình tiêu tán năng lượng do sóng vỡ sinh ra chuyển động rối dị hướng, được mô tả qua phương trình tổng năng lượng rối $k = \frac{1}{2}(\overline{u'^2} + \overline{v'^2} + \overline{w'^2})$ kết hợp định luật phổ năng lượng Kolmogorov-Obukhov. Hệ phương trình nước nông 1 chiều Saint-Venant (dạng ma trận hyperbolic phi tuyến) được thiết lập dựa trên cơ chế tương đồng vật lý giữa sự sụp đổ của bore sóng tràn và bài toán vỡ đập thủy lực (Dam-Break).
Phương pháp nghiên cứu
Nghiên cứu tích hợp phương pháp đo đạc thực nghiệm trường phân giải cao và mô hình hóa toán học số trị hiện đại. Dữ liệu thực nghiệm được thu thập trực tiếp tại bãi biển Nha Trang qua 2 chiến dịch đo đạc vào tháng 5 và tháng 12 năm 2013 thuộc đề tài Khoa học Công nghệ cấp Nhà nước. Tác giả sử dụng thiết bị đo vận tốc âm học Doppler đa điểm Vectrino II ADV (Nortek) với tần số lấy mẫu cực cao 0,015 s (~66,7 Hz), cấu hình 35 cell đo phân bố trên dải vi trắc diện 3,5 cm (độ phân giải 1 mm/cell) đặt cách đáy biển tối thiểu 7,5 cm đến 8,0 cm.
Nhằm khắc phục tình trạng phơi cảm biến trong pha nước rút, phương pháp đồng bộ dữ liệu Vectrino với máy quay video tốc độ cao theo thời gian thực được ứng dụng, lọc bỏ hoàn toàn các chuỗi dữ liệu có tỷ số tín hiệu trên nhiễu (SNR) dưới 20 dB. Ngoài khơi tại độ sâu 10 m, trạm quan trắc sóng AWAC ghi nhận liên tục các thông số phổ sóng, kết hợp bộ dữ liệu tái phân tích 24 năm từ mô hình sóng WaveWatch-III (1988–2012). Về mặt số trị, mô hình thủy lực 1D giải hệ phương trình Saint-Venant bằng phương pháp thể tích hữu hạn Godunov, sử dụng bộ giải Riemann Harten-Lax-Van Leer (HLL), phép tái cấu trúc biến không gian bậc hai MUSCL với hàm giới hạn độ dốc minmod và tích phân thời gian Runge-Kutta bậc hai theo điều kiện ổn định hội tụ Courant-Friedrichs-Lewy ($C = 0,75$).
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Thứ nhất, trường dòng chảy trong vùng sóng tràn bộc lộ tính chất bất đối xứng rõ rệt giữa pha nước dồn lên (uprush) và pha nước rửa ngược (downrush). Điển hình trong đợt quan trắc lúc 08h49 ngày 29/5/2013 (chuỗi 6,22 s), vận tốc pha dồn lên đạt cực đại 1,153 m/s, trong khi vận tốc pha nước rút ra biển lên tới 1,760 m/s (chênh lệch hơn 52,6%). Ngược lại, vào lúc 16h35 cùng ngày (chuỗi 5,10 s), pha nước lên đạt vận tốc 0,622 m/s, cao hơn 24,9% so với pha nước xuống (0,498 m/s).
Thứ hai, thời gian diễn thế của các chu kỳ sóng tràn diễn ra cực nhanh, biến động từ 5,10 s đến 21,44 s, phản ánh sự chiếm ưu thế tuyệt đối của sóng trọng lực bề mặt tần số cao (0,07–0,20 Hz). Trong chuỗi đo 21,44 s lúc 09h40 ngày 30/5/2013, ghi nhận 5 đợt nước dồn lên với vận tốc đỉnh 0,0705 m/s và 9 đợt nước rút với vận tốc đỉnh 0,0985 m/s, minh chứng cho hiện tượng bắt giữ bore sóng (bore capture) khi các đợt sóng liên tiếp tương tác và triệt tiêu lẫn nhau.
Thứ ba, cấu trúc năng lượng rối (TKE) phân tầng mạnh mẽ trong lớp biên sát đáy 3,5 cm. Năng lượng rối đạt giá trị cực đại tại hai thời điểm: khi mũi bore chạm tới giới hạn dồn lên cao nhất trên mái dốc và khi pha nước rút chạm mép nước tĩnh, chứng minh cơ chế khuấy động trầm tích đáy cát diễn ra mãnh liệt nhất ở các điểm đảo chiều pha dòng chảy.
Thứ tư, mô hình Dam-Break 1D cho độ tương thích xuất sắc khi đối sánh dữ liệu thực đo với kết quả mô phỏng tại các trường hợp 1, 3 và 8, tái hiện chính xác dạng truyền sóng bore, phân bố biên độ vận tốc và mức độ suy giảm năng lượng rối qua đới sóng vỡ.
Thảo luận kết quả
Sự bất đối xứng trong chu kỳ sóng tràn bắt nguồn từ sự kết hợp của lực ma sát đáy bãi cát, độ dốc địa hình và tương tác ngấm - thoát nước ngầm (infiltration-exfiltration). Trong pha dồn lên, một phần thể tích nước ngấm vào bãi cát chưa bão hòa, làm giảm lưu lượng nhưng gia tăng độ dốc mặt nước của bore; trong pha rửa ngược, nước ngầm thoát ra kết hợp lực hấp dẫn gia tốc dòng chảy hướng ra biển. Các chuỗi dữ liệu được trực quan hóa chi tiết qua biểu đồ vi trắc diện thời gian thực (micro-profile) và đồ thị phân bố năng lượng rối.
Đồng thời, biểu đồ hoa sóng từ dữ liệu WaveWatch-III 24 năm chỉ rõ hướng sóng Đông Đông Bắc (ENE) chiếm tỷ trọng áp đảo 38,32% với chiều cao cực đại lên tới 5,75 m (trung bình 2,5 m) trong mùa đông (tháng X–IV), trong khi hướng Nam Đông Nam (SSE) chiếm 18,29% với chiều cao sóng từ 0 m đến 1,5 m trong mùa hè (tháng V–IX). Nhóm sóng có độ cao từ 0,5 m đến 1,5 m chiếm tần suất lớn nhất với 58,93%, tiếp theo là nhóm 1,5 m đến 2,5 m chiếm 17,27%. Kết quả đối sánh hoàn toàn tương đồng với các công bố quốc tế của Hughes et al. (1997) và Masselink & Puleo (2006), đồng thời khẳng định việc ứng dụng mô hình vỡ đập là giải pháp tính toán hiệu quả, chính xác cho các bãi biển khuếch tán ven bờ miền Trung Việt Nam.
Đề xuất và khuyến nghị
- Thiết lập mạng lưới quan trắc thủy động lực học bờ biển tự động: Lắp đặt 3 trạm quan trắc tích hợp thiết bị đo sóng đa tầng AWAC và cảm biến Vectrino ADV tần số cao tại vịnh Nha Trang giai đoạn 2026–2028, do Sở Tài nguyên và Môi trường tỉnh Khánh Hòa chủ trì phối hợp cùng Viện Hải dương học nhằm giám sát liên tục diễn biến xói lở đường bờ.
- Ứng dụng chuẩn hóa mô hình toán Dam-Break và Saint-Venant trong quy hoạch đới bờ: Triển khai tích hợp mô hình số trị vào bài toán dự báo bồi xói tại các cửa sông ven biển, giảm sai số ước tính độ xói lở bùn cát xuống dưới 12% trong giai đoạn 2026–2027, do Viện Khoa học Thủy lợi miền Trung và Trường Đại học Khoa học Tự nhiên phối hợp thực hiện.
- Thiết kế công trình bảo vệ bờ biển mềm thích ứng động lực sóng: Ứng dụng giải pháp nuôi bãi nhân tạo kết hợp đê ngầm giảm sóng ngoài đới sóng vỡ tại bãi biển Nha Trang nhằm giảm thiểu 35–40% năng lượng sóng tác động trực tiếp lên bờ cát trước năm 2030, do Ban Quản lý Dự án Đầu tư Xây dựng các Công trình Nông nghiệp và Phát triển Nông thôn tỉnh Khánh Hòa thực hiện.
- Xây dựng cơ sở dữ liệu số vi cấu trúc dòng chảy và năng lượng rối quốc gia: Hoàn thiện hệ thống lưu trữ và chia sẻ dữ liệu thủy động lực học sóng tần số cao (>50 Hz) trước quý IV/2027 dưới sự quản lý của Bộ Khoa học và Công nghệ, phục vụ công tác nghiên cứu và đào tạo chuyên sâu.
Đối tượng nên tham khảo luận văn
- Nhà khoa học và nghiên cứu sinh ngành Hải dương học, Thủy văn biển: Nắm bắt phương pháp xử lý tín hiệu Doppler Vectrino tần số cao (66,7 Hz), kỹ thuật lọc nhiễu SNR và quy trình tính toán phân bố năng lượng rối thực nghiệm tại đới bờ.
- Kỹ sư công trình thủy công và đê điều ven biển: Sử dụng các thông số phân bố vận tốc đỉnh (từ 0,50 m/s đến 1,76 m/s) và mô hình Dam-Break để tính toán tải trọng va đập sóng, tối ưu hóa thiết kế mặt cắt đê kè giảm sóng và bảo vệ bờ biển.
- Cơ quan quản lý tài nguyên và môi trường đới bờ (Sở TN&MT Khánh Hòa, Chi cục Biển và Hải đảo): Áp dụng số liệu thống kê chế độ sóng 24 năm và cơ chế xói lở bãi biển để xây dựng đồ án quy hoạch không gian biển và bảo tồn cảnh quan du lịch vịnh Nha Trang.
- Giảng viên và học viên cao học chuyên ngành Kỹ thuật Biển, Động lực học Thủy khí: Sử dụng tài liệu như một cẩm nang thực hành mẫu về giải thuật thể tích hữu hạn Godunov, sơ đồ HLL Riemann và kỹ thuật đồng bộ video trong nghiên cứu hải văn hiện đại.
Câu hỏi thường gặp
Vùng sóng tràn (swash zone) là gì và đóng vai trò gì trong biến đổi bờ biển?
Vùng sóng tràn là dải bờ biển nằm giữa mép nước triều kiệt và điểm cao nhất mà sóng có thể dồn lên. Đây là khu vực tiêu tán năng lượng sóng mạnh nhất, nơi dòng chảy đạt vận tốc từ 0,01 m/s đến 1,76 m/s, trực tiếp khuấy đảo và vận chuyển bùn cát gây biến đổi địa mạo bãi tắm Nha Trang.
Thiết bị Vectrino II ADV mang lại ưu thế gì trong đo đạc sóng tràn?
Vectrino II ADV của hãng Nortek cho phép đo đạc 3 chiều với tần số siêu cao 0,015 s (~66,7 Hz) trên dải vi trắc diện 3,5 cm gồm 35 cell (1 mm/cell). Thiết bị giúp thu thập chi tiết cấu trúc vận tốc dòng chảy và năng lượng rối sát đáy với độ chính xác vượt trội khi chỉ số SNR lớn hơn 20 dB.
Tại sao mô hình vỡ đập (Dam-Break) lại được áp dụng hiệu quả cho vùng sóng tràn?
Khi sóng biển đổ vỡ ở độ sâu 5–10 m, khối nước sụp đổ và dồn lên bãi biển hình thành bore sóng có động học tương đồng với hiện tượng chất lỏng tràn qua đập vỡ. Mô hình 1D giải phương trình Saint-Venant bằng sơ đồ Godunov giúp mô phỏng chuẩn xác bước sóng và vận tốc lan truyền bore nước.
Hiện tượng bất đối xứng sóng tràn diễn ra như thế nào tại bãi biển Nha Trang?
Dữ liệu thực đo ngày 29/5/2013 cho thấy vận tốc pha nước rút (1,76 m/s) lớn hơn đáng kể so với pha nước dồn lên (1,153 m/s), chu kỳ dao động chỉ từ 5,1 s đến 9,5 s. Sự chênh lệch lưu lượng kết hợp tương tác ngấm nước ngầm tạo ra sự bất đối xứng động lực học chi phối quá trình xói lở.
Chế độ sóng tại vịnh Nha Trang có đặc trưng mùa vụ như thế nào?
Theo số liệu WaveWatch-III 24 năm (1988–2012), hướng sóng thịnh hành mùa đông (tháng X–IV) là Đông Đông Bắc (chiếm 38,32%) với chiều cao trung bình 2,5 m, cực đại 5,75 m; trong khi mùa hè (tháng V–IX) hướng sóng chủ đạo là Nam Đông Nam (chiếm 18,29%) với chiều cao sóng dưới 1,5 m.
Kết luận
- Phân tích toàn diện trường dòng chảy vi mô vùng sóng tràn tại bãi biển Nha Trang với dải vận tốc thực đo biến thiên từ 0,0138 m/s đến 1,760 m/s.
- Xác lập quy trình xử lý và đồng bộ dữ liệu Vectrino II ADV tần số 66,7 Hz với video độ phân giải cao, lọc bỏ chính xác các đoạn tín hiệu có SNR dưới 20 dB.
- Định lượng cấu trúc năng lượng rối (TKE) phân tầng trong lớp biên đáy 3,5 cm, làm sáng tỏ cơ chế phát sinh xoáy cuộn sau đới sóng đổ vỡ.
- Đánh giá chi tiết chế độ sóng 24 năm từ WaveWatch-III với hướng ENE chủ đạo (38,32%) và chiều cao sóng cực đại 5,75 m tác động đến bờ biển Khánh Hòa.
- Kiểm chứng thành công mô hình số trị Dam-Break 1D giải bằng sơ đồ thể tích hữu hạn Godunov và HLL Riemann solver trong mô phỏng bore sóng tràn.
Định hướng tiếp theo tập trung mở rộng mô hình hóa thủy động lực học 2D/3D kết hợp vận chuyển bùn cát trong giai đoạn 2026–2028. Quý độc giả, nhà nghiên cứu và kỹ sư công trình biển hãy tham khảo ngay toàn văn luận văn thạc sĩ này để khai thác trọn vẹn bộ số liệu thực nghiệm và thuật toán mô phỏng giá trị cho các dự án bảo vệ đới bờ!