Tổng quan nghiên cứu

Biển Đông là một trong những vùng biển nhiệt đới lớn nhất trên thế giới với diện tích khoảng 3,5 triệu km², giữ vị trí chiến lược đặc biệt quan trọng về quốc phòng, an ninh và phát triển kinh tế biển. Để phục vụ công tác dự báo hạn ngắn các trường khí tượng thủy văn, các nhà khoa học biển phụ thuộc rất lớn vào việc xây dựng và phát triển các mô hình số trị thủy động lực học 3D. Tuy nhiên, một hạn chế lớn của các mô hình hiện nay như Delft 3D hay ROMS là thường xem điều kiện tại các biên lỏng là hằng số điều hòa mực nước hoặc chỉ cập nhật số liệu nhiệt - muối bề mặt, chưa xét đến cấu trúc dòng chảy và trao đổi nước phân tầng theo chiều thẳng đứng qua các eo biển.

Nghiên cứu được thực hiện nhằm mục tiêu xác định quy luật và định lượng thông lượng trao đổi nước tại các biên lỏng của Biển Đông thông qua trường nhiệt độ và độ muối trung bình nhiều năm. Phạm vi không gian nghiên cứu bao phủ toàn bộ vùng biển từ vĩ độ 0°15'B đến 32°15'B và từ kinh độ 100°15'Đ đến 123°15'Đ, tập trung vào 5 eo biển xung yếu: eo Đài Loan, eo Bashi, eo Mindoro, eo Balabac và eo Singapore. Chuỗi số liệu phân tích được thu thập liên tục trong 50 năm từ năm 1958 đến năm 2008. Kết quả nghiên cứu cung cấp bộ dữ liệu biên lỏng chuẩn xác, nâng cao độ tin cậy của các mô hình số trị dự báo hoàn lưu và các trường thủy văn biển phục vụ kinh tế và an ninh hàng hải.

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 ứng dụng lý thuyết hoàn lưu đại dương và động lực học chất lỏng địa vật lý, trong đó trọng tâm là lý thuyết dòng chảy địa chuyển và dòng chảy mật độ trong môi trường chất lỏng phân tầng nghiêng áp. Sự phân bố không đồng đều của bức xạ mặt trời, lượng bốc hơi, lượng giáng thủy và trao đổi nhiệt biển - khí quyển tạo ra trường nhiệt độ và độ muối biến thiên, kéo theo sự nghiêng của các mặt đẳng mật độ và mặt đẳng áp. Dưới tác dụng của trường trọng lực và lực Coriolis, sự cân bằng địa chuyển được hình thành khi gradien áp suất theo phương ngang cân bằng hoàn toàn với lực Coriolis tác động vuông góc với hướng chuyển động.

Cơ sở toán học của phương pháp động lực dựa trên việc tích phân phương trình thủy tĩnh từ mặt đẳng áp tham chiếu lên mặt biển để xác định độ cao động lực của từng trạm hải văn. Nghiên cứu sử dụng công thức Bjerknes để tính toán thể tích riêng quy ước của nước biển, kết hợp các số hạng hiệu chỉnh theo nhiệt độ, độ muối và áp suất thủy tĩnh. Dòng chảy gradien và dòng địa chuyển được mô tả qua mối quan hệ vi phân giữa gradien độ cao động lực và vận tốc dòng chảy trực giao với mặt cắt.

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

Nguồn dữ liệu đầu vào sử dụng bộ số liệu tái phân tích và đồng hóa đại dương SODA được trích xuất từ máy chủ OPeNDAP thuộc Trung tâm Đồng hóa Dữ liệu Hải dương của Đại học Texas A&M. Bộ dữ liệu có độ phân giải không gian lưới 0,5° x 0,5° kinh vĩ, gồm đầy đủ thông số nhiệt độ và độ muối từ tầng mặt đến đáy theo 12 tháng trong giai đoạn 50 năm từ 1958 đến 2008. Lý do chọn bộ dữ liệu này xuất phát từ thực tế tổng số 40.678 trạm đo đạc mặt rộng trong nước và quốc tế qua 4.812 chuyến khảo sát tại Biển Đông phân bố không đồng đều, thiếu hẳn chuỗi đo liên tục đồng thời trên cùng mặt cắt.

Quy trình tính toán vận tốc và lưu lượng dòng chảy qua mặt cắt được thực hiện qua các bước:

  1. Tính toán giá trị nhiệt độ và độ muối trung bình tháng nhiều năm cho từng điểm lưới không gian.
  2. Xác định tỷ trọng chuẩn ở 0°C và tính mật độ quy ước theo công thức hải dương học tiêu chuẩn.
  3. Tính thể tích riêng quy ước tại các tầng sâu và áp dụng thuật toán bổ sung độ cao động lực cho các cặp trạm có độ sâu chênh lệch.
  4. Tích phân độ cao động lực từ mặt chuẩn hoặc từ đáy lên mặt biển để dựng bản đồ đẳng vận tốc trực giao qua 5 mặt cắt eo biển đại diện cho 12 tháng và 2 mùa gió chính.

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

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

Nghiên cứu đã định lượng hóa bức tranh phân bố nhiệt muối và cấu trúc trao đổi nước qua 5 biên lỏng của Biển Đông với những kết quả nổi bật:

  • Tại eo Bashi (độ sâu từ 2.341 m đến 2.600 m): Đây là cửa ngõ trao đổi nước mạnh nhất của Biển Đông với Tây Thái Bình Dương. Vào mùa đông, khối nước độ muối cao trên 34‰ xâm nhập rất mạnh vào Biển Đông qua tầng mặt và tầng trung gian. Vận tốc dòng chảy địa chuyển qua mặt cắt đạt cực đại tại tầng mặt và suy giảm dần xuống độ sâu 1.000 m.
  • Tại eo Đài Loan (độ sâu trung bình 69 m đến 70 m, chiều rộng khoảng 127 km): Cấu trúc dòng chảy thể hiện tính mùa tương phản sâu sắc. Mùa đông, dưới tác động của gió mùa Đông Bắc, dòng nước lạnh chảy từ biển Hoa Đông xuống phía Nam; mùa hè, dòng chảy đổi hướng hoàn toàn chảy về phía Bắc thoát ra biển Hoa Đông, độ muối có sự phân tầng mạnh mẽ theo phương thẳng đứng.
  • Tại eo Mindoro (độ sâu 329 m đến 450 m) và eo Balabac (độ sâu 49 m đến 100 m): Các đường đẳng nhiệt và đẳng muối phân bố nghiêng rõ rệt từ Đông sang Tây, phản ánh sự lưu chuyển nước liên tục giữa Biển Đông và các biển nội địa Philippines và biển Sulu.
  • Tại eo Singapore (độ sâu nông dưới 70 m): Mùa đông toàn bộ tầng nước mang tính đồng nhất cao với nhiệt độ ổn định, trong khi mùa hè xuất hiện sự phân tầng nhiệt yếu từ 29°C ở tầng mặt xuống 27°C ở sát đáy.
+------------------+------------------+-------------------+-------------------+
|     Eo biển      |  Độ sâu đáy (m)  | Nhiệt độ TB hè/đông| Độ muối TB (‰)   |
+------------------+------------------+-------------------+-------------------+
| Eo Bashi         |   2.341 - 2.600  |    28°C / 23°C    |    33,5 - 34,5    |
| Eo Đài Loan      |      69 - 70     |    27°C / 18°C    |    31,0 - 33,0    |
| Eo Mindoro       |     329 - 450    |    29°C / 25°C    |    33,0 - 34,2    |
| Eo Balabac       |      49 - 100    |    29°C / 26°C    |    32,5 - 33,8    |
| Eo Singapore     |      35 - 70     |    29°C / 27°C    |    30,5 - 32,5    |
+------------------+------------------+-------------------+-------------------+

Thảo luận kết quả

Các kết quả tính toán dòng chảy mật độ phù hợp chặt chẽ với quy luật hoàn lưu gió mùa trên Biển Đông. Mùa đông, sự chênh lệch nhiệt độ tầng mặt giữa khu vực Tây Bắc và Đông Nam Biển Đông lên tới 8°C tạo ra một gradien mật độ ngang cực lớn, kết hợp với trường ứng suất gió Đông Bắc hình thành nên hoàn lưu xoáy thuận quy mô lớn trên toàn biển và làm gia tăng vận tốc dòng chảy ven bờ miền Trung Việt Nam lên mức 0,75 m/s đến 0,85 m/s, thậm chí vượt quá 1,0 m/s tại một số thời điểm.

Khi thể hiện dữ liệu qua các biểu đồ mặt cắt thẳng đứng 2 chiều, các đường đẳng vận tốc cho thấy rõ sự phân hóa giữa dòng chảy tầng mặt do gió chi phối và dòng chảy ngầm sâu do gradien mật độ dẫn dắt. Tại eo Bashi, dòng nhánh Kuroshio uốn lượn xâm nhập tạo nên các hoàn lưu xoáy thứ cấp ở Tây Bắc đảo Luzon. Sự tương đồng cao giữa kết quả tính toán động lực từ chuỗi số liệu nhiệt muối SODA với các kết quả khảo sát thực tế trong các Atlas Khí tượng Thủy văn Biển Đông khẳng định tính khả thi của phương pháp trong việc cung cấp thông số đầu vào cho mô hình toán học.

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

Dựa trên kết quả nghiên cứu, tác giả đưa ra 4 nhóm giải pháp cụ thể:

  1. Tích hợp ma trận vận tốc dòng chảy biên lỏng 3D vào các hệ thống mô hình dự báo hải văn quốc gia (như ROMS, Delft 3D, POM) do Trung tâm Dự báo Khí tượng Thủy văn Quốc gia vận hành trong giai đoạn 2026-2028, đặt mục tiêu giảm sai số dự báo mực nước và dòng chảy xuống dưới 10%.
  2. Thiết lập chương trình quan trắc thủy văn định kỳ 6 tháng/lần bằng các thiết bị hiện đại (như phao ngầm tự hành Argo, đầu đo CTD sâu) dọc theo 2 mặt cắt trọng yếu là eo Bashi và eo Đài Loan nhằm cập nhật cơ sở dữ liệu kiểm chứng thực tế cho Viện Hải dương học và Tổng cục Biển và Hải đảo Việt Nam.
  3. Xây dựng nền tảng cơ sở dữ liệu số hóa khí tượng - hải văn chuẩn hóa OpenDAP cấp quốc gia trước năm 2027, cho phép các cơ sở nghiên cứu kết nối, đồng hóa và chia sẻ dữ liệu nhiệt - muối đa tầng cho toàn bộ vùng biển đặc quyền kinh tế 200 hải lý của Việt Nam.
  4. Ứng dụng ngay số liệu thông lượng trao đổi nước qua các eo biển vào bài toán mô phỏng lan truyền chất ô nhiễm và sự cố tràn dầu trên các tuyến hàng hải quốc tế, do Ủy ban Quốc gia Ứng phó Sự cố Thiên tai và Tìm kiếm Cứu nạn chủ trì triển khai.

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

  • Chuyên gia mô hình hóa khí tượng thủy văn và hải dương học: Sử dụng thuật toán và bộ giá trị phân bố vận tốc qua 5 eo biển để làm điều kiện biên chính xác cho các bài toán mô phỏng 3D.
  • Kỹ sư dự báo tác nghiệp biển và an toàn hàng hải: Ứng dụng các quy luật dòng chảy mật độ mùa để tối ưu hóa hải trình tàu thuyền và phục vụ công tác cứu hộ, cứu nạn trên Biển Đông.
  • Cơ quan quản lý tài nguyên môi trường biển: Tham khảo số liệu trao đổi nước để đánh giá sức tải môi trường, dự báo lan truyền ô nhiễm và hoạch định chiến lược bảo vệ hệ sinh thái biển ven bờ.
  • Giảng viên, học viên cao học và nghiên cứu sinh chuyên ngành Khoa học Trái đất: Sử dụng tài liệu như một cẩm nang học thuật mẫu mực về phương pháp động lực tính dòng chảy biển từ phương trình thủy tĩnh và số liệu nhiệt muối thực nghiệm.

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

1. Phương pháp động lực tính dòng chảy biển có ưu điểm gì nổi bật? Phương pháp động lực cho phép gián tiếp xác định trường dòng chảy địa chuyển và thông lượng trao đổi nước trên quy mô không gian rộng lớn chỉ từ số liệu phân bố nhiệt độ và độ muối, giải quyết bài toán thiếu hụt thiết bị đo dòng chảy trực tiếp ngoài biển sâu vốn rất tốn kém và khó duy trì liên tục.

2. Tại sao điều kiện biên lỏng lại quyết định độ chính xác của mô hình số trị? Biên lỏng là nơi tiếp giáp giữa miền tính toán với đại dương mở. Nếu điều kiện biên chỉ gán giá trị mực nước tĩnh hoặc bỏ qua dòng chảy phân tầng theo độ sâu, mô hình sẽ không phản ánh đúng động năng truyền tải vào lưu vực, gây sai số tích lũy lớn đối với hoàn lưu 3D bên trong nội miền.

3. Sự khác biệt lớn nhất của dòng chảy qua eo Đài Loan giữa hai mùa là gì? Vào mùa đông, dưới áp lực gió Đông Bắc, nước lạnh chảy từ phía Bắc xuống phía Nam dọc theo eo biển; sang mùa hè, hướng dòng chảy đảo ngược hoàn toàn từ Nam lên Bắc, đưa nước ấm từ Biển Đông thoát ra biển Hoa Đông với sự phân tầng độ muối sâu sắc.

4. Chuỗi số liệu SODA 50 năm có ý nghĩa gì đối với nghiên cứu này? Số liệu SODA giai đoạn 1958-2008 cung cấp trường nhiệt muối 3D đồng hóa khách quan với bước lưới 0,5° x 0,5°, giúp khắc phục hoàn toàn sự gián đoạn và thiếu đồng bộ của các trạm đo thực nghiệm, tạo ra chuỗi giá trị trung bình tháng có độ tin cậy thống kê cao.

5. Cách xử lý sai số độ cao động lực giữa hai trạm hải văn có độ sâu khác nhau? Khi hai trạm có độ sâu đáy chênh lệch, tác giả áp dụng công thức hiệu chỉnh bổ sung thể tích riêng trung bình ở lớp đáy cho trạm nông hơn, đảm bảo hai trạm có cùng một mặt đẳng áp chuẩn để so sánh và tính toán chính xác độ chênh động lực.

Kết luận

  • Luận văn đã thiết lập thành công quy trình chuẩn hóa tính toán dòng chảy mật độ qua các biên lỏng Biển Đông dựa trên lý thuyết động lực học chất lỏng địa vật lý.
  • Xây dựng hoàn chỉnh bộ bản đồ phân bố nhiệt độ, độ muối và trường vận tốc dòng chảy thẳng đứng qua 5 eo biển xung yếu cho 12 tháng và 2 mùa gió đặc trưng.
  • Làm sáng tỏ vai trò kiểm soát vượt trội của eo Bashi và eo Đài Loan trong cơ chế trao đổi khối lượng và năng lượng giữa Biển Đông và Tây Thái Bình Dương.
  • Cung cấp bộ dữ liệu điều kiện biên 3D chuẩn mực cho các mô hình thủy động lực học phục vụ công tác dự báo tác nghiệp khí tượng thủy văn biển tại Việt Nam.
  • Trong giai đoạn 2026-2030, việc mở rộng tích hợp dữ liệu vệ tinh độ cao mặt biển và phao trôi tự hành sẽ tiếp tục nâng tầm hoàn thiện các mô hình dự báo biển quốc gia.