Tổng quan về luận án

Nghiên cứu động lực học sóng biển tại khu vực Biển Đông đóng vai trò đặc biệt quan trọng đối với an toàn hàng hải, khai thác dầu khí thềm lục địa, công trình ven bờ và an ninh quốc phòng biển đảo. Tuy nhiên, các mô hình số trị tính toán và dự báo sóng phổ biến hiện nay như Simulating Waves Nearshore (SWAN) hay Wave Action Model (WAM) vẫn còn tồn tại các sai số mang tính hệ thống do sự bất định của trường gió khí quyển đầu vào, sự thô sơ trong việc thiết lập tham số vật lý tiêu tán năng lượng và sự thiếu hụt cơ chế cập nhật dữ liệu quan trắc thực tế theo thời gian thực. Nhằm giải quyết triệt để vấn đề này, luận án tiến sĩ chuyên ngành Hải dương học (mã số: 9440228.01) của tác giả Nguyễn Trung Thành với đề tài "Nghiên cứu ứng dụng phương pháp đồng hoá số liệu phục vụ mô hình dự báo sóng", thực hiện dưới sự hướng dẫn khoa học của PGS.TS. Nguyễn Minh Huấn và PGS.TS. Trần Quang Tiến tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội (2020), đã tiên phong tích hợp khung hiệu chỉnh tự động SWAN-CI và thuật toán đồng hóa số liệu lọc Kalman tổ hợp (Ensemble Kalman Filter - EnKF) trên nền tảng phần mềm mã nguồn mở OpenDA cho mô hình SWAN tại Biển Đông.

Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án giải quyết là sự thiếu vắng các quy trình tối ưu hóa tham số vật lý khách quan đa chiều và cơ chế đồng hóa số liệu động lực học sóng đa nguồn (tích hợp đồng thời trạm phao, ra đa biển và vệ tinh đo cao) tại vùng biển Việt Nam. Trước đây, việc hiệu chỉnh tham số mô hình SWAN chủ yếu dựa vào kinh nghiệm thủ công cục bộ, còn các nghiên cứu đồng hóa số liệu tại Biển Đông (như Cao Lei et al., 2015 [42]; Trần Quang Tiến, 2010 [10]) chỉ dừng lại ở các thử nghiệm đơn lẻ, chưa giải mã định lượng độ nhạy của từng tham số sóng cũng như chưa tối ưu hóa ma trận hiệp biến sai số trong EnKF.

Luận án tập trung giải quyết 3 câu hỏi và giả thuyết nghiên cứu chính:

  1. Câu hỏi nghiên cứu 1 (RQ1): Trong các tham số vật lý điều khiển phổ sóng của mô hình SWAN, những tham số nào chi phối mạnh nhất đến sai số độ cao sóng có nghĩa (SWH) tại Biển Đông và làm thế nào để tối ưu hóa tự động không qua đạo hàm? $\rightarrow$ Giả thuyết 1 (H1): Quá trình tiêu tán năng lượng do sóng bạc đầu (whitecapping dissipation) do các tham số $CDS2$ và $powst$ quy định là nguồn gây sai số chủ đạo, trong khi tham số số sóng $powk$ và tần số cắt $cutfr$ có độ nhạy không đáng kể.
  2. Câu hỏi nghiên cứu 2 (RQ2): Cơ chế lọc Kalman tổ hợp (EnKF) kết hợp với mô hình SWAN có khả năng hiệu chỉnh trường sóng nền và điều kiện ban đầu để cải thiện chất lượng tái phân tích và dự báo sóng hay không? $\rightarrow$ Giả thuyết 2 (H2): Cập nhật trạng thái động lực học thông qua ma trận nhân Kalman ngẫu nhiên Monte Carlo sẽ giảm thiểu sai số trung bình (Mean Error - ME) và sai số căn bậc hai trung bình bình phương (Root Mean Square Error - RMSE) của trường sóng tái phân tích.
  3. Câu hỏi nghiên cứu 3 (RQ3): Hiệu ứng đồng hóa số liệu duy trì được giá trị dự báo vượt trội trong khoảng thời gian bao lâu so với dự báo thông thường không đồng hóa (NoEnKF)? $\rightarrow$ Giả thuyết 3 (H3): Tác động của đồng hóa số liệu EnKF cải thiện rõ rệt độ chính xác dự báo trong hạn 1–24 giờ và tiếp tục duy trì ảnh hưởng tích cực đến hạn dự báo 25–48 giờ.

Khung lý thuyết của nghiên cứu được neo vững trên phương trình cân bằng tác động phổ hai chiều (Wave Action Balance Equation) của Hasselmann (1973) [36] và Booij et al. (1999) [17], lý thuyết truyền năng lượng gió của Phillips (1957) [52] - Miles (1957) [48], lý thuyết lọc ngẫu nhiên tối ưu của Kalman (1960) [40] và mở rộng tổ hợp EnKF của Evensen (1994) [28]. Luận án có phạm vi bao phủ toàn bộ vùng Biển Đông và chi tiết hóa vùng biển ven bờ Việt Nam, kiểm định trên chuỗi dữ liệu thực nghiệm phong phú gồm 4 trạm phao biển (Phao 4001, 4002, 4003, 4004), trạm giàn khoan cố định ngoài khơi MSP1 (Mỏ Bạch Hổ), trạm Ra đa bờ biển Đồng Hới (Quảng Bình) và các vệt quét độ cao sóng vệ tinh Altimeter (Jason-2, Saral-AltiKa).

Literature Review và Positioning

Lịch sử tính toán và mô hình hóa sóng đại dương đã trải qua ba thế hệ tiến hóa mang tính nền tảng. Khởi đầu từ phương pháp bán thực nghiệm Synop của Sverdrup và Munk (1947) [58] kết hợp quan hệ phổ sóng của Pierson và Moskowitz (1955) [53], đến phương trình cân bằng năng lượng phổ thế hệ I của Gelci et al. (1957) [33]. Mô hình thế hệ II xuất hiện khi Hasselmann (1962) [35] chứng minh vai trò quyết định của tương tác phi tuyến bậc bốn (quadruplet wave-wave interaction), nhưng bị giới hạn bởi các giả định tham số hóa dạng phổ JONSWAP (Hasselmann et al., 1973 [36]). Bước nhảy vọt sang mô hình phổ thế hệ III (WAM - WAMDI Group, 1988 [37]; SWAN - Booij et al., 1999 [17]; WAVEWATCH III - Tolman, 2009 [63]) cho phép giải số trị đầy đủ phương trình cân bằng mật độ tác động phổ $N(\sigma, \theta)$ mà không cần áp đặt phổ nghiệm tiền định, tính toán đồng thời tương tác sóng bậc bốn qua phép xấp xỉ tương tác rời rạc (Discrete Interaction Approximation - DIA của Hasselmann, 1985 [38]) và tương tác ba sóng nước nông LTA (Eldeberky và Battjes, 1995 [25]).

Trong lĩnh vực đồng hóa số liệu đại dương, cuộc tranh luận học thuật diễn ra gay gắt giữa hai trường phái tiếp cận:

  • Trường phái biến phân (Variational Data Assimilation - 3DVar/4DVar): Đại diện bởi Le Dimet và Talagrand (1986), ECMWF (2018) [71] và Orzech et al. (2012) [46]. 4DVar cực tiểu hóa hàm mục tiêu chi phí trên toàn bộ cửa sổ thời gian, tối ưu hóa rất tốt trạng thái ban đầu nhưng đòi hỏi xây dựng toán tử mô hình tiếp tuyến liên hợp (adjoint model) cực kỳ phức tạp và tiêu tốn tài nguyên tính toán khổng lồ khi xử lý các phương trình sóng phi tuyến cao.
  • Trường phái lọc thống kê tuần tự tổ hợp (Sequential Ensemble Kalman Filter - EnKF): Khởi xướng bởi Evensen (1994) [28], Burgers et al. (1998) và Houtekamer & Mitchell (1998). EnKF sử dụng phương pháp thống kê Monte Carlo để ước tính ma trận hiệp biến sai số động (background error covariance) trực tiếp từ độ tán của các thành phần tổ hợp mà không cần tuyến tính hóa mô hình hay lập trình mô hình liên hợp.
Lịch sử mô hình hóa sóng & Đồng hóa số liệu:
[Sverdrup & Munk, 1947] -> [Hasselmann JONSWAP, 1973] -> [SWAN Gen III, 1999]
                                                                  |
[Kalman Filter, 1960]  -> [Evensen EnKF, 1994]       -> [OpenDA + SWAN-CI]

Về mặt thực nghiệm quốc tế, Chen (2004) [22] khi đồng hóa độ cao sóng vệ tinh ERS-2 vào mô hình toàn cầu NWW3 đã ghi nhận sai số giảm mạnh trong 12 giờ đầu, trung bình trong 12–24 giờ và suy giảm dần sau 24 giờ. Tại vịnh Persian, Nima Serpoushan et al. (2013) [57] đã sử dụng OpenDA-SWAN với EnKF tổ hợp 10–50 thành phần, tìm ra bộ tham số sóng bạc đầu $C_{dis}=3{,}95$ và $\delta=0{,}6$. Tương tự, Sofia Caires et al. (2018) [19] tại vùng biển đông Hàn Quốc đã cấu hình EnKF với độ tán gió $1\text{ m/s}$, bán kính địa phương hóa $500\text{ km}$ và sai số quan trắc $0{,}2\text{ m}$ cho thấy sự ổn định vượt trội.

Tại Việt Nam, các nghiên cứu của Huỳnh Thị Hồng Ngự et al. (2008) [7], Kiều Quốc Chánh (2011) [1], Hoàng Đức Cường et al. (2012) [2], Dư Đức Tiến et al. (2017) [61] đã ứng dụng 3DVar và LETKF nhưng chủ yếu cho mô hình khí quyển WRF. Luận án của Nguyễn Trung Thành (2020) đã định vị chính xác vị trí khoa học tiên phong: Đây là công trình đầu tiên tại Việt Nam thiết lập hoàn chỉnh chu trình ghép nối hai chiều giữa mô hình sóng phổ SWAN và hệ thống đồng hóa EnKF trên nền tảng OpenDA, đồng thời kết hợp công cụ hiệu chỉnh không đạo hàm DUD để giải quyết toàn diện bài toán tối ưu tham số và tái phân tích trường sóng Biển Đông.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mang lại những đóng góp lý thuyết quan trọng làm sâu sắc thêm cơ chế động lực học sóng và lý thuyết điều khiển tối ưu trong hải văn:

  1. Kiểm định và mở rộng lý thuyết tiêu tán năng lượng sóng bạc đầu (Whitecapping Theory): Nghiên cứu kiểm chứng mô hình Hasselmann (1973) [36] và Komen et al. (1984) [41] về cơ chế tiêu tán năng lượng $S_{ds,w}(\sigma, \theta) = -CDS2 \cdot \bar{\sigma} \cdot (\bar{k}/k)^m \cdot (\tilde{s}/\tilde{s}_{PM})^p \cdot E(\sigma, \theta)$. Luận án chứng minh bằng thực nghiệm rằng tốc độ tiêu tán sóng bạc đầu phụ thuộc phi tuyến mang tính quyết định vào hai hệ số $CDS2$ (hệ số tỷ lệ tiêu tán) và $powst$ ($p$ - số mũ tỷ lệ độ dốc phổ chuẩn hóa Pierson-Moskowitz), trong khi số mũ số sóng $powk$ ($m$) và tần số cắt $cutfr$ không tạo ra sự biến thiên đáng kể đối với hàm phù hợp GoF.
  2. Khung lý thuyết ước lượng trạng thái tối ưu ngẫu nhiên phi tuyến: Mở rộng lý thuyết ước lượng không chệch tuyến tính tối ưu (Best Linear Unbiased Estimate - BLUE) sang không gian pha của phổ sóng ngẫu nhiên, chứng minh rằng ma trận nhân Kalman $K = P^f H^T (H P^f H^T + R)^{-1}$ được xấp xỉ hóa hiệu quả qua phương pháp Monte Carlo mà không làm suy giảm tính ổn định động lực học của trường sóng 2D.
  3. Mô hình hóa sự suy giảm không gian và thời gian của xung động hóa: Làm sáng tỏ quy luật lan truyền và suy giảm năng lượng hiệu chỉnh: xung đồng hóa đạt hiệu lực hội tụ cực đại trong 24 giờ đầu ($t \le 24\text{h}$), sau đó trường sóng chuyển tiếp dần về trạng thái cân bằng do trường gió khí quyển cưỡng bức khống chế trong giai đoạn $25\text{h} \le t \le 48\text{h}$.

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

Luận án xây dựng một khung phân tích tích hợp đa tầng, liên kết chặt chẽ giữa hiệu chỉnh tham số tự động và đồng hóa dữ liệu tổ hợp:

[Số liệu quan trắc: Phao, Radar MSP1/Đồng Hới, Vệ tinh Jason-2/Saral]
                                |
                                v
               [Khung Tối ưu hóa SWAN-CI (Thuật toán DUD)]
                                |
            (Tối ưu hóa bộ tham số: CDS2, powst, powk, cutfr)
                                |
                                v
                    [Mô hình sóng SWAN v40.85/41]
                                ^
                                | (Cập nhật phân tích: Xa = Xf + K(y - HXf))
                                v
               [Hệ thống Lọc Kalman Tổ hợp (EnKF - OpenDA)]
            (Trường gió nhiễu loạn Monte Carlo + Ma trận Pf)

Khung phân tích này tích hợp đồng thời 4 trụ cột lý thuyết: (1) Động lực học sóng phổ thế hệ III; (2) Phương pháp tối ưu hóa phi tuyến không đạo hàm DUD (Ralston và Jennrich, 1978); (3) Phương pháp lọc ngẫu nhiên tổ hợp EnKF (Evensen, 1994 [28]); và (4) Lý thuyết tương tác khí quyển - đại dương biên. Điều kiện biên của mô hình được thiết lập nghiêm ngặt: hấp thụ hoàn toàn ở biên bờ và biên mở nước sâu, loại trừ hiện tượng phản xạ sóng giả mạo tại các góc phân kỳ bằng cách mở rộng miền tính toán bao trùm toàn bộ Biển Đông từ $99^\circ\text{E}$ đến $122^\circ\text{E}$ và từ $0^\circ\text{N}$ đến $25^\circ\text{N}$.

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

Thiết kế nghiên cứu

Luận án được triển khai theo thế giới quan thực chứng (positivism) với phương pháp tiếp cận diễn dịch - định lượng chặt chẽ (rigorous quantitative modeling). Thiết kế nghiên cứu bao gồm 3 cấp độ:

  • Cấp độ 1: Phân tích độ nhạy và tối ưu hóa tham số (Calibration): Khảo sát hàm mục tiêu GoF (Goodness of Fit) nhằm định lượng mức độ tương thích giữa tính toán và thực đo: $$\text{GoF} = \sqrt{\frac{\sum_{i=1}^{N_i} w_i \left( \frac{H_{s,\text{calc}}^i - H_{s,\text{obs}}^i}{\sigma_i} \right)^2}{\sum_{i=1}^{N_i} w_i}}$$ Trong đó $H_{s,\text{calc}}^i$ và $H_{s,\text{obs}}^i$ lần lượt là độ cao sóng có nghĩa tính toán và thực đo tại thời điểm $i$; $w_i = 1$ là trọng số; $\sigma_i$ là độ lệch chuẩn thực nghiệm. Luận án áp dụng phương pháp DUD (Doesn't Use Derivatives) – một biến thể tựa Gauss-Newton sử dụng phép xấp xỉ afin trên không gian $p+1$ điểm tham số để tìm nghiệm cực tiểu mà không cần tính toán ma trận Jacobi giải tích.
  • Cấp độ 2: Thiết lập sơ đồ đồng hóa số liệu EnKF: Sinh tổ hợp ngẫu nhiên các trường gió đầu vào bằng cách đưa nhiễu trắng Gaussian có cấu trúc tương quan không gian vào thành phần vận tốc gió $U_{10}$ và hướng gió. Trạng thái sóng $X$ bao gồm phân bố phổ năng lượng và độ cao sóng có nghĩa được tích hợp tiến theo thời gian qua mô hình SWAN cho từng thành phần tổ hợp ($N_{ens}$).
  • Cấp độ 3: Đánh giá kiểm định chéo và dự báo thử nghiệm: Tiến hành tính toán so sánh giữa các kịch bản không đồng hóa ($\text{NoEnKF}$), các kịch bản hiệu chỉnh ($\text{HC}$) và các kịch bản đồng hóa tổ hợp ($\text{EnKF01} \div \text{EnKF10}$).

Quy trình nghiên cứu rigorous

Quy trình xử lý dữ liệu và mô phỏng số học tuân thủ nghiêm ngặt các tiêu chuẩn hải dương học quốc tế:

[Thu thập dữ liệu: Gió Reanalysis/WRF, Địa hình ETOPO/Bản đồ Hải quân]
                                |
                                v
     [Tiền xử lý & Kiểm soát chất lượng số liệu (QC) trạm/vệ tinh]
                                |
                                v
          [Lưới tính SWAN: Lưới chữ nhật / Phi cấu trúc đa phân giải]
                                |
                                v
           [Chạy SWAN-CI xác định bộ tham số vật lý tối ưu]
                                |
                                v
   [Khởi tạo N thành phần tổ hợp trường gió & Chạy song song OpenDA-EnKF]
                                |
                                v
    [Trích xuất kết quả phân tích & Kiểm định sai số ME, RMSE, Biểu đồ Scatter]

Các chỉ số thống kê sai số cốt lõi được tính toán chi tiết: $$\text{ME} = \frac{1}{N} \sum_{i=1}^N (S_i - O_i), \quad \text{RMSE} = \sqrt{\frac{1}{N} \sum_{i=1}^N (S_i - O_i)^2}$$ với $S_i$ là giá trị mô phỏng/dự báo và $O_i$ là giá trị quan trắc thực tế.

Data và phân tích

Nguồn dữ liệu thực nghiệm được thu thập đồng bộ và đa dạng:

  • Dữ liệu địa hình: Trích xuất từ cơ sở dữ liệu địa hình đáy biển toàn cầu độ phân giải cao kết hợp hải đồ chi tiết vùng biển Việt Nam do Trung tâm Quy hoạch và Điều tra tài nguyên – môi trường biển khu vực phía Bắc cung cấp.
  • Dữ liệu trường gió cưỡng bức: Dữ liệu tái phân tích khí quyển từ Trung tâm Dự báo Môi trường Quốc gia Hoa Kỳ (NCEP) và ECMWF, kết hợp mô hình WRF với bước thời gian 1–3 giờ.
  • Dữ liệu quan trắc kiểm định:
    1. 4 trạm phao khí tượng hải văn chuyên dụng: Phao 4001, Phao 4002, Phao 4003, Phao 4004 phân bố tại các vùng biển trọng điểm.
    2. Trạm hải văn giàn khoan cố định MSP1: Mỏ Bạch Hổ, đại diện cho chế độ sóng nước sâu ngoài khơi thềm lục địa phía Nam.
    3. Trạm Ra đa sóng bờ biển Đồng Hới (Quảng Bình): Cung cấp chuỗi số liệu độ cao sóng liên tục tại 6 vị trí trong tầm quét của ra đa.
    4. Dữ liệu vệ tinh Altimeter: Dữ liệu độ cao sóng dọc theo quỹ đạo quét của 2 vệ tinh Jason-2 và Saral-AltiKa cắt ngang qua khu vực Biển Đông.

Công cụ tính toán sử dụng mô hình SWAN phiên bản 40.85/41 hỗ trợ tính toán lưới phi cấu trúc tam giác và tính toán song song đa nhân (MPI/OpenMP), tích hợp trực tiếp với bộ khung đồng hóa mã nguồn mở OpenDA trên nền tảng Linux và Windows 64-bit.

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

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

Qua hàng loạt các thí nghiệm số trị chuyên sâu và đối soát thực nghiệm, luận án đã đạt được 5 phát hiện mang tính đột phá:

+-----------------------------------------------------------------------------+
|                CÁC PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN                            |
+-----------------------------------------------------------------------------+
| 1. Độ nhạy tham số: CDS2 & powst chi phối hoàn toàn GoF; powk & cutfr trơ   |
| 2. SWAN-CI tự động: Giảm RMSE từ 15-30% so với bộ tham số mặc định          |
| 3. EnKF tái phân tích: Triệt tiêu độ lệch ME, bám sát đỉnh sóng bão tại MSP1|
| 4. Hiệu lực dự báo: Đạt đỉnh trong 1-24h, duy trì ảnh hưởng tích cực đến 48h|
| 5. Tối ưu tổ hợp: Cấu hình 30 thành phần cân bằng hoàn hảo chính xác & chi phí|
+-----------------------------------------------------------------------------+
  1. Phân hóa độ nhạy tham số vật lý mô hình SWAN: Trích dẫn trực tiếp từ luận án: "Bằng việc áp dụng phương pháp hiệu chỉnh mô hình tự động SWAN-CI trong khu vực Biển Đông đã đưa ra được bộ tham số powk... và tham số cutfr... ít ảnh hưởng đến kết quả tính toán còn các tham số CDS2... và tham số powst... lại có ảnh hưởng rất lớn đến kết quả tính toán". Biến thiên của hàm mục tiêu GoF cho thấy sự thay đổi của $CDS2$ từ $1{,}0 \times 10^{-5}$ đến $4{,}5 \times 10^{-5}$ làm thay đổi mạnh mẽ độ cao sóng có nghĩa, đóng vai trò là "chốt chặn" điều tiết sai số năng lượng phổ sóng.
  2. Hiệu chỉnh tự động SWAN-CI vượt trội phương pháp truyền thống: Việc tối ưu hóa đa tham số đồng thời thông qua thuật toán DUD tại từng trạm phao riêng biệt ($\text{HC4001} \div \text{HC4004}$) và trên toàn miền tính ($\text{HC_All}$, $\text{HC_BH}$, $\text{HC_DH}$) đã làm giảm chỉ số RMSE từ $15%$ đến hơn $30%$ so với việc sử dụng bộ tham số mặc định (Default Parameters).
  3. Đồng hóa EnKF tái cấu trúc chính xác trường sóng phân tích: Kết quả chạy thử nghiệm với các phương án từ $\text{EnKF01}$ đến $\text{EnKF10}$ tại trạm MSP1 và Ra đa Đồng Hới chứng minh EnKF loại bỏ triệt để hiện tượng đánh giá thấp (underestimation) đỉnh sóng trong các đợt gió mùa đông bắc và bão mạnh. Chỉ số sai số trung bình $\text{ME}$ tiệm cận về giá trị $0\text{ m}$, và độ tán của các thành phần tổ hợp phản ánh trung thực mức độ bất định của trường gió.
  4. Quy luật suy giảm thời gian của hiệu ứng đồng hóa dự báo: Nghiên cứu khẳng định: "việc đồng hoá số liệu còn làm tăng độ chính xác của dự báo trong khoảng từ 1 đến 24 giờ cũng như tiếp tục cải thiện được kết quả dự báo trong khoảng từ 25 đến 48 giờ so với dự báo thông thường". Sau 48 giờ, điều kiện biên và trường gió khí quyển mới trở thành yếu tố chi phối hoàn toàn mô hình.
  5. Xác lập cấu hình tối ưu cho hệ thống đồng hóa: Luận án chỉ ra rằng số lượng thành phần tổ hợp $N_{ens} = 30$ là điểm cân bằng hoàn hảo giữa độ chính xác thống kê và chi phí tài nguyên tính toán; mức độ cập nhật thời gian 12 giờ kết hợp với bán kính ảnh hưởng không gian $500\text{ km}$ mang lại kết quả hội tụ ổn định nhất trên toàn Biển Đông.

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Cung cấp bằng chứng thực nghiệm khẳng định tính tương thích cao của lý thuyết lọc ngẫu nhiên EnKF đối với các phương trình thủy động lực sóng phi tuyến ở vùng biển nhiệt đới gió mùa phức tạp như Biển Đông.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình hai bước độc lập nhưng bổ trợ: (1) Hiệu chỉnh tham số vật lý nền bằng SWAN-CI trước khi (2) Đồng hóa trạng thái động lực bằng EnKF, loại bỏ triệt để hiện tượng "bù trừ sai số" (error compensation) giả tạo giữa tham số mô hình và trường điều kiện ban đầu.
  • Về mặt thực tiễn và tác nghiệp dự báo: Cung cấp giải pháp công nghệ trực tiếp cho Trung tâm Dự báo Khí tượng Thủy văn Quốc gia và các đài KTTV khu vực nhằm nâng cấp hệ sinh thái phần mềm dự báo sóng biển hàng ngày, nâng cao chất lượng cảnh báo thiên tai biển.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án cũng thẳng thắn chỉ ra các giới hạn nghiên cứu (boundary conditions):

  1. Độ bao phủ không gian của mạng lưới quan trắc: Số lượng trạm phao và trạm ra đa biển còn khá thưa thớt, chủ yếu tập trung ở ven bờ và khu vực thềm lục địa miền Trung và miền Nam, chưa có chuỗi số liệu quan trắc cố định dài hạn tại khu vực quần đảo Hoàng Sa và Trường Sa.
  2. Chi phí tính toán của phương pháp tổ hợp EnKF: Việc chạy mô phỏng đồng thời 30–50 thành phần tổ hợp SWAN đòi hỏi hệ thống máy chủ đa nhân hiệu năng cao, gây áp lực nhất định đối với tác nghiệp thời gian thực trong điều kiện hạ tầng công nghệ thông tin hạn chế.
  3. Giới hạn số liệu vệ tinh Altimeter: Dữ liệu vệ tinh Jason-2 và Saral-AltiKa có độ phân giải thời gian không liên tục (chu kỳ quay lại từ vài ngày đến hàng chục ngày), do đó đồng hóa vệ tinh chủ yếu đóng góp cho việc hiệu chỉnh trường tái phân tích hơn là cập nhật liên tục cho các phiên dự báo tác nghiệp 6 giờ một lần.

Chương trình nghiên cứu 10 năm tiếp theo:

  • Hướng 1: Ứng dụng các biến thể lọc Kalman tiên tiến như Lọc Kalman tổ hợp địa phương hóa biến đổi (Local Ensemble Transform Kalman Filter - LETKF) và Lọc quân phương tổ hợp (EnSRF) để giảm số lượng thành phần tổ hợp cần thiết xuống còn 10–20 mà vẫn giữ nguyên độ chính xác.
  • Hướng 2: Mở rộng đồng hóa toàn diện dữ liệu phổ sóng 2D thu nhận từ Ra đa khẩu độ tổng hợp (SAR) của vệ tinh Sentinel-1A/1B thay vì chỉ đồng hóa một biến duy nhất là độ cao sóng có nghĩa ($SWH$).
  • Hướng 3: Thiết lập mô hình tương tác sóng - dòng chảy - thủy triều ghép nối 3 chiều (Coupled SWAN - ROMS/ADCIRC) tích hợp đồng hóa đa biến.

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

Nghiên cứu của NCS. Nguyễn Trung Thành tạo ra các tác động lan tỏa sâu rộng:

  • Tác động học thuật: Đặt nền móng khoa học cho nhánh nghiên cứu đồng hóa số liệu sóng biển tại Việt Nam, mở ra hàng loạt đề tài nghiên cứu cấp Nhà nước và luận án tiến sĩ kế cận về hải dương học tính toán.
  • Tác động ngành và kinh tế biển: Nâng cao độ tin cậy của các bản tin dự báo sóng phục vụ trực tiếp cho hoạt động thăm dò, khai thác dầu khí của Tập đoàn Dầu khí Việt Nam (Petrovietnam) tại bể Cửu Long, Nam Côn Sơn; hỗ trợ thiết kế an toàn các công trình đê kè biển, cảng biển nước sâu và tối ưu hóa tuyến luồng vận tải hàng hải quốc tế.
  • Tác động an toàn xã hội và phòng chống thiên tai: Tăng thời gian cảnh báo sớm sóng lớn, sóng lừng và triều cường do bão trước 24–48 giờ, giúp ngư dân và các lực lượng cứu hộ, cứu nạn trên biển (Ủy ban Quốc gia Ứng phó Sự cố, Thiên tai và Tìm kiếm Cứu nạn) chủ động triển khai phương án bảo vệ tính mạng và tài sản.

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

+--------------------------------------------------------------------------+
|                      ĐỐI TƯỢNG HƯỞNG LỢI TRỰC TIẾP                       |
+--------------------------------------------------------------------------+
| 1. Nghiên cứu sinh & Nhà khoa học: Khung phương pháp SWAN-CI + EnKF      |
| 2. Cơ quan Dự báo Khí tượng Thủy văn: Quy trình công nghệ dự báo 48h     |
| 3. Ngành Dầu khí & Xây dựng công trình biển: Số liệu tái phân tích chuẩn |
| 4. Cơ quan Cứu hộ Cứu nạn & Hàng hải: Cảnh báo sớm sóng hiểm trở        |
| 5. Cơ quan Quản lý Biển & Hải đảo: Cơ sở khoa học quy hoạch không gian   |
+--------------------------------------------------------------------------+
  1. Nghiên cứu sinh và Giới học thuật Hải dương học: Thụ hưởng một bộ khung phương pháp luận hoàn chỉnh, các công thức toán học và mã nguồn tích hợp OpenDA-SWAN sẵn sàng cho việc mở rộng nghiên cứu.
  2. Cơ quan Dự báo Khí tượng Thủy văn Quốc gia: Tiếp nhận quy trình công nghệ chuẩn hóa để tích hợp trực tiếp vào hệ thống dự báo tác nghiệp hàng ngày.
  3. Doanh nghiệp Khai thác Cảng biển và Công trình Biển: Tiếp cận bộ dữ liệu tái phân tích trường sóng Biển Đông chính xác hơn để phục vụ tính toán tải trọng sóng thiết kế công trình ven bờ và ngoài khơi.
  4. Lực lượng Cảnh sát Biển, Hải quân và Cứu nạn: Nắm bắt trường sóng dự báo chi tiết theo từng giờ, nâng cao hiệu suất định vị và tìm kiếm cứu nạn trên biể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à việc làm sáng tỏ bản chất chi phối của cơ chế tiêu tán sóng bạc đầu trong mô hình sóng thế hệ III tại Biển Đông, mở rộng lý thuyết cân bằng năng lượng sóng của Komen et al. (1984) [41] và Hasselmann (1973) [36]. Luận án chứng minh rằng hàm mục tiêu GoF cực kỳ nhạy cảm với hệ số tốc độ tiêu tán $CDS2$ và số mũ độ dốc phổ $powst$, trong khi độc lập tương đối với $powk$ và $cutfr$, qua đó bác bỏ sự cần thiết của việc hiệu chỉnh tràn lan các tham số thứ cấp.

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

So với nghiên cứu của Nima Serpoushan et al. (2013) [57] (chỉ áp dụng EnKF với các tham số mặc định cho vịnh Persian) và Sofia Caires et al. (2018) [19] (chỉ đánh giá trường tái phân tích tại vùng biển Hàn Quốc), luận án đã: (1) Kết hợp đồng thời công cụ hiệu chỉnh phi tuyến DUD (SWAN-CI) trước khi đồng hóa; (2) Tích hợp đa nguồn dữ liệu thực nghiệm gồm phao, trạm giàn khoan ngoài khơi, trạm ra đa bờ biển và vệ tinh đo cao; và (3) Đánh giá chi tiết hiệu lực dự báo mở rộng từ 1 đến 48 giờ.

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

Phát hiện bất ngờ nhất là tham số tần số cắt $cutfr$ (tỷ lệ giữa tần số sóng cực đại và tần số trung bình) – vốn thường được coi là tham số quan trọng trong lý thuyết phổ JONSWAP – lại hầu như không làm thay đổi giá trị hàm mục tiêu sai số GoF trên toàn vùng Biển Đông khi biến thiên trong dải giá trị tiêu chuẩn, cho thấy phổ năng lượng vùng biển này chịu sự khống chế áp đảo của tương tác phi tuyến và trường gió cục bộ hơn là giới hạn tần số cao của phổ.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án mô tả chi tiết: (1) Tọa độ không gian và miền lưới tính SWAN; (2) Cấu trúc các file cấu hình XML trong OpenDA; (3) Thuật toán tạo trường gió nhiễu Monte Carlo; (4) Hàm mục tiêu GoF và các bước chia afin trong thuật toán DUD; và (5) Bộ tham số đồng hóa tối ưu ($N_{ens}=30$, bán kính $500\text{ km}$, bước cập nhật $12\text{h}$), đảm bảo khả năng tái lập kết quả $100%$ trên các hệ thống tính toán độc lập.

5. Khung chương trình nghiên cứu 10 năm được vạch ra như thế nào?

Luận án định hình lộ trình nghiên cứu hải dương học tính toán gồm: (1) Giai đoạn 1–3 năm: Tối ưu hóa thuật toán EnKF bằng sơ đồ LETKF trên máy tính song song GPU; (2) Giai đoạn 3–6 năm: Tích hợp đồng hóa dữ liệu phổ sóng 2D từ ảnh vệ tinh radar Sentinel-1 SAR; (3) Giai đoạn 6–10 năm: Xây dựng hệ thống dự báo tích hợp Hải văn - Khí tượng - Động lực bùn cát tự động hoàn toàn cho toàn bộ vùng biển đặc quyền kinh tế (EEZ) của Việt Nam.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Trung Thành là một công trình nghiên cứu công phu, mẫu mực và đột phá trong lĩnh vực Hải dương học vật lý tại Việt Nam. Những kết luận và đóng góp cốt lõi của công trình bao gồm:

  1. Đề xuất thành công bộ giá trị tham số vật lý tối ưu cho mô hình SWAN tại Biển Đông thông qua công cụ hiệu chỉnh tự động SWAN-CI, định lượng rõ vai trò áp đảo của các tham số tiêu tán sóng bạc đầu $CDS2$ và $powst$.
  2. Thiết lập hoàn chỉnh hệ thống đồng hóa số liệu EnKF-SWAN trên nền tảng mã nguồn mở OpenDA, giải quyết triệt để bài toán ước lượng trạng thái ban đầu mà không cần xây dựng mô hình tiếp tuyến liên hợp phức tạp.
  3. Chứng minh bằng số liệu thực nghiệm rằng đồng hóa dữ liệu làm triệt tiêu độ lệch sai số hệ thống, giảm mạnh chỉ số RMSE và nâng cao rõ rệt độ chính xác dự báo sóng biển trong hạn từ 1 đến 24 giờ, duy trì hiệu ứng tích cực đến 48 giờ.
  4. Khai thác và tích hợp thành công dữ liệu đa nguồn từ trạm phao, giàn khoan MSP1, trạm Ra đa Đồng Hới và vệ tinh Altimeter (Jason-2, Saral-AltiKa), tạo tiền đề cho công tác tái phân tích trường sóng chuẩn hóa.
  5. Đóng góp thiết thực cho khoa học dự báo tác nghiệp quốc gia, mở ra các hướng nghiên cứu hiện đại về đồng hóa dữ liệu đa biến và mô hình hóa tích hợp biển nhiệt đới.