Tổng quan nghiên cứu

Vịnh Bắc Bộ là vùng biển chiến lược có diện tích khoảng 160.000 km², chu vi xấp xỉ 1.950 km với chiều dài bờ biển phía Việt Nam trải dài hơn 740 km. Đây là khu vực tập trung các tuyến hàng hải trọng điểm, ngư trường lớn và giữ vai trò tối quan trọng đối với an ninh quốc phòng. Tuy nhiên, việc giám sát động lực học tầng mặt tại vùng biển này trước đây gặp nhiều thách thức do thiếu hụt công nghệ quan trắc liên tục với độ phân giải cao theo thời gian thực. Đề tài luận văn thạc sĩ chuyên ngành Hải dương học của tác giả Nguyễn Thanh Trang tập trung giải quyết bài toán cốt lõi: khai thác, xử lý và ứng dụng chuỗi số liệu radar biển tần số cao (HF) để phân tích chi tiết chế độ dòng chảy tầng mặt và mô phỏng quỹ đạo vật thể trôi tại Vịnh Bắc Bộ.

Nghiên cứu được triển khai dựa trên hệ thống radar Seasonde thuộc Giai đoạn 1 của Đề án xây dựng hệ thống radar biển Việt Nam do Trung tâm Hải văn (Tổng cục Biển và Hải đảo Việt Nam) quản lý. Phạm vi không gian của nghiên cứu bao phủ vùng biển giới hạn từ 17,7°N đến 21,13°N và từ 105,0°E đến 108,56°E, tiếp nhận dữ liệu từ 3 trạm quan trắc ven bờ gồm Hòn Dáu, Nghi Xuân và Đồng Hới với bán kính quét tối đa lên tới 300 km. Kết quả nghiên cứu mang lại ý nghĩa học thuật và thực tiễn sâu sắc: không chỉ làm sáng tỏ các cấu trúc hoàn lưu quy mô vừa và nhỏ biến động theo mùa với độ chính xác sai số vận tốc dưới 7 cm/s, mà còn xây dựng cơ sở khoa học tin cậy để nâng cao hiệu suất tìm kiếm cứu nạn và ứng phó sự cố tràn dầu trên biển lên trên 90%.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Cơ sở lý thuyết của luận văn dựa trên sự kết hợp liên ngành giữa vật lý sóng vô tuyến, động lực học biển và mô hình toán chuyển dịch chất điểm:

  • Lý thuyết tán xạ sóng Bragg và hiệu ứng Doppler: Hệ thống radar HF Seasonde vận hành ở dải tần số cao từ 4 đến 50 MHz (bước sóng điện từ từ 6 đến 75 m). Tín hiệu phát xạ tương tác với sóng biển bề mặt tạo hiện tượng tán xạ cộng hưởng Bragg. Vận tốc pha của sóng biển và dòng chảy tầng mặt gây ra sự dịch chuyển tần số Doppler (Doppler shift), cho phép trích xuất vận tốc dòng chảy hướng tâm với độ chính xác cao.
  • Lý thuyết tổng hợp trường vector dòng chảy hai chiều: Vận tốc dòng chảy thực tế tại mỗi nút lưới được hợp thành từ ít nhất hai thành phần vector hướng tâm từ các trạm radar khác nhau, sử dụng phương pháp bình phương tối thiểu (Least Squares Method) và thuật toán nội suy ranh giới hình học.
  • Mô hình quỹ đạo hạt Lagrangian hai chiều (Lagrangian Particle Tracking Model): Quá trình chuyển dịch của vật thể hoặc vệt dầu được giải theo phương trình vi phân chuyển động hạt Euler-Lagrange, kết hợp giữa thành phần bình lưu do dòng chảy bề mặt và tác động của trường khuếch tán rối với hệ số khuếch tán chuẩn hóa $D = 56\text{ m}^2/\text{s}$.

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

Nguồn dữ liệu thực nghiệm được thu thập liên tục từ ngày 05/06/2012 đến ngày 15/11/2014 với bước thời gian 1 giờ/lần tại 3 trạm radar ven bờ. Tổng quy mô mẫu quan trắc ghi nhận gồm 11.434 lượt đo tại trạm Hòn Dáu (đạt mật độ bao phủ 55 - 60%), 15.163 lượt đo tại trạm Nghi Xuân (đạt mật độ bao phủ 70 - 75%) và khoảng 12.000 lượt đo tại trạm Đồng Hới (đạt mật độ 55 - 60%). Phương pháp chọn mẫu là quan trắc tự động định kỳ đồng bộ toàn mạng lưới, giúp loại bỏ tính chủ quan và ghi nhận trọn vẹn các dao động thủy triều cũng như chu kỳ thời tiết.

Phần mềm chuyên dụng CODAR SeaSonde Combine Suite Release 7 được áp dụng để xử lý dữ liệu và thiết lập miền lưới tính toán gồm 55 điểm theo phương vĩ tuyến và 63 điểm theo phương kinh tuyến với bước lưới không gian 6 km, tâm lưới tại vị trí 19°24,6'N và 107°00'E. Phương pháp bình phương tối thiểu được lựa chọn vì tính ưu việt trong việc khử nhiễu ngẫu nhiên và giải quyết bài toán giao thoa đa nguồn. Đối với các vùng ranh giới có góc cắt giữa hai tia radar nhỏ hơn 30°, mô hình áp dụng thuật toán nội suy tuyến tính đối với thành phần trực giao và bình phương tối thiểu cho thành phần song song để đảm bảo tính liên tục của trường vận tốc. Cuối cùng, mô hình CODAR Leeway được ứng dụng với hệ số dòng chảy $C_c = 0,787$ để mô phỏng quỹ đạo chuyển động của vật thể trong 30 ngày (kịch bản trung bình mùa) và 7 ngày (kịch bản thời gian thực).

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Phân tích chuỗi số liệu trường dòng chảy trung bình tháng tại Vịnh Bắc Bộ đã mang lại những phát hiện quan trọng về cấu trúc thủy động lực tầng mặt:

  • Chế độ hoàn lưu mùa đông (tháng 12 đến tháng 2): Toàn vịnh hình thành một xoáy thuận khép kín ở vùng trung tâm với vận tốc trung bình từ 10 đến 15 cm/s. Dòng chảy ven bờ phía tây (từ nam Thanh Hóa đến bắc Quảng Bình) bị nén mạnh và chảy xuôi về phía nam với vận tốc cao đạt từ 25 đến 30 cm/s. Đáng chú ý, ở phía bắc vịnh xuất hiện một xoáy nghịch quy mô nhỏ có đường kính từ 50 đến 70 km, vận tốc duy trì từ 5 đến 8 cm/s.
  • Chế độ hoàn lưu mùa xuân (tháng 3 đến tháng 5): Cấu trúc xoáy thuận trung tâm vẫn duy trì với vận tốc từ 8 đến 12 cm/s. Vận tốc dải ven bờ hạ nhiệt nhẹ xuống mức 18 đến 25 cm/s. Xoáy nghịch phía bắc mở rộng xuống phía nam và suy yếu dần rồi biến mất hoàn toàn vào tháng 5.
  • Chế độ hoàn lưu mùa hè (tháng 6 đến tháng 8): Bức tranh dòng chảy đảo chiều sâu sắc. Khu vực trung tâm vịnh chuyển đổi thành một hoàn lưu xoáy nghịch với vận tốc từ 8 đến 12 cm/s do khối nước phía nam đẩy lên qua cửa vịnh. Ngược lại, phía bắc vịnh xuất hiện một xoáy thuận nhỏ.
  • Chế độ hoàn lưu mùa thu (tháng 9 đến tháng 11): Hoàn lưu xoáy thuận trung tâm tái lập mạnh mẽ với vận tốc từ 15 đến 20 cm/s. Dòng chảy ven bờ Thanh Hóa - Quảng Bình đạt cường độ cực đại trong năm từ 30 đến 40 cm/s (nhiều điểm đạt 45 cm/s, tăng hơn 50% so với mùa xuân). Đặc biệt, nghiên cứu phát hiện cấu trúc xoáy đôi (dipole) tồn tại ổn định ở phía bắc vịnh gồm một xoáy thuận trên đỉnh vịnh và một xoáy nghịch ngay bên dưới với quy mô từ 70 đến 90 km, vận tốc từ 8 đến 10 cm/s.

Thảo luận kết quả

Các kết quả phân tích trường dòng chảy tổng hợp được trực quan hóa thông qua hệ thống bản đồ vector hai chiều và biểu đồ phân bố vận tốc trung bình tháng. Những cấu trúc này khẳng định sự chi phối mạnh mẽ của trường gió mùa (gió mùa Đông Bắc trong mùa đông - thu và gió mùa Tây Nam trong mùa hè) kết hợp với hiệu ứng đường bờ cong khúc khuỷu và độ sâu nông dưới 100 m của Vịnh Bắc Bộ.

Khi so sánh với các khảo sát lịch sử trước đây như Báo cáo Điều tra Tổng hợp Vịnh Bắc Bộ (1964) hay nghiên cứu của Ding và cộng sự (1960 - 1962), kết quả từ số liệu radar biển thể hiện sự đồng nhất cao về các hoàn lưu quy mô lớn. Tuy nhiên, nhờ độ phân giải không gian mịn (6 km) và tần suất đo liên tục (1 giờ/lần), dữ liệu radar biển đã bộc lộ rõ các xoáy quy mô trung bình (50 - 90 km) ở phía bắc vịnh – cấu trúc mà các nghiên cứu trước đây dùng tàu khảo sát hoặc trạm phao rời rạc không thể nắm bắt được. Thử nghiệm mô phỏng quỹ đạo vật thể trôi (kích thước mẫu 5m x 20m) bằng mô hình Lagrangian trong 30 ngày cũng cho thấy độ tương thích vượt trội: trong mùa đông, vật thể ở Nghệ An trôi dạt nhanh dọc theo bờ biển hướng đông nam về phía Đồng Hới, trong khi vật thể tại Hải Phòng di chuyển chậm và xoay tròn theo xoáy nghịch phía bắc.

Đề xuất và khuyến nghị

Nhằm tối ưu hóa giá trị ứng dụng của hệ thống radar biển trong nghiên cứu hải dương học và quản lý nhà nước, nghiên cứu đưa ra 4 nhóm khuyến nghị trọng tâm:

  • Tích hợp trường gió và sóng dự báo vào mô hình trôi dạt: Cần nâng cấp mô hình CODAR Leeway bằng cách ghép nối dữ liệu gió bề mặt 10m từ các mô hình khí tượng số trị như WRF hoặc ECMWF để tính toán đầy đủ hệ số cản của gió ($C_d$ từ 2,5% đến 4,4%), giúp giảm thiểu sai số dự báo vị trí vật thể trôi xuống dưới 5% trong vòng 24 - 48 giờ tới.
  • Mở rộng mật độ trạm radar giai đoạn 2: Kiến nghị Tổng cục Biển và Hải đảo Việt Nam đẩy nhanh tiến độ lắp đặt thêm 2 trạm radar tại phía bắc Vịnh Bắc Bộ và 3 trạm tại vùng biển Nam Bộ trong kế hoạch hoàn thiện mạng lưới 17 trạm toàn quốc, nâng tỷ lệ bao phủ dữ liệu đồng thời lên trên 85%.
  • Đồng hóa số liệu radar vào các mô hình thủy động lực học 3D: Khuyến nghị các viện nghiên cứu chuyên ngành sử dụng trường dòng chảy và sóng radar làm điều kiện biên ban đầu (Data Assimilation) cho các mô hình số trị ROMS, POM hay Delft3D nhằm tăng độ chính xác của các bản tin dự báo biển hạn ngắn lên 20 - 30%.
  • Xây dựng quy trình phối hợp tìm kiếm cứu nạn thời gian thực: Cục Hàng hải Việt Nam và Ủy ban Quốc gia Tìm kiếm Cứu nạn cần chuẩn hóa việc tiếp nhận trực tiếp dữ liệu dòng chảy từ Trung tâm Hải văn để lập phương án cứu hộ cho tàu cá và thuyền viên bị nạn ngay trong 6 giờ đầu tiên sau khi nhận tín hiệu cấp cứu.

Đối tượng nên tham khảo luận văn

Luận văn là tài liệu tham khảo học thuật và cẩm nang kỹ thuật giá trị cho 4 nhóm đối tượng cụ thể:

  • Cơ quan quản lý nhà nước và lực lượng Tìm kiếm Cứu nạn (SAR): Giúp Trung tâm Phối hợp tìm kiếm cứu nạn hàng hải Việt Nam (VMRCC) và Cảnh sát Biển định vị nhanh chóng tọa độ trôi dạt của tàu thuyền tai nạn, tối ưu hóa vùng tìm kiếm và tiết kiệm chi phí nhiên liệu cứu nạn.
  • Các viện nghiên cứu và chuyên gia Khí tượng Thủy văn - Hải dương học: Cung cấp nguồn tài liệu nền tảng về phương pháp xử lý dữ liệu radar HF Seasonde, thuật toán tổng hợp trường vector 2D và cơ chế hoàn lưu phức tạp theo mùa tại Vịnh Bắc Bộ.
  • Doanh nghiệp vận tải biển, dầu khí và khai thác thủy sản: Hỗ trợ các thuyền trưởng tối ưu hóa lộ trình hàng hải dựa trên dòng chảy thực tế để tiết kiệm từ 5 đến 10% nhiên liệu, đồng thời giúp ngư dân nhận diện các tâm xoáy nước hội tụ nguồn lợi hải sản cao gấp 7 - 8 lần vùng biển thông thường.
  • Giảng viên, học viên cao học và sinh viên các trường đại học chuyên ngành: Là học liệu mẫu mực phục vụ nghiên cứu phương pháp luận xử lý dữ liệu viễn thám biển, mô hình hóa toán học Lagrangian và phân tích động lực học vùng nước nông ven bờ.

Câu hỏi thường gặp

1. Hệ thống radar biển HF Seasonde đo đạc dòng chảy tầng mặt dựa trên nguyên lý vật lý nào?
Radar biển HF hoạt động ở tần số 4 - 50 MHz dựa trên nguyên lý tán xạ cộng hưởng Bragg và hiệu ứng Doppler. Sóng radar phát xạ chạm vào mặt biển chuyển động sẽ bị lệch pha tần số tỉ lệ thuận với vận tốc dòng chảy tầng mặt, từ đó tính toán được vận tốc dòng chảy với sai số dưới 7 cm/s.

2. Tại sao cần tối thiểu hai trạm radar ven bờ để tổng hợp trường dòng chảy hai chiều?
Mỗi trạm radar đơn lẻ chỉ đo được thành phần vận tốc hướng tâm (radial velocity) dọc theo tia quét. Để xác định vector dòng chảy thực tế với 2 thành phần kinh hướng và vĩ hướng, cần có sự giao thoa từ ít nhất hai trạm radar với góc cắt chéo lý tưởng lớn hơn 30° nhằm tránh suy biến hình học.

3. Hoàn lưu dòng chảy Vịnh Bắc Bộ trong mùa hè có điểm gì khác biệt cốt lõi so với mùa đông?
Trong mùa đông, trung tâm vịnh tồn tại một xoáy thuận ổn định với vận tốc 10 - 15 cm/s và dòng chảy ven bờ chảy mạnh về phía nam đạt 25 - 30 cm/s. Ngược lại, vào mùa hè, hoàn lưu trung tâm đảo chiều thành một xoáy nghịch với vận tốc 8 - 12 cm/s do khối nước phía nam đẩy lên.

4. Mô hình quỹ đạo Lagrangian áp dụng số liệu radar đạt độ tin cậy bao nhiêu trong mô phỏng vật thể trôi?
Mô hình quỹ đạo Lagrangian kết hợp phương pháp Monte Carlo và trường dòng chảy radar đạt khoảng tin cậy 95% sau các bước hiệu chỉnh thông số. Khi áp dụng hệ số dòng chảy chuẩn 0,787 và hệ số khuếch tán rối 56 $\text{m}^2/\text{s}$, mô hình tái hiện chính xác quỹ đạo trôi của tàu cá trong 30 ngày.

5. Ứng dụng thực tế nổi bật nhất của số liệu radar biển trong bảo vệ môi trường biển là gì?
Số liệu radar thời gian thực cho phép mô phỏng chính xác hướng lan truyền của dầu tràn và nước thải độc hại. Điển hình như các ứng dụng quốc tế đã giúp ngăn chặn việc xả thải vào các rạn san hô bảo tồn hoặc dự báo kịp thời vệt dầu tràn dạt vào bờ biển.

Kết luận

  • Làm chủ công nghệ xử lý dữ liệu radar HF: Luận văn là công trình đầu tiên tại Việt Nam khai thác thành công chuỗi dữ liệu thực nghiệm liên tục hơn 38.000 giờ đo từ mạng lưới 3 trạm radar Seasonde tại Vịnh Bắc Bộ.
  • Phát hiện mới về động lực học biển: Khám phá chi tiết sự đảo chiều hoàn lưu theo mùa và lần đầu tiên ghi nhận cấu trúc xoáy đôi (dipole) quy mô 70 - 90 km tồn tại ổn định tại phía bắc vịnh trong các tháng mùa thu.
  • Đóng góp giá trị mô phỏng thực tiễn: Ứng dụng thành công mô hình Lagrangian hai chiều để mô phỏng chính xác đường đi của tàu thuyền trôi dạt trong cả kịch bản trung bình mùa và kịch bản thời gian thực.
  • Định hướng phát triển tương lai: Cần tiếp tục tích hợp trường gió động lực học và mở rộng hệ thống trạm quan trắc trên toàn vùng biển Việt Nam theo lộ trình đề án quốc gia.
  • Lời kêu gọi hành động: Các cơ quan quản lý hàng hải, cứu hộ cứu nạn và nghiên cứu biển cần nhanh chóng ứng dụng cơ sở dữ liệu radar biển vào hệ thống điều hành tác nghiệp nhằm nâng cao hiệu quả bảo vệ chủ quyền và phát triển bền vững kinh tế biển Việt Nam.