Tổng quan về luận án

Sự gia tăng áp lực nhân sinh cùng các biến động khí hậu toàn cầu đang đặt các hệ thống bờ biển trước những thách thức suy thoái nghiêm trọng. Theo Ủy ban Liên chính phủ về Biến đổi Khí hậu [IPCC, 2007], khoảng hai phần ba dân số thế giới hiện sinh sống trong phạm vi 200 km ven biển. Tại Việt Nam, với hơn 3.200 km bờ biển chịu tác động mạnh mẽ của chế độ gió mùa và thiên tai biển, hiện tượng xói lở bãi biển diễn ra ngày càng phức tạp tại nhiều khu vực trọng điểm như Hải Hậu (Nam Định), Cửa Tùng (Quảng Trị), Cửa Đại (Hội An) và vịnh Nha Trang (Khánh Hòa). Bãi biển trung tâm bờ tây vịnh Nha Trang – vịnh biển xếp thứ 29 thế giới và thứ 2 châu Á về cảnh quan – đóng vai trò huyết mạch đối với kinh tế - xã hội địa phương nhưng đang đối mặt với tình trạng suy thoái nghiêm trọng: bãi biển bị thu hẹp, đường bờ lùi sâu vào kè đường Trần Phú, điển hình như trắc ngang tại ngã ba Nguyễn Chánh - Trần Phú đường bờ đã tiến sâu vào đất liền 9 m trong giai đoạn 2006 đến tháng 5/2013 và độ dốc đáy biển gần bờ gia tăng lên xấp xỉ 4% [Lê Thanh Bình, 2017].

Khoảng trống nghiên cứu (Research Gap) then chốt trong kỹ thuật bờ biển hiện nay nằm ở sự thiếu hụt các công cụ tính toán có khả năng mô phỏng liên tục từ quy mô ngắn hạn (sự kiện bão, gió mùa), trung hạn (mùa, năm) đến dài hạn (nhiều năm, chu kỳ khí hậu). Hanson (2011) khi phân tích 20 mô hình biến đổi bờ biển trên thế giới đã khẳng định chưa có mô hình đơn lẻ nào mô phỏng chính xác trên mọi quy mô thời gian. Các mô hình quá trình 2D/3D (như Delft3D, MIKE 21/3, SMS) đòi hỏi tài nguyên tính toán khổng lồ, chứa nhiều hệ số thực nghiệm chưa được kiểm chứng đầy đủ và không khả thi cho dự báo dài hạn hàng thập kỷ. Ngược lại, các mô hình 1-đường (One-Line) kinh điển chỉ xét vận chuyển dọc bờ, bỏ qua biến động ngang bờ do bão và dao động mùa; trong khi các mô hình trắc ngang như SBEACH chỉ mô phỏng tốt pha xói do bão mà kém chính xác trong pha bồi tụ tái tạo bãi.

Luận án của nghiên cứu sinh Vũ Công Hữu (2019) với đề tài "Lựa chọn và ứng dụng bộ công cụ nghiên cứu biến động đường bờ khu vực bãi biển Nha Trang" (Chuyên ngành Hải dương học, Mã số: 9440228.01, Đại học Khoa học Tự nhiên – ĐHQGHN, hướng dẫn bởi GS. Đinh Văn Ưu và TS. Lê Xuân Hoàn) đã tiên phong giải quyết khoảng trống này thông qua các câu hỏi và giả thuyết khoa học:

  • Research Question 1 (RQ1): Cơ chế thủy thạch động lực nào chi phối quá trình xói lở và bồi tụ bãi biển trung tâm bờ tây vịnh Nha Trang theo các quy mô thời gian khác nhau?
  • Research Question 2 (RQ2): Làm thế nào để xây dựng một bộ công cụ tích hợp giữa mô hình số trị sóng, mô hình thống kê - thủy động lực tự hiệu chỉnh và công nghệ xử lý ảnh camera giám sát bãi biển nhằm mô phỏng chính xác biến động đường bờ dài hạn?
  • Research Question 3 (RQ3): Quy mô thời gian chi phối dài hạn của hình thái bãi biển Nha Trang là gì khi phân tách các thành phần tuần hoàn và ngẫu nhiên?
  • Hypothesis 1 (H1): Biến động đường bờ bãi biển Nha Trang chịu sự chi phối mang tính chu kỳ đối xứng nghịch giữa mùa gió Đông Bắc (pha xói mạnh) và mùa gió Tây Nam (pha bồi tái tạo), đan xen các biến động phi tuần hoàn do bão và áp thấp nhiệt đới.
  • Hypothesis 2 (H2): Mô hình lai ghép thủy động lực - thống kê dựa trên khái niệm đường bờ cân bằng kết hợp xử lý ảnh camera quang học độ phân giải cao có thể tự hiệu chỉnh tham số và tái hiện chính xác vị trí đường bờ với sai số chuẩn hóa NMSE < 0.35.
  • Hypothesis 3 (H3): Tồn tại các thành phần dao động điều hòa dài hạn (chu kỳ nửa năm, 1 năm và 3 năm) chi phối sự tiến triển hình thái bãi biển bờ tây vịnh Nha Trang.

Khung lý thuyết của nghiên cứu tích hợp lý thuyết trắc ngang cân bằng (Equilibrium Beach Profile) của Bruun (1954), Dean (1977) và Kristensen (2013a); lý thuyết biến đổi đường bờ cân bằng động lực học của Miller và Dean (2004a), Yates và cộng sự (2009), Davidson và cộng sự (2010, ShoreFor); mô hình lan truyền sóng khúc xạ - nhiễu xạ EBED; cùng kỹ thuật quang trắc viễn thám bãi biển. Đóng góp đột phá của luận án là việc xây dựng thành công bộ công cụ phần mềm viết bằng ngôn ngữ Fortran, kết hợp chuỗi dữ liệu sóng NOAA 20 năm (1997-2016), ảnh vệ tinh Landsat 26 năm (1988-2014) và hệ thống video camera giám sát thực địa (2013-2014) tại bờ tây vịnh Nha Trang (phạm vi từ cửa sông Cái đến quân cảng/Vinpearl, chiều dài ~7 km).

Literature Review và Positioning

Nghiên cứu hình thái động lực ven bờ trên thế giới đã phát triển qua nhiều trường phái học thuật với những bước chuyển biến quan trọng:

Trường phái mô hình quá trình dọc bờ bắt đầu từ lý thuyết 1-đường kinh điển của Pelnard-Considere (1956), mô tả biến động đường bờ dựa trên phương trình bảo toàn khối lượng trầm tích dọc bờ $\frac{\partial y}{\partial t} + \frac{1}{D}\frac{\partial Q_{ls}}{\partial x} = 0$. Hanson và Kraus (1989) phát triển thành mô hình số GENESIS, áp dụng hiệu quả cho các bờ biển mở chịu chi phối bởi gradient dòng vận chuyển bùn cát dọc bờ $Q_{ls}$ trong quy mô thập kỷ. Lê Xuân Hoàn (2011) đã nâng cấp mô hình này bằng cách tích hợp mô hình phát triển lạch triều của Kraus (2000, 2002), bổ sung công thức tính dòng bùn cát trong vùng sóng tràn (swash zone) dựa trên nguyên lý Bagnold (1963, 1966) và Watanabe (1991), kiểm nghiệm thành công tại Long Island (Mỹ) và Badreveln Spit (Thụy Điển). Tuy nhiên, các mô hình này hoàn toàn loại bỏ quá trình vận chuyển ngang bờ, dẫn đến sai số lớn khi áp dụng cho các bờ biển vịnh hình cung có sự trao đổi bùn cát ngang bãi chiếm ưu thế.

Trường phái mô hình quá trình ngang bờ và tiêu tán năng lượng khởi nguồn từ khái niệm trắc ngang cân bằng $h = A x^{2/3}$ của Bruun (1954) và Dean (1977), dựa trên giả thiết năng lượng sóng tiêu tán đồng đều trên một đơn vị thể tích trong vùng sóng vỡ. Kriebel và Dean (1985) đề xuất mô hình EDUNE dựa trên độ lệch tiêu tán năng lượng $D_{cs} = K(D - D^*)$, sau đó được Zheng và Dean (1997) hoàn thiện. Larson và Kraus (1989) phát triển mô hình SBEACH chia trắc ngang thành 4 vùng động lực để tính toán xói lở do bão. Dù vậy, các mô hình này gặp hạn chế lớn vì chỉ mô phỏng tốt pha xói ngắn hạn trong bão mà không thể mô phỏng chính xác sự bồi tụ phục hồi bãi sau bão và sự dịch chuyển của bar cát ven bờ.

Trường phái mô hình lai ghép (Hybrid) và thống kê động lực học hiện đại xuất hiện nhằm khắc phục các nhược điểm trên. Miller và Dean (2004a) xây dựng mô hình biến động đường bờ bán thực nghiệm phụ thuộc vào độ cao sóng và mực nước tổng cộng, đạt sai số chuẩn hóa $NMSE = 0.30$ tại Wildwood (New Jersey) và $NMSE = 0.33$ tại Long Beach (Washington). Yates và cộng sự (2009, 2011) đề xuất mô hình cân bằng năng lượng sóng tức thời, hiệu chỉnh thành công trên 5 bãi biển California với hệ số tương quan $R^2 > 0.5$. Davidson, Lewis và Turner (2010, 2013) phát triển mô hình ShoreFor tại Gold Coast và Collaroy-Narrabeen (Úc), chứng minh rằng ngưỡng độ dốc sóng phân định bồi - xói thay đổi động theo tiền sử trường sóng và phân loại hình thái bãi. Karunarathna và Reeve (2009, 2011, 2013) thiết lập mô hình khuếch tán 1D/2D dạng lai ghép $\frac{\partial h}{\partial t} = \frac{\partial}{\partial x}(K\frac{\partial h}{\partial x}) + S(x,t)$ để mô phỏng hình thái cửa sông Humber và trắc ngang bãi biển.

Trong y văn tồn tại hai cuộc tranh luận học thuật sâu sắc:

  1. Bản chất của tiêu tán năng lượng: Dean (1977) giả thiết tiêu tán năng lượng sóng trên một đơn vị thể tích là động lực chính duy trì trắc ngang cân bằng. Ngược lại, Lee (1994) chỉ ra rằng năng lượng là đại lượng vô hướng nên giả thiết của Dean không phản ánh được tính định hướng của chuyển động hạt trầm tích dưới tác động của trường vận tốc sóng phi đối xứng.
  2. Quy mô chi phối biến động dài hạn: Các tác giả theo trường phái One-Line cổ điển (Hanson & Kraus, 1989) cho rằng biến động dài hạn hàng thập kỷ hoàn toàn do gradient vận chuyển dọc bờ quyết định. Tuy nhiên, Ruggiero và cộng sự (2010), Yates và cộng sự (2009) và Splinter và cộng sự (2013) đã chứng minh tại các vùng bờ biển có dao động mùa rõ rệt, sự mất cân bằng trong vận chuyển ngang bờ dưới tác động của chuỗi sóng bão liên tiếp có thể tích lũy thành biến động dài hạn không thể đảo ngược nếu thiếu nguồn bùn cát bổ cập.

Positioning của luận án nằm ở điểm giao thoa giữa mô hình sóng khúc xạ - nhiễu xạ vùng ven bờ phức tạp (EBED), mô hình biến đổi đường bờ cân bằng tự hiệu chỉnh tham số và công nghệ quan trắc video camera bãi biển thời gian thực. So với nghiên cứu tại bãi biển Gold Coast (Úc) của Davidson và cộng sự (2010) và bãi biển Torrey Pines (Mỹ) của Yates và cộng sự (2009), nghiên cứu của Vũ Công Hữu đã mở rộng áp dụng cho vùng vịnh nhiệt đới gió mùa có điều kiện địa hình bán kín, chịu ảnh hưởng phức tạp của 19 hòn đảo chắn sóng (tiêu biểu là đảo Hòn Tre) và nguồn bồi tích sông Cái bị kiểm soát bởi công trình kè chỉnh trị ven bờ.

Đó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 trắc ngang cân bằng của Dean (1977) và dạng cải tiến của Kristensen (2013a): $$h(x) = \begin{cases} h_0 & \text{khi } x \le x_0 \ h_0 + A(x - x_0)^m & \text{khi } x_0 < x < x_{cl} \ h_{cl} & \text{khi } x \ge x_{cl} \end{cases}$$ trong đó $h_0$ là độ cao gờ bãi (berm) tại vị trí $x_0$, $h_{cl}$ là độ sâu tới hạn vận chuyển trầm tích tại vị trí $x_{cl}$, $A$ là tham số hình dạng phụ thuộc đường kính hạt cát $D_{50}$, và $m \approx 2/3$.

Luận án phát triển khung lý thuyết động lực học đường bờ cân bằng:

  • Proposition 1: Tốc độ dịch chuyển ngang của đường bờ $\frac{dy}{dt}$ tỷ lệ thuận với mức độ mất cân bằng giữa vị trí đường bờ hiện tại $y(t)$ và vị trí đường bờ cân bằng $y_{eq}(t)$, được điều khiển bởi năng lượng sóng tức thời và mực nước dâng do sóng (NDS): $$\frac{dy}{dt} = k [y_{eq}(t) - y(t)]$$ trong đó hệ số tốc độ $k$ nhận giá trị $k_a$ trong pha bồi ($y < y_{eq}$) và $k_e$ trong pha xói ($y > y_{eq}$).
  • Proposition 2: Vị trí đường bờ cân bằng $y_{eq}(t)$ là hàm phi tuyến của độ cao sóng có nghĩa $H_s$, chu kỳ sóng $T_p$, hướng sóng vỡ $\alpha_b$ và mực nước tổng cộng bao gồm thủy triều và nước dâng do sóng $S_w$: $$y_{eq}(t) = \bar{Y} + C_w \frac{H_s^2 L_0}{h_{cl}} - \frac{S_w(t) + \Delta \eta(t)}{\tan \beta}$$ với $\tan \beta$ là độ dốc bãi trước, $L_0$ là bước sóng nước sâu, $C_w$ là hệ số không thứ nguyên.

Minh chứng cho bước chuyển Paradigm: Luận án thúc đẩy sự chuyển dịch từ các mô hình tiền định thuần túy (Deterministic Process-based) nặng nề sang mô hình lai ghép Động lực - Thống kê tự hiệu chỉnh tham số (Data-driven Hydrodynamic Hybrid Paradigm), kết hợp dữ liệu giám sát quang học độ phân giải cao để liên tục cập nhật trạng thái cân bằng động của hệ thống bờ bãi.

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

Khung phân tích của luận án tích hợp đồng bộ 3 hợp phần khoa học:

  1. Hợp phần thủy động lực học sóng ven bờ: Ứng dụng mô hình sai phân hữu hạn giải phương trình cân bằng năng lượng sóng có xét đến hiện tượng khúc xạ và nhiễu xạ (EBED - Energy Balance Equation with Diffraction) viết bằng Fortran để tính toán trường sóng chi tiết từ ngoài khơi truyền vào vùng nước nông vịnh Nha Trang.
  2. Hợp phần xử lý ảnh quang học (Video Camera Photogrammetry): Thiết lập hệ phương trình chuyển đổi hình học từ tọa độ mặt phẳng ảnh $(u, v)$ sang tọa độ thực địa $(X, Y, Z)$ dựa trên các điểm khống chế mặt đất GCP (Ground Control Points) và độ cao đặt camera $Z_0$: $$u = f \frac{r_{11}(X - X_0) + r_{12}(Y - Y_0) + r_{13}(Z - Z_0)}{r_{31}(X - X_0) + r_{32}(Y - Y_0) + r_{33}(Z - Z_0)}$$ $$v = f \frac{r_{21}(X - X_0) + r_{22}(Y - Y_0) + r_{23}(Z - Z_0)}{r_{31}(X - X_0) + r_{32}(Y - Y_0) + r_{33}(Z - Z_0)}$$ cho phép tự động giải đoán đường mép nước (đường bờ thực tế) với tần suất liên tục hàng giờ.
  3. Hợp phần mô hình số dự báo biến động đường bờ: Module thuật toán tự hiệu chỉnh tham số (Auto-calibration module) sử dụng phương pháp bình phương tối thiểu (Least Squares Optimization) để xác định bộ tham số tối ưu ($k_a, k_e, \Delta Y_0$) từ chuỗi số liệu camera, từ đó vận hành phương trình vi phân mô phỏng diễn biến đường bờ theo thời gian.

Điều kiện biên (Boundary Conditions) được xác định rõ ràng: Biên ngoài khơi lấy từ số liệu tái phân tích sóng NOAA (1997-2016); biên phía bờ khống chế bởi cao trình kè cứng đường Trần Phú; độ sâu tới hạn khống chế vận chuyển bùn cát $h_{cl}$ được tính toán theo công thức Hallermeier (1981) liên hệ với chiều cao sóng bão hàng năm.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường triết học Thực chứng (Positivism) kết hợp Hiện thực phê phán (Critical Realism) trong hải dương học vật lý và địa mạo động lực bờ biển. Thiết kế nghiên cứu đa tầng bậc (Multi-level Research Design) bao gồm:

  • Tầng vi mô (Micro-level / Event scale): Mô phỏng chi tiết các quá trình thủy động lực sóng, dòng chảy ven bờ, nước dâng do sóng trong các sự kiện bão điển hình (như bão Nari, bão Haiyan năm 2013) và các đợt không khí lạnh gió mùa Đông Bắc mạnh.
  • Tầng trung mô (Meso-level / Seasonal scale): Khôi phục và đánh giá biến động trắc ngang, độ rộng bãi biển và vị trí đường bờ giữa hai mùa gió mùa Đông Bắc (NE) và Tây Nam (SW) thông qua dữ liệu camera giám sát liên tục.
  • Tầng vĩ mô (Macro-level / Engineering & Long-term scale): Phân tích xu thế biến động đường bờ trong 26 năm (1988-2014) từ ảnh vệ tinh Landsat và chuỗi số liệu sóng NOAA 20 năm (1997-2016), tách lọc các dao động điều hòa dài hạn.

Quy trình nghiên cứu rigorous

Chiến lược thu thập và xử lý số liệu được chuẩn hóa nghiêm ngặt:

  • Hệ thống camera giám sát: Lắp đặt trên bờ tây vịnh Nha Trang với phạm vi quan trắc bao quát bãi biển trung tâm. Tọa độ thực địa được hiệu chỉnh thông qua mạng lưới điểm khống chế GCP đo bằng máy định vị GPS vi sai độ chính xác cao. Đã tiến hành khảo sát độ nhạy của sai số vị trí theo tiêu cự camera $f$ để triệt tiêu biến dạng quang học.
  • Dữ liệu viễn thám và trắc đạc: Sử dụng ảnh vệ tinh Landsat đa thời gian (tháng 6/1988, 4/1995, 2/1996, 6/1996, 8/1996, 3/1997, 1999, 2014) kết hợp số liệu khảo sát mặt cắt địa hình đáy biển do Viện Hải dương học và Đề tài Nghị định thư Việt - Pháp (MOST2:NDT.FRA/16) thực hiện trong các năm 2013-2014.
  • Tam giác hóa dữ liệu (Data Triangulation): Kiểm chứng chéo độc lập giữa 3 nguồn dữ liệu: (1) Số liệu giải đoán đường bờ từ video camera; (2) Số liệu đường bờ trích xuất từ ảnh vệ tinh Landsat; (3) Kết quả tính toán từ mô hình số trị EBED và mô hình biến đổi đường bờ.
  • Độ tin cậy và chuẩn hóa mô hình: Đánh giá thông qua chỉ số sai số bình phương quân phương chuẩn (Normalized Mean Square Error - NMSE) và sai số bình phương quân phương (Root Mean Square Error - RMSE): $$NMSE = \frac{\sum (y_{obs} - y_{sim})^2}{\sum (y_{obs} - \bar{y}{obs})^2}$$ $$RMSE = \sqrt{\frac{1}{N}\sum{i=1}^N (y_{obs,i} - y_{sim,i})^2}$$

Data và phân tích

Dữ liệu khí tượng hải văn đầu vào:

  • Chế độ gió (Trạm Khí tượng Nha Trang 1995-2014): Gió mùa Đông Bắc (tháng 10 đến tháng 3 năm sau) có tốc độ trung bình 2.7 - 4.1 m/s, hướng thịnh hành N-NE (tần suất 19.4 - 35.8%), tốc độ cực đại khi có bão/không khí lạnh đạt 15 - 24 m/s (thậm chí 30 m/s như cơn bão ngày 10/11/1998). Gió mùa Tây Nam (tháng 5 đến tháng 8) có tốc độ trung bình 1.5 - 1.8 m/s, hướng thịnh hành SW (tần suất 19.8 - 24.4%), tần suất lặng gió cao (40.8 - 44.8%).
  • Chế độ sóng (NOAA 1997-2016): Mùa Đông Bắc sóng ngoài khơi có hướng chủ đạo NE với độ cao $H_s$ từ 2.0 - 4.0 m (tháng 11, 12, 1); mùa Tây Nam sóng có độ cao nhỏ hơn đáng kể ($H_s < 1.2\text{ m}$). Các cơn bão mạnh lịch sử như Kyle (11/1993) và Lola (12/1993) tạo sóng ngoài khơi Khánh Hòa cực đại $> 13\text{ m}$, chu kỳ $T_p \approx 15\text{ s}$.
  • Phân tích điều hòa chuỗi thời gian (Harmonic Analysis): Áp dụng biến đổi Fourier và hàm điều hòa phân rã chuỗi vị trí đường bờ $y(t)$ thành tổng các thành phần tần số: $$y(t) = \bar{y} + \sum_{k=1}^M A_k \cos(\omega_k t - \varphi_k)$$ để xác định chính xác biên độ $A_k$ và pha $\varphi_k$ tương ứng với chu kỳ nửa năm ($T = 0.5\text{ năm}$), chu kỳ năm ($T = 1.0\text{ năm}$) và chu kỳ liên năm ($T = 3.0\text{ năm}$).

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

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

  1. Tính bất đối xứng hình thái động lực theo mùa sâu sắc: Kết quả mô phỏng và số liệu camera chỉ rõ bãi biển Nha Trang trải qua hai pha đối lập hoàn toàn trong năm. Trong mùa gió Đông Bắc (tháng 10 - tháng 3), trường sóng mạnh kết hợp nước dâng do sóng làm bãi biển bị xói lở dữ dội, đường bờ lùi sâu vào đất liền từ 10 - 25 m tùy vị trí, trắc ngang bãi bị hạ thấp và độ dốc đáy gia tăng. Ngược lại, trong mùa gió Tây Nam (tháng 5 - tháng 8), sóng êm hơn tạo điều kiện cho bùn cát di chuyển ngược từ ngoài khơi vào bờ, bãi biển được tái tạo bồi tụ trở lại, đường bờ tiến ra biển.
  2. Độ chính xác cao của bộ công cụ mô phỏng: Tại các vị trí kiểm định dọc bờ tây vịnh Nha Trang ($x = 50\text{ m}, 100\text{ m}, 150\text{ m}, 200\text{ m}, 250\text{ m}$), mô hình biến đổi đường bờ tích hợp thuật toán tự hiệu chỉnh cho kết quả khớp hoàn hảo với vị trí đường bờ giải đoán từ camera. Chỉ số sai số chuẩn hóa NMSE đạt mức rất thấp: tại vị trí $x = 200\text{ m}$, $NMSE \approx 0.18 - 0.25$, vượt trội so với các mô hình kinh nghiệm quốc tế như mô hình Miller & Dean (2004a) tại bãi biển Duck ($NMSE = 0.95$) hay Gold Coast ($NMSE = 0.52$).
  3. Tác động xói lở đột biến trong các sự kiện bão cực đoan: Trong thời đoạn chịu ảnh hưởng liên tiếp của hai cơn bão Nari (10/2013) và siêu bão Haiyan (11/2013), năng lượng sóng tăng vọt làm đường bờ bị xói lùi tức thời từ 5 - 12 m chỉ trong vòng 48 đến 72 giờ. Mô hình đã tái hiện chính xác pha sụt giảm bãi và quá trình phục hồi chậm sau bão kéo dài nhiều tuần.
  4. Phát hiện chu kỳ biến động 3 năm chi phối hình thái dài hạn: Phân tích hàm điều hòa chuỗi số liệu khôi phục chỉ ra rằng bên cạnh chu kỳ nửa năm và 1 năm gắn liền với gió mùa, biến động đường bờ bãi biển Nha Trang chứa đựng một thành phần điều hòa rõ rệt với chu kỳ xấp xỉ 3 năm. Đây là phát hiện mới, phản ánh sự liên kết giữa hình thái bãi vịnh Nha Trang với các hiện tượng hoàn lưu khí quyển - đại dương quy mô lớn khu vực Tây Thái Bình Dương (ENSO) và chu kỳ biến thiên lượng bồi tích từ lưu vực sông Cái.
  5. Hệ quả tiêu cực từ các công trình nhân sinh: Việc xây dựng kè bảo vệ bờ đường Trần Phú và cầu Trần Phú tại cửa sông Cái đã làm thay đổi trường thủy động lực cục bộ, chặn dòng cát bùn bồi bổ tự nhiên dọc bờ, dẫn đến hiện tượng bãi biển phía trước UBND tỉnh bị xói mất hoàn toàn bãi cát; trắc ngang ngã ba Nguyễn Chánh - Trần Phú bị xói lùi 9 m (2006-2013), đáy gần bờ dốc đứng (~4%), tiềm ẩn hiểm họa đối với hoạt động du lịch tắm biển.

Implications đa chiều

  • Về mặt lý thuyết: Luận án đã khẳng định tính đúng đắn của lý thuyết đường bờ cân bằng trong điều kiện bờ biển vịnh nhiệt đới gió mùa; chứng minh rằng sự kết hợp giữa mô hình khuếch tán hình thái và dữ liệu viễn thám thực địa giải quyết triệt để bài toán thiếu hụt số liệu đo đạc địa hình đáy biển truyền thống.
  • Về mặt phương pháp luận: Thiết lập quy trình công nghệ hoàn chỉnh từ thu nhận ảnh camera quang học $\rightarrow$ tự động giải đoán đường bờ $\rightarrow$ khôi phục trường sóng ven bờ bằng EBED $\rightarrow$ tự hiệu chỉnh tham số $\rightarrow$ dự báo diễn biến đường bờ. Bộ công cụ viết bằng Fortran có mã nguồn mở, tốc độ tính toán nhanh, chi phí vận hành thấp, dễ dàng chuyển giao áp dụng cho các dải bờ biển khác tại Việt Nam và các quốc gia đang phát triển.
  • Về mặt ứng dụng thực tiễn: Cung cấp cơ sở khoa học định lượng cho các cơ quan quản lý đô thị và du lịch tỉnh Khánh Hòa trong việc thiết kế các giải pháp công trình mềm: nuôi bãi nhân tạo (beach nourishment), xây dựng đê ngầm giảm sóng tạo bãi (submerged breakwater), và quy hoạch hành lang bảo vệ bờ biển an toàn.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận các giới hạn khoa học:

  1. Độ dài chuỗi số liệu camera thực địa: Hệ thống video camera quan trắc liên tục mới chỉ được thu thập và khai thác chi tiết trong giai đoạn 2013-2014 do hạn chế về thời gian thực hiện các đề tài nghị định thư; cần chuỗi dữ liệu camera liên tục 5-10 năm để nâng cao độ tin cậy của thuật toán tự hiệu chỉnh tham số dài hạn.
  2. Giả định về tính đồng nhất của trắc ngang: Mô hình biến đổi đường bờ 1-đường mở rộng vẫn dựa trên giả thiết hình dạng trắc ngang duy trì trạng thái cân bằng dạng Dean ($h = Ax^{2/3}$), chưa mô tả chi tiết sự phân bố đa chiều của các thanh bar cát ngầm di chuyển phức tạp trong đới sóng vỡ.
  3. Hiệu ứng che chắn 3D của hệ thống đảo: Mô hình sóng EBED dù đã xét hiện tượng nhiễu xạ nhưng chưa mô phỏng hết tương tác 3D phi tuyến phức tạp giữa sóng, triều và dòng chảy hoàn lưu vịnh qua các luồng lạch hẹp giữa 19 đảo (đặc biệt là hòn Tre, hòn Mun, hòn Tằm).

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

  • Phát triển mạng lưới camera quan trắc thông minh ứng dụng trí tuệ nhân tạo (Deep Learning) để tự động nhận diện mép nước thời gian thực dưới mọi điều kiện chiếu sáng và thời tiết.
  • Tích hợp mô hình biến động đường bờ với các kịch bản biến đổi khí hậu và nước biển dâng của Bộ Tài nguyên và Môi trường (kịch bản RCP 4.5 và RCP 8.5) đến năm 2050 và 2100.
  • Kết nối mô hình đường bờ với mô hình dòng bùn cát lưu vực sông Cái để lượng hóa tác động của các hồ chứa thượng nguồn đến cán cân trầm tích ven vịnh.
  • Mở rộng mô hình sang dạng tích hợp 2D/3D phi cấu trúc kết hợp học máy (Physics-Informed Machine Learning) nhằm tối ưu hóa các phương án nuôi bãi nhân tạo tại vịnh Nha Trang.

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

  • Tác động học thuật (Academic Impact): Công trình là luận án tiến sĩ hải dương học tiêu biểu trong nước kết hợp thành công công nghệ xử lý ảnh camera viễn thám độ phân giải cao với mô hình toán thủy thạch động lực; mở ra hướng nghiên cứu mới về mô hình lai ghép hình thái ven bờ tại Việt Nam; có tiềm năng trích dẫn cao trong các tạp chí chuyên ngành như Coastal Engineering, Ocean Engineering, và Tạp chí Khoa học ĐHQGHN.
  • Chuyển đổi ngành & Ứng dụng kỹ thuật (Industry & Engineering): Cung cấp bộ công cụ tính toán thủy lực bờ biển hỗ trợ các viện nghiên cứu, doanh nghiệp tư vấn xây dựng công trình thủy (như Vinpearl, Ban Quản lý Dự án tỉnh Khánh Hòa) trong việc tính toán ổn định công trình kè, đê biển, và bến cảng.
  • Tác động chính sách (Policy Influence): Là luận cứ khoa học trực tiếp cho UBND tỉnh Khánh Hòa, Sở Nông nghiệp & Phát triển Nông thôn, Sở Tài nguyên & Môi trường trong việc ban hành Quy chế quản lý không gian bờ biển Nha Trang, thiết lập ranh giới hành lang bảo vệ bờ biển theo Luật Tài nguyên, môi trường biển và hải đảo.
  • Lợi ích xã hội và môi trường (Societal & Environmental Benefits): Giúp bảo tồn tài nguyên bãi cát tự nhiên phục vụ cộng đồng và du lịch; cảnh báo sớm các vị trí bờ biển có nguy cơ xói lở sâu gây mất an toàn cho người dân và du khách tắm biển; giảm thiểu thiệt hại kinh tế do sạt lở hạ tầng đường Trần Phú.
  • Ý nghĩa quốc tế (International Relevance): Đóng góp dữ liệu thực nghiệm quý giá về động lực học hình thái bãi biển dạng vịnh nhiệt đới chịu ảnh hưởng sâu sắc của gió mùa Đông Nam Á vào mạng lưới quan trắc bờ biển toàn cầu (CoastSnap/Argus network).

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

  • Nghiên cứu sinh & Nhà nghiên cứu trẻ: Tiếp cận mã nguồn mở Fortran, phương pháp luận tích hợp viễn thám - mô hình số trị, và các khoảng trống nghiên cứu về tương tác thủy thạch động lực biển nhiệt đới.
  • Các nhà khoa học & Giáo sư chuyên ngành: Khai thác khung lý thuyết mở rộng về đường bờ cân bằng và phương pháp phân tích điều hòa chuỗi hình thái bãi biển đa quy mô thời gian.
  • Kỹ sư R&D & Tư vấn công trình thủy: Ứng dụng bộ công cụ để tính toán trường sóng, dự báo xói lở phục vụ thiết kế các công trình bảo vệ bờ biển, công trình nuôi bãi nhân tạo với chi phí đầu tư tối ưu.
  • Cơ quan quản lý nhà nước & Hoạch định chính sách: Sử dụng các bản đồ dự báo biến động đường bờ và phân tích nguyên nhân xói lở để đưa ra các quyết định quy hoạch phát triển kinh tế biển bền vững.

Câu hỏi chuyên sâu

  1. Đóng góp lý thuyết độc đáo nhất của luận án là gì? Luận án đã mở rộng thành công lý thuyết động lực học đường bờ cân bằng của Miller & Dean (2004a) và Yates và cộng sự (2009) sang môi trường vịnh bán kín nhiệt đới gió mùa. Nghiên cứu đã chứng minh rằng trạng thái cân bằng động của đường bờ không chỉ phụ thuộc đơn thuần vào độ dốc sóng ngoài khơi mà chịu sự chi phối quyết định bởi tương tác phi tuyến giữa nước dâng do sóng (NDS) và sự chuyển đổi hướng sóng thịnh hành giữa hai mùa gió đối nghịch NE-SW.
  2. Phương pháp nghiên cứu có điểm đổi mới nào khi so sánh với các nghiên cứu quốc tế trước đây? So với mô hình GENESIS (Hanson & Kraus, 1989) chỉ xét vận chuyển dọc bờ và mô hình SBEACH (Larson & Kraus, 1989) chỉ xét vận chuyển ngang bờ ngắn hạn, bộ công cụ của luận án tích hợp trực tiếp mô hình sóng khúc xạ - nhiễu xạ EBED với module tự hiệu chỉnh tham số dựa trên dữ liệu giải đoán từ camera quang học mặt đất. Điểm đột phá so với nghiên cứu của Davidson và cộng sự (2010) tại Úc là việc giải quyết bài toán biến dạng quang học camera bằng ma trận chuyển đổi tọa độ GCP thực địa kết hợp kiểm định tiêu cự $f$, cho phép thu nhận vị trí mép nước độ chính xác dưới mét tại vùng biển có biên độ triều nhỏ.
  3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì? Phát hiện bất ngờ nhất là sự tồn tại của chu kỳ biến động hình thái 3 năm rõ rệt thông qua phân tích hàm điều hòa chuỗi số liệu vị trí đường bờ. Điều này chứng minh rằng bên cạnh tính tuần hoàn mùa vụ (NE-SW) làm bãi xói lở và tái bồi hàng năm, bãi biển Nha Trang còn chịu sự điều tiết bởi một chu kỳ tích lũy/suy giảm trầm tích 3 năm gắn liền với biến động khí hậu quy mô khu vực và chế độ xả lũ - bồi tích của sông Cái.
  4. Luận án có cung cấp quy trình tái lập (Replication Protocol) không? Luận án cung cấp đầy đủ và chi tiết quy trình tính toán: từ hệ phương trình toán học, sơ đồ khối thuật toán tự hiệu chỉnh tham số, bảng danh mục tọa độ điểm khống chế GCP, đến cấu trúc chương trình tính toán viết bằng ngôn ngữ Fortran trong phần phụ lục, đảm bảo khả năng tái lập hoàn toàn trên các tập dữ liệu bãi biển khác.
  5. Chương trình nghiên cứu 10 năm tới được định hình ra sao? Chương trình nghiên cứu 10 năm tới tập trung vào: (1) Tự động hóa hệ thống quan trắc bờ biển bằng AI Camera thời gian thực; (2) Phát triển mô hình ghép thủy thạch động lực - hình thái 3D phi cấu trúc tích hợp biến đổi khí hậu; (3) Hoàn thiện cẩm nang kỹ thuật hướng dẫn công nghệ nuôi bãi nhân tạo thích ứng cho toàn bộ dải ven biển Nam Trung Bộ.

Kết luận

  1. Xây dựng thành công bộ công cụ tích hợp hoàn chỉnh: Luận án đã lựa chọn, phát triển và tích hợp thành công bộ công cụ mô hình số trị bao gồm mô hình sóng EBED, mô đun xử lý ảnh camera giám sát bãi biển và mô hình biến đổi đường bờ tự hiệu chỉnh tham số viết bằng ngôn ngữ Fortran.
  2. Làm rõ quy luật động lực học hình thái theo mùa: Xác lập định lượng quy luật biến động đối lập giữa mùa gió Đông Bắc (pha xói lở sâu, đường bờ lùi từ 10 - 25 m, đáy dốc) và mùa gió Tây Nam (pha bồi tụ tái tạo, bãi thoai thoải mở rộng ra biển).
  3. Khám phá các quy mô thời gian chi phối dài hạn: Lần đầu tiên tại bãi biển Nha Trang, nghiên cứu đã phân tách và chỉ rõ 3 quy mô thời gian chi phối: quy mô sự kiện tức thời (bão, không khí lạnh), quy mô mùa/năm (chu kỳ gió mùa NE-SW), và quy mô dài hạn 3 năm (dao động điều hòa liên năm).
  4. Chỉ ra nguyên nhân suy thoái bãi do nhân sinh: Cung cấp bằng chứng khoa học khẳng định công trình kè bảo vệ bờ đường Trần Phú và cầu Trần Phú làm biến dạng chế độ bùn cát cửa sông Cái, gây mất hoàn toàn bãi trước UBND tỉnh và làm xói lùi 9 m tại ngã ba Nguyễn Chánh (2006-2013).
  5. Đạt độ chính xác cao và khả thi trong dự báo: Mô hình đạt sai số chuẩn hóa $NMSE \approx 0.18 - 0.25$ tại các trắc ngang bãi trung tâm, mở ra giải pháp công nghệ chi phí thấp, hiệu quả cao phục vụ quy hoạch, bảo vệ và tôn tạo bền vững bãi biển Nha Trang cũng như các vùng bờ biển tương tự trên cả nước.