Tổng quan về luận án

Nghiên cứu tính toán và dự báo sóng biển đóng vai trò then chốt trong bảo đảm an toàn hàng hải, khai thác tài nguyên dầu khí, công trình thủy công ven bờ và phòng chống thiên tai bão lũ. Trong hải dương học hiện đại, các mô hình số trị phổ sóng thế hệ ba như WAM, WAVEWATCH III và SWAN (Simulating Waves Nearshore) đã thay thế hoàn toàn các phương pháp kinh nghiệm cổ điển và đồ giải Synop. Tuy nhiên, tại khu vực Biển Đông – một vùng biển nhiệt đới gió mùa có chế độ khí động lực phức tạp và địa hình phân cắt mạnh – độ chính xác của các mô hình dự báo sóng đơn thuần thường suy giảm do sự bất định của trường gió khí quyển đầu vào, sai số tham số vật lý và việc thiếu hụt điều kiện ban đầu tối ưu.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: các nghiên cứu dự báo sóng tại Việt Nam trước đây chủ yếu dừng lại ở việc hiệu chỉnh thủ công mô hình số trị dựa trên kinh nghiệm chuyên gia cho từng chuỗi quan trắc ngắn hạn, chưa khai thác phương pháp hiệu chỉnh tự động đa tham số và chưa áp dụng kỹ thuật đồng hóa số liệu (data assimilation) đa nguồn. Trong khi đó, các trung tâm dự báo tiên tiến trên thế giới như ECMWF, NOAA/NCEP hay DWD đã tích hợp sâu rộng các sơ đồ đồng hóa số liệu vệ tinh và trạm phao. Để giải quyết khoảng trống này, luận án tập trung vào hai câu hỏi nghiên cứu then chốt:

  • RQ1: Làm thế nào để xác định khách quan và tối ưu bộ tham số vật lý của mô hình phổ sóng SWAN trên khu vực Biển Đông thông qua công cụ hiệu chỉnh tự động?
  • RQ2: Phương pháp đồng hóa số liệu lọc Kalman tổ hợp (EnKF) cải thiện độ chính xác của trường phân tích sóng và nâng cao chất lượng dự báo sóng ở các hạn dự báo từ 1 đến 48 giờ như thế nào khi tích hợp đồng thời số liệu phao biển, trạm ra đa và máy đo cao vệ tinh?

Hai giả thuyết khoa học tương ứng được xác lập:

  • H1: Áp dụng thuật toán tối ưu hóa tự động SWAN-CI sẽ bóc tách được độ nhạy cục bộ của các tham số tiêu tán năng lượng, làm giảm thiểu hàm mục tiêu sai số toàn miền so với bộ tham số mặc định.
  • H2: Sơ đồ lọc Kalman tổ hợp kết hợp mô hình sóng SWAN trong môi trường OpenDA cho phép cập nhật tối ưu ma trận hiệp phương sai sai số trường nền, tạo ra điều kiện ban đầu chính xác giúp giảm đáng kể sai số dự báo sóng trong 24 giờ đầu và duy trì hiệu quả tích cực đến 48 giờ.

Khung lý thuyết của nghiên cứu dựa trên phương trình cân bằng mật độ tác động sóng hai chiều (Wave Action Balance Equation) kết hợp với lý thuyết ước lượng trạng thái ngẫu nhiên phi tuyến (Stochastic Dynamic State Estimation). Luận án mang lại đóng góp đột phá khi thiết lập quy trình tích hợp hoàn chỉnh giữa công cụ hiệu chỉnh SWAN-CI và thuật toán EnKF, định lượng cụ thể mức độ suy giảm sai số RMSE (Root Mean Square Error) và ME (Mean Error) trên toàn miền tính Biển Đông. Phạm vi nghiên cứu bao quát toàn bộ lưu vực Biển Đông (tọa độ từ 0°N đến 25°N và 99°E đến 122°E) với hệ số liệu thực nghiệm đa dạng gồm 4 trạm phao (4001, 4002, 4003, 4004), trạm giàn khoan cố định MSP1 (mỏ Bạch Hổ), trạm Ra đa biển Đồng Hới (Quảng Bình) và các đường quét đo cao vệ tinh Altimeter (Jason-2, Saral-AltiKa, ENVISAT).

Literature Review và Positioning

Lịch sử phát triển mô hình tính sóng trải qua ba thế hệ gắn liền với các mốc lý thuyết nền tảng: từ quan hệ bán thực nghiệm của Sverdrup và Munk (1947), phân tích phổ ngẫu nhiên của Pierson và nnk (1955), khái niệm phương trình cân bằng năng lượng sóng thế hệ I của Gelci và nnk (1957), đến các cơ chế truyền năng lượng gió của Phillips (1957) và Miles (1957). Bước nhảy vọt tiếp theo là việc Hasselmann (1962, 1973) xây dựng hàm tương tác sóng-sóng phi tuyến bốn sóng (quadruplet) qua dự án phổ JONSWAP, định hình nên các mô hình thế hệ II và sau đó là xấp xỉ tương tác rời rạc DIA (Discrete Interaction Approximation) của Hasselmann và nnk (1985) mở đường cho mô hình phổ thế hệ III (WAMDI Group, 1988; Booij và nnk, 1999).

Song song với động lực học sóng, lý thuyết đồng hóa số liệu phát triển từ phương pháp nội suy tối ưu OI (Optimal Interpolation), lọc Kalman cổ điển (Kalman, 1960), biến phân 3D-Var/4D-Var (Le Dimet và Talagrand, 1986) cho đến lọc Kalman tổ hợp EnKF (Ensemble Kalman Filter) do Evensen (1994) đề xuất dựa trên mô phỏng ngẫu nhiên Monte Carlo. Trong hải dương học quốc tế, các cuộc tranh luận học thuật tập trung vào hai luồng quan điểm:

  • Nhóm phương pháp biến phân (4D-Var): Ưu điểm là tối ưu hóa quỹ đạo động lực liên tục trong cửa sổ đồng hóa nhưng đòi hỏi xây dựng mô hình tiếp tuyến liên hợp (adjoint model) vô cùng phức tạp, chi phí tính toán cực lớn đối với các quá trình phi tuyến cao như sóng vỡ và sóng bạc đầu (Orzech và nnk, 2012).
  • Nhóm phương pháp tổ hợp (EnKF, DEnKF, EnOI): Tránh được việc đạo hàm giải tích mô hình liên hợp bằng cách ước lượng ma trận hiệp phương sai sai số từ một tập hợp các trạng thái mô hình phân tán, rất linh hoạt khi tích hợp với các mô hình phổ sóng phi tuyến (Evensen, 1994; Burgers và nnk, 1998).

So sánh với các công trình quốc tế tiêu biểu:

  • Nghiên cứu của Chen (2004) trên mô hình toàn cầu NWW3 chỉ ra rằng đồng hóa số liệu phao và vệ tinh ERS-2 cải thiện rõ rệt dự báo độ cao sóng có nghĩa (SWH) trong 12 giờ đầu, duy trì mức trung bình trong 12–24 giờ và suy giảm mạnh sau 24 giờ.
  • Nghiên cứu của Nima Serpoushan và nnk (2013) tại vịnh Persian sử dụng SWAN kết hợp EnKF trong OpenDA, chỉ hiệu chỉnh hai tham số sóng bạc đầu $C_{dis} = 3.95$ và $\delta = 0.6$, khuyến nghị tổ hợp từ 50 đến 100 thành phần nhưng chỉ dừng lại ở bài toán tái phân tích (hindcast) mà chưa đánh giá hạn dự báo (forecast).
  • Nghiên cứu của Sofia Caires và nnk (2018) tại vùng biển phía đông Hàn Quốc áp dụng SWAN-EnKF với 30 thành phần tổ hợp, độ tán gió 1 m/s, bán kính tương quan không gian 500 km, nhưng chỉ sử dụng số liệu gió WRF và trạm ven bờ.
  • Nghiên cứu của Almeida và nnk (2015) tại bờ biển Bồ Đào Nha thử nghiệm SWAN-EnKF 20 thành phần cho 4 trạm phao ven bờ, khẳng định hiệu quả vượt trội so với mô phỏng nền.

Tại Việt Nam, các công trình của Huỳnh Thị Hồng Ngự và nnk (2008), Kiều Quốc Chánh (2011), Trần Tân Tiến và nnk (2011), Hoàng Đức Cường và nnk (2012), Dư Đức Tiến và nnk (2017) đã ứng dụng 3D-Var, EnKF, LETKF chủ yếu cho mô hình khí quyển WRF dự báo bão và mưa lớn. Trong lĩnh vực sóng biển, Nguyễn Minh Huấn (2010) và Trần Quang Tiến (2014) đã bước đầu thử nghiệm lọc Kalman cho mô hình sóng ven bờ nhưng chỉ giới hạn ở 01 chuỗi số liệu ra đa cục bộ, chưa có công cụ hiệu chỉnh tham số tự động và chưa đánh giá hệ thống tác động của đồng hóa số liệu đa nguồn trên toàn Biển Đông. Luận án này định vị chính xác điểm giao thoa giữa hiệu chỉnh tham số tự động SWAN-CI và đồng hóa số liệu tổ hợp EnKF, tạo bước tiến mới trong mô hình hóa hải dương học tại khu vực 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 mở rộng và kiểm chứng lý thuyết cân bằng tác động phổ sóng trong môi trường biển nhiệt đới thông qua việc phân định độ nhạy vật lý của các cơ chế tiêu tán năng lượng. Nghiên cứu kiểm chứng trực tiếp công thức tiêu tán do sóng bạc đầu của Komen và nnk (1984) được tích hợp trong mô hình SWAN:

$$S_{ds,w}(\sigma, \theta) = - C_{DS2} \left( \frac{\alpha}{\alpha_{PM}} \right)^{powst} \left( \frac{k}{\tilde{k}} \right)^{powk} \left( \frac{\sigma}{\tilde{\sigma}} \right)^{cutfr} \tilde{\sigma} E(\sigma, \theta)$$

Trong đó, luận án làm sáng tỏ hai mệnh đề lý thuyết quan trọng:

  • Mệnh đề 1: Độ nhạy của trường phổ sóng Biển Đông bị chi phối áp đảo bởi hệ số tốc độ tiêu tán sóng bạc đầu ($CDS2$) và tỉ lệ độ dốc phổ chuẩn hóa Pierson-Moskowitz ($powst$), trong khi các số mũ chuẩn hóa số sóng ($powk$) và tỉ lệ tần số cắt ($cutfr$) có mức độ ảnh hưởng không đáng kể đến độ cao sóng có nghĩa.
  • Mệnh đề 2: Lý thuyết ước lượng tuyến tính không thiên lệch tối ưu (BLUE) trong không gian trạng thái tổ hợp EnKF cho phép hiệu chỉnh đồng thời trường sóng và trường gió cưỡng bức, bù đắp sự thiếu hụt thông tin động lực học tại các vùng biển mở không có trạm quan trắc cố định.

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

Khung phân tích của luận án tích hợp ba trụ cột lý thuyết: Động lực học sóng phổ thế hệ III (Booij và nnk, 1999), Lý thuyết tối ưu hóa toán học không đạo hàm DUD (Ralston và Jennrich, 1978) và Lý thuyết lọc ngẫu nhiên EnKF (Evensen, 1994). Tính độc đáo thể hiện ở quy trình hai bước tuần tự có kiểm soát:

  1. Hiệu chuẩn tham số vật lý nền tảng: Sử dụng công cụ SWAN-CI tự động điều chỉnh các hệ số nguồn/tiêu tán dựa trên cực tiểu hóa hàm sai số $GoF$ (Goodness of Fit) trước khi đưa vào vận hành mô hình.
  2. Đồng hóa số liệu động lực học: Sử dụng nền tảng phần mềm mã nguồn mở OpenDA liên kết với SWAN để cập nhật trường sóng và điều kiện biên tại từng chu kỳ phân tích, xử lý triệt để bài toán phi tuyến mà không làm mất tính ổn định số của sơ đồ sai phân ẩn BSBT (Backward Space - Backward Time).

Điều kiện biên giới hạn (boundary conditions) của khung phân tích được xác lập rõ: áp dụng cho vùng biển sâu và vùng ven bờ chịu ảnh hưởng của gió mùa Đông Bắc và Tây Nam; các quá trình tương tác sóng - dòng chảy mạnh tại cửa sông và biến động triều quy mô cực nhỏ được xem là nằm ngoài phạm vi chính của bài toán phổ quy mô Biển Đông.

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng (positivism) kết hợp chủ nghĩa hiện thực phản biện (critical realism), dựa trên việc kiểm chứng định lượng các mô hình toán học - vật lý bằng chuỗi dữ liệu thực nghiệm quan trắc độc lập. Thiết kế nghiên cứu kết hợp giữa mô phỏng số trị đa quy mô và phân tích thống kê tổ hợp Monte Carlo.

Miền tính toán bao phủ toàn bộ Biển Đông từ vĩ độ 0°N đến 25°N, kinh độ 99°E đến 122°E với độ phân giải không gian $0.1^\circ \times 0.1^\circ$ (khoảng $11 \times 11\text{ km}$). Lưới địa hình đáy biển được trích xuất từ cơ sở dữ liệu toàn cầu ETOPO1/GEBCO kết hợp bình đồ độ sâu ven bờ chi tiết của Việt Nam. Trường gió cưỡng bức mặt biển được thu thập từ mô hình tái phân tích khí quyển ERA-Interim của ECMWF và mô hình dự báo hạn vừa NCEP/GFS với độ phân giải thời gian 3–6 giờ.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được thực hiện qua các giao thức chuẩn hóa nghiêm ngặt:

  • Sơ đồ số trong SWAN: Sử dụng sơ đồ sai phân ẩn BSBT trong không gian địa lý $(x, y)$ và phép xấp xỉ trung tâm trong không gian phổ tần số - hướng sóng $(\sigma, \theta)$, với 24 hướng phổ ($15^\circ/\text{hướng}$) và 25 dải tần số logarit từ 0.04 Hz đến 1.0 Hz.
  • Giao thức hiệu chỉnh tự động SWAN-CI: Ứng dụng thuật toán DUD (Doesn't Use Derivative) để tìm kiếm cực tiểu của hàm sai số $GoF$:

$$GoF = \frac{\sum_{i=1}^{N_i} w_i \left( \frac{SWH_{sim,i} - SWH_{obs,i}}{\sigma_{obs}} \right)^2}{\sum_{i=1}^{N_i} w_i}$$

Thuật toán DUD xấp xỉ ma trận Jacobian thông qua phép nội suy afin giữa $N+1$ đánh giá trong bộ nhớ, loại bỏ nhu cầu tính đạo hàm giải tích phức tạp của phương trình sóng.

  • Giao thức đồng hóa số liệu EnKF: Hệ thống lọc Kalman tổ hợp được thiết lập trong OpenDA. Trạng thái phân tích $\mathbf{x}^a$ tại mỗi bước thời gian được tính toán thông qua công thức BLUE:

$$\mathbf{x}^a = \mathbf{x}^f + \mathbf{K} (\mathbf{y} - \mathbf{H} \mathbf{x}^f)$$

$$\mathbf{K} = \mathbf{P}^f \mathbf{H}^T (\mathbf{H} \mathbf{P}^f \mathbf{H}^T + \mathbf{R})^{-1}$$

Trong đó $\mathbf{x}^f$ là véctơ trạng thái dự báo của các thành phần tổ hợp, $\mathbf{y}$ là véctơ quan trắc, $\mathbf{H}$ là toán tử quan trắc ánh xạ từ không gian lưới mô hình sang vị trí đo, $\mathbf{R}$ là ma trận hiệp phương sai sai số quan trắc, $\mathbf{P}^f$ là ma trận hiệp phương sai sai số trường nền ước lượng từ $N_{ens}$ thành phần tổ hợp Monte Carlo.

Tam giác đạc dữ liệu (data triangulation) được bảo đảm thông qua 4 nguồn số liệu độc lập:

  1. Trạm phao hải văn quốc tế: Các trạm 4001, 4002, 4003, 4004 đo đạc liên tục độ cao sóng có nghĩa và chu kỳ sóng.
  2. Trạm hải văn cố định: Trạm giàn khoan MSP1 (Bạch Hổ) tại thềm lục địa phía Nam.
  3. Trạm Ra đa sóng ven bờ: Ra đa biển Đồng Hới (Quảng Bình) đo trường sóng mặt rộng ven bờ miền Trung.
  4. Đo cao vệ tinh Altimeter: Vệ tinh Jason-2 và Saral-AltiKa cung cấp các vết quét độ cao sóng dọc quỹ đạo qua Biển Đông.

Data và phân tích

Các chỉ số thống kê định lượng nghiêm ngặt được sử dụng để đánh giá sai số giữa kết quả mô hình và thực đo:

  • Sai số trung bình (Mean Error - ME):

$$ME = \frac{1}{N} \sum_{i=1}^{N} (SWH_{sim,i} - SWH_{obs,i})$$

  • Sai số căn bậc hai trung bình bình phương (Root Mean Square Error - RMSE):

$$RMSE = \sqrt{\frac{1}{N} \sum_{i=1}^{N} (SWH_{sim,i} - SWH_{obs,i})^2}$$

  • Hệ số tương quan ($R$), chỉ số tán xạ ($SI$) và hàm phù hợp $GoF$.

Các thử nghiệm độ nhạy về cấu hình EnKF được tiến hành với quy mô tổ hợp biến thiên $N_{ens} \in {10, 20, 30, 50, 100}$, mức độ tán gió (wind spread) từ 0.5 m/s đến 2.0 m/s, chu kỳ cập nhật đồng hóa từ 3 giờ đến 12 giờ, và bán kính ảnh hưởng không gian (localization radius) từ 100 km đến 500 km.

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

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

Kết quả mô phỏng và kiểm định trong luận án mang lại 4 phát hiện mang tính đột phá:

Thứ nhất, phân tích độ nhạy đa tham số bằng SWAN-CI đã chứng minh sự phân hóa rõ nét giữa hai nhóm tham số trong phương trình sóng bạc đầu:

  • Nhóm chi phối mạnh: Tham số tốc độ tiêu tán sóng bạc đầu $CDS2$ và số mũ độ dốc phổ $powst$. Giá trị $GoF$ biến thiên mạnh theo dạng đường cong parabol khi thay đổi $CDS2$ (từ $1.0 \times 10^{-5}$ đến $4.5 \times 10^{-5}$) và $powst$ (từ 1.0 đến 3.5), cho thấy đây là hai biến điều khiển cốt lõi quyết định biên độ phổ sóng Biển Đông.
  • Nhóm ít nhạy cảm: Tham số chuẩn hóa số sóng $powk$ và tỉ lệ tần số cắt $cutfr$ hầu như không làm thay đổi giá trị hàm mục tiêu $GoF$ khi quét trên toàn dải giá trị vật lý khả dĩ.

Thứ hai, công cụ hiệu chỉnh tự động SWAN-CI xác lập được bộ tham số tối ưu đặc thù cho từng tiểu vùng và bộ tham số tổng thể toàn Biển Đông, vượt trội hoàn toàn so với tham số mặc định của mô hình (MD: $CDS2 = 2.36 \times 10^{-5}, powst = 2.0$):

  • Tại trạm phao 4001: Bộ tham số hiệu chỉnh tối ưu đạt $CDS2 = 3.82 \times 10^{-5}, powst = 1.45$.
  • Tại trạm phao 4002: $CDS2 = 3.15 \times 10^{-5}, powst = 1.62$.
  • Tại trạm phao 4003: $CDS2 = 2.94 \times 10^{-5}, powst = 1.78$.
  • Tại trạm phao 4004: $CDS2 = 3.41 \times 10^{-5}, powst = 1.55$.
  • Phương án tổng hợp toàn miền (HC_All): Xác định được dải tối ưu chung với $CDS2 \approx 3.20 \times 10^{-5}$ và $powst \approx 1.60$.

Thứ ba, việc áp dụng đồng hóa số liệu EnKF mang lại sự cải thiện vượt bậc về chỉ số thống kê sai số tại tất cả các trạm quan trắc thực nghiệm:

  • Tại trạm phao 4001: Phương án mô phỏng mặc định (NoEnKF) có $RMSE = 0.58\text{ m}, ME = -0.32\text{ m}$; khi áp dụng hiệu chỉnh tự động (HC4001) $RMSE$ giảm xuống $0.44\text{ m}$; khi tích hợp đồng hóa EnKF (EnKF03), $RMSE$ giảm mạnh chỉ còn $0.28\text{ m}$ (giảm trên 51% sai số so với mặc định) và $ME$ tiến sát về $0.03\text{ m}$.
  • Tại trạm giàn khoan MSP1: $RMSE$ trong giai đoạn đồng hóa giảm từ $0.62\text{ m}$ (mô hình gốc) xuống còn $0.31\text{ m}$ (phương án EnKF), triệt tiêu hiện tượng đánh giá thấp (underestimation) độ cao sóng trong các đợt gió mùa mạnh.
  • Tại trạm Ra đa Đồng Hới: Sai số $RMSE$ độ cao sóng giảm từ $0.52\text{ m}$ xuống $0.29\text{ m}$, trường hướng sóng và chu kỳ sóng tái hiện trơn tru và trùng khớp với phổ quan trắc ra đa.
  • Kiểm định dọc vết quét vệ tinh Jason-2 và Saral-AltiKa: Sai số $RMSE$ trên toàn tuyến vệ tinh giảm trung bình từ $0.65\text{ m}$ xuống $0.34\text{ m}$.

Thứ tư, xác định chính xác quy luật suy giảm theo thời gian (temporal decay) của hiệu ứng đồng hóa số liệu trong dự báo sóng hạn ngắn:

  • Trong khoảng thời gian từ 1 đến 24 giờ đầu ($T+1$ đến $T+24$): Đồng hóa số liệu EnKF tạo ra tác động mạnh mẽ nhất, biên độ giảm sai số $RMSE$ đạt từ 35% đến 52% so với dự báo thông thường không đồng hóa.
  • Trong khoảng thời gian từ 25 đến 48 giờ ($T+25$ đến $T+48$): Mặc dù mô hình dần bị chi phối trở lại bởi sai số của trường gió khí quyển dự báo, hiệu ứng của trường điều kiện ban đầu tối ưu từ EnKF vẫn tiếp tục duy trì tác động tích cực, giúp giảm từ 12% đến 20% sai số $RMSE$ so với phương án NoEnKF.

Implications đa chiều

  • Đóng góp học thuật: Cung cấp cơ sở lý luận và phương pháp luận hoàn chỉnh để ghép nối mô hình phổ sóng phi tuyến với hệ thống đồng hóa ngẫu nhiên OpenDA; làm sáng tỏ cơ chế tiêu tán năng lượng sóng trong điều kiện khí hậu hải văn Biển Đông.
  • Ứng dụng thực tiễn: Hoàn thiện quy trình công nghệ phục vụ dự báo sóng nghiệp vụ tại Trung tâm Dự báo Khí tượng Thủy văn Quốc gia và Viện Hải dương học; nâng cao năng lực cảnh báo sớm hiểm họa sóng lớn cho tàu thuyền, giàn khoan dầu khí (Vietsovpetro) và các hoạt động kinh tế biển.
  • Đề xuất chính sách: Khuyến nghị cơ quan quản lý nhà nước (Tổng cục Khí tượng Thủy văn, Cục Biển và Hải đảo Việt Nam) đầu tư mở rộng mạng lưới trạm phao hải văn ngầm định kỳ và tích hợp luồng dữ liệu vệ tinh thời gian thực vào quy chuẩn dự báo hải văn quốc gia.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận một số hạn chế:

  • Hạn chế về tương tác sóng - dòng chảy: Mô hình chưa xét đến ảnh hưởng hai chiều (two-way coupling) giữa sóng và hoàn lưu dòng chảy quy mô lớn do thiếu mô hình thủy động lực học ghép nối đồng thời trong hệ thống đồng hóa.
  • Giới hạn về mật độ dữ liệu: Tần suất quét của vệ tinh đo cao Altimeter qua Biển Đông còn thưa theo không gian và thời gian; số lượng trạm phao cố định ở vùng biển sâu còn hạn chế.
  • Tài nguyên tính toán: Phương pháp EnKF với số thành phần tổ hợp lớn ($N_{ens} \ge 50$) đòi hỏi năng lực tính toán hiệu năng cao (HPC), gây thách thức cho việc vận hành tác nghiệp thời gian thực nếu không tối ưu hóa mã nguồn song song.

Chương trình nghiên cứu tiếp theo (4 hướng đi cụ thể):

  1. Phát triển hệ thống ghép nối hoàn lưu - sóng - khí quyển ba chiều (WRF - ROMS - SWAN) tích hợp sơ đồ đồng hóa tổ hợp biến phân lai (Hybrid EnVar/LETKF).
  2. Đồng hóa trực tiếp phổ sóng hai chiều 2D trích xuất từ Ra đa khẩu độ tổng hợp SAR (Synthetic Aperture Radar) của vệ tinh Sentinel-1 thay vì chỉ đồng hóa độ cao sóng có nghĩa dạng vô hướng (SWH).
  3. Ứng dụng trí tuệ nhân tạo (AI/Physics-Informed Neural Networks - PINNs) để tăng tốc độ tính toán ma trận hiệp phương sai sai số trường nền trong EnKF.
  4. Mở rộng thử nghiệm đồng hóa trong các tình huống thời tiết cực đoan như siêu bão, áp thấp nhiệt đới và không khí lạnh tăng cường mạnh.

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

  • Tác động học thuật: Luận án tạo lập một bộ số liệu tái phân tích trường sóng chuẩn hóa cho Biển Đông, cung cấp tài liệu tham khảo cho các nghiên cứu sinh chuyên ngành Hải dương học, Khí tượng học và Thủy khí động lực học biển.
  • Chuyển đổi ngành công nghiệp: Cung cấp giải pháp tính toán tải trọng sóng chính xác hơn cho thiết kế công trình biển, tuyến luồng hàng hải quốc tế, tối ưu hóa tuyến hành trình tàu biển nhằm tiết kiệm nhiên liệu và giảm thiểu rủi ro tai nạn hàng hải.
  • Tác động xã hội và môi trường: Nâng cao độ chính xác cảnh báo sóng phục vụ công tác tìm kiếm cứu nạn hàng hải (SAR), bảo vệ an toàn sinh mạng ngư dân và hỗ trợ quản lý đới bờ bền vững trước tác động của biến đổi khí hậu và nước dâng do bão.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận hiệu chuẩn tự động SWAN-CI và thuật toán lọc Kalman tổ hợp mã nguồn mở, khai thác các khoảng trống nghiên cứu về đồng hóa số liệu hải dương.
  • 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ệ đồng hóa số liệu đã được 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.
  • Doanh nghiệp Khai thác Dầu khí và Năng lượng Tái tạo Biển: Tối ưu hóa kế hoạch vận hành giàn khoan, lắp đặt tuabin gió ngoài khơi dựa trên các bản tin dự báo sóng có độ tin cậy cao từ 1 đến 48 giờ.
  • Cơ quan Ứng phó Sự cố và Cứu nạn Hàng hải: Khai thác trường sóng tái phân tích và dự báo độ chính xác cao để mô phỏng quỹ đạo trôi dạt của nạn nhân và phương tiện chìm đắm 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ì?

Luận án đã mở rộng lý thuyết cân bằng tác động phổ sóng của Hasselmann (1962) và Booij và nnk (1999) bằng cách chứng minh trên thực nghiệm Biển Đông rằng cơ chế tiêu tán năng lượng sóng bạc đầu bị kiểm soát chủ yếu bởi cặp tham số ${CDS2, powst}$, đồng thời bác bỏ sự cần thiết của việc tinh chỉnh phức tạp các tham số ${powk, cutfr}$. Đây là bằng chứng vật lý quan trọng giúp tinh giản không gian tìm kiếm tham số tối ưu trong mô hình hóa phổ sóng nhiệt đới.

2. Tính đổi mới trong phương pháp luận so với các nghiên cứu trước đây?

So với nghiên cứu của Nima Serpoushan và nnk (2013) hay Sofia Caires và nnk (2018) vốn chỉ sử dụng tham số mặc định của EnKF cho bài toán tái phân tích cục bộ, luận án đã: (1) Kết hợp tuần tự hai bước giữa hiệu chuẩn tự động đa tham số SWAN-CI (dùng thuật toán DUD) và đồng hóa động lực học EnKF trong OpenDA; (2) Thực hiện tam giác đạc số liệu đồng thời từ 4 nguồn độc lập (phao, giàn khoan MSP1, ra đa Đồng Hới, vệ tinh Jason-2/Saral-AltiKa); (3) Đánh giá định lượng tác động của đồng hóa đến chất lượng dự báo mở rộng từ 1 đến 48 giờ.

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à việc đồng hóa dữ liệu đo cao vệ tinh dọc quỹ đạo hẹp vẫn có khả năng lan tỏa hiệu chỉnh ra toàn trường sóng xung quanh trong bán kính lên tới hàng trăm kilomet nhờ cấu trúc ma trận hiệp phương sai sai số động $\mathbf{P}^f$ trong EnKF, giúp duy trì mức giảm sai số $RMSE$ từ 12% đến 20% ở hạn dự báo 25–48 giờ, vượt xa giả định trước đây cho rằng hiệu ứng đồng hóa sóng sẽ tiêu tán hoàn toàn sau 24 giờ.

4. Luận án có cung cấp giao thức tái lập nghiên cứu (Replication Protocol) không?

Có. Luận án mô tả chi tiết toàn bộ cấu hình lưới tính SWAN ($0.1^\circ \times 0.1^\circ$), 24 hướng phổ, 25 dải tần số, các tham số thuật toán DUD trong OpenDA, hàm chi phí $GoF$, các thiết lập nhiễu ngẫu nhiên trường gió Monte Carlo ($wind\ spread = 1.0\text{ m/s}$), bán kính định vị không gian và ma trận phương sai sai số quan trắc $\mathbf{R}$, bảo đảm khả năng tái lập hoàn toàn trên các hệ thống tính toán độc lập.

5. Chương trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Lộ trình 10 năm tập trung vào: (1) Xây dựng hệ thống đồng hóa đại dương tích hợp đa thành phần (khí quyển - sóng - dòng chảy - triều) quy mô toàn diện Biển Đông; (2) Tích hợp luồng dữ liệu lớn từ vệ tinh SWOT và máy bay không người lái hải văn; (3) Thay thế toán tử phân tích truyền thống bằng các mô hình học sâu biến phân (Deep Variational Data Assimilation) nhằm tiến tới dự báo sóng thời gian thực với độ phân giải siêu cao 1 km.

Kết luận

  1. Luận án đã giải quyết trọn vẹn bài toán nâng cao độ chính xác mô phỏng và dự báo sóng biển khu vực Biển Đông thông qua việc tích hợp phương pháp hiệu chỉnh tự động SWAN-CI và phương pháp đồng hóa số liệu lọc Kalman tổ hợp (EnKF) trên nền tảng phần mềm OpenDA.
  2. Xác lập được bản chất vật lý và độ nhạy của các tham số mô hình phổ sóng SWAN, chứng minh cặp tham số $CDS2$ và $powst$ có ý nghĩa quyết định trong khi $powk$ và $cutfr$ ít ảnh hưởng, từ đó đề xuất thành công bộ tham số tối ưu toàn miền ($CDS2 \approx 3.20 \times 10^{-5}, powst \approx 1.60$).
  3. Ứng dụng thành công sơ đồ EnKF đa nguồn (kết hợp đồng thời số liệu phao 4001–4004, trạm giàn khoan MSP1, trạm Ra đa Đồng Hới và vệ tinh đo cao Jason-2, Saral-AltiKa), giúp giảm sai số $RMSE$ độ cao sóng có nghĩa từ 35% đến trên 51% trong giai đoạn đồng hóa và tái phân tích.
  4. Định lượng chính xác hiệu năng cải thiện của đồng hóa số liệu trong chu trình dự báo: nâng cao vượt bậc độ chính xác dự báo trong hạn ngắ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ờ.
  5. Đặt nền móng phương pháp luận vững chắc cho công nghệ dự báo hải văn số trị hiện đại tại Việt Nam, mở ra hướng ứng dụng mở rộng cho các bài toán thủy động lực học, lan truyền ô nhiễm và tìm kiếm cứu nạn trên biển.
  6. Kết quả nghiên cứu có giá trị khoa học và thực tiễn cao, sẵn sàng chuyển giao phục vụ công tác nghiệp vụ dự báo khí tượng thủy văn quốc gia và chiến lược phát triển kinh tế biển bền vững của đất nước.