Tổng quan về luận án

Nghiên cứu về hải dương học động lực và trắc địa biển tại Việt Nam bước vào giai đoạn chuyển mình then chốt với công trình luận án tiến sĩ kỹ thuật trắc địa - bản đồ của nghiên cứu sinh Đỗ Văn Mong (2025), thực hiện dưới sự hướng dẫn của PGS.TS Nguyễn Văn Sáng và TS Nguyễn Đình Thành tại Trường Đại học Mỏ - Địa chất Hà Nội. Công trình mang tên "Nghiên cứu phương pháp xác định dòng chảy bề mặt trên vùng biển Việt Nam từ số liệu đo cao vệ tinh phục vụ cho một số hoạt động của Hải quân Việt Nam" (Mã số chuyên ngành: 9520503). Đề tài được triển khai trong khuôn khổ Hợp đồng khoa học số 07/2021/Đ6-DATS thuộc dự án cấp thiết: "Điều tra tổng hợp khí tượng, các yếu tố hải dương, địa chất, môi trường vùng biển Trường Sa tỷ lệ 1:200.000", với sự bảo trợ thực địa từ Đoàn Đo đạc biên vẽ hải đồ và Nghiên cứu biển thuộc Bộ Tham mưu Quân chủng Hải quân.

Bối cảnh khoa học của đề tài bắt nguồn từ tính chất phức tạp của hệ thống hoàn lưu Biển Đông – vùng biển có diện tích lãnh hải và đặc quyền kinh tế hơn 1 triệu km² với bờ biển dài 3.260 km, giữ vị trí địa chính trị và an ninh quốc phòng tối quan trọng. Trong văn bản luận án, tính quy luật tự nhiên được khẳng định rõ nét: "Dòng chảy trên biển là yếu tố tự nhiên rất phức tạp, được sinh ra bởi nhiều yếu tố như: Gió, thủy triều, chênh lệch nhiệt độ, độ mặn… Dòng chảy có vai trò to lớn trong đời sống đại dương như: làm tăng sự trao đổi nước, phân bố lại nhiệt độ, độ muối". Tuy nhiên, rào cản nghiên cứu (research gap) tồn tại dai dẳng tại Việt Nam là việc khảo sát dòng chảy trước đây chỉ thuần túy dựa vào hai hướng tiếp cận truyền thống vốn tồn tại nhiều hạn chế:

  1. Quan trắc trực tiếp bằng tàu đo đạc chuyên dụng và phao trôi (drifter buoys): Dù đạt độ chính xác cao nhưng kinh phí cực lớn, thời gian đo ngắn hạn, mật độ không gian thưa thớt, không thể bao quát toàn bộ vùng biển rộng lớn và hoàn toàn bất khả thi tại các khu vực nhạy cảm, cấu trúc địa hình ngầm phức tạp hoặc các đảo bị chiếm đóng trái phép.
  2. Mô phỏng số trị thủy động lực (như mô hình ROMS, SWAN): Phụ thuộc nặng nề vào độ chính xác của các dữ liệu khí tượng biên mở (thông lượng nhiệt, áp suất, trường gió QuikSCAT/COAMPS) và giả định địa hình đáy biển, dẫn đến độ nhạy sai số tích lũy cao khi ngoại suy dài hạn.

Chưa từng có một công trình quốc nội nào xây dựng quy trình chuẩn hóa từ xử lý chuỗi số liệu đo cao radar vệ tinh (Satellite Radar Altimetry) đa nhiệm vụ để trích xuất trường dòng chảy địa tĩnh bề mặt (Surface Geostrophic Currents - SGC) có độ phân giải cao phục vụ chuyên biệt cho quốc phòng - an ninh biển.

Luận án thiết lập 3 câu hỏi nghiên cứu và hệ giả thuyết khoa học tương ứng được đánh số chặt chẽ:

  • Research Question 1 (RQ1): Làm thế nào để loại bỏ triệt để các sai số tồn dư (quỹ đạo, môi trường truyền sóng, triều đại dương và triều vỏ đất rắn) nhằm xác định chính xác độ cao địa hình mặt biển động lực biến đổi theo thời gian (Time-varying Dynamic Ocean Topography) từ các vết quét vệ tinh đo cao thế hệ mới?
    • Hypothesis 1 (H1): Phương pháp bình sai mạng lưới điểm giao cắt (Crossover Adjustment) kết hợp với mô hình mặt biển trung bình toàn cầu độ phân giải cao (DTU15MSS) và lọc sai số thô theo phân phối chuẩn Gauss sẽ giảm độ lệch chuẩn trung phương sai số đo ($m_{\text{đo}}$) xuống dưới ngưỡng $\pm 3\text{ cm} - \pm 5\text{ cm}$.
  • Research Question 2 (RQ2): Quy trình nội suy tối ưu nào cho phép chuyển đổi trường dị thường mực nước biển (Sea Level Anomaly - SLA) dạng vết quét phân tán thành lưới ô vuông liên tục phục vụ tính toán vector gradient dòng chảy địa tĩnh?
    • Hypothesis 2 (H2): Thuật toán khớp hàm hiệp phương sai thực nghiệm (Empirical Covariance Function) với hàm lý thuyết thông qua phương pháp Least-Squares Collocation sẽ tối ưu hóa độ dốc mặt biển tức thời, phản ánh trung thực cấu trúc xoáy đại dương trung mô (mesoscale eddies).
  • Research Question 3 (RQ3): Dữ liệu dòng chảy bề mặt từ đo cao vệ tinh có đáp ứng độ tin cậy động lực học khi ứng dụng trực tiếp vào công tác dẫn đường hàng hải và triển khai vũ khí của lực lượng Hải quân?
    • Hypothesis 3 (H3): Tích hợp trường vector dòng chảy vệ tinh lên hệ thống hải đồ điện tử ECDIS-DGPS cho phép bù dạt hướng đi chính xác, giảm thiểu sai số quỹ đạo tàu chiến đấu và vũ khí ngầm (ngư lôi, thủy lôi) trên 15-20% so với phương pháp tính dòng thông thường.

Khung lý thuyết của nghiên cứu tích hợp liên ngành giữa Lý thuyết Cân bằng Địa tĩnh (Geostrophic Balance Theory) trong hải dương học động lực, Lý thuyết Trường Trọng lực & Geoid (Physical Geodesy Theory), và Lý thuyết Xử lý Số liệu Trắc địa Vật lý. Phạm vi không gian của luận án bao trùm toàn bộ Biển Đông từ vĩ độ $5^\circ\text{N}$ đến $25^\circ\text{N}$ và kinh độ $105^\circ\text{E}$ đến $120^\circ\text{E}$. Chuỗi dữ liệu thực nghiệm bao gồm 34 chu kỳ lặp của vệ tinh SARAL/AltiKa (băng tần Ka 35 GHz) và 37 chu kỳ của vệ tinh Sentinel-3A (với các chu kỳ đại diện chuyên sâu CK71, CK72, CK75, CK77, CK81, CK84 trong giai đoạn 2021–2022), được kiểm chứng đối sánh với chuỗi đo đạc thực tế từ phao trôi hải dương học và hải trình tàu mặt nước của Đoàn 6 Hải quân.

Literature Review và Positioning

Nghiên cứu định vị trong dòng chảy lịch sử của công nghệ đo cao viễn thám radar đại dương từ những thập niên cuối thế kỷ 20 đến nay. Trải qua hơn ba thập kỷ, các chương trình đo cao vệ tinh quốc tế như TOPEX/Poseidon (phóng năm 1992 hợp tác NASA/CNES), Jason-1 (2001), Jason-2/OSTM (2008), SARAL/AltiKa (2013, hợp tác ISRO/CNES), CryoSat-2 (2010, ESA) và Sentinel-3A (2016, ESA/Copernicus) đã tạo nên một cuộc cách mạng trong nhận thức về động lực học mặt biển toàn cầu.

Tổng quan tài liệu quốc tế cho thấy các bước tiến lớn về phương pháp luận:

  • Cheinway Hwang và Sung-An Chen (2000) đã tiên phong sử dụng dữ liệu đo cao TOPEX/Poseidon giai đoạn 1993–1999 để phát hiện các hoàn lưu xoáy trung mô quy mô lớn trên Biển Đông, ghi nhận sự phân bố dày đặc của các xoáy nghịch mùa hè và xoáy thuận mùa đông dọc theo các vết đo cao.
  • Jayne (2003) nghiên cứu hoàn lưu Bắc Đại Tây Dương bằng việc kết hợp số liệu đo cao radar với mô hình Geoid GGM02 từ vệ tinh trọng lực GRACE, chứng minh bước tiến nhảy vọt khi phân tách thành phần địa hình động trung bình (Mean Dynamic Topography - MDT).
  • Pail và các cộng sự (2011) đã tạo đột phá khi công bố ước tính dòng chảy địa tĩnh Bắc Đại Tây Dương từ sứ mệnh vệ tinh đo gradient trọng trường GOCE. Nghiên cứu của Pail khẳng định: chỉ với 2 tháng dữ liệu GOCE, trường MDT thu được có độ chính xác vượt trội hơn hẳn so với chuỗi quan sát 8 năm tích lũy từ GRACE, do tín hiệu trọng trường bước sóng ngắn của GOCE ít bị nhiễu hơn, giảm thiểu sai số do làm mịn (smoothing filter), bảo toàn cường độ tín hiệu hải văn.
  • Chao (2010) tổng hợp dữ liệu đo cao đa vệ tinh (T/P, Jason-1/2, ERS-1/2, GEOSAT) tích hợp mô hình trọng lực GRACE và GOCE trên quy mô toàn cầu, khẳng định đo cao vệ tinh xác định chuẩn xác các dòng chảy có bước sóng 100 km và hướng tới độ phân giải 10 km.
  • Martina Idžanović, Vegard Ophaug và Ole Baltazar Andersen (2017) tại vùng ven biển Na Uy đã khai thác dữ liệu CryoSat-2 và GOCE, chỉ ra rằng đo cao vệ tinh đạt độ chính xác MDT từ 3 cm đến 5 cm, song kết quả phụ thuộc nghiêm ngặt vào mô hình Geoid tham chiếu địa phương.
  • Mulero-Martínez và J. Bruno (2018) đã chứng minh tính khả thi khi ứng dụng đo cao radar độ phân giải cao tại Vịnh Cádiz (Tây Ban Nha), đồng thời nhấn mạnh tại dải ven bờ (3–25 km), hiệu ứng ma sát đáy và ứng suất gió ven bờ cần được tính toán bổ sung để bù trừ độ lệch địa tĩnh.

Tại Việt Nam, các nghiên cứu hoàn lưu trước đây chủ yếu tập trung vào mô hình hóa khí tượng - thủy văn hoặc quan trắc cục bộ:

  • Đinh Văn Ưu và cộng sự (2015) sử dụng mô hình 3D ROMS đánh giá tác động biến đổi khí hậu lên trường dòng chảy mặt Biển Đông theo các kịch bản RCP.
  • Đỗ Như Kiều và Lê Đình Mầu (2013) nghiên cứu dòng chảy dọc bờ do sóng đổ nhào tại Tuy An (Phú Yên) bằng mô hình SWAN kết hợp phao trôi và máy đo sóng ALEC AWH-16M.
  • Phạm Xuân Dương (2014) khảo sát hoàn lưu vịnh Thái Lan bằng mô hình ROMS lưới $90 \times 100$ điểm với 5 lớp sigma, chỉ ra dòng chảy ven bờ mùa gió Đông Bắc từ Bình Thuận xuống Tây Nam.
  • Trần Mạnh Cường và Nguyễn Kim Cương (2016) phân tích dòng chảy tầng mặt vịnh Bắc Bộ qua hệ thống trạm Radar sóng ngắn bờ biển, phát hiện dòng ven bờ đạt 25–30 cm/s tại Thanh Hóa – Quảng Bình vào mùa đông.
  • Lê Đình Mầu và cộng sự (2017) phân tích trường dòng chảy ven biển Bắc Trung Bộ (Hà Tĩnh đến Thừa Thiên Huế), ghi nhận tác động chặn dòng của mũi Lăng Cô.

Khoảng trống học thuật và sự tranh luận cốt lõi: Tồn tại một cuộc tranh luận học thuật sâu sắc giữa hai trường phái: Trường phái Thủy động lực thuần túy (dựa vào mô hình số trị Navier-Stokes định hướng dòng chảy bằng trường gió và bức xạ nhiệt) và Trường phái Trắc địa Hải dương học (dựa vào đo đạc hình học bề mặt thực tế từ vệ tinh phản ánh trạng thái cân bằng áp suất). Trường phái thứ nhất dễ bị trôi dạt sai số mô hình khi thiếu dữ liệu viễn thám thực thời, trong khi trường phái thứ hai từng gặp hạn chế do độ phân giải vết quét vệ tinh thưa và sai số vùng nước nông ven bờ.

Luận án của Đỗ Văn Mong định vị chính xác tại điểm giao thoa này: Tận dụng ưu thế đo đạc trực tiếp bề mặt hình học từ các máy đo cao radar băng tần Ka (SARAL/AltiKa) và SAR Altimeter (Sentinel-3A), kết hợp với mô hình mặt biển trung bình toàn cầu tối tân DTU15MSS và mô hình trọng trường EGM2008/GOCE, thiết lập một phương pháp luận trắc địa chuẩn hóa nhằm giải quyết triệt để bài toán khôi phục trường dòng chảy mặt tức thời và trung bình trên toàn bộ vùng biển Việt Nam.

Đó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à làm sâu sắc thêm Lý thuyết Cân bằng Địa tĩnh (Geostrophic Equilibrium) do Laplace và Coriolis đặt nền móng, được hiện đại hóa trong bối cảnh Trắc địa Hải dương học (Physical Geodesy & Ocean Dynamics) theo các mô hình của Heiskanen & Moritz. Trong chuyển động đại dương quy mô lớn, lực gradien áp suất nằm ngang cân bằng với lực Coriolis sinh ra do chuyển động tự quay của Trái đất.

Mô hình lý thuyết được thiết lập qua hệ phương trình vi phân đạo hàm riêng cấp một xác định hai thành phần vector vận tốc dòng chảy địa tĩnh bề mặt:

  • Thành phần vận tốc theo trục Đông – Tây ($u$, dương về hướng Đông): $$u = -\frac{g}{f R} \frac{\partial (MDT)}{\partial \varphi} = -\frac{g}{f} \frac{\partial (MDT)}{\partial y}$$
  • Thành phần vận tốc theo trục Nam – Bắc ($v$, dương về hướng Bắc): $$v = \frac{g}{f R \cos\varphi} \frac{\partial (MDT)}{\partial \lambda} = \frac{g}{f} \frac{\partial (MDT)}{\partial x}$$

Trong đó:

  • $g$ là gia tốc trọng trường cục bộ tại mặt biển ($m/s^2$).
  • $f = 2\Omega \sin\varphi$ là tham số Coriolis ($\Omega = 7.292115 \times 10^{-5}\text{ rad/s}$ là vận tốc góc tự quay của Trái đất; $\varphi$ là vĩ độ địa lý; $\lambda$ là kinh độ địa lý).
  • $R$ là bán kính Trái đất chuẩn ($R \approx 6378137\text{ m}$).
  • $MDT$ (Mean Dynamic Topography) là độ cao địa hình động trung bình, hoặc thay thế bằng tổng hợp $MDT = MSS - N = SLA + MDT_{\text{ref}}$, với $MSS$ là mặt biển trung bình, $N$ là độ cao Geoid, $SLA$ là dị thường mực nước biển (Sea Level Anomaly).

Cường độ vận tốc tổng hợp ($V$) và hướng dòng chảy ($\theta$) được xác định bởi: $$V = \sqrt{u^2 + v^2}, \quad \theta = \arctan\left(\frac{u}{v}\right) \quad (\text{chuẩn hóa theo góc phương vị hải dương học } 0^\circ - 360^\circ)$$

Hệ luận đề lý thuyết (Theoretical Propositions) được luận án kiểm chứng:

  • Proposition 1: Độ dốc bề mặt địa hình động lực học mặt biển tức thời được trích xuất từ dữ liệu viễn thám đo cao radar (sau khi triệt tiêu toàn bộ thành phần triều và hiệu ứng khí quyển nghịch đảo) phản ánh tương đương 90–95% động năng của hoàn lưu bề mặt tại các vùng biển khơi có độ sâu trên 200m.
  • Proposition 2: Sự biến đổi không gian - thời gian của trường dị thường mực nước biển (SLA) biểu diễn chính xác vị trí tâm xoáy, bán kính ảnh hưởng và cường độ dòng xiết của các cấu trúc xoáy đại dương trung mô tại vùng thềm lục địa miền Trung và khu vực quần đảo Trường Sa.

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

Khung phân tích của luận án tích hợp một cách chặt chẽ 3 trụ cột lý thuyết:

  1. Lý thuyết Quỹ đạo và Tín hiệu Đo cao Radar (Radar Altimetry Signal Processing).
  2. Lý thuyết Trường Trọng lực Toàn cầu và Mặt Chuẩn Geoid/MSS (Geodetic Surface Modeling).
  3. Lý thuyết Ước lượng Tối ưu và Xử lý Thống kê Không gian (Least-Squares Collocation & Spatial Covariance Estimation).
                          [Tín hiệu Radar Đo cao Vệ tinh]
                       (SARAL/AltiKa Băng Ka, Sentinel-3A SRAL)
                      [Xử lý Hiệu chỉnh Sai số Vật lý]
                 (Ion, Đối lưu Khô/Ướt, Khí áp, Triều, Sea State Bias)
                     [Tính Độ cao Mặt biển Tức thời: SSH]
                  [Loại bỏ Mặt biển Trung bình: DTU15MSS]
                   [Trích xuất Dị thường Mực nước biển: SLA]
                  [Lọc Dữ liệu Bất thường theo Chuẩn Gauss]
                            (Ngưỡng sai số 3σ)
                  [Bình sai Mạng lưới Giao cắt Crossover]
                 (Xác định sai số hệ thống & Triệt tiêu nhiễu)
                   [Ước lượng Hàm Hiệp phương sai Không gian]
                       (Tính EmpCov & Khớp Covfit lý thuyết)
                     [Xây dựng Lưới Ô vuông SLA Liên tục]
                            (Nội suy Collocation)
                  [Tính Toán Vector Gradient Dòng Địa tĩnh]
                  (Thành phần Vận tốc u, v, Độ lớn V, Hướng θ)
               [Ứng dụng Quân sự & Tác chiến Hải quân Việt Nam]
             (Điều hướng Tàu, Dẫn đường ECDIS, SAR, Ngư lôi/Thủy lôi)

Điều kiện biên và giới hạn lý thuyết (Boundary Conditions) được luận án làm rõ: Lý thuyết dòng chảy địa tĩnh suy biến tại vùng xích đạo ($\varphi \to 0$ dẫn đến $f \to 0$, mẫu số tiến tới 0). Trong phạm vi nghiên cứu của luận án từ $5^\circ\text{N}$ đến $25^\circ\text{N}$, tham số Coriolis $f$ luôn đảm bảo tính ổn định toán học, tránh hiện tượng kỳ dị số học. Tại dải ven bờ cách đất liền dưới 10 km, các yếu tố ma sát đáy biển nông và nhiễu sóng phản xạ địa hình đất liền đòi hỏi việc ứng dụng các thuật toán lọc cửa sổ trượt đặc thù.

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

Thiết kế nghiên cứu

Nghiên cứu đứng vững trên lập trường triết học thực chứng (Positivism) với định lượng toán học - trắc địa nghiêm ngặt, dữ liệu được kiểm nghiệm qua các thuật toán bình sai xác suất và đối chứng thực tế.

Thiết kế nghiên cứu đa tầng bao gồm:

  • Tầng dữ liệu viễn thám: Thu thập từ Trung tâm Dữ liệu Hải dương học AVISO (Pháp) và Cơ quan Vũ trụ Châu Âu (ESA). Sử dụng dữ liệu đo cao từ 2 vệ tinh có công nghệ vượt trội:
    • SARAL/AltiKa: Vệ tinh của ISRO/CNES, sử dụng máy đo cao AltiKa hoạt động ở băng tần Ka (35 GHz), độ cao quỹ đạo ~800 km, góc nghiêng $98.55^\circ$, chu kỳ lặp 35 ngày. Tần số sóng Ka mang lại kích thước vết quét nhỏ hơn, cải thiện đáng kể độ chính xác tại vùng ven bờ và giảm thiểu hiện tượng tán xạ do mưa.
    • Sentinel-3A: Vệ tinh thuộc chương trình Copernicus (ESA), trang bị máy đo cao SAR Altimeter (SRAL) hoạt động ở chế độ khẩu độ tổng hợp SAR (độ phân giải dọc vết đo cực cao), máy đo bức xạ vi sóng MWR, hệ thống định vị quỹ đạo chính xác POD tích hợp DORIS, GNSS đa kênh và phản xạ laser LRR. Chu kỳ lặp 27 ngày, góc nghiêng $98.65^\circ$, độ cao 814.5 km.
  • Tầng dữ liệu trọng lực và mặt biển chuẩn: Thu thập từ ICGEM (International Centre for Global Earth Models) bao gồm các mô hình trọng trường Trái đất EGM2008, XGM2019, mô hình Geoid vệ tinh GOCE/GRACE, và mô hình mặt biển trung bình toàn cầu DTU15MSS của Đại học Kỹ thuật Đan Mạch với độ phân giải $1' \times 1'$.

Quy trình nghiên cứu rigorous

Quy trình công nghệ gồm 6 bước chuẩn hóa:

Bước 1: Hiệu chỉnh sai số đường truyền và tính độ cao mặt biển tức thời ($SSH$). Độ cao từ vệ tinh đến mặt biển ($h$) được đo qua thời gian lan truyền tín hiệu hai chiều $\Delta t$ theo công thức cơ bản: $h = c \cdot \frac{\Delta t}{2}$ (với $c$ là tốc độ ánh sáng). Sau khi xác định độ cao quỹ đạo vệ tinh ($H$) so với Ellipsoid tham chiếu WGS-84 thông qua hệ thống POD/DORIS, độ cao mặt biển tức thời được tính: $$SSH = H - h - h_{cor}$$ Trong đó, số hiệu chỉnh tổng hợp $h_{cor}$ bao gồm các thành phần vật lý nghiêm ngặt:

  • Hiệu chỉnh tầng đối lưu khô (Dry Troposphere): Đạt giá trị lớn nhất (khoảng 2.3 m), ít biến thiên thời gian, tính theo mô hình giải tích Saastamoinen kết hợp dữ liệu áp suất bề mặt từ ECMWF/NCEP.
  • Hiệu chỉnh tầng đối lưu ướt (Wet Troposphere): Biến thiên từ 0 đến 50 cm, trích xuất từ dữ liệu đo bức xạ vi sóng MWR đồng thời trên vệ tinh.
  • Hiệu chỉnh tầng điện ly (Ionosphere): Từ 0 đến 50 cm, tính toán dựa trên trị đo đa tần số (Ku-band và C-band hoặc đo Doppler đa tần DORIS).
  • Hiệu chỉnh nghịch đảo khí áp (Dynamic Atmospheric Correction): Bù trừ phản ứng mực nước biển do chênh lệch áp suất khí quyển ($\approx 1\text{ cm}$ trên mỗi $1\text{ mbar}$).
  • Hiệu chỉnh độ gồ ghề mặt biển (Sea State Bias - SSB): Bù trừ hiện tượng chân sóng phản xạ tín hiệu radar tốt hơn đỉnh sóng.
  • Hiệu chỉnh thủy triều (Tidal Correction): Loại bỏ triều đại dương thông qua mô hình FES99/FES2014 và triều vỏ đất rắn (Solid Earth Tide) với độ chính xác dưới 1 cm.

Bước 2: Phân tích độ chính xác tại điểm giao cắt (Crossover Analysis). Trong văn bản luận án, cơ sở toán học đánh giá chất lượng số liệu đo cao được công thức hóa chi tiết qua độ chênh lệch độ cao $dH$ tại điểm giao cắt giữa cung bay thăng (ascending track) và cung bay giáng (descending track): $$dH = SSH_C^t - SSH_C^g$$ Tại giao điểm $C$, kỳ vọng toán học lý thuyết $E(dH) = 0$. Luận án thiết lập hai trường hợp đánh giá sai số:

  • Khi dãy trị đo không có sai số hệ thống ($E(dH) = 0$), sai số trung phương của trị đo tại điểm giao cắt ($m_{dH}$) và sai số một lần đo ($m_{\text{đo}}$) tính theo công thức Gauss: $$m_{dH} = \pm \sqrt{\frac{[dH^2]}{m}} \implies m_{\text{đo}} = \frac{m_{dH}}{\sqrt{2}} = \pm \sqrt{\frac{[dH^2]}{2m}}$$
  • Khi tồn tại sai số hệ thống ($E(dH) \neq 0$), áp dụng công thức Bessel với số hiệu chỉnh $v = dH - \overline{dH}$: $$m_{dH} = \pm \sqrt{\frac{[vv]}{m-1}} \implies m_{\text{đo}} = \pm \sqrt{\frac{[vv]}{2(m-1)}}$$

Bước 3: Tách chiết Dị thường Mực nước biển (SLA) và Lọc thô thống kê. Dị thường mực nước biển được tách chiết bằng cách loại bỏ mô hình mặt biển trung bình: $$SLA = SSH - DTU15MSS$$ Dữ liệu SLA dọc theo vết quét được lọc sạch các giá trị bất thường (outliers) bằng tiêu chuẩn xác suất thống kê $3\sigma$. Các giá trị vượt ngưỡng $|SLA_i - \overline{SLA}| > 3\sigma$ bị loại bỏ hoàn toàn để bảo toàn tính chuẩn xác của trường dữ liệu.

Bước 4: Ước lượng hàm hiệp phương sai không gian (Covariance Modeling). Tính toán hiệp phương sai thực nghiệm (Empirical Covariance - EmpCov) giữa các cặp điểm đo cao cách nhau khoảng cách khoảng cách $d_{ij}$: $$C(d_k) = \frac{1}{N_k} \sum_{i,j} (SLA_i \cdot SLA_j)$$ Sử dụng công cụ Covfit để làm khớp các giá trị thực nghiệm với hàm hiệp phương sai lý thuyết dạng Gauss hoặc hàm mũ Markov cấp hai, xác định chính xác bán kính tương quan không gian ($D_0$) và phương sai tín hiệu ($C_0$).

Bước 5: Xây dựng lưới ô vuông SLA và MDT bằng thuật toán nội suy Least-Squares Collocation. Nội suy trường độ cao mặt biển liên tục trên lưới chính quy với kích thước ô lưới phù hợp ($0.25^\circ \times 0.25^\circ$ hoặc $0.1^\circ \times 0.1^\circ$).

Bước 6: Xác định trường vector dòng chảy địa tĩnh và nội suy điểm bất kỳ. Tính đạo hàm riêng số trị của bề mặt độ cao theo hai phương kinh - vĩ độ để giải hệ phương trình Laplace - Coriolis, trích xuất hai thành phần $u, v$, tính độ lớn $V$ và hướng dòng $\theta$. Áp dụng thuật toán nội suy song tuyến tính (Bilinear Interpolation) để xác định vector dòng chảy tại tọa độ hải trình bất kỳ.

Data và phân tích

Hệ thống số liệu xử lý trong luận án có khối lượng cực lớn:

  • Đánh giá chất lượng 34 chu kỳ vệ tinh SARAL/AltiKa và 37 chu kỳ vệ tinh Sentinel-3A trên Biển Đông.
  • Độ lệch chuẩn sai số đo cao tại điểm giao cắt của Sentinel-3A đạt mức $m_{\text{đo}} = \pm 2.8\text{ cm} - \pm 4.2\text{ cm}$, của SARAL/AltiKa đạt $\pm 3.1\text{ cm} - \pm 4.8\text{ cm}$, hoàn toàn đáp ứng tiêu chuẩn quốc tế về độ chính xác nghiên cứu dòng chảy bề mặt đại dương.
  • Bộ công cụ và phần mềm xử lý chuyên dụng: Basic Radar Altimetry Toolbox (BRAT) của ESA/CNES, phần mềm Covfit, hệ thống thuật toán tự phát triển trên môi trường Matlab/Python, kết hợp phần mềm chuyên dụng xử lý hải đồ của Quân chủng Hải quân.
  • Dữ liệu đối chiếu thực địa: Chuỗi số liệu từ hệ thống phao đo vệt trôi (drifter buoys) thả trên vùng biển miền Trung và Trường Sa trong năm 2022 vào các ngày: 30/5/2022, 31/5/2022, 03-04/6/2022, 11/8/2022, 19/8/2022, 22/8/2022, 02-05/10/2022, 02/11/2022. Cùng với đó là dữ liệu hải trình thực tế của tàu khảo sát trên mặt cắt ngang vĩ tuyến $11^\circ 40' 00''\text{N}$, kinh độ từ $109^\circ 15' 00''\text{E}$ đến $112^\circ 00' 00''\text{E}$.

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

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

  1. Khẳng định độ chính xác vượt trội của công nghệ đo cao băng tần Ka và SAR Altimetry trên Biển Đông: Kết quả bình sai giao cắt chứng minh sai số trung phương của dữ liệu Sentinel-3A (vết quét SAR) đạt độ ổn định cao với $m_{\text{đo}} < \pm 3.5\text{ cm}$, triệt tiêu hoàn toàn hiện tượng lệch dải tín hiệu giữa cung thăng và cung giáng sau bình sai mạng lưới.
  2. Khôi phục thành công bức tranh hoàn lưu động lực học theo mùa tại vùng biển Việt Nam: Bản đồ trường dòng chảy tỷ lệ 1/2.000.000 thành lập từ số liệu đo cao Sentinel-3A (các chu kỳ CK71, CK72, CK75, CK77, CK81, CK84) phản ánh rõ nét tính chu kỳ gió mùa:
    • Trong thời kỳ gió mùa mùa hè (tháng 6 – tháng 8, điển hình chu kỳ CK72, CK75): Xuất hiện dòng tia mạnh ven bờ Nam Trung Bộ chảy hướng Đông Bắc với vận tốc lớn đạt từ 0.6 m/s đến trên 1.1 m/s, tách bờ tại khoảng vĩ độ $12^\circ - 13^\circ\text{N}$ tạo thành cấu trúc lưỡi dòng hướng ra giữa Biển Đông, kèm theo các tâm xoáy nghịch (anticyclonic eddies) quy mô lớn ở phía Đông bán đảo Cam Ranh.
    • Trong thời kỳ chuyển tiếp và mùa đông (chu kỳ CK77, CK81): Trường dòng chuyển hướng chủ đạo sang Tây Nam dọc thềm lục địa miền Trung với vận tốc ổn định 0.4 – 0.7 m/s.
  3. Độ tương đồng cao giữa số liệu vệ tinh và phao đo vệt trôi: Kết quả so sánh đối ứng giữa vector dòng chảy tính từ đo cao vệ tinh với quỹ đạo thực tế của phao trôi hải dương trong các đợt thả ngày 30/5, 31/5, 03-04/6, 11/8, 19/8, 22/8, 02/10, 05/10 và 02/11/2022 cho thấy độ lệch hướng dòng chảy trung bình dưới $\pm 12^\circ - \pm 18^\circ$, độ lệch vận tốc tuyệt đối dao động trong khoảng $0.05 - 0.12\text{ m/s}$. Đây là minh chứng thực nghiệm vững chắc xác nhận dòng chảy địa tĩnh chiếm ưu thế tuyệt đối trong động lực học tầng mặt khu vực nghiên cứu.
  4. Hiệu ứng dòng chảy cục bộ và tác động dạt tàu thực tế: Thực nghiệm hành trình thực tế trên mặt cắt vĩ độ $11^\circ 40' 00''\text{N}$ (từ $109^\circ 15' 00''\text{E}$ đến $112^\circ 00' 00''\text{E}$) chia thành 3 phân đoạn tuyến A-B, B-C, C-D đã định lượng chính xác góc dạt tàu do dòng chảy ngang mặt cắt gây ra, làm trôi dạt thân tàu từ $1.5^\circ$ đến $4.8^\circ$ so với hướng lái la bàn khi không có hiệu chỉnh.

Implications đa chiều

  • Về mặt Lý thuyết: Công trình đã hoàn thiện cơ sở lý luận về trắc địa hải dương học tại Việt Nam, thiết lập mối liên kết toán học hữu cơ giữa thế trọng trường Trái đất, mô hình Geoid/MDT và động lực học dòng chảy địa tĩnh, mở ra hướng nghiên cứu mới trong việc sử dụng dữ liệu viễn thám không gian để giải các bài toán thủy văn biển sâu.
  • Về mặt Phương pháp luận: Cung cấp một quy trình công nghệ khép kín từ tiền xử lý dữ liệu thô, lọc bỏ nhiễu, khớp hàm tương quan không gian đến tự động hóa trích xuất bản đồ dòng chảy diện rộng với chi phí vận hành tiệm cận bằng 0 so với phương pháp khảo sát tàu biển truyền thống.
  • Về mặt Thực tiễn và Quốc phòng - An ninh Hải quân:
    • Tối ưu hóa hành trình tàu chiến đấu: Luận án khẳng định nắm bắt quy luật dòng chảy cho phép sĩ quan điều hải tính toán tốc độ tối ưu và góc bẻ lái bù dạt, tiết kiệm từ 8% đến 15% nhiên liệu hành trình, đồng thời rút ngắn thời gian cơ động tác chiến trên biển.
    • Dẫn đường tích hợp: Toàn bộ dữ liệu trường dòng chảy dạng số được tích hợp trực tiếp vào hệ thống thông tin hiển thị hải đồ điện tử (ECDIS) kết hợp hệ thống định vị DGPS trên các tàu Hải quân hiện đại, cung cấp thông tin trợ giúp quyết định thời gian thực cho thuyền trưởng.
    • Tác chiến vũ khí dưới nước (Ngư lôi và Thủy lôi): Dòng chảy tác động trực tiếp làm biến dạng quỹ đạo dẫn hướng của ngư lôi tự dẫn và khí tài ngầm. Nắm chắc trường vector vận tốc cho phép hiệu chỉnh tham số bắn chính xác. Đặc biệt trong tác chiến phong tỏa thủy lôi, văn bản nêu rõ: "người chỉ huy sẽ thả những quả thủy lôi và lợi dụng hướng của dòng chảy để đưa những quả thủy lôi vào đội hình tàu địch hoặc các tuyến luồng nơi mà tàu địch thường xuyên phải qua lại".
    • Tìm kiếm cứu nạn (SAR) và Ứng phó sự cố tràn dầu: Cung cấp trường vector dòng chảy làm đầu vào mô hình hóa bài toán Monte Carlo dự báo vết trôi dạt của nạn nhân, phao cứu sinh, hoặc vệt dầu loang, thu hẹp diện tích tìm kiếm cứu nạn khẩn cấp từ hàng ngàn km² xuống khu vực trọng tâm.
    • Khảo sát trắc địa công trình biển: Tối ưu hóa hướng tuyến quét của thiết bị Sonar quét sườn (Side-scan Sonar) và thả thiết bị đo sâu đa tia đáy biển dọc theo hướng dòng chảy để loại bỏ rung lắc và nhiễu sóng thân tàu.

Limitations và Future Research

Luận án thẳng thắn chỉ ra các giới hạn kỹ thuật nội tại:

  • Giới hạn vùng dải ven bờ (Coastal Blind Zone): Tại khoảng cách dưới 5–10 km cách bờ, tín hiệu radar đo cao bị ô nhiễm nặng do phản xạ hỗn hợp từ địa hình đất liền và công trình nhân tạo, dẫn đến sai số độ cao mặt biển gia tăng, đòi hỏi phải áp dụng các kỹ thuật Re-tracking chuyên biệt.
  • Độ phân giải thời gian của vệ tinh đơn lẻ: Chu kỳ lặp 27 ngày của Sentinel-3A hay 35 ngày của SARAL/AltiKa tạo khoảng trống thời gian không quan sát được các biến động thủy triều tức thời hoặc các hiện tượng nước trồi ngắn hạn dưới 10 ngày.
  • Thành phần dòng phi địa tĩnh: Phương pháp đo cao chỉ trích xuất thành phần dòng cân bằng địa tĩnh (Geostrophic), chưa tích hợp đầy đủ thành phần dòng trôi Ekman do gió bão tức thời tạo ra tại tầng nước sát mặt biển dày 0–5 m.

Chương trình nghiên cứu tiếp nối trong 5–10 năm tới:

  1. Đồng hóa dữ liệu đo cao đa vệ tinh (Sentinel-3A/3B, Sentinel-6 Michael Freilich, và vệ tinh giao thoa đo cao góc rộng SWOT) vào các mô hình số trị thủy động lực 3D (ROMS, FVCOM) nhằm tái tạo trường dòng chảy 3 chiều theo thời gian thực với độ phân giải siêu cao 1–2 km.
  2. Ứng dụng trí tuệ nhân tạo (Deep Learning / Physics-Informed Neural Networks - PINNs) để tái cấu trúc tín hiệu đo cao radar ven bờ bị nhiễu và khôi phục trường dòng chảy tầng sâu hỗ trợ hoạt động của tàu ngầm Kilo 636.
  3. Mở rộng nghiên cứu tương tác giữa dòng chảy đáy đại dương và địa hình ngầm phục vụ xây dựng mạng lưới định vị thủy âm ngầm dưới biển.

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

  • Ảnh hưởng học thuật: Thiết lập quy chuẩn tính toán trắc địa biển mới cho các cơ sở nghiên cứu và đào tạo sau đại học tại Việt Nam, đóng góp chuỗi bài báo chất lượng cao trên các tạp chí chuyên ngành trắc địa - địa vật lý uy tín.
  • Hiện đại hóa Quân chủng Hải quân: Chuyển đổi toàn diện phương thức bảo đảm hải văn từ thủ công, kinh nghiệm sang số hóa tự động tích hợp hải đồ điện tử ECDIS, nâng cao rõ rệt hiệu suất sẵn sàng chiến đấu, tuần tra bảo vệ chủ quyền tại hai quần đảo Hoàng Sa và Trường Sa.
  • Giá trị kinh tế - xã hội: Cung cấp dữ liệu nền tảng cho Cục Hàng hải Việt Nam, các đội tàu vận tải biển viễn dương, ngành đánh bắt hải sản xa bờ, quy hoạch điện gió ngoài khơi và các trung tâm cứu nạn hàng hải quốc gia (MRCC).
  • Khẳng định chủ quyền quốc gia: Các ấn phẩm khoa học và công trình nghiên cứu sử dụng dữ liệu đo đạc chi tiết trên Biển Đông là cơ sở pháp lý và khoa học vững chắc khẳng định chủ quyền, quyền chủ quyền và quyền tài phán của Việt Nam trên trường quốc tế.

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

  • Nghiên cứu sinh và Giảng viên ngành Trắc địa - Hải dương học: Tiếp cận quy trình tính toán Least-Squares Collocation, xử lý số liệu đo cao radar và mô hình Geoid hiện đại.
  • Bộ Tư lệnh Quân chủng Hải quân & Đoàn Đo đạc Biên vẽ Hải đồ (Đoàn 6): Sở hữu giải pháp công nghệ trực tiếp phục vụ công tác dẫn đường hàng hải, tác chiến tàu ngầm, bố trí bãi thủy lôi, phóng ngư lôi và tìm kiếm cứu nạn trên biển.
  • Ngành Hàng hải & Khai thác Dầu khí Biển: Ứng dụng trường dòng chảy tối ưu hóa lộ trình tàu chở dầu, tàu hàng container, giảm phát thải carbon và bảo đảm an toàn cho các giàn khoan ngoài khơi.
  • Cơ quan Hoạch định Chính sách & Khí tượng Thủy văn: Nâng cao năng lực cảnh báo sớm sự cố môi trường biển, nước trồi sinh thái và giám sát tác động của biến đổi khí hậu lên hải lưu khu vực.

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ì? Đó là việc mở rộng thành công Lý thuyết Cân bằng Địa tĩnh (Geostrophic Equilibrium) trong hệ quy chiếu trắc địa vật lý cục bộ Biển Đông, thiết lập giải pháp toán học kết hợp mô hình mặt biển trung bình DTU15MSS với thuật toán bình sai giao cắt để triệt tiêu sai số quỹ đạo vệ tinh mà không làm suy giảm biên độ tín hiệu hải văn trung mô.

2. Điểm mới về phương pháp luận so với các nghiên cứu tiền nhiệm quốc tế là gì? So với nghiên cứu của Pail et al. (2011) ở Bắc Đại Tây Dương hay Martina Idžanović et al. (2017) tại Na Uy, luận án đã xây dựng quy trình kết hợp đa nguồn linh hoạt giữa dữ liệu SAR Altimetry băng tần Ku (Sentinel-3A) và sóng siêu ngắn băng tần Ka (SARAL/AltiKa), áp dụng hàm hiệp phương sai lý thuyết tự thích nghi theo vĩ độ nhiệt đới Biển Đông, giải quyết triệt để bài toán thiếu hụt trạm đo kiểm chuẩn trọng lực ngoài khơi xa.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì? Sự tương đồng chặt chẽ đến kinh ngạc giữa vector dòng chảy tính từ độ dốc mặt biển của vệ tinh Sentinel-3A với quỹ đạo thực tế của phao trôi hải dương học (độ lệch hướng chỉ $\pm 12^\circ - \pm 18^\circ$, độ lệch vận tốc dưới $0.12\text{ m/s}$) ngay cả tại vùng chuyển tiếp thủy văn phức tạp ngoài khơi Ninh Thuận – Bình Thuận trong mùa hè 2022, chứng minh trường dòng địa tĩnh chi phối áp đảo so với các thành phần dòng phi địa tĩnh tại thềm lục địa sâu.

4. Nghiên cứu có cung cấp giao thức tái lặp (Replication Protocol) hoàn chỉnh không? Có. Luận án mô tả chi tiết từng bước thuật toán: từ mã hóa dữ liệu nhị phân của AVISO/ESA, công thức tính các số hiệu chỉnh vật lý Saastamoinen/MWR/DORIS, thuật toán bình sai giao cắt Gauss-Bessel, phương pháp tính hàm hiệp phương sai thực nghiệm EmpCov đến thuật toán nội suy vi phân giải hệ phương trình Laplace - Coriolis.

5. Lộ trình nghiên cứu 10 năm tới được phác thảo như thế nào? Tập trung vào 3 trọng tâm: (1) Tích hợp dữ liệu siêu phân tích 2D từ vệ tinh thế hệ mới SWOT; (2) Xây dựng hệ thống đồng hóa dữ liệu viễn thám thời gian thực vào mô hình số trị đại dương đa tầng; (3) Hoàn thiện phần mềm chuyên dụng phục vụ tác chiến ngầm dưới đáy biển sâu cho lực lượng tàu ngầm Hải quân Việt Nam.

Kết luận

  1. Luận án đã xây dựng thành công phương pháp luận và quy trình công nghệ chuẩn hóa, có cơ sở khoa học chặt chẽ và tính khả thi tuyệt đối trong việc trích xuất trường dòng chảy bề mặt từ số liệu đo cao vệ tinh trên toàn bộ vùng biển Việt Nam.
  2. Ứng dụng thành công công nghệ đo cao tiên tiến từ vệ tinh SARAL/AltiKa và Sentinel-3A, chứng minh độ chính xác đo cao tại điểm giao cắt đạt mức sai số trung phương $m_{\text{đo}} < \pm 3.5\text{ cm}$, đáp ứng hoàn hảo yêu cầu tính toán hải văn độ chính xác cao.
  3. Thành lập hệ thống bản đồ dòng chảy bề mặt Biển Đông theo các chu kỳ mùa với tỷ lệ 1/2.000.000, phản ánh trung thực quy luật hoàn lưu và các cấu trúc xoáy đại dương trung mô, được kiểm chứng chuẩn xác qua chuỗi dữ liệu thực địa từ phao đo vệt trôi và tàu khảo sát của Quân chủng Hải quân.
  4. Đề xuất các giải pháp kỹ thuật mang tính đột phá ứng dụng trực tiếp vào công tác dẫn đường hàng hải trên hệ thống ECDIS, tối ưu hóa cơ động tàu chiến đấu, nâng cao độ chính xác triển khai vũ khí ngầm (thủy lôi, ngư lôi), phục vụ tìm kiếm cứu nạn và khảo sát địa hình đáy biển.
  5. Mở ra một kỷ nguyên mới cho chuyên ngành Kỹ thuật Trắc địa - Bản đồ biển tại Việt Nam, khẳng định chủ quyền biển đảo thiêng liêng của Tổ quốc thông qua các thành tựu nghiên cứu khoa học công nghệ vũ trụ đỉnh cao.