Tổng quan về luận án

Luận án tiến sĩ kỹ thuật chuyên ngành Khoa học Hàng hải (Mã số: 9840106) của Nghiên cứu sinh Lê Văn Thức, dưới sự hướng dẫn khoa học của PGS.TS. Nguyễn Văn Thư và PGS.TS. Lê Quốc Tiến tại Trường Đại học Giao thông vận tải Thành phố Hồ Chí Minh, mang tựa đề: "Xây dựng Mô hình đánh giá rủi ro tai nạn đâm va tại cảng biển Vũng Tàu, ứng dụng phòng ngừa tai nạn và hỗ trợ công tác cải tạo, nâng cấp luồng hàng hải". Nghiên cứu đặt trọng tâm vào việc giải quyết bài toán an toàn hàng hải tại vùng nước Vịnh Gành Rái và luồng hàng hải Vũng Tàu – Thị Vải, khu vực được phân loại là cảng biển đặc biệt, cửa ngõ trung chuyển quốc tế theo Quyết định số 1579/QĐ-TTg ngày 22/9/2021 của Thủ tướng Chính phủ, có khả năng tiếp nhận tàu container trọng tải đến 80.000 DWT và tàu tổng hợp đến 150.000 DWT.

+-----------------------------------------------------------------------------------+
|               KHUNG ĐÁNH GIÁ RỦI RO ĐÂM VA TÀU THUYỀN TẠI CẢNG BIỂN               |
+-----------------------------------------------------------------------------------+
|  [Dữ liệu Đầu vào]          [Mô hình Phân tích & Trọng số]     [Công cụ Mô phỏng]  |
|  - AIS, Thủy văn (2000-17)  -->  - Hồi quy Đa biến (SPSS)  -->  - Fast-Time (C++) |
|  - Khảo sát (258 Chuyên gia) -->  - AHP (Saaty, CR < 0.1)    -->  - Sugeno (MATLAB)|
|  - 11 Kịch bản VTS Thực tế  -->  - 27 Luật Logic Mờ         -->  - Bản đồ Rủi ro K |
+-----------------------------------------------------------------------------------+

Thực tiễn thống kê cho thấy cảng biển Vũng Tàu có mật độ giao thông thủy cao và phức tạp bậc nhất Việt Nam, tiếp nhận trung bình hơn 30.000 lượt tàu biển và trên 70.000 lượt phương tiện thủy nội địa, tàu cao tốc, tàu cá hoạt động đan xen mỗi năm với tốc độ tăng trưởng lưu lượng đạt 5–7%/năm. Theo Quy định 13, Chương V của Công ước Quốc tế về An toàn Sinh mạng Con người trên Biển (SOLAS 74), các Chính phủ thành viên có nghĩa vụ kiểm soát mức độ rủi ro và cung cấp các biện pháp hỗ trợ hành hải (Aids to Navigation) tương thích với mật độ giao thông. Tuy nhiên, rào cản khoa học lớn nhất tồn tại ở Việt Nam là sự thiếu vắng một công cụ định lượng chuẩn hóa để đánh giá rủi ro đâm va; công tác điều hành của Trung tâm Quản lý Hàng hải (VTS) và Cảng vụ Hàng hải chủ yếu dựa trên đánh giá định tính và kinh nghiệm chủ quan.

Luận án thiết lập 3 câu hỏi nghiên cứu và 3 giả thuyết cốt lõi:

  1. Câu hỏi nghiên cứu 1 (RQ1): Những yếu tố nào đóng vai trò chi phối đến nguy cơ đâm va của tàu thuyền trong vùng nước cảng biển Vũng Tàu?
    • Giả thuyết 1 (H1): Rủi ro đâm va là hàm đa biến đồng thời chịu sự tương tác của 5 nhóm yếu tố: Dòng chảy ($V_c, \beta_c$), Gió ($V_w, Y_r$), Tình huống đối hướng (Head-on), Tình huống cắt hướng (Crossing) và Nhân tố con người (Human Factor).
  2. Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để định lượng hóa tác động phi tuyến tính và tính chất mờ của nhân tố con người trong quá trình ra quyết định điều động tránh va?
    • Giả thuyết 2 (H2): Nhân tố con người có thể được định lượng chính xác thông qua hệ suy mờ mờ Sugeno bậc 0 với 3 biến đầu vào: Trình độ chuyên môn, Kinh nghiệm đi biển và Yếu tố tâm lý.
  3. Câu hỏi nghiên cứu 3 (RQ3): Mô hình toán và phần mềm mô phỏng có thể lượng hóa mức độ rủi ro theo thời gian thực và hỗ trợ công tác quy hoạch, mở rộng luồng hàng hải đạt chuẩn PIANC hay không?
    • Giả thuyết 3 (H3): Chỉ số rủi ro tích hợp $K$ cho phép tái hiện chân thực áp lực điều động trên từng phân đoạn luồng và dự báo chính xác mức suy giảm rủi ro khi điều chỉnh bình đồ tuyến luồng.

Phạm vi nghiên cứu bao quát dữ liệu khí tượng thủy văn 17 năm (2000–2017), hồ sơ điều tra tai nạn hàng hải 15 năm (2007–2021) từ Cảng vụ Hàng hải Vũng Tàu, mẫu khảo sát 258 chuyên gia (Thuyền trưởng, Đại phó, Hoa tiêu ngoại hạng, Điều hành viên VTS) và dữ liệu hành trình AIS của 11 tình huống điều động tránh va đối kháng trực tiếp trên luồng Vũng Tàu – Thị Vải.

Literature Review và Positioning

Nghiên cứu an toàn hàng hải quốc tế được phân hóa thành các trường phái tiếp cận riêng biệt với những ưu thế và hạn chế nội tại:

+---------------------------------------------------------------------------------------+
|                SO SÁNH CÁC MÔ HÌNH ĐÁNH GIÁ RỦI RO HÀNG HẢI QUỐC TẾ                   |
+---------------+-------------------+----------------------+----------------------------+
| Mô hình       | Tác giả / Tổ chức | Nguyên lý Tiếp cận   | Hạn chế khi ứng dụng VN    |
+---------------+-------------------+----------------------+----------------------------+
| PAWSA         | IALA / IMO (2010) | Định tính (Delphi)   | Thiếu định lượng, tốn kém  |
| IWRAP Mk2     | DAMSA / Fujii     | Xác suất Hình học    | Sai lệch vùng luồng phức tạp|
| ES Model      | Kinzo Inoue (2002)| Áp lực buồng lái 180°| Chỉ chuẩn hóa cho Nhật Bản |
| PARK Model    | Young-Soo Park    | Hồi quy Môi trường   | Phụ thuộc cấu trúc Busan   |
| Luận án (VN)  | Lê Văn Thức (2022)| AHP + Fuzzy Sugeno   | Tích hợp chuyên sâu Con người|
+---------------+-------------------+----------------------+----------------------------+

Trường phái thứ nhất đại diện bởi mô hình PAWSA (Port and Waterway Safety Assessment) do Hiệp hội Báo hiệu Hàng hải và Hải đăng Quốc tế (IALA) khởi xướng và được IMO ban hành theo Thông tư MSC/Circ.296 (2010). PAWSA vận hành dựa trên phương pháp chuyên gia Delphi (Estimate-Talk-Estimate - ETE) thông qua các hội thảo 2 ngày. Nhược điểm căn bản của PAWSA là mang bản chất "định tính chủ quan", chi phí tổ chức cao, khó cập nhật liên tục và theo khuyến cáo của IMO, PAWSA chỉ đóng vai trò nhận diện sơ bộ trước khi chuyển sang mô hình định lượng.

Trường phái thứ hai là mô hình IWRAP Mk2 (IALA Waterway Risk Assessment Program) do Cục An toàn Hàng hải Đan Mạch (DAMSA) và Đại học Kỹ thuật Đan Mạch phát triển dựa trên mô hình va chạm hình học kinh điển của Fujii và cộng sự (1974, 1984) kết hợp hệ số nguyên nhân ($P_c = 0.5 \times 10^{-4}$ đối với tàu đối hướng và $1.3 \times 10^{-4}$ đối với tàu cắt hướng/chuyển hướng). Công thức tổng quát $N_{coll} = N_G \times P_c$ của IWRAP Mk2 tỏ ra hữu hiệu trên các tuyến phân luồng hàng hải rộng mở (như eo biển Malacca hay eo biển Dover), nhưng bộc lộ sai số nghiêm trọng khi áp dụng vào luồng hẹp và luồng uốn khúc. Điển hình, nghiên cứu ứng dụng IWRAP Mk2 tại Vịnh Gành Rái của nhóm tác giả Nguyễn Xuân Thành, Park Young-Soo và Park Jin-Soo (2013) kết luận khu vực cửa sông Dinh (từ phao 8–9 đến phao 1–2) có nguy cơ đâm va cao nhất. Tuy nhiên, trích dẫn dữ liệu điều tra thực tế của Cảng vụ Hàng hải Vũng Tàu từ năm 2007 đến 2021: "kết quả thống kê tai nạn thực tế trong 15 năm qua cho thấy cửa sông Dinh không xảy ra bất kỳ vụ tai nạn đâm va nào", chứng minh rằng các hàm phân bố giao thông chuẩn tắc trong IWRAP không phản ánh đúng hành vi bám luồng thực tế của phương tiện.

Trường phái thứ ba là mô hình áp lực môi trường ES Model (Environmental Stress Model) của Giáo sư Kinzo Inoue (2002) công bố tại Hội nghị PIANC Sydney lần thứ 30 và mô hình PARK (Potential Risk Assessment Model) của Young-Soo Park và cộng sự (Busan, Hàn Quốc). ES Model đánh giá mức độ an toàn qua chỉ số áp lực môi trường tổng hợp ($ESA = ES_L + ES_S$), tương quan với biến thiên sinh lý nhịp tim (R-R interval spectral analysis). Dù vậy, ES Model chỉ quét góc nhìn $180^\circ$ phía trước mũi tàu, bỏ qua các nguy cơ tiếp cận từ mạn và phía sau lái, đồng thời bộ trọng số thực nghiệm được chuẩn hóa trên mẫu Thuyền trưởng người Nhật Bản, không tương thích với văn hóa điều động và nhận thức rủi ro tại vùng nước hỗn tạp của Việt Nam.

Luận án định vị khoảng trống học thuật (research gap): Chưa có một công trình nào tích hợp đồng thời cơ chế động lực học thủy khí, hình học va chạm cục bộ và định lượng hóa nhân tố con người mờ (Human Factor) cho một vùng nước luồng hẹp có mật độ tàu cá, tàu nội địa và tàu biển đan xen như tại Vũng Tàu.

Đó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 lý thuyết an toàn giao thông thủy thông qua việc tái định nghĩa rủi ro đâm va: Rủi ro không thuần túy là xác suất hình học khách quan ($N_G$) như cách tiếp cận của MacDuff (1974) hay Pedersen (1995), mà là hàm số phản ánh "mức độ khó khăn của người điều khiển phương tiện khi dẫn tàu trong một môi trường giao thông và ngoại cảnh xác định".

+-----------------------------------------------------------------------------------+
|                 CẤU TRÚC MÔ HÌNH SUY MỜ SUGENO ĐỊNH LƯỢNG NHÂN TỐ CON NGƯỜI        |
+-----------------------------------------------------------------------------------+
|  [Biến Đầu Vào]                 [Cơ Sở 27 Luật Mờ Sugeno]       [Đầu Ra Định Lượng]|
|  - Trình độ ($H_1 \in [-2, +2]$)   -->  IF ($H_1$ is High &       -->  Hệ số Con người|
|  - Kinh nghiệm ($H_2 \in [-2, +2]$)-->      $H_2$ is Low &              $H_{con\,nguoi}$  |
|  - Tâm lý ($H_3 \in [-2, +2]$)     -->      $H_3$ is Stressed)          $\in [-2, +2]$    |
|                                         THEN $z = f(H_1, H_2, H_3)$                       |
+-----------------------------------------------------------------------------------+

Khung lý thuyết xây dựng biểu thức toán học tổng quát xác định hệ số an toàn rủi ro tổng hợp $K$: $$K = f(F_{dong}, F_{gio}, F_{doi_huong}, F_{cat_huong}, F_{con_nguoi})$$

Trong đó, đóng góp đột phá về mặt lý thuyết là việc mô hình hóa cấu trúc nhân tố con người bằng Logic mờ (Fuzzy Logic) của Lotfi Zadeh (1965) kết hợp thuật toán suy luận Sugeno (1985). Nhân tố con người chiếm tỷ lệ thống kê trên 80% nguyên nhân tai nạn được phân rã thành 3 biến ngôn ngữ mờ:

  • $H_1$: Trình độ chuyên môn (Kiến thức COLREGs 72, quy chế cảng biển, kỹ năng radar/ECDIS).
  • $H_2$: Kinh nghiệm công tác (Thâm niên đi biển, số lượt dẫn tàu qua luồng Vũng Tàu).
  • $H_3$: Trạng thái tâm lý (Áp lực điều động, mức độ căng thẳng, sự mệt mỏi buồng lái).

Mỗi biến đầu vào được gán 3 hàm liên thuộc (Membership Functions: Thấp/Kém, Trung bình, Cao/Tốt) trên miền biến thiên $[-2, +2]$, thiết lập nên cơ sở 27 luật mờ dạng khối $3 \times 3 \times 3$ đồng nhất: $$\text{Rule } i: \text{IF } (H_1 \text{ is } A_i) \text{ AND } (H_2 \text{ is } B_i) \text{ AND } (H_3 \text{ is } C_i) \text{ THEN } z_i = a_i H_1 + b_i H_2 + c_i H_3 + k_i$$

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

Khung phân tích của luận án kết hợp 3 trụ cột khoa học liên ngành: Thủy khí động lực học tàu thủy, Phân tích đa tiêu chí AHP và Hệ suy mờ thông minh:

  1. Mô đun Khí động lực học gió: Tính toán lực tác động và mô-men quay do gió gây dạt tàu dựa trên lý thuyết Isherwood (1972): $$X_{wind} = \frac{1}{2} C_X(Y_r) \rho_a V_r^2 A_T$$ $$Y_{wind} = \frac{1}{2} C_Y(Y_r) \rho_a V_r^2 A_L$$ $$N_{wind} = \frac{1}{2} C_N(Y_r) \rho_a V_r^2 A_L L$$ (Trong đó $C_X, C_Y, C_N$ là các hệ số thực nghiệm; $\rho_a$ là mật độ không khí; $A_T, A_L$ là diện tích cản gió chính diện và mạn; $V_r, Y_r$ là tốc độ và góc tới của gió tương đối).
  2. Mô đun Thủy động lực học dòng chảy: Xác định ảnh hưởng của trường dòng triều Vịnh Gành Rái phân rã theo hệ tọa độ tàu ($V_b = V_c \cos \beta_c$; $U_b = V_c \sin \beta_c$).
  3. Mô đun Tương tác đối kháng hai tàu: Ứng dụng lý thuyết nhiễu loạn hai chiều cho tình huống đối hướng (Head-on) $F = \frac{1}{2} C_F \rho B_0 T_0 V_p^2$ và chỉ số động học tiếp cận gần nhất đối với tình huống cắt hướng (Crossing) qua tọa độ $CPA$ (Closest Point of Approach) và $TCPA$ (Time to Closest Point of Approach): $$CPA = \sqrt{R^2 (C_{ot} - B_t)}$$ $$TCPA = \frac{R \cos(C_{ot} - B_t)}{V_{0t}}$$

Điều kiện biên xác lập: Khung phân tích áp dụng tối ưu cho tàu tự hành có trang bị AIS hành hải trong luồng hạn chế độ sâu và luồng có chiều rộng đáy từ $150\text{ m}$ đến $350\text{ m}$.

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

+-----------------------------------------------------------------------------------+
|               TIẾN TRÌNH PHƯƠNG PHÁP NGHIÊN CỨU LIÊN HOÀN (MIXED METHODS)         |
+-----------------------------------------------------------------------------------+
|  [Thu thập Dữ liệu Đa nguồn]                                                      |
|   ├── AIS VTS Vũng Tàu + Thủy triều/Gió (2000-2017)                               |
|   └── Hồ sơ Điều tra 15 năm Tai nạn Hàng hải                                      |
|                                                                                   |
|  [Khảo sát Chuyên gia Ngữ cảnh Hóa]                                               |
|   └── N = 258 Chuyên gia (Thuyền trưởng, Hoa tiêu, VTS) Thang đo Likert 5 mức     |
|                                                                                   |
|  [Xử lý Thống kê & Chuẩn hóa Trọng số]                                            |
|   ├── SPSS: Phân tích Hồi quy Đa biến & ANOVA (p < 0.001)                         |
|   └── AHP (Saaty): Ma trận Nhất quán Cặp (CR < 0.10)                              |
|                                                                                   |
|  [Phát triển & Kiểm định Mô phỏng]                                                |
|   ├── C++: Fast-Time Simulation (Giao diện Quản lý Thời gian thực)                |
|   └── MATLAB: 27 Luật Mờ Sugeno -> Kiểm chứng 11 Kịch bản Thực tế                 |
+-----------------------------------------------------------------------------------+

Thiết kế nghiên cứu

Nghiên cứu áp dụng thế giới quan thực chứng (Positivism) kết hợp chủ nghĩa thực tế phản biện (Critical Realism). Thiết kế nghiên cứu định lượng phức hợp (Multi-method Quantitative Design) kết nối quy nạp dữ liệu điều tra thực địa với diễn dịch toán học mô phỏng. Nghiên cứu phân tầng thành 3 cấp độ:

  • Cấp độ vĩ mô: Đánh giá phân bố lưu lượng tổng thể từ hải đồ điện tử và AIS.
  • Cấp độ trung mô: Đánh giá hình học tương tác giữa các luồng hàng hải (Vũng Tàu – Thị Vải, Sài Gòn – Vũng Tàu, luồng Sông Dinh).
  • Cấp độ vi mô: Mô hình hóa từng pha điều động tránh va cục bộ giữa từng cặp tàu ($L > 200\text{ m}$, $150\text{–}200\text{ m}$, $100\text{–}150\text{ m}$, $<100\text{ m}$).

Quy trình nghiên cứu rigorous

Quy trình khảo sát thực địa được chuẩn hóa qua các bước:

  1. Chiến lược chọn mẫu: Mẫu thuận tiện kết hợp phân tầng chuyên môn cao gồm $N = 258$ chuyên gia có thâm niên công tác thực tế từ 3 đến trên 13 năm (gồm Thuyền trưởng tàu biển viễn dương, Hoa tiêu hàng hải Vũng Tàu, Cán bộ điều phối VTS).
  2. Công cụ đo lường: Phiếu khảo sát cấu trúc đóng với thang đo 5 mức độ an toàn chuẩn hóa: $-2$ (Rất nguy hiểm/Rất mất an toàn), $-1$ (Nguy hiểm/Mất an toàn), $0$ (Bình thường/Không ảnh hưởng), $+1$ (An toàn/Thuận lợi), $+2$ (Rất an toàn/Rất thuận lợi).
  3. Quy trình Tam giác giác đạc (Triangulation): Kết hợp chéo giữa dữ liệu hồi cứu tai nạn thực tế 2007–2021, dữ liệu cảm biến AIS/Radar VTS và kết quả đánh giá chuyên gia độc lập.
  4. Kiểm định thang đo: Dữ liệu khảo sát chuyên gia sau khi mã hóa được đưa vào phần mềm IBM SPSS Statistics để kiểm định độ tin cậy Cronbach's Alpha ($\alpha > 0.85$), phân tích ANOVA và xác lập bảng ma trận hệ số quan trọng tương đối AHP (Analytic Hierarchy Process) của Thomas L. Saaty. Tỷ số nhất quán (Consistency Ratio - $CR$) của tất cả các ma trận so sánh cặp đều đạt giá trị $CR < 0.10$, bảo đảm tính nhất quán toán học.

Data và phân tích

Phân tích hồi quy tuyến tính đa biến trên SPSS xác định được phương trình toán học lượng hóa rủi ro đâm va trên luồng Vũng Tàu: $$K = 0.285 \times F_{dong} + 0.198 \times F_{gio} + 0.246 \times F_{doi_huong} + 0.312 \times F_{cat_huong} + 0.384 \times F_{con_nguoi} + \epsilon$$

Kết quả kiểm định mô hình thống kê cho thấy: Hệ số xác định $R^2 = 0.782$, kiểm định ANOVA có mức ý nghĩa $p\text{-value} < 0.001$, khẳng định mô hình có độ tương thích xuất sắc với dữ liệu thực nghiệm.

Dựa trên thuật toán hồi quy và trọng số AHP, nhóm nghiên cứu đã xây dựng phần mềm Mô phỏng nhanh (Fast Time Simulation) độc lập bằng ngôn ngữ C++. Phần mềm nạp trực tiếp tọa độ AIS, dữ liệu mớn nước ($D$), chiều dài lớn nhất ($LOA$), tốc độ thực tế ($SOG$), hướng đi ($COG$), kết hợp trường gió và dòng triều theo thời gian thực để xuất ra bản đồ nhiệt hiển thị hệ số an toàn $K$ trên từng lý trình luồng.

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

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

+-----------------------------------------------------------------------------------+
|               KẾT QUẢ KIỂM ĐỊNH THỰC NGHIỆM TRÊN 11 KỊCH BẢN VTS VŨNG TÀU          |
+-----+-----------------------+-------------------------+-----------+---------------+
| Case| Tàu Chủ               | Tàu Mục Tiêu            | Tình huống| Kết quả K     |
+-----+-----------------------+-------------------------+-----------+---------------+
| 1   | Wan Hai 302 (199.9m)  | Medi Valencia (190m)    | Đối hướng | Nguy hiểm ($K<0$)|
| 2   | Stever (199.9m)       | Starship Ursa (172.2m)  | Cắt hướng | Cảnh báo vừa  |
| 3   | Stever (199.9m)       | CSC Xin Hai (125m)      | Cắt hướng | An toàn ($K>0$)|
| 4   | Tsing Tao Exp. (335m) | Delaperouse (99m)       | Cắt hướng | Rất an toàn   |
| 5   | Buxhansa (208m)       | APL Cairo (207.4m)      | Đối hướng | Nguy hiểm cao |
| 6   | KMTC Gwangyang (172m) | APL Cairo (207.4m)      | Đối hướng | Áp lực lớn    |
| 7   | Tsing Tao Exp. (335m) | APL Cairo (207.4m)      | Cắt hướng | Thao diễn gấp |
| 8   | YM Intelligent (172m) | East Bangkok (177m)     | Đối hướng | Ổn định       |
| 9   | YM Intelligent (172m) | Zhong Xin Pearl (225m)  | Đối hướng | Căng thẳng cao|
| 10  | YM Intelligent (172m) | Jinyunhe (182.9m)       | Đối hướng | Nguy hiểm ($K<0$)|
| 11  | Kiểm chứng tổng hợp   | Mật độ tàu cá cao điểm  | Hỗn hợp   | Mất an toàn   |
+-----+-----------------------+-------------------------+-----------+---------------+
  1. Xác lập định lượng vị trí điểm nóng rủi ro: Nghiên cứu bác bỏ kết luận của các mô hình trước đây và chứng minh rằng đoạn luồng từ Phao số 0 đến Phao số 3 và khu vực Phao 7–9 luồng Vũng Tàu – Thị Vải (khu vực Vịnh Gành Rái giao cắt luồng Sài Gòn – Vũng Tàu và tuyến hành trình tàu cá) là điểm hội tụ rủi ro đâm va cao nhất toàn tuyến ($p < 0.01$).
  2. Phát hiện quy luật phi đối xứng trong tình huống đối hướng (Head-on): Dữ liệu quan trắc từ 11 kịch bản thực tế (như Case 1: Wan Hai 302 gặp Medi Valencia; Case 5: Buxhansa gặp APL Cairo) chứng minh khi hai tàu có $LOA > 200\text{ m}$ đi đối hướng trong điều kiện dòng triều chảy xiết ($V_c > 1.8\text{ m/s} \approx 3.5\text{ knots}$), khoảng cách an toàn chính ngang bắt buộc phải duy trì tối thiểu bằng 1 lần chiều dài tàu chủ ($1 \times LOA \approx 200\text{ m}$). Mọi tình huống khoảng cách giảm xuống dưới $0.5 \times LOA$ (như Case 9: YM Intelligent vs Zhong Xin Pearl, khoảng cách $86\text{ m}$) đều đẩy hệ số $K$ vào ngưỡng nguy hiểm nghiêm trọng ($K < -1.5$).
  3. Lượng hóa tương quan nhường đường trong tình huống cắt hướng (Crossing): Trong Case 7 (Tsing Tao Express, $LOA = 335\text{ m}, DWT = 103.800$ cắt hướng APL Cairo, $LOA = 207.4\text{ m}$), phần mềm mô phỏng chỉ ra tàu được nhường đường chịu áp lực tâm lý tương đương tàu phải nhường đường nếu hành động chuyển hướng tránh va của tàu đối phương diễn ra muộn ở cự ly dưới $1.2\text{ hải lý}$.
  4. Khớp nối hoàn hảo giữa hàm mờ MATLAB và quan sát thực địa: Chạy thử nghiệm 11 kịch bản trên bộ mờ Sugeno MATLAB (với các thông số dòng $V_c = 1.7\text{ hl/h}, \beta_c = 25^\circ$, gió cấp 4 $V_w = 7.2\text{ m/s}$) cho kết quả đầu ra hệ số rủi ro hoàn toàn tương thích với các biên bản điều động tránh va thực tế được VTS ghi nhận.

Implications đa chiều

  • Về mặt học thuật: Cung cấp một phương pháp luận lai ghép (Hybrid Methodology) kết hợp AHP, Hồi quy và Fuzzy Sugeno, có khả năng mở rộng để nghiên cứu rủi ro hàng hải tại các vùng nước luồng hẹp phức tạp tại Đông Nam Á (như luồng sông Soài Rạp, luồng Lạch Huyện, luồng eo biển Malacca).
  • Về mặt thực tiễn khai thác: Phần mềm mô phỏng nhanh bằng C++ là công cụ trực quan hóa hỗ trợ trực tiếp cho điều hành viên VTS phát hiện sớm các kịch bản đâm va tiềm ẩn trước 15–20 phút, đưa ra khuyến cáo kịp thời cho tàu thuyền thay vì can thiệp thụ động.
  • Về mặt quy hoạch chính sách: Cung cấp cơ sở khoa học trực tiếp cho Cục Hàng hải Việt Nam và Bộ Giao thông Vận tải trong việc lập dự án cải tạo, nâng cấp mở rộng luồng Vũng Tàu – Thị Vải (đoạn từ phao số 0 vào khu bến cảng container Cái Mép), chuẩn hóa bề rộng luồng chạy tàu hai chiều đáp ứng chuẩn an toàn hàng hải quốc tế PIANC.

Limitations và Future Research

Nghiên cứu tồn tại một số giới hạn học thuật cần được nhìn nhận minh bạch:

  1. Độ bao phủ của dữ liệu mờ: 27 luật mờ Sugeno trong mô hình được thiết lập dựa trên chuyên gia chủ yếu là người Việt Nam và một số thuyền trưởng quốc tế hoạt động thường xuyên tại khu vực; do đó mức độ nhạy cảm tâm lý có thể biến thiên đối với các thuyền bộ đa quốc gia lần đầu cập cảng.
  2. Mô phỏng 2D: Các lực tương tác thủy động tàu - tàu trong mô hình được tính toán gần đúng thông qua phương pháp nhiễu loạn hai chiều, chưa mô phỏng đầy đủ hiệu ứng cạn (Shallow Water Effect), hiệu ứng hút mạn đồi bờ (Bank Effect) và tương tác 3 chiều phức tạp dưới đáy luồng (Squat effect) trong điều kiện mớn nước tới hạn ($Draft \approx Depth$).
  3. Yếu tố tầm nhìn hạn chế: Mô hình chủ yếu kiểm chứng trong điều kiện khí hậu nhiệt đới gió mùa điển hình (gió Đông Nam/Đông Bắc cấp 4–6), chưa khảo sát sâu tình huống sương mù dày đặc hoặc mưa dông cực đoan làm suy giảm tầm nhìn thị giác dưới 0.5 hải lý.

Chương trình nghiên cứu tương lai cần hướng tới:

  • Mở rộng thuật toán bằng việc tích hợp Mạng nơ-ron nhân tạo (Artificial Neural Networks - ANN) hoặc Mạng thích nghi Thần kinh - Mờ (ANFIS) để mô hình có khả năng tự học trực tiếp từ "Big Data" dữ liệu AIS hàng triệu lượt tàu.
  • Nâng cấp mô phỏng thủy động lực học 3D (CFD) để đánh giá chi tiết lực hút đẩy thủy động giữa hai tàu siêu lớn ($ULCV > 18.000\text{ TEU}$) khi chạy đối hướng trong luồng hẹp.
  • Tích hợp thêm các kịch bản va chạm của phương tiện nổi không người lái (Maritime Autonomous Surface Ships - MASS) trong tương lai.

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

Luận án tạo ra những tác động lan tỏa mạnh mẽ trên nhiều phương diện:

  • Giá trị học thuật: Là công trình tiên phong tại Việt Nam giải quyết trọn vẹn bài toán lượng hóa rủi ro đâm va bằng công cụ toán học và mô phỏng số, mở ra hướng trích dẫn quan trọng cho các nghiên cứu sinh chuyên ngành Khoa học Hàng hải, Quản lý Cảng và Logistics.
  • Chuyển đổi công nghệ ngành hàng hải: Chuyển đổi mô hình quản lý an toàn hàng hải của các Cảng vụ từ phương thức "quản lý dựa trên sự vụ" (Reactive management) sang "quản lý và dự báo rủi ro chủ động" (Proactive Risk Assessment).
  • Tác động kinh tế - xã hội: Đảm bảo an toàn thông suốt cho tuyến luồng hàng hải huyết mạch trung chuyển hàng hóa quốc tế Cái Mép – Thị Vải, giảm thiểu nguy cơ tai nạn đâm va có thể gây tắc nghẽn luồng biển, tổn thất hàng trăm triệu USD/ngày cho nền kinh tế, đồng thời phòng ngừa triệt để nguy cơ tràn dầu và ô nhiễm môi trường biển khu sinh quyển Cần Giờ và Vũng Tàu.

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

  • Nghiên cứu sinh và Giới học thuật: Thụ hưởng bộ cơ sở dữ liệu thực nghiệm chuẩn hóa về luồng Vũng Tàu, hệ thống 27 luật mờ Sugeno và khung lý thuyết đánh giá an toàn giao thông đường thủy.
  • Hoa tiêu hàng hải và Thuyền trưởng: Tiếp cận công cụ tính toán và hiểu rõ các ranh giới an toàn tới hạn ($CPA, TCPA$, khoảng cách chính ngang $1 \times LOA$) khi điều động trên luồng hẹp.
  • Cơ quan quản lý VTS và Cảng vụ Hàng hải: Sở hữu phần mềm mô phỏng nhanh hỗ trợ kiểm soát không gian luồng, phân bổ vị trí neo đậu, quy định khoảng cách vượt tàu và tốc độ tối ưu trên từng phân đoạn luồng.
  • Các Viện Thiết kế, Tư vấn Công trình Hàng hải: Sử dụng phần mềm làm công cụ kiểm chứng tĩnh và động khi thiết kế bình đồ, mở rộng luồng hàng hải và tính toán kích thước vũng quay tàu.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của nghiên cứu là gì và đã mở rộng lý thuyết nào? Trả lời: Đóng góp lý thuyết độc đáo nhất là việc mô hình hóa thành công nhân tố con người mờ (Human Factor) - yếu tố chiếm trên 80% nguyên nhân tai nạn nhưng thường bị bỏ qua trong các mô hình xác suất cổ điển - vào biểu thức rủi ro hàng hải thông qua Hệ suy mờ Sugeno (Sugeno Fuzzy Inference System). Nghiên cứu đã mở rộng Lý thuyết Tập mờ của Lotfi Zadeh và Lý thuyết Áp lực Môi trường của Kinzo Inoue bằng việc tích hợp 3 biến tâm lý học hành vi (Trình độ, Kinh nghiệm, Tâm lý) vào hệ tọa độ điều động thực tế.

2. Điểm cải tiến phương pháp luận của luận án khi so sánh với 2 mô hình quốc tế PAWSA và IWRAP Mk2? Trả lời: So với PAWSA (thuần túy định tính, phụ thuộc vào hội thảo Delphi), luận án xây dựng được thuật toán định lượng chạy trên phần mềm máy tính (C++ và MATLAB). So với IWRAP Mk2 (tính toán hình học xác suất thuần túy, cho kết quả sai lệch tại luồng hẹp cong như trường hợp Vịnh Gành Rái), luận án tính toán rủi ro theo bản chất động lực học thực tế của tàu, trường thủy triều/gió và khả năng thao diễn tránh va thực tế, loại bỏ giả định tĩnh phi thực tế của IWRAP.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì? Trả lời: Phát hiện bất ngờ nhất là đoạn cửa sông Dinh – khu vực từng bị mô hình quốc tế IWRAP Mk2 cảnh báo có nguy cơ đâm va cao nhất – trên thực tế suốt 15 năm (2007–2021) không hề xảy ra vụ đâm va nào. Ngược lại, đoạn luồng tưởng chừng an toàn hơn từ Phao 0 đến Phao 3 và Phao 7–9 mới chính là "điểm đen" đâm va thực sự do sự giao cắt phức tạp của các tuyến tàu cá tự phát và dòng triều lệch trục luồng $35^\circ$.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) không? Trả lời: Có. Luận án cung cấp trọn vẹn quy trình tái lập nghiên cứu bao gồm: Cấu trúc bảng câu hỏi khảo sát 5 mức độ Likert, phương pháp mã hóa và xử lý ma trận trên IBM SPSS Statistics, bảng trọng số AHP đã kiểm định nhất quán ($CR < 0.1$), toàn bộ mã nguồn mô phỏng trên nền tảng MATLAB (Fuzzy Logic Toolbox) tại Phụ lục 5 và mã nguồn phần mềm Fast-Time Simulation viết bằng C++ tại Phụ lục 7.

5. Lộ trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào? Trả lời: Lộ trình 10 năm tập trung vào 3 trọng tâm: (1) Chuẩn hóa phần mềm mô phỏng thành mô-đun tích hợp trực tiếp vào hệ thống ECDIS và Radar VTS quốc gia; (2) Mở rộng xây dựng bản đồ số rủi ro cho toàn bộ hệ thống cảng biển loại 1 và loại đặc biệt của Việt Nam (Hải Phòng, Đà Nẵng, Quy Nhơn, TP.HCM); (3) Nghiên cứu tích hợp trí tuệ nhân tạo (Deep Learning) để tự động nhận dạng hành vi bất thường của tàu cá và phương tiện thủy nội địa nhằm điều khiển tự động tàu không người lái MASS.

Kết luận

  1. Nghiên cứu đã hoàn thành xuất sắc việc xây dựng cơ sở dữ liệu khoa học toàn diện, có độ tin cậy cao về điều kiện khí tượng thủy văn, lưu lượng tàu thuyền và lịch sử tai nạn hàng hải tại cảng biển Vũng Tàu.
  2. Xác lập luận điểm khoa học đột phá: Lượng hóa rủi ro đâm va trên cơ sở đánh giá mức độ khó khăn của người điều khiển phương tiện dưới tác động tổng hợp của 5 nhóm yếu tố ngoại cảnh và nội tại.
  3. Thiết lập thành công hệ suy mờ Sugeno gồm 27 luật mờ, giải quyết triệt để bài toán định lượng hóa nhân tố con người trong an toàn hàng hải.
  4. Phát triển thành công 02 công cụ ứng dụng thực tiễn: Phần mềm Mô phỏng nhanh (Fast Time Simulation) trên nền tảng C++ và Mô hình tính toán hàm mờ trên MATLAB, được kiểm chứng chính xác tuyệt đối qua 11 kịch bản điều động thực tế.
  5. Cung cấp luận cứ khoa học vững chắc phục vụ trực tiếp cho việc hoạch định chính sách, hoàn thiện văn bản pháp quy hàng hải và dự án cải tạo, nâng cấp luồng hàng hải Vũng Tàu – Thị Vải, đáp ứng chiến lược phát triển kinh tế biển bền vững của Việt Nam.