Tổng quan về luận án

Trong kỷ nguyên phát triển của kinh tế hàng hải và sự gia tăng không ngừng về kích thước đội tàu thương mại quốc tế, công tác bảo đảm an toàn hàng hải tại các tuyến luồng tiếp cận cảng biển đối mặt với những thách thức khắt khe về mặt thủy động lực học. Tại Việt Nam, hệ thống luồng hàng hải Hải Phòng giữ vai trò huyết mạch của cửa ngõ quốc tế phía Bắc, tiếp nhận gần 20.000 lượt tàu biển mỗi năm (năm 2019) với mật độ lưu thông qua luồng Lạch Huyện đạt trung bình 70 lượt tàu/ngày. Tuy nhiên, luồng hàng hải Hải Phòng có cấu trúc hình học phức tạp, bao gồm luồng một chiều có ga tránh, bề ngang luồng hạn chế (chỉ từ 60 m đến 140 m), nhiều khúc cua gấp và chịu tác động mạnh mẽ của chế độ nhật triều thuần nhất với biên độ dao động lên tới 2,0 - 4,0 m cùng tốc độ biến đổi mực nước đạt 0,5 m/giờ.

Thực tiễn điều động tàu biển trên các tuyến luồng hẹp đòi hỏi hoa tiêu và thuyền trưởng phải tuân thủ nguyên tắc hàng hải căn bản: điều động tàu ngược dòng chảy triều nhằm duy trì tối đa khả năng kiểm soát hướng đi của tàu và góc bẻ lái. Dẫu vậy, khoảng trống nghiên cứu (research gap) lớn nhất hiện nay xuất phát từ sự hạn chế mang tính cố hữu của các công cụ tra cứu thủy văn truyền thống. Bảng thủy triều Việt Nam (gồm 3 tập do Trung tâm Hải văn xuất bản) cũng như Bịch thủy triều Anh (Admiralty Tide Tables - ATT) do Văn phòng Thủy đạc Hoàng gia Anh (UKHO) ấn hành vốn chỉ cung cấp dữ liệu rời rạc về thời gian và độ cao nước lớn/nước ròng. Trích dẫn nguyên văn luận án: "Bảng thuỷ triều chỉ cho biết thời gian và độ cao nước lớn và nước ròng trong một ngày tại một khu vực địa lý, nghĩa là mang tính tức thời. Vì vậy, việc dự đoán thuỷ triều tại các thời điểm khác cần phải tính toán nội suy. Việc này ít nhiều dẫn đến suy giảm độ chính xác, không có tính liên tục và không trực quan, gây khó khăn cho Hoa tiêu và Thuyền trưởng trong quá trình dẫn tàu, đặc biệt khi dẫn tàu hành trình trong tuyến luồng hàng hải".

Đồng thời, trên các hải đồ cảng (Habour charts) và hải đồ điện tử (ENC) hiện hành: "Trên các bình đồ hay còn gọi hải đồ cảng (Habour chart), kể cả hải đồ điện tử của khu vực hàng hải, đặc biệt là tuyến luồng hàng hải chưa thể hiện hướng và tốc độ dòng chảy mang tính liên tục".

Để giải quyết triệt để hạn chế này, luận án của NCS. Nguyễn Văn Cang đặt ra hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:

  • RQ1: Làm thế nào để mô hình hóa liên tục trường vận tốc và hướng dòng chảy triều theo thời gian thực trên toàn bộ tuyến luồng Hải Phòng với cấu trúc địa hình đáy phức tạp?
  • RQ2: Phương pháp xác lập điều kiện biên thủy động lực học nào có khả năng tích hợp chính xác dữ liệu từ trạm quan trắc thực tế, bảng thủy triều thiên văn và số liệu đo đạc thực địa?
  • RQ3: Mẫu bản đồ dòng chảy dạng giấy và số hóa cần được thiết kế và chuẩn hóa như thế nào để tích hợp tối ưu vào quy trình dẫn tàu an toàn của hoa tiêu?

Tương ứng với đó là các giả thuyết khoa học:

  • H1: Ứng dụng động lực học chất lưu tính toán (Computational Fluid Dynamics - CFD) giải hệ phương trình Navier-Stokes trung bình Reynolds (RANS) kết hợp mô hình dòng chảy rối $(k-\varepsilon)$ trên mô hình tuyến luồng phân đoạn cho phép tái tạo trường dòng triều chi tiết với sai số vận tốc dưới 10% so với thực địa.
  • H2: Việc thiết lập các hệ số hiệu chỉnh thực nghiệm ($k_{lh}, k_{sc}, k_{ar}$) kết hợp hệ số hiệu chỉnh bảng triều mới ($k_{xd}$) sẽ tối ưu hóa điều kiện biên đầu vào ($V_{dc}$), phản ánh chính xác tương tác thủy lực giữa các nhánh sông đan xen.
  • H3: Bản đồ dòng chảy trực quan hóa liên tục 24 giờ/ngày trong 365 ngày sẽ giúp loại bỏ hoàn toàn các sai số chủ quan do tính toán nội suy thủ công, làm giảm nguy cơ đâm va và mắc cạn tại các điểm đen hàng hải.

Khung lý thuyết của luận án được xây dựng dựa trên sự giao thoa giữa Thủy động lực học công trình thủy (Navier-Stokes Equations, RANS), Lý thuyết sóng nước nông (Shallow Water Wave Theory), Lý thuyết phân tích điều hòa thủy triều (Harmonic Tidal Analysis của Darwin và Doodson) và Động lực học điều động tàu (Maneuvering Hydrodynamics). Phạm vi không gian nghiên cứu bao quát 10 đoạn luồng trọng yếu thuộc tuyến luồng Hải Phòng từ km 44+400 đến km 14+000 (tổng chiều dài nghiên cứu 44,4 km), kết hợp dữ liệu chuỗi thời gian từ 3 trạm quan trắc tự động (Arroyo, Sông Cấm, Lạch Huyện) và khảo sát thực nghiệm trực tiếp bằng máy đo dòng chảy XZ-2 trong năm 2019 - 2020. Đóng góp đột phá của công trình là việc lần đầu tiên tại Việt Nam thiết lập thành công quy trình công nghệ tính toán và xuất bản bộ bản đồ dòng chảy tuyến luồng hàng hải chi tiết 365 ngày dạng giấy và số hóa.

Literature Review và Positioning

Lịch sử nghiên cứu thủy triều và dòng chảy biển trên thế giới đã đạt được những bước tiến vượt bậc, dẫn đầu bởi các cơ quan thủy đạc quốc tế như Văn phòng Thủy đạc Hoàng gia Anh (UKHO), Cơ quan Quản trị Khí quyển và Đại dương Quốc gia Hoa Kỳ (NOAA) và Cơ quan Thủy đạc Hải quân Pháp (SHOM). Các ấn phẩm hàng hải kinh điển như China Sea Pilot (NP 30) hay hệ thống hải đồ chuẩn BA charts (№ 1965, 3875, 3881, 3882) của UKHO cung cấp hệ số điều hòa, độ cao triều và các mũi tên dòng chảy mang tính đại diện tĩnh tại một số điểm đo cô lập. Đối với mô phỏng số dòng chảy quy mô lớn vùng Vịnh Bắc Bộ, các nghiên cứu quốc tế tiêu biểu như Manh et al. (1998, 2000) với mô hình số trị ba chiều (3D numerical model of tide and tidal current in the Gulf of Tongking), Zu et al. (2008), hay Nguyen et al. (2014) (Tidal characteristics of the Gulf of Tonkin, Continental Shelf Research) đã giải quyết thành công bài toán hoàn lưu đại dương, dòng chảy dư và sóng triều quy mô khu vực. Tuy nhiên, các mô hình quy mô vùng biển mở (oceanic scale) sử dụng lưới phân giải thô (hàng trăm mét đến vài kilômét), không thể tái hiện được đặc tính thủy động lực học vi mô trong lòng các tuyến luồng hàng hải hạn chế có biên bờ và địa hình đáy thay đổi đột ngột.

Tại Việt Nam, các hướng nghiên cứu về hải dương học và cơ học thủy khí phát triển qua các công trình của GS. Lương Công Nhớ, PGS. Phạm Kỳ Quang, PGS. Nguyễn Cảnh Sơn, TS. Bùi Văn Cường (Trường ĐH Hàng hải Việt Nam) về động lực học tàu và xâm thực chân vịt; PGS. Nguyễn Xuân Phương, PGS. Nguyễn Phùng Hưng (Trường ĐH Giao thông vận tải TP.HCM) về dòng chảy bao quanh thân tàu và mô phỏng CFD vùng neo đậu; GS. Vũ Duy Quang, PGS. Trương Việt Anh (Trường ĐH Bách khoa Hà Nội) về thủy khí động lực học cánh dẫn; GS. Trần Tân Tiến, PGS. Phạm Văn Huấn (Trường ĐH Khoa học Tự nhiên - ĐHQG Hà Nội) về hải dương học tính toán; và các công trình của Trần Văn Chung, Bùi Hồng Long (2007, 2015) trên Tạp chí Khoa học và Công nghệ Biển ứng dụng mô hình 3D phi tuyến mô phỏng dòng triều Vịnh Bắc Bộ.

Trong y văn tồn tại cuộc tranh luận học thuật sâu sắc giữa hai trường phái tiếp cận:

  • Trường phái 1 (Mô hình tích hợp 2D vùng rộng): Sử dụng các bộ công cụ thương mại giải phương trình nước nông 2 chiều tích phân theo chiều sâu (như MIKE 21, Delft3D, TELEMAC-2D). Ưu điểm là tính toán nhanh trên diện rộng nhưng nhược điểm lớn là bỏ qua cấu trúc biến thiên vận tốc theo phương thẳng đứng (profil vận tốc dòng chảy), đánh giá thấp hiệu ứng ma sát biên thành luồng hẹp và xoáy cục bộ tại các ngã ba sông.
  • Trường phái 2 (Mô phỏng CFD 3D trường dòng vi mô): Ứng dụng phương pháp thể tích hữu hạn (FVM) với các phần mềm chuyên sâu như ANSYS Fluent, Flow-3D để mô phỏng chi tiết dòng chảy rối 3 chiều. Trường phái này thường bị chỉ trích là tiêu tốn tài nguyên tính toán và gặp khó khăn khi mở rộng miền biên ngoài cửa biển.

Luận án của NCS. Nguyễn Văn Cang định vị chính xác tại điểm giao thoa giữa hai trường phái: kế thừa dữ liệu mực nước thiên văn vĩ mô từ bảng thủy triều và trạm quan trắc, phân đoạn thủy lực tuyến luồng (fairway hydraulic segmentation) thành 10 module độc lập, sau đó ứng dụng công cụ CFD ANSYS Fluent để giải bài toán thủy động lực học cục bộ với độ phân giải lưới cực cao. So với các nghiên cứu quốc tế của Manh et al. (1998) hay nghiên cứu trong nước của Trần Văn Chung và Bùi Hồng Long (2015), luận án tạo ra bước tiến vượt bậc khi chuyển hóa kết quả mô phỏng số học thuần túy thành sản phẩm hải đồ chuyên dụng phục vụ trực tiếp cho tác nghiệp hàng hải của hoa tiêu và thuyền trưởng theo thời gian thực.

Đó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 phạm vi ứng dụng của Lý thuyết lớp biên thủy động lực học (Hydrodynamic Boundary Layer Theory) và Phương trình dòng chảy rối Navier-Stokes trung bình Reynolds (RANS) vào miền không gian luồng hàng hải cửa sông ven biển chịu chi phối của nhật triều. Công trình đóng góp hệ thống mệnh đề lý thuyết (Theoretical Propositions) liên kết trực tiếp động lực học dòng chảy với động lực học chuyển động tàu:

  • Mệnh đề P1: Trong luồng hàng hải hạn chế, góc dạt của tàu ($\beta$) tỷ lệ thuận với thành phần vận tốc dòng chảy vuông góc với trục luồng ($V_{dc} \cdot \sin\theta$) và tỷ lệ nghịch với vận tốc thật của tàu ($V_t$), được mô tả qua phương trình: $$\beta = \arcsin\left(\frac{V_{dc} \cdot \sin(\theta - HT)}{V_{HTT}}\right)$$ trong đó hướng đi thực tế của tàu dưới tác động của dòng chảy được xác định bởi $HTT_\beta = HT \pm \beta$.
  • Mệnh đề P2: Sự phân bố năng lượng dòng chảy dọc theo tuyến luồng tuân theo định luật bảo toàn năng lượng Bernoulli mở rộng cho chất lưu thực có tổn thất ma sát đáy luồng (bed shear stress) và tổn thất cục bộ tại các đoạn uốn cong: $$z_1 + \frac{p_1}{\rho g} + \frac{\alpha_1 V_1^2}{2g} = z_2 + \frac{p_2}{\rho g} + \frac{\alpha_2 V_2^2}{2g} + h_f + h_j$$
  • Mệnh đề P3: Chuyển đổi trạng thái dòng chảy từ chảy êm sang chảy xiết được kiểm soát thông qua tỷ số không thứ nguyên số Froude ($Fr = \frac{V}{\sqrt{g \cdot h_{tb}}}$) và số Reynolds ($Re = \frac{\rho V R_h}{\mu}$), quyết định cấu trúc xoáy thủy lực tại các vùng hợp lưu.

Luận án đánh dấu một bước chuyển dịch mô thức (paradigm shift) từ phương pháp ước lượng bán thực nghiệm tĩnh dựa trên nội suy bảng số sang mô hình động lực học trường vector dòng chảy liên tục (continuous vector field dynamics).

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

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

  1. Lý thuyết Thủy động lực học tính toán (CFD): Giải hệ phương trình RANS kết hợp mô hình dòng chảy rối $(k-\varepsilon)$ tiêu chuẩn (Standard $k-\varepsilon$ model) gồm phương trình động năng rối $k$ và tỷ lệ tiêu tán năng lượng rối $\varepsilon$.
  2. Lý thuyết Điều hòa Thủy triều (Tidal Harmonic Analysis): Tích hợp 4 thành phần triều chính ($M_2, S_2, K_1, O_1$) để xác định dao động mực nước biên.
  3. Lý thuyết Động lực học Tàu thủy (Ship Hydrodynamics): Xác định lực đẩy và mômen dạt tác động lên vỏ tàu do phân bố vận tốc dòng chảy không đồng đều.

Trích dẫn nguyên văn đóng góp mới của luận án: "Đưa ra phương pháp luận và tính toán hệ số hiệu chỉnh bảng thủy triều mới xây dựng (kxd). Đặc biệt phương pháp xác định điều kiện biên là vận tốc dòng chảy (Vdc) đã tính đến hai hệ số hiệu chỉnh thực nghiệm klh và kdc".

Cơ chế xác lập điều kiện biên vận tốc $V_{dc}$ được chuẩn hóa qua hàm toán học: $$V_{dc}(t) = f(h_b(t), k_{xd}, k_{station}) = \frac{\Delta h_b(t)}{\Delta t} \cdot \frac{A_{surface}}{A_{cross}} \cdot k_{xd} \cdot k_{station}$$ trong đó $k_{station} \in {k_{ar}, k_{sc}, k_{lh}}$ là hệ số hiệu chỉnh thực nghiệm tại các trạm Arroyo, Sông Cấm, Lạch Huyện; $k_{xd}$ là hệ số hiệu chỉnh bảng thủy triều mới; $A_{surface}$ và $A_{cross}$ lần lượt là diện tích mặt thoáng và diện tích mặt cắt ướt của đoạn luồng khảo sát.

Điều kiện biên biên giới hạn (boundary conditions) được xác lập rõ ràng:

  • Biên vào (Inlet): Vận tốc dòng chảy $V_{dc}$ xác định theo giờ dựa trên độ dốc đường cong triều dâng/triều rút.
  • Biên ra (Outlet): Áp suất thủy tĩnh tự do (Pressure-outlet) tương ứng với mực nước tức thời tại đầu hạ lưu.
  • Biên thành bờ và đáy luồng (Wall): Điều kiện không trượt (No-slip condition) kết hợp hàm vách chuẩn (Standard Wall Functions) với độ nhám đáy tương đương $k_s$.
  • Mặt thoáng tự do: Biên đối xứng (Symmetry) hoặc mô hình thể tích chất lưu (VOF - Volume of Fluid) đối với các pha nước - không khí.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên nền tảng triết lý thực chứng (positivism) kết hợp kiểm chứng thực nghiệm chặt chẽ, dựa trên nguyên lý xác định luận của vật lý thủy động. Thiết kế nghiên cứu đa cấp độ (multi-level design) được triển khai qua việc phân chia toàn bộ tuyến luồng Hải Phòng thành 10 đoạn luồng tính toán độc lập nhằm đảm bảo tính đồng nhất về mặt hình học thủy lực và kiểm soát tối ưu các điều kiện biên:

  • Đoạn 1 (km 44+400 đến km 43+200): Dài 1,2 km; Chiều rộng mặt thoáng/luồng: 310/60 m; Độ sâu trung bình: 7,5 m.
  • Đoạn 2 (km 43+200 đến km 40+200): Dài 3,0 km; Chiều rộng mặt thoáng/luồng: 315/80 m; Độ sâu trung bình: 8,0 m (riêng khu vực Phà Bính sâu 13,3 m).
  • Đoạn 3 (km 40+200 đến km 38+800): Dài 1,4 km; Chiều rộng mặt thoáng/luồng: 300/120 m; Độ sâu trung bình: 8,1 m (khu vực ngã ba Sông Cấm - Sông Ruột Lợn, trạm Sông Cấm).
  • Đoạn 4 (km 38+800 đến km 33+800): Dài 5,0 km; Chiều rộng mặt thoáng/luồng: 300/80 m; Độ sâu trung bình: 7,4 m.
  • Đoạn 5 (km 33+800 đến km 30+200): Dài 3,6 km; Chiều rộng mặt thoáng/luồng: 315/80 m; Độ sâu trung bình: 9,1 m (khu vực ngã ba Sông Cấm - Sông Bạch Đằng, trạm Arroyo).
  • Đoạn 6 (km 30+200 đến km 26+600): Dài 3,6 km; Chiều rộng mặt thoáng/luồng: 600/80 m; Độ sâu trung bình: 6,5 m.
  • Đoạn 7 (km 26+600 đến km 24+600): Dài 2,0 km; Chiều rộng mặt thoáng/luồng: 1200/140 m; Độ sâu trung bình: 7,1 m.
  • Đoạn 8 (km 24+600 đến km 20+200): Dài 4,4 km; Chiều rộng mặt thoáng/luồng: 242/80 m; Độ sâu trung bình: 6,9 m (giao tuyến Sông Cấm, Nam Triệu, Cái Tráp, kênh Hà Nam; trạm Lạch Huyện).
  • Đoạn 9 (km 20+200 đến km 15+800): Dài 4,4 km; Chiều rộng mặt thoáng/luồng: 1500/100 m; Độ sâu trung bình: 8,5 m (kênh Hà Nam - Cái Tráp, luồng Lạch Huyện, cặp phao 19 - 21 - 23).
  • Đoạn 10 (km 15+800 đến km 14+000): Dài 1,8 km; Chiều rộng mặt thoáng/luồng: 1500/100 m; Độ sâu trung bình: 11,0 m (cửa biển tiếp cận luồng Lạch Huyện).

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được chuẩn hóa qua 5 giai đoạn liên hoàn:

  1. Thu thập dữ liệu trắc địa và hải văn: Số hóa bình đồ độ sâu chuẩn VN4HP003, VN4HP005, VN4HP007, VN4HP008 (Tổng công ty Bảo đảm An toàn Hàng hải Miền Bắc) và bình đồ I-25-022, I-25-031 (Bộ Tư lệnh Hải quân) theo chuẩn tọa độ WGS-84.
  2. Khai thác chuỗi dữ liệu trạm quan trắc tự động: Trích xuất dữ liệu mực nước tức thời với tần suất quét 10 phút/lần tại 3 trạm:
    • Trạm Arroyo: $\varphi = 20^\circ 49'54''\text{N}; \lambda = 106^\circ 49'16''\text{E}$ (hoạt động từ 10/2013).
    • Trạm Lạch Huyện: $\varphi = 20^\circ 49'19''\text{N}; \lambda = 106^\circ 54'28''\text{E}$ (hoạt động từ 09/2018).
    • Trạm Sông Cấm: $\varphi = 20^\circ 52'43''\text{N}; \lambda = 106^\circ 42'16''\text{E}$ (hoạt động từ 09/2018).
  3. Khảo sát thực địa: Tiến hành đo đạc trực tiếp vận tốc thực tế ($V_{tt}$) và hướng dòng chảy tại các tọa độ tim luồng và mép luồng bằng thiết bị đo dòng chảy chuyên dụng XZ-2.
  4. Mô phỏng số học CFD: Xây dựng hình học 3D, chia lưới phi cấu trúc (unstructured mesh) với phần tử tam giác và tứ diện kết hợp lớp biên lăng trụ (prism boundary layers), thiết lập điều kiện biên giải bằng ANSYS Fluent.
  5. Tam giác đạc và kiểm định: Đối chuẩn chéo 3 nguồn dữ liệu (Bảng triều lý thuyết - Trạm quan trắc tự động - Đo đạc thực địa XZ-2).

Data và phân tích

Mô hình số trị sử dụng thuật toán phân tách áp suất - vận tốc SIMPLE (Semi-Implicit Method for Pressure Linked Equations) với sơ đồ sai phân ngược bậc hai (Second Order Upwind) cho các phương trình động lượng, động năng rối ($k$) và tỷ lệ tiêu tán ($\varepsilon$). Tiêu chuẩn hội tụ (residual convergence criteria) được đặt ở mức $10^{-5}$ đối với tất cả các biến số liên tục và động lượng. Phân tích độ nhạy lưới (mesh independence test) được thực hiện với 3 mức mật độ lưới (thô, trung bình, mịn) để loại bỏ sai số do chia lưới. Sai số tương đối giữa vận tốc tính toán mô phỏng ($V_{dc}$) và vận tốc đo thực tế ($V_{tt}$) được đánh giá qua tỷ lệ phần trăm chênh lệch: $$\text{Error} (%) = \frac{|V_{tt} - V_{dc}|}{V_{tt}} \times 100%$$

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

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

  1. Xác thực độ chính xác của mô hình CFD: Kết quả so sánh giữa vận tốc mô phỏng ($V_{dc}$) và vận tốc thực nghiệm ($V_{tt}$) đo bằng máy XZ-2 trong ngày 01/01/2020 trên các đoạn luồng cho thấy độ tương đồng cao. Sai số trung bình duy trì dưới 8,5%, chứng minh tính đúng đắn của việc lựa chọn mô hình dòng rối $(k-\varepsilon)$ và phương pháp xác lập điều kiện biên vận tốc qua hệ số hiệu chỉnh $k_{xd}$ và $k_{lh}$.
  2. Phát hiện dòng chảy ngang nguy hiểm tại Đoạn luồng số 5 (Ngã ba Sông Cấm - Sông Bạch Đằng): Dữ liệu mô phỏng và đo đạc chứng minh rằng khi triều dâng, vận tốc dòng chảy đạt từ 2,0 knots đến gần 3,0 knots (~1,0 - 1,5 m/s) với vector dòng chảy tạo góc xiên lớn so với tim luồng, luôn tạo ra lực dạt đẩy thân tàu mạnh về phía bờ Thủy Nguyên. Đây là nguyên nhân trực tiếp lý giải các vụ va chạm và mắc cạn tại khu vực này trong thống kê tai nạn hàng hải giai đoạn 2015 - 2019.
  3. Xác định điểm nút động lực học phức tạp tại Đoạn luồng số 8 (Khu vực Bạch Đằng - Hà Nam - Nam Triệu): Đây là giao điểm của 4 nhánh thủy đạo (Sông Cấm, Nam Triệu, Cái Tráp, kênh Hà Nam). Tại đây xuất hiện hiện tượng xoáy thủy lực cục bộ và phân kỳ dòng chảy với tốc độ dòng chảy thay đổi đột ngột giữa hai mép luồng (chênh lệch vận tốc ngang luồng lên tới 0,4 m/s), đe dọa nghiêm trọng tính ổn định hướng đi của tàu khi hành trình qua phao 30 - 32.
  4. Đặc trưng dòng chảy mạnh tại Đoạn luồng số 9 (Khu vực kênh Hà Nam - Lạch Huyện): Tốc độ dòng triều dâng và rút tại khu vực các cặp phao 19 - 21 - 23 đạt xấp xỉ 3,0 knots. Tại thời điểm triều rút cường, dòng chảy xiết kết hợp mật độ phương tiện thủy nội địa cắt ngang luồng làm tăng vọt nguy cơ đâm va đối với các tàu container trọng tải lớn tiếp cận cảng Lạch Huyện.
  5. Quy luật phi tuyến tính giữa đỉnh triều và vận tốc dòng chảy: Nghiên cứu chỉ ra kết quả phi trực giác: thời điểm vận tốc dòng chảy đạt giá trị cực đại ($V_{max}$) không trùng với thời điểm nước lớn nhất (High Water) trên bảng triều, mà xuất hiện trước đó từ 1,5 đến 2,5 giờ (thời kỳ triều dâng nhanh nhất với tốc độ biến đổi mực nước $\Delta h / \Delta t$ đạt cực đại).

Implications đa chiều

  • Về mặt lý thuyết: Mở rộng cơ sở dữ liệu thủy động lực học luồng hàng hải hạn chế vùng nhiệt đới gió mùa chịu tác động nhật triều; cung cấp bằng chứng thực nghiệm củng cố lý thuyết điều khiển học tàu biển trong điều kiện dòng chảy xiên.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn 5 bước ứng dụng CFD mô phỏng dòng triều theo module phân đoạn, có khả năng nhân rộng cho tất cả các tuyến luồng hàng hải khác tại Việt Nam (như luồng Soài Rạp, luồng Vũng Tàu - Thị Vải, luồng Quy Nhơn, luồng Đà Nẵng).
  • Về mặt thực tiễn hàng hải: Cung cấp công cụ trực quan hóa cho hoa tiêu và thuyền trưởng, cho phép tra cứu chính xác hướng và tốc độ dòng chảy theo từng giờ trong 365 ngày của năm, tối ưu hóa thời điểm điều động tàu cập cầu ngược nước an toàn, giảm thời gian chờ đợi con nước và loại bỏ sự phụ thuộc vào trực giác cá nhân.
  • Về mặt chính sách và quản lý: Cung cấp cơ sở khoa học cho Cục Hàng hải Việt Nam, Cảng vụ Hàng hải Hải Phòng và Tổng công ty Bảo đảm An toàn Hàng hải Miền Bắc trong việc quy hoạch lại hệ thống báo hiệu hàng hải, phân luồng điều tiết giao thông và cấp phép mớn nước cho tàu biển hành trình qua luồng.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận 4 giới hạn khoa học:

  1. Bỏ qua các hiện tượng khí tượng thủy văn cực đoan: Mô hình chưa tính toán đến sự biến đổi bất thường do bão lớn, lũ lụt từ thượng nguồn sông Hồng - sông Thái Bình, sóng biển quy mô lớn do gió mùa Đông Bắc cực mạnh gây nước dâng do bão (storm surges).
  2. Giới hạn tính toán 3D toàn tuyến: Do giới hạn về tài nguyên tính toán, mô hình CFD tập trung mô phỏng chi tiết 3D cho các đoạn luồng xung yếu (đoạn 3, 5, 8, 9), trong khi các đoạn luồng thẳng, ổn định được chuẩn hóa theo điều kiện biên giả định 2D/3D đồng nhất.
  3. Sự cố gián đoạn số liệu trạm đo: Chuỗi dữ liệu trạm Arroyo bị gián đoạn từ tháng 5/2018 do tai nạn hàng hải và trạm Sông Cấm bị lỗi thiết bị từ tháng 10/2019, dẫn đến việc phải sử dụng các thuật toán nội suy và hệ số tương quan để bù đắp dữ liệu biên.
  4. Biến động địa hình do sa bồi: Địa hình đáy luồng được cố định tại thời điểm đo đạc chuẩn hóa (tháng 9/2020), chưa tích hợp mô hình bồi lắng - xói lở động lực hình thái học đáy (morphodynamic sediment transport).

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Phát triển mô hình ghép nối thời gian thực (Real-time Coupled Hydrodynamic Model) tích hợp dữ liệu khí tượng - thủy triều - dòng chảy từ hệ thống IoT quan trắc trực tiếp.
  • Mở rộng thuật toán mô phỏng 3D tương tác thân tàu - dòng chảy - đáy luồng (Ship-Bank-Current Interaction) ứng dụng lưới chuyển động (Overset Grid / Dynamic Mesh).
  • Ứng dụng trí tuệ nhân tạo (Machine Learning / Deep Neural Networks) dự báo tức thời trường dòng chảy dựa trên dữ liệu lớn quan trắc khí quyển.
  • Tích hợp lớp bản đồ dòng chảy động (Dynamic Current Layer) trực tiếp vào hệ thống thông tin và hiển thị hải đồ điện tử (ECDIS) và dịch vụ điều phối tàu thuyền (VTS).

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

  • Tác động học thuật: Đặt nền móng khoa học cho chuyên ngành Khoa học Hàng hải (Mã số: 9840106) tại Việt Nam trong lĩnh vực số hóa thủy văn luồng tàu; công bố 12 công trình khoa học chuyên ngành, trong đó có 01 sách chuyên khảo chuẩn mực (Ứng dụng CFD xây dựng bản đồ dòng chảy tuyến luồng hàng hải, NXB Hàng hải, ISBN: 978-604-937-217-9) và chủ trì/tham gia các đề tài cấp Bộ GTVT (Mã số: DT183035).
  • Chuyển đổi công nghiệp hàng hải: Trực tiếp hỗ trợ hơn 50 hoa tiêu thuộc Công ty Hoa tiêu hàng hải khu vực II và hàng trăm thuyền trưởng điều động tàu an toàn ra vào cụm cảng Hải Phòng và Lạch Huyện; giảm thiểu rủi ro tai nạn hàng hải (vốn ghi nhận trung bình 5-11 vụ/năm giai đoạn 2015 - 2019), tiết kiệm chi phí thuê tàu lai dắt và giảm phát thải khí nhà kính nhờ tối ưu hóa tốc độ hành trình.
  • Tác động chính sách: Cung cấp luận cứ kỹ thuật cho Bộ Giao thông vận tải ban hành các quy chuẩn kỹ thuật quốc gia về ấn phẩm hàng hải điện tử và tiêu chuẩn an toàn luồng tàu hẹp.
  • Lợi ích xã hội và tầm vóc quốc tế: Góp phần bảo vệ môi trường biển chống nguy cơ tràn dầu do đâm va, nâng cao năng lực cạnh tranh của cụm cảng cửa ngõ quốc tế Hải Phòng trong chuỗi logistics toàn cầu.

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

  • Nghiên cứu sinh và học viên sau đại học: Tiếp cận phương pháp luận mẫu mực về kết hợp mô phỏng số CFD với thực nghiệm hiện trường trong lĩnh vực hàng hải và kỹ thuật công trình thủy.
  • Các nhà khoa học và giảng viên chuyên ngành: Kế thừa bộ khung lý thuyết và dữ liệu thủy động lực học phục vụ giảng dạy các môn học Điều động tàu, Thủy văn hàng hải và Thủy khí động lực học.
  • Hoa tiêu hàng hải và Thuyền trưởng: Sở hữu công cụ tra cứu trực quan hướng và tốc độ dòng chảy 24/7, giúp lập kế hoạch chuyến đi (Passage Planning) và điều động cập/rời cầu cảng chuẩn xác tuyệt đối.
  • Các cơ quan quản lý nhà nước (Cục Hàng hải, Cảng vụ, VMS North): Có công cụ hỗ trợ ra quyết định phân luồng, điều tiết giao thông hàng hải, phê duyệt phương án dẫn tàu chở hàng siêu trường, siêu trọng.

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ì và đã mở rộng lý thuyết nào?
Luận án mở rộng Lý thuyết lớp biên thủy động lực học và mô hình dòng rối RANS $(k-\varepsilon)$ cho môi trường luồng hàng hải tự nhiên cửa sông chịu chi phối của chế độ nhật triều biên độ lớn. Đóng góp độc đáo nhất là xây dựng thành công hàm liên kết giải tích giữa độ dốc mực nước thiên văn $\Delta h_b / \Delta t$ với vector vận tốc dòng chảy biên $V_{dc}$ thông qua cặp hệ số hiệu chỉnh thực nghiệm ($k_{xd}, k_{station}$), giải quyết bài toán thiếu hụt dữ liệu dòng chảy liên tục mà các mô hình lý thuyết trước đây chưa lượng hóa được.

2. Đổi mới phương pháp luận của nghiên cứu khi so sánh với ít nhất 2 công trình quốc tế tiêu biểu?
So với mô hình số trị 3D vùng Vịnh Bắc Bộ của Manh et al. (1998, 2000) vốn sử dụng lưới hữu hạn quy mô đại dương không thể phản ánh cấu trúc luồng hẹp, và so với hệ thống hải đồ tĩnh truyền thống của Văn phòng Thủy đạc Anh (UKHO NP 30 / BA charts) chỉ cung cấp thông tin điểm rời rạc, luận án đã đổi mới vượt bậc bằng kỹ thuật "Phân đoạn thủy lực 10 module kết hợp chia lưới mịn cục bộ trong ANSYS Fluent". Phương pháp này vừa kiểm soát chính xác điều kiện biên thủy lực từng đoạn, vừa tái tạo hoàn hảo trường dòng chảy 3 chiều tại các khúc cua hiểm trở.

3. Phát hiện thực nghiệm gây bất ngờ nhất có sự hỗ trợ của dữ liệu là gì?
Phát hiện bất ngờ nhất là tại Đoạn luồng số 5 và Đoạn luồng số 8, đỉnh vận tốc dòng chảy ($V_{max} \approx 3,0\text{ knots}$) xuất hiện sớm hơn đỉnh triều thiên văn từ 1,5 đến 2,5 giờ. Hiện tượng này xảy ra do hiệu ứng co thắt thủy lực cửa sông (estuarine hydraulic contraction) khi lưu lượng nước tràn qua mặt cắt hẹp đạt cực đại trong giai đoạn triều dâng nhanh nhất, thay vì tại thời điểm nước đứng đỉnh triều như quan niệm phổ biến của nhiều hoa tiêu.

4. Giao thức tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?
Luận án cung cấp tường minh 100% giao thức tái lập nghiên cứu: toàn bộ thông số hình học đáy của 10 đoạn luồng (bảng 1.4), tọa độ địa lý chuẩn WGS-84 của 3 trạm quan trắc (Arroyo, Sông Cấm, Lạch Huyện), cấu hình chi tiết phần mềm ANSYS Fluent (mô hình $(k-\varepsilon)$, sơ đồ SIMPLE, điều kiện biên áp suất/vận tốc), thông số kỹ thuật máy đo XZ-2 và bảng tổng hợp thuật toán tính toán hệ số hiệu chỉnh $k_{xd}$ tại 6 phụ lục chi tiết.

5. Chương trình nghiên cứu 10 năm tiếp theo được phác thảo như thế nào?
Chương trình 10 năm tập trung vào 3 trụ cột: (1) Số hóa toàn diện và tích hợp bản đồ dòng chảy động 3D vào hệ thống Hải đồ điện tử quốc gia S-100/S-111 (Dynamic Water Current Standard); (2) Xây dựng hệ thống bản sao số (Digital Twin) mô phỏng tương tác thời gian thực giữa tàu biển và luồng hàng hải Hải Phòng; (3) Mở rộng xây dựng bản đồ dòng chảy cho toàn bộ 45 tuyến luồng hàng hải trên phạm vi toàn quốc.

Kết luận

  1. Chuẩn hóa cơ sở khoa học phân đoạn luồng: Thiết lập thành công mô hình phân đoạn thủy lực 10 đoạn luồng cho 44,4 km tuyến luồng Hải Phòng bám sát điều kiện biên thực địa, chỉ rõ các đặc trưng thủy động lực học và các điểm đen mất an toàn hàng hải tại các đoạn luồng số 3, 5, 8 và 9.
  2. Quy trình công nghệ mô phỏng CFD tiên phong: Xây dựng hoàn chỉnh quy trình ứng dụng ANSYS Fluent giải hệ phương trình RANS kết hợp mô hình dòng rối $(k-\varepsilon)$ để tính toán trường vector dòng triều chi tiết với sai số kiểm định thực nghiệm dưới 8,5%.
  3. Phương pháp luận xác định điều kiện biên đột phá: Đề xuất và chứng minh thành công phương pháp xác định vận tốc dòng chảy biên ($V_{dc}$) thông qua hệ số hiệu chỉnh bảng thủy triều mới ($k_{xd}$) và các hệ số hiệu chỉnh trạm quan trắc ($k_{lh}, k_{sc}, k_{ar}$).
  4. Xuất bản bộ sản phẩm bản đồ dòng chảy 365 ngày: Xây dựng thử nghiệm thành công bộ bản đồ dòng chảy tuyến luồng Hải Phòng dạng giấy và dạng số hóa hiển thị liên tục hướng và tốc độ dòng chảy theo từng giờ trong ngày cho 365 ngày trong năm, khắc phục hoàn toàn nhược điểm gián đoạn của bảng thủy triều truyền thống.
  5. Mở ra 3 hướng nghiên cứu học thuật mới: Mở đường cho các nghiên cứu về tích hợp dữ liệu dòng chảy động vào chuẩn hải đồ điện tử tương lai S-111, mô phỏng số tương tác thủy động đa vật thể (tàu - dòng chảy - bờ luồng) và dự báo thủy văn hàng hải thời gian thực ứng dụng trí tuệ nhân tạo.
  6. Di sản thực tiễn lâu dài: Đóng góp trực tiếp vào mục tiêu nâng cao năng lực an toàn hàng hải quốc gia, phòng ngừa tai nạn đâm va, mắc cạn, bảo vệ luồng hàng hải và thúc đẩy phát triển kinh tế biển bền vững tại cụm cảng cửa ngõ quốc tế Hải Phòng.