Tổng quan nghiên cứu

Vùng Vịnh Bắc Bộ có diện tích khoảng 160.950 km² với chiều dài bờ biển phía Việt Nam trải dài hơn 740 km, giữ vị trí chiến lược đặc biệt quan trọng trong phát triển kinh tế biển và an ninh quốc phòng. Tuy nhiên, công tác điều tra cơ bản và giám sát động lực học biển truyền thống trước đây chủ yếu dựa vào các chuyến khảo sát ngắn ngày bằng tàu thuyền hoặc các trạm hải văn ven bờ đơn lẻ, dẫn đến sự thiếu hụt nghiêm trọng nguồn dữ liệu liên tục theo thời gian thực và đồng bộ trên diện rộng. Nhằm giải quyết điểm nghẽn này, dự án xây dựng hệ thống radar biển tần số cao (HF) đã được triển khai với trạm đầu tiên vận hành từ năm 2011 tại Đồng Hới, tiếp nối là giai đoạn 1 hoàn thành 3 trạm quan trắc ven bờ và 1 trạm điều hành trung tâm tại Hà Nội.

Mục tiêu cốt lõi của nghiên cứu là khai thác và đánh giá toàn diện nguồn dữ liệu radar biển thu thập trong giai đoạn từ tháng 6 năm 2012 đến tháng 11 năm 2014 tại Vịnh Bắc Bộ. Nghiên cứu tập trung giải quyết hai bài toán học thuật và thực tiễn lớn: tính toán các đặc trưng trường dòng chảy tầng mặt trung bình tháng trên lưới tính chi tiết 55x63 điểm với bước lưới 6 km, và ứng dụng thử nghiệm mô hình động lực học để phân tích, dự báo quỹ đạo vật thể trôi dạt trên biển trong thời gian 7 ngày liên tục.

Kết quả nghiên cứu mang ý nghĩa khoa học và thực tiễn vượt bậc, chứng minh độ chính xác của công nghệ radar HF với sai số vận tốc dòng chảy kiểm soát dưới 7 cm/s và sai số hướng dưới 10 độ. Nguồn dữ liệu này tạo tiền đề vững chắc cho công tác tìm kiếm cứu nạn hàng hải, định vị vệt dầu tràn, hỗ trợ khai thác thủy hải sản và cung cấp dữ liệu đầu vào chuẩn xác để đồng hóa cho các mô hình số trị dự báo khí tượng thủy văn biển tại Việt Nam.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng kết hợp giữa lý thuyết tán xạ sóng điện từ cao tần và cơ học chất lưu hải dương học. Nguyên lý cốt lõi thứ nhất là hiện tượng tán xạ Bragg và hiệu ứng Doppler của sóng vô tuyến tần số cao HF trong dải tần 4-50 MHz (tương ứng bước sóng 6-75 m). Tín hiệu vô tuyến phát ra từ anten trạm radar khi chạm vào các sóng biển chuyển động sẽ phản hồi lại với độ lệch pha tần số Doppler Shift tỷ lệ thuận với vận tốc pha của sóng biển và thành phần vận tốc dòng chảy tầng mặt hướng tâm.

Nguyên lý cốt lõi thứ hai là mô hình vận chuyển Lagrangian mô tả chuyển động của các hạt phần tử nước hoặc vật thể trôi dạt trong trường dòng chảy bề mặt và trường gió khí quyển. Phương trình dịch chuyển hạt kết hợp vận tốc bình lưu và hệ số khuếch tán rối được biểu diễn qua quan hệ tuyến tính giữa vận tốc dòng chảy tầng mặt và vận tốc gió tại độ cao 10 m, trong đó hệ số kéo khí động học của gió dao động trong khoảng 2,5% đến 4,4% (trung bình đạt 3,5%) và hệ số khuếch tán rối đạt 56 m²/s.

Các khái niệm chuyên ngành then chốt được chuẩn hóa trong nghiên cứu bao gồm:

  1. Sóng cao tần HF (High Frequency Radar): Dải sóng điện từ phản xạ bề mặt biển phục vụ đo đạc hải văn tầm xa tới 300 km.
  2. Vận tốc dòng chảy hướng tâm (Radial Velocity): Thành phần vận tốc dòng chảy chiếu theo phương xuyên tâm nối từ vị trí trạm radar đến điểm quan trắc.
  3. Trường vectơ dòng chảy 2 chiều (2D Total Vector): Trường vectơ tổng hợp gồm hai thành phần trực giao kinh hướng và vĩ hướng (U, V) được tái tạo tại từng nút lưới tính.
  4. Hệ thống dự báo dòng chảy ngắn hạn STPS (Short-Term Predictive System): Công cụ mô hình hóa số trị kết hợp dữ liệu quá khứ để dự báo trường vận tốc dòng chảy trong tương lai gần.

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ừ hệ thống 3 trạm radar biển HF đặt tại Hòn Dáu (Hải Phòng), Nghi Xuân (Hà Tĩnh) và Đồng Hới (Quảng Bình) trong khung thời gian từ ngày 05/06/2012 đến ngày 15/11/2014 với bước thời gian 1 giờ/lần đo. Cỡ mẫu nghiên cứu rất lớn với hơn 11.163 lượt quan trắc tại trạm Hòn Dáu và hơn 12.000 lượt quan trắc tại trạm Đồng Hới, phản ánh đầy đủ các biến động hải văn qua các mùa trong năm.

Phương pháp chọn mẫu và xây dựng miền tính được thực hiện trên hệ tọa độ chữ nhật giới hạn từ vĩ tuyến 17,7°N đến 21,13°N và kinh tuyến từ 105,0°E đến 108,56°E với tâm lưới đặt tại 19°24,6'N và 107°00'E. Miền tính bao gồm 55 điểm theo phương kinh tuyến và 63 điểm theo phương vĩ tuyến với độ phân giải không gian đạt 6 km.

Lý do lựa chọn phương pháp phân tích:

  1. Phương pháp bình phương tối thiểu (Least Squares Method): Lựa chọn phần mềm chuyên dụng SeaSonde Combine Suite Release 7 để giải bài toán tổng hợp vectơ 2D từ ít nhất hai vectơ vận tốc hướng tâm giao thoa, giúp triệt tiêu sai số ngẫu nhiên tại các vùng nước mở.
  2. Phương pháp nội suy đường ranh giới tuyến tính: Được áp dụng đối với các nút lưới có góc giao thoa giữa hai tia xuyên tâm nhỏ hơn 30 độ (alpha, beta = 30°), khắc phục hiện tượng suy biến ma trận và lấp đầy khoảng trống dữ liệu sát bờ cũng như ngoài khơi xa.
  3. Phương pháp tích phân Euler và mô phỏng Monte Carlo: Ứng dụng trong mô hình TESEO và CODAR Leeway nhằm tính toán quỹ đạo vật thể trôi dạt với khoảng tin cậy thống kê 95%, đảm bảo độ tin cậy cao nhất khi so sánh với dữ liệu thực tế.

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

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

Thứ nhất, nghiên cứu đã định lượng hóa chất lượng và độ phủ của dữ liệu quan trắc radar biển tại Vịnh Bắc Bộ. Trạm Nghi Xuân đạt mật độ số liệu phân bố đồng đều nhất với tỷ lệ hoàn chỉnh từ 70% đến 75%, trong khi trạm Hòn Dáu và Đồng Hới đạt mật độ từ 55% đến 60%. Sai số trung bình của vận tốc dòng chảy đo bằng radar HF trong điều kiện bình thường duy trì ở mức dưới 7 cm/s và sai số hướng dưới 10 độ, hoàn toàn tương đồng với các nghiên cứu quốc tế sử dụng thiết bị đo dòng tự ghi ADCP có sai số quân phương dưới 5 cm/s.

Thứ hai, cấu trúc hoàn lưu tầng mặt Vịnh Bắc Bộ được tái hiện chi tiết theo quy luật mùa rõ rệt. Trong các tháng mùa đông (tháng 12, tháng 1, tháng 2), dòng chảy ven bờ phía Tây Vịnh Bắc Bộ chảy mạnh theo hướng Bắc - Nam dọc bờ biển Việt Nam với vận tốc trung bình phổ biến từ 20 cm/s đến 45 cm/s, nhận thêm khối nước bổ sung từ phía Nam đảo Hải Nam và eo Quỳnh Châu đi vào. Ngược lại, vào mùa hè, hoàn lưu biển đổi hướng với sự xuất hiện của cấu trúc xoáy đôi đối lập dấu gồm xoáy nghịch ở phía Bắc và xoáy thuận ở khu vực trung tâm vịnh.

Thứ ba, các thử nghiệm dự báo quỹ đạo vật thể trôi dạt 7 ngày trong tháng 7 và tháng 11 năm 2014 cho thấy sự phân hóa mạnh mẽ do trường dòng chảy chi phối. Trong tháng 7/2014, hạt trôi di chuyển quãng đường hơn 35 km về hướng Tây Bắc dưới tác động của hoàn lưu mùa hè. Trong tháng 11/2014, dưới tác động của gió mùa Đông Bắc và dòng ven bờ chảy xiết, vật thể di chuyển nhanh về hướng Nam Tây Nam với quãng đường trôi vượt 120 km chỉ sau 168 giờ tính toán, với dải phân tán được bao bọc trọn vẹn trong elip tin cậy 95%.

Thảo luận kết quả

Cơ chế hình thành bức tranh hoàn lưu mùa được giải thích bởi sự kết hợp giữa trường gió mùa thịnh hành (gió mùa Đông Bắc trong mùa đông và gió mùa Tây Nam/Đông Nam trong mùa hè), tương tác nhiệt biển khí và quá trình trao đổi nước qua cửa vịnh phía Nam rộng khoảng 230 km (độ sâu trên 100 m) cùng eo biển Quỳnh Châu hẹp 18 km (độ sâu 20 m). Kết quả tính toán từ số liệu radar đã kiểm chứng và làm mịn bức tranh hoàn lưu trong các báo cáo điều tra tổng hợp Vịnh Bắc Bộ năm 1964-1965 và Atlas Quốc gia Việt Nam năm 1999, đồng thời giải thích rõ cấu trúc xoáy quy mô nhỏ mà các nghiên cứu trước đây của các tác giả quốc tế chưa thể hiện được do độ phân giải dữ liệu thô.

Về phương thức biểu diễn trực quan học thuật:

  1. Trường vectơ dòng chảy 2 chiều được biểu diễn tối ưu qua bản đồ phân bố vectơ mũi tên (Quiver Plot), trong đó độ dài và màu sắc mũi tên thể hiện độ lớn vận tốc (cm/s), góc mũi tên chỉ hướng dòng chảy tại từng nút lưới 6x6 km.
  2. Mật độ khả dụng của trạm radar được minh họa trực quan qua biểu đồ đẳng trị dạng nhiệt (Heatmap), thể hiện dải phần trăm bao phủ từ 50% đến 80% theo cự ly quét từ trạm ra xa đến 300 km.
  3. Quỹ đạo vật thể trôi dạt được mô tả qua biểu đồ phân tán đường đi nhiều nhánh (Spaghetti Diagram) kết hợp đường bao elip xác suất 95%, giúp các đơn vị tìm kiếm trực quan hóa ngay vùng quét ưu tiên theo từng mốc thời gian 24h, 48h, 72h và 168h.

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

  1. Mở rộng và hoàn thiện mạng lưới quan trắc radar biển quốc gia: Đẩy nhanh tiến độ xây dựng 5 trạm radar biển trong giai đoạn II (gồm 2 trạm tại khu vực Vịnh Bắc Bộ và 3 trạm tại vùng biển Nam Bộ), hướng tới mục tiêu hoàn thành toàn bộ 17 trạm radar trải dài dọc bờ biển và hải đảo Việt Nam trước năm 2030, nâng tầm kiểm soát hải văn liên tục trên vùng đặc quyền kinh tế.
  2. Tích hợp trực tiếp dữ liệu radar thời gian thực vào quy trình cứu nạn hàng hải: Trung tâm Hải văn phối hợp với Ủy ban Quốc gia Ứng phó Sự cố, Thiên tai và Tìm kiếm Cứu nạn cài đặt module phần mềm tính toán trôi dạt tự động CODAR Leeway, đặt mục tiêu rút ngắn 50% thời gian xác định vị trí tàu thuyền tai nạn và giảm thiểu 40% diện tích vùng tìm kiếm trên biển trong 24 giờ đầu tiên.
  3. Chuẩn hóa quy trình đồng hóa số liệu dòng chảy bề mặt vào các mô hình số trị 3D: Viện Hải dương học và các trường đại học chuyên ngành thiết lập khung thuật toán đồng hóa số liệu radar vào các mô hình động lực học biển phân giải cao (như ROMS, POM), phấn đấu hoàn thành bộ quy chuẩn kỹ thuật trong vòng 2 năm tới nhằm nâng độ chính xác dự báo hải văn 72 giờ lên trên 85%.
  4. Phát triển bản đồ cảnh báo môi trường và ngư trường khai thác thủy sản: Tổng cục Thủy sản cùng Tổng cục Biển và Hải đảo định kỳ xử lý dữ liệu radar để công bố bản đồ vị trí tâm xoáy nước và dải hội tụ dòng chảy hàng tuần, hỗ trợ ngư dân nâng cao năng suất đánh bắt hải sản gấp 5 đến 7 lần và chủ động khoanh vùng ứng phó nếu xảy ra sự cố tràn dầu trên diện tích 160.950 km² của vịnh.

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

  1. Cơ quan quản lý nhà nước và lực lượng cứu nạn hàng hải: Điển hình là Cục Hàng hải Việt Nam, Tổng cục Biển và Hải đảo, Trung tâm Phối hợp Tìm kiếm Cứu nạn Hàng hải (MRCC). Ứng dụng nguồn dữ liệu dòng chảy thực tế để điều hành phân luồng giao thông tàu biển, chỉ đạo phương án cứu hộ cứu nạn và xây dựng kịch bản ứng phó sự cố tràn dầu ven bờ.
  2. Viện nghiên cứu và nhà khoa học chuyên ngành Hải dương - Khí tượng: Bao gồm các nhà nghiên cứu về vật lý biển, động lực học cửa sông ven biển và mô hình hóa số trị. Luận văn cung cấp phương pháp luận chặt chẽ về xử lý tín hiệu radar HF, kỹ thuật nội suy biên và thuật toán tổng hợp vectơ 2D.
  3. Doanh nghiệp vận tải biển, khai thác dầu khí và năng lượng tái tạo: Các tổng công ty vận tải thủy, ban quản lý dự án điện gió ngoài khơi khai thác thông số sóng và dòng chảy tại 10 cell quan trắc (cự ly 5,8 km đến 58,3 km) để thiết kế kết cấu công trình biển an toàn và tối ưu hóa hải trình tiết kiệm nhiên liệu.
  4. Giảng viên, nghiên cứu sinh và học viên cao học: Thuộc các chuyên ngành Hải dương học, Khí tượng - Thủy văn, Kỹ thuật Bờ biển và Địa tin học. Luận văn là tài liệu tham khảo mẫu mực về kỹ thuật phân tích chuỗi thời gian hải văn và mô hình hóa chuyển động hạt Lagrangian.

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

  1. Công nghệ Radar sóng cao tần (HF) đo đạc dòng chảy tầng mặt biển dựa trên nguyên lý nào? Hệ thống radar HF (như SeaSonde) hoạt động trên nguyên lý tán xạ Bragg và dịch pha tần số Doppler Shift. Sóng điện từ phát ra ở tần số 4-50 MHz tương tác với sóng biển có bước sóng bằng một nửa bước sóng radar (6-75 m). Tốc độ dịch chuyển của sóng biển do dòng chảy mang lại làm biến đổi tần số sóng phản hồi, từ đó phần mềm giải mã chính xác vận tốc và hướng dòng chảy tầng mặt.

  2. Hệ thống radar biển Vịnh Bắc Bộ có phạm vi quan trắc và độ phân giải không gian ra sao? Mỗi trạm radar HF trong hệ thống có tầm quét quan trắc dòng chảy tối đa lên tới 300 km từ bờ. Dữ liệu sóng được đo đạc chi tiết tại 10 cell trên tia trung tuyến với cự ly từ 5,8 km đến 58,3 km. Trên toàn miền tính Vịnh Bắc Bộ, lưới tính được thiết lập với 55x63 điểm nút, tương ứng độ phân giải không gian là 6 km x 6 km và tần suất cập nhật dữ liệu liên tục 1 giờ/lần.

  3. Độ chính xác của dữ liệu dòng chảy đo bằng radar biển có đáp ứng yêu cầu thực tế? Dữ liệu quan trắc từ radar biển có độ tin cậy rất cao. Trong điều kiện môi trường bình thường, sai số vận tốc dòng chảy trung bình đạt dưới 7 cm/s và sai số hướng dưới 10 độ. Khi kiểm chứng chéo với các thiết bị đo dòng chuyên dụng tầng sâu ADCP, sai số trung bình quân phương của vận tốc dòng chảy tầng mặt đo bằng radar được xác nhận dưới 5 cm/s.

  4. Tại sao cần áp dụng phương pháp nội suy đường ranh giới khi tổng hợp trường dòng chảy 2 chiều? Khi hai tia vận tốc hướng tâm từ 2 trạm radar tạo với nhau một góc nhỏ hơn 30 độ (vùng ranh giới giữa 2 trạm hoặc vùng rất xa bờ), phép tính bình phương tối thiểu sẽ phát sinh sai số lớn do ma trận suy biến. Việc áp dụng phương pháp nội suy tuyến tính đường ranh giới với các tham số góc alpha, beta bằng 30 độ giúp tính toán thành phần trực giao chính xác và lấp đầy khoảng trống dữ liệu.

  5. Mô hình Lagrangian tích hợp dữ liệu radar hỗ trợ công tác tìm kiếm cứu nạn như thế nào? Mô hình Lagrangian kết hợp trường dòng chảy thời gian thực từ radar và trường gió 10 m (với hệ số kéo khí động học trung bình 3,5% và hệ số khuếch tán 56 m²/s) để tích phân phương trình dịch chuyển hạt theo thời gian. Thuật toán cung cấp vị trí dự báo kèm đường bao xác suất 95%, giúp lực lượng cứu hộ khoanh vùng chính xác mục tiêu trôi dạt sau 24 đến 168 giờ.

Kết luận

  • Luận văn đã xử lý và đánh giá thành công tập dữ liệu hải văn đồ sộ từ 3 trạm radar biển HF đầu tiên tại Việt Nam với hơn 35.000 giờ quan trắc thực tế trong giai đoạn 2012-2014.
  • Xây dựng hoàn chỉnh lưới tính 55x63 điểm với bước lưới 6 km, làm chủ quy trình công nghệ tổng hợp trường vectơ dòng chảy 2 chiều và kỹ thuật xử lý ranh giới góc hẹp dưới 30 độ.
  • Xác lập bức tranh hoàn lưu chi tiết theo mùa tại Vịnh Bắc Bộ, chứng minh sự chiếm ưu thế của dòng chảy ven bờ Tây hướng Nam vào mùa đông với vận tốc 20-45 cm/s và hoàn lưu xoáy đôi phân hóa trong mùa hè.
  • Ứng dụng thành công mô hình quỹ đạo hạt Lagrangian mô phỏng chính xác đường đi của vật thể trôi dạt trong 7 ngày với độ tin cậy thống kê 95%, kiểm soát sai số vận tốc dưới 7 cm/s.
  • Khẳng định giá trị chiến lược của công nghệ radar biển trong việc hiện đại hóa công tác cảnh báo thiên tai, tìm kiếm cứu nạn hàng hải và bảo vệ an ninh chủ quyền biển đảo.

Kế hoạch tiếp theo cần tập trung hoàn thiện 5 trạm radar giai đoạn II và chuẩn hóa quy trình đồng hóa số liệu radar vào các mô hình dự báo đại dương quốc gia trong vòng 24 tháng tới. Hãy liên hệ các đơn vị nghiên cứu hải dương học chuyên sâu để tiếp cận trọn vẹn thuật toán và nguồn dữ liệu quý giá này ngay hôm nay.