Tổng quan nghiên cứu
Biển Đông là một trong những biển rìa lớn nhất thế giới với diện tích khoảng 3,5 triệu km², sở hữu địa hình phức tạp gồm khoảng 80% diện tích là thềm lục địa nông dưới 200m và 20% diện tích là bồn trũng sâu từ 2000m đến 5500m. Để phục vụ phát triển kinh tế biển và bảo đảm an ninh quốc phòng, công tác dự báo các trường khí tượng thủy văn biển đóng vai trò đặc biệt quan trọng. Trong điều kiện hiện nay, việc ứng dụng hệ thống mô hình số trị ba chiều như Delft3D hay ROMS là giải pháp tối ưu nhằm cung cấp thông tin hải văn chi tiết. Tuy nhiên, độ chính xác của các mô hình thủy lực phụ thuộc rất lớn vào dữ liệu trường ban đầu và điều kiện tại các biên lỏng hở.
Vấn đề nghiên cứu then chốt đặt ra là hầu hết các nghiên cứu mô hình hóa trước đây chỉ thiết lập điều kiện biên dựa trên hằng số điều hòa mực nước hoặc trường nhiệt muối tầng mặt, chưa tính toán chi tiết sự trao đổi nước phân tầng theo độ sâu qua các eo biển. Mục tiêu cụ thể của luận văn là xây dựng quy trình tính toán cấu trúc dòng chảy và lưu lượng trao đổi nước qua các biên lỏng trọng yếu ở Biển Đông dựa trên trường phân bố nhiệt độ và độ mặn. Phạm vi nghiên cứu bao quát các eo biển chiến lược gồm eo Bashi, Đài Loan, Mindoro, Balabac và Singapore trong chu kỳ các tháng đại diện cho mùa đông và mùa hè. Kết quả nghiên cứu có ý nghĩa quyết định trong việc giảm thiểu sai số của các mô hình hoàn lưu xuống dưới 15%, đồng thời cung cấp bộ dữ liệu biên chuẩn xác cho bài toán lan truyền chất ô nhiễm và dự báo cứu nạn trên biển.
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 nền tảng lý thuyết hoàn lưu hải dương học vật lý kết hợp phương pháp động lực học biển cổ điển. Hai khung lý thuyết chính được tích hợp bao gồm: lý thuyết cân bằng địa chuyển mô tả trạng thái cân bằng giữa lực Coriolis và lực gradient áp suất thủy tĩnh; và lý thuyết dòng chảy mật độ trong chất lỏng phân tầng nghiêng áp. Khái niệm độ cao động lực biểu thị công dịch chuyển một đơn vị khối lượng nước biển chống lại trọng lực giữa hai mặt đẳng áp chuẩn.
Bên cạnh đó, khái niệm mặt không động lực được áp dụng để xác định tầng nước sâu từ 1000m đến 1500m ở đại dương nơi vận tốc dòng chảy mật độ tiệm cận giá trị 0 m/s. Trọng lượng riêng quy ước và thể tích riêng quy ước của nước biển được tính toán dựa trên phương trình trạng thái phi tuyến của nước biển, phụ thuộc trực tiếp vào nhiệt độ thực tế, độ muối thực tế và áp suất thủy tĩnh tại từng độ sâu quan trắc.
Phương pháp nghiên cứu
Nguồn dữ liệu của luận văn được tổng hợp từ cơ sở dữ liệu hải văn nhiệt muối nhiều năm với hơn 50 trạm quan trắc chuẩn phân bố trên 8 mặt cắt thẳng đứng qua các eo biển lớn của Biển Đông. Cỡ mẫu nghiên cứu bao gồm hàng nghìn điểm đo nhiệt độ và độ mặn từ tầng mặt xuống đến độ sâu hơn 2000m. Phương pháp chọn mẫu tập trung vào các mặt cắt đại diện có tọa độ địa lý cố định: tuyến dọc 23°15'N qua eo Đài Loan, tuyến 120°45'E qua eo Bashi, tuyến 120°15'E qua eo Mindoro, tuyến 117°45'E qua eo Balabac và tuyến 2°15'N qua eo Singapore.
Lý do lựa chọn phương pháp động lực học kết hợp giải tích số thay vì đo đạc trực tiếp là bởi phương pháp này cho phép tái tạo trường vận tốc góc vuông với mặt cắt từ trường mật độ sẵn có với chi phí tối ưu nhưng vẫn đạt độ chính xác trên 85% so với dữ liệu thả phao trôi thực tế. Quá trình tính toán tích phân phương trình thủy tĩnh để xác định độ cao động lực và dựng bản đồ đẳng vận tốc được thực hiện theo quy trình chuẩn hóa số liệu qua các tháng mùa đông (tháng 1) và mùa hè (tháng 7).
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 chỉ ra rằng eo Bashi (Luzon) với độ sâu đáy từ 2341m đến 2600m là cửa ngõ trao đổi nước quan trọng nhất giữa Biển Đông và Thái Bình Dương. Vào mùa đông, nhánh xâm nhập của dòng hải lưu Kuroshio qua eo Bashi mang khối nước có độ mặn cao trên 34,5‰ vào sâu bồn trũng phía bắc, hình thành nên cấu trúc dòng chảy thẳng đứng nhiều tầng phức tạp.
Thứ hai, tại dải ven biển miền Trung Việt Nam trong mùa gió Đông Bắc, sự kết hợp giữa dòng chảy do gió và dòng chảy nhiệt muối tạo nên dòng chảy ven bờ mạnh đột biến với vận tốc trung bình đạt 0,75 đến 0,85 m/s, nhiều thời điểm vượt ngưỡng cực đại 1,0 m/s. Dòng hải lưu này làm dịch chuyển đường đẳng nhiệt 25°C sâu xuống phía nam tới tận vĩ tuyến 4°N - 5°N.
Thứ ba, vào mùa hè dưới tác động của gió mùa Tây Nam, hoàn lưu Biển Đông đảo chiều thành một hệ thống xoáy nghịch quy mô lớn. Vận tốc dòng chảy trung bình tầng mặt đạt khoảng 0,25 m/s và cực đại đạt 0,5 m/s ngoài khơi Đông Nam Bộ. Dòng nước ấm với nhiệt độ trung bình trên 29°C và độ mặn từ 33‰ đến 34‰ chiếm ưu thế trên toàn bộ vùng biển khơi.
Thứ tư, tại các eo biển phía nam như eo Balabac (độ sâu 49m đến 100m) và eo Singapore (độ sâu dưới 40m), sự trao đổi nước với biển Java và biển Sulu diễn ra chủ yếu ở tầng nông dưới 50m, với biên độ dao động lưu lượng theo mùa đạt mức chênh lệch hơn 40% giữa hai mùa gió chính.
+-------------------------------------------------------------------------+
| ĐẶC TRƯNG TRAO ĐỔI NƯỚC TẠI CÁC EO BIỂN |
+-------------------+----------------+---------------+--------------------+
| Tên eo biển | Chiều rộng (km)| Độ sâu tb (m) | Vận tốc max (m/s) |
+-------------------+----------------+---------------+--------------------+
| Eo Đài Loan | 69 - 70 | 127 | 0,60 - 0,80 |
| Eo Luzon (Bashi) | 2341 - 2600 | 372 | 0,85 - 1,10 |
| Eo Mindoro | 329 - 450 | 78 | 0,40 - 0,55 |
| Eo Balabac | 49 - 100 | 49 | 0,30 - 0,45 |
| Eo Singapore | 12 - 30 | 35 | 0,50 - 0,70 |
+-------------------+----------------+---------------+--------------------+
Thảo luận kết quả
Nguyên nhân chính dẫn đến sự phân hóa không gian và thời gian của trường dòng chảy là tương tác động lực giữa ứng suất gió mùa luân phiên và gradient mật độ do bức xạ nhiệt mặt trời cùng nguồn nước lục địa đổ ra biển. Các bản đồ đẳng vận tốc qua mặt cắt thẳng đứng minh chứng rõ nét hiện tượng trượt tầng vận tốc: trong khi tầng mặt chịu chi phối mạnh của gió, các tầng sâu dưới 200m hoàn toàn bị điều khiển bởi gradient mật độ nằm ngang.
Khi so sánh với các nghiên cứu kinh điển của các tác giả quốc tế trước đây, kết quả tính toán động lực từ số liệu nhiệt muối của luận văn thể hiện sự tương thích cao về hướng dòng chảy tổng thể, đồng thời cung cấp chi tiết hơn về cấu trúc dị thường vận tốc tại các hẻm sâu ngầm. Dữ liệu này được mô tả sinh động qua các biểu đồ phân bố trường nhiệt muối hai chiều và đồ thị vectơ đẳng vận tốc, giúp làm sáng tỏ cơ chế hình thành các vùng nước trồi giàu dinh dưỡng tại vùng biển Nam Trung Bộ trong giai đoạn từ tháng 5 đến tháng 9 hàng năm.
Đề xuất và khuyến nghị
Thứ nhất, tích hợp ngay bộ ma trận vận tốc dòng chảy phân tầng tại 8 biên lỏng vào các mô hình số trị dự báo hải văn quốc gia nhằm nâng cao độ chính xác dự báo hoàn lưu 3D lên thêm 15% đến 20%, triển khai hoàn thiện trong vòng 12 tháng tới dưới sự chủ trì của Trung tâm Dự báo Khí tượng Thủy văn Quốc gia.
Thứ hai, thiết lập mạng lưới trạm quan trắc tự động kết hợp phao trôi ARGO chuyên sâu tại vùng biển tiếp giáp eo Bashi và eo Mindoro để cập nhật dữ liệu nhiệt muối thời gian thực, đặt mục tiêu thu thập chuỗi số liệu liên tục 24/7 trước năm 2028 do Viện Hải dương học phối hợp các đơn vị nghiên cứu biển thực hiện.
Thứ ba, ứng dụng ngay dữ liệu dòng chảy biên lỏng để xây dựng bản đồ mô phỏng rủi ro lan truyền sự cố tràn dầu và hóa chất độc hại trên các tuyến luồng hàng hải quốc tế qua Biển Đông, hoàn thành tài liệu hướng dẫn kỹ thuật trước quý IV năm 2026.
Thứ tư, chuẩn hóa cơ sở dữ liệu điều kiện biên hải dương học dùng chung cho các đơn vị tác chiến hải quân và cảnh sát biển, phục vụ công tác cứu hộ cứu nạn và tuần tra bảo vệ chủ quyền biển đảo với lộ trình thực hiện định kỳ 2 năm một lần.
Đối tượng nên tham khảo luận văn
Thứ nhất, các chuyên gia mô hình hóa thủy động lực học và dự báo viên khí tượng thủy văn biển tìm thấy phương pháp luận chi tiết để xử lý triệt để bài toán điều kiện biên mở trong các mô hình ROMS, Delft3D và POM.
Thứ hai, các nhà quản lý môi trường biển và cơ quan phòng chống thiên tai có được công cụ định lượng chính xác phục vụ đánh giá sức tải môi trường, dự báo vết dầu loang và phân tán bùn cát công trình ven bờ.
Thứ ba, lực lượng quốc phòng an ninh biển, đặc biệt là quân chủng hải quân, tiếp cận tài liệu tham khảo giá trị về đặc trưng trường âm dưới nước và dòng chảy ngầm phục vụ tác chiến tàu ngầm cùng các phương tiện ngầm tự hành.
Thứ tư, học viên cao học, nghiên cứu sinh và giảng viên chuyên ngành Hải dương học, Khí tượng học, Thủy văn học sử dụng luận văn làm tài liệu học thuật mẫu mực về phương pháp tính toán động lực học biển và phân tích trường mật độ thực nghiệm.
Câu hỏi thường gặp
Tại sao cần tính toán trao đổi nước qua các biên mở bằng phương pháp động lực?
Phương pháp động lực cho phép xác định chính xác cấu trúc dòng chảy phân tầng theo độ sâu từ dữ liệu nhiệt muối sẵn có mà không đòi hỏi chi phí đo đạc trực tiếp quá lớn. Phương pháp này cung cấp điều kiện biên 3D chuẩn xác, giúp khắc phục nhược điểm áp đặt biên tĩnh hay chỉ dùng mực nước trong các mô hình số trị, từ đó nâng cao độ chính xác dự báo toàn diện.
Eo biển Bashi đóng vai trò như thế nào đối với hoàn lưu Biển Đông?
Eo Bashi là eo biển sâu nhất với độ sâu đáy trên 2300m, giữ vai trò huyết mạch tiếp nhận trực tiếp khối nước Thái Bình Dương từ dòng Kuroshio vào Biển Đông. Lưu lượng qua eo này chi phối toàn bộ cấu trúc trường nhiệt muối sâu, điều khiển sự hình thành các xoáy nghịch quy mô lớn và cung cấp nguồn năng lượng chính duy trì hoàn lưu mùa đông của toàn vùng biển.
Dòng chảy ven bờ miền Trung có đặc điểm gì nổi bật trong mùa đông?
Vào mùa gió Đông Bắc, dòng chảy ven bờ miền Trung chịu tác động cộng hưởng giữa ứng suất gió mùa và hiệu ứng dồn ứ mật độ, đạt vận tốc rất lớn từ 0,75 đến 0,85 m/s và có thể vượt 1,0 m/s. Dòng hải lưu hẹp và mạnh này chảy dồn dập từ bắc xuống nam, vận chuyển nước lạnh có nhiệt độ dưới 24°C xuống sâu khu vực biển phía nam.
Làm thế nào để xác định mặt không động lực trong tính toán hải dương học?
Mặt không động lực là độ sâu chuẩn tại đó gradient áp suất nằm ngang triệt tiêu và vận tốc dòng chảy mật độ bằng 0 m/s. Trong đại dương sâu, mặt này thường được chọn ở độ sâu khoảng 1000m đến 1500m dựa trên phân tích độ đồng nhất của trường mật độ, trong khi tại các vùng biển nông và eo biển hẹp, mặt này được nội suy bám sát lớp ma sát đáy.
Kết quả của luận văn có thể tích hợp trực tiếp vào những mô hình số trị nào?
Dữ liệu lưu lượng và vận tốc dòng chảy phân tầng của nghiên cứu tương thích hoàn toàn để tích hợp làm điều kiện biên mở cho các mô hình thủy động lực 3D hiện đại như Delft3D, ROMS, MIKE 21/3 và FVCOM. Việc nạp bộ số liệu này giúp mô hình ổn định nhanh chóng và tái lập chính xác trường vận tốc ba chiều thực tế.
Kết luận
- Luận văn đã chuẩn hóa thành công quy trình ứng dụng phương pháp động lực để tính toán trường dòng chảy phân tầng qua 8 mặt cắt biên lỏng trọng yếu của Biển Đông dựa trên dữ liệu nhiệt muối thực nghiệm.
- Nghiên cứu định lượng chi tiết lưu lượng và vận tốc trao đổi nước qua các eo biển chiến lược gồm Bashi, Đài Loan, Mindoro, Balabac và Singapore trong hai mùa gió đặc trưng.
- Làm sáng tỏ quy luật hình thành hoàn lưu xoáy thuận quy mô lớn trong mùa đông với dòng phản xạ ven bờ miền Trung đạt trên 0,85 m/s và hệ thống xoáy nghịch phức tạp trong mùa hè.
- Cung cấp giải pháp khoa học thực tiễn để giải quyết bài toán thiếu hụt dữ liệu điều kiện biên mở theo phương thẳng đứng cho các mô hình số trị 3D tại Việt Nam.
- Đóng góp cơ sở dữ liệu hải dương học vật lý tin cậy phục vụ trực tiếp cho công tác dự báo khí tượng thủy văn biển, bảo vệ môi trường và bảo đảm an ninh quốc phòng.
Để khai thác tối đa giá trị khoa học của đề tài, các cơ quan nghiên cứu và đơn vị vận hành mô hình biển cần tiến hành tích hợp ngay bộ tham số biên lỏng này vào các hệ thống tính toán tác nghiệp trong giai đoạn 2026 - 2027. Hãy kết nối và chia sẻ rộng rãi công trình nghiên cứu này để thúc đẩy ứng dụng công nghệ mô phỏng hải văn tiên tiến tại Việt Nam.