Tổng quan luận án

Trong hoạt động vận tải biển, tàu luôn chịu tác động liên tục từ các ngoại lực môi trường như sóng và gió. Khả năng tự quay trở về trạng thái cân bằng ban đầu khi các lực kích động bên ngoài chấm dứt tác động là điều kiện kỹ thuật cốt lõi để duy trì sự an toàn của con tàu trong hành trình hàng hải. Khả năng cơ học này được định nghĩa là tính ổn định của tàu. Nhằm đảm bảo an toàn, các quy định hàng hải quốc tế yêu cầu sỹ quan quản lý và khai thác tàu phải kiểm tra, đánh giá và duy trì các thông số ổn định tuân thủ tiêu chuẩn của Tổ chức Hàng hải Quốc tế (IMO), tiêu biểu là Bộ luật quốc tế về Ổn định nguyên vẹn (IS Code) cùng các quy định bổ sung cho tàu chở gỗ, tàu chở chất lỏng và tàu chở hàng rời. Bên cạnh đó, tiêu chuẩn ổn định khi tàu bị hư hỏng được chuẩn hóa trong Công ước quốc tế về an toàn sinh mạng con người trên biển (SOLAS 74 sửa đổi, bổ sung năm 2014).

Dữ liệu thống kê giai đoạn 2009–2018 của tổ chức bảo hiểm Allianz ghi nhận 1.036 vụ tai nạn hàng hải có mức độ thiệt hại nghiêm trọng; trong đó tai nạn chìm tàu chiếm số lượng cao nhất với 551 vụ (dao động từ 46 đến 132 vụ mỗi năm). Riêng đối với phân khúc tàu chở hàng rời, ghi nhận 93 vụ tai nạn nghiêm trọng và khoảng 80% nguyên nhân chìm đắm xuất phát từ việc tàu không đảm bảo các điều kiện về ổn định nguyên vẹn và ổn định tai nạn.

Về mặt cơ học lý thuyết, tính ổn định của tàu biển nói chung và tàu hàng rời nói riêng được đặc trưng bởi hai thông số cơ bản: chiều cao thế vững ($G_0M$) và cánh tay đòn ổn định tĩnh ($G_0Z$). Thực tế khai thác hiện nay cho thấy quy trình tính toán và kiểm tra độ ổn định chỉ được thực hiện tại các mốc thời gian cố định: lập kế hoạch xếp hàng, trong tiến trình xếp dỡ, kết thúc xếp hàng và kiểm tra lần cuối trước khi rời bến. Khi đã thỏa mãn các tiêu chuẩn của IMO ở thời điểm xuất phát, tính ổn định này thường được mặc định là duy trì an toàn trong suốt chuyến đi.

Tuy nhiên, trong quá trình hành trình trên biển, nhiều yếu tố bất thường có thể phát sinh: hàng hóa bị xô lệch cơ học, hiện tượng hóa lỏng của hàng rời, sự hình thành mặt thoáng chất lỏng do điều chỉnh nước dằn (ballast) hoặc các khiếm khuyết ẩn tì của thân vỏ. Các yếu tố này làm suy giảm nghiêm trọng chiều cao thế vững và cánh tay đòn hồi phục, có thể dẫn đến nguy cơ lật tàu. Sĩ quan hàng hải trên thực tế gặp nhiều trở ngại trong việc tính toán lại các giá trị ổn định hiện thời khi tàu đang hành trình. Do đó, khoảng trống nghiên cứu cấp thiết đặt ra là phải phát triển một giải pháp cho phép theo dõi, cập nhật liên tục các thông số ổn định của tàu hàng rời theo thời gian thực để kịp thời đưa ra cảnh báo nguy hiểm. Đề tài "Nghiên cứu thiết kế hệ thống thông báo ổn định theo thời gian thực cho tàu hàng rời" được xây dựng nhằm giải quyết trực tiếp yêu cầu này.

Mục tiêu nghiên cứu của luận án:

  • Thiết kế hệ thống thông báo ổn định của tàu hàng rời có khả năng hiển thị liên tục các thông số đặc trưng cho ổn định nguyên vẹn theo thời gian thực, xác định xu hướng thay đổi giá trị và phát tín hiệu cảnh báo khi tính ổn định của tàu không thỏa mãn các tiêu chuẩn quy định trong Bộ luật quốc tế về ổn định nguyên vẹn (IS Code 2008) của Tổ chức Hàng hải Quốc tế (IMO).

Đối tượng nghiên cứu:

  • Tàu hàng rời;
  • Tính ổn định của tàu hàng rời;
  • Phương pháp đánh giá tính ổn định của tàu hàng rời.

Phạm vi nghiên cứu:

  • Đánh giá ổn định nguyên vẹn của tàu hàng rời khi hành trình trên biển;
  • Điều kiện khảo sát: Tàu nổi ở trạng thái cân bằng, không bị nghiêng ngang tĩnh ban đầu.

Tổng quan tài liệu và vị trí của luận án

Tổng quan các công trình nghiên cứu trên thế giới cho thấy nhiều hướng tiếp cận khác nhau đối với bài toán ổn định tàu:

  • Năm 2012, tác giả Apostolos Papanikolaou (Đại học Kỹ thuật Quốc gia Athens, Hy Lạp) công bố bài báo "Mục đích của tiêu chuẩn an toàn và ổn định tàu biển", tổng kết các kết quả thu được từ dự án GOALDS của Ủy ban châu Âu. Nghiên cứu này tập trung phân tích và đánh giá các công thức tính toán xác suất ổn định hư hỏng theo quy định tại Chương II-1 của Công ước SOLAS 74 đối với phân khúc tàu hàng khô và tàu khách.
  • Năm 2015, nhóm tác giả Mohamed Abdelkader Djebli, Benameur Hamoudi, Omar Imine và Lahouari Adjlout (Đại học Khoa học Công nghệ Oran Mohamed, Algeria) thực hiện nghiên cứu "Ứng dụng điện thoại thông minh trong việc thử nghiệm ổn định của tàu". Nhóm tác giả khai thác cụm cảm biến gia tốc tích hợp sẵn trên điện thoại thông minh nhằm ghi nhận chu kỳ lắc ngang, từ đó tính toán chiều cao thế vững thông qua công thức thực nghiệm của IMO. Tuy nhiên, thử nghiệm thực tế của công trình mới chỉ dừng lại trên mô hình tàu đánh cá tỷ lệ thu nhỏ 1/25 trong điều kiện thực nghiệm mô phỏng.
  • Năm 2015, tác giả Shadman Sakib (Đại học Kỹ thuật và Công nghệ Bangladesh) công bố công trình "Thiết bị mới để xác định ổn định động và xếp tải động cho tàu nội địa dựa trên chuyển động lắc ngang của nó". Tác giả sử dụng cảm biến đo gia tốc góc để xác định chu kỳ dao động lắc ngang, song nghiên cứu chỉ được triển khai thực nghiệm trong quy mô phòng thí nghiệm thông qua các cụm thiết bị tự tạo, áp dụng riêng cho các phương tiện thủy cỡ nhỏ hoạt động tại vùng nước nội địa Bangladesh.
  • Năm 2017, các tác giả Przemysław Krata và Wojciech Wawrzyński (Đại học Hàng hải Gdynia, Ba Lan) công bố nghiên cứu "Dự đoán vùng nguy hiểm thông qua mối quan hệ giữa chu kỳ lắc riêng và chiều cao thế vững tương đương". Công trình khảo sát tác động của chu kỳ dao động riêng lên tính chất phi tuyến của đường cong cánh tay đòn tĩnh $GZ$, qua đó đề xuất sử dụng thông số chiều cao thế vững tương đương ($GM_{eq}$) thay cho giá trị $GM$ ban đầu khi phân tích an toàn ổn định.
  • Về trang thiết bị thương mại: Trên thị trường quốc tế đã xuất hiện các thiết bị đo nghiêng điện tử thỏa mãn Nghị quyết MSC.363(92) của IMO do các đơn vị như Công ty TNHH Hàng hải UC (UC Marine Company Limited) hay Công ty TNHH Kỹ thuật Nanjing Ninglu sản xuất. Các thiết bị này thực hiện chức năng ghi nhận góc nghiêng ngang, góc nghiêng dọc và chu kỳ lắc ngang, nhưng hoàn toàn chưa tích hợp tính năng tự động tính toán các chỉ số ổn định theo thời gian thực cũng như chưa có cơ chế dự báo xu hướng suy giảm ổn định của tàu.

Tình hình nghiên cứu trong nước:

  • Năm 2016, tác giả Nguyễn Kim Phương và Bùi Văn Hưng công bố bài báo "Một số khuyến nghị trong công tác vận hành tàu hàng rời có nguy cơ hóa lỏng" trên Tạp chí Khoa học Công nghệ Hàng hải (số 46-3/2016). Nghiên cứu phân tích bản chất cơ học và vật lý của hiện tượng hóa lỏng hàng rời, từ đó đưa ra các hướng dẫn nghiệp vụ hỗ trợ sĩ quan hàng hải trong khai thác an toàn.
  • Năm 2018, các tác giả Nguyễn Kim Phương, Nguyễn Công Vịnh và Từ Mạnh Chiến công bố bài báo "Mô phỏng tiến trình xếp dỡ hàng trên tàu hàng rời phục vụ đào tạo và huấn luyện hàng hải" trên Tạp chí Khoa học Công nghệ Hàng hải (số 54-4/2018). Nghiên cứu xây dựng thuật toán mô phỏng diễn tiến ổn định của một tàu hàng rời cụ thể trong thời gian làm hàng tại cảng nhằm phục vụ công tác huấn luyện sĩ quan và sinh viên chuyên ngành điều khiển tàu biển.

Khoảng trống nghiên cứu được luận án xác định: Các công trình công bố trong nước và quốc tế trước đây chưa xây dựng được phương pháp luận hoàn chỉnh để tính toán và đánh giá các thông số ổn định của tàu hàng rời một cách liên tục theo thời gian thực trong điều kiện hành hải thực tế. Đồng thời, chưa có một giải pháp tích hợp đồng bộ giữa thiết bị đo phần cứng và thuật toán phần mềm để tự động xác định, hiển thị thông số và đưa ra các cảnh báo chủ động khi độ ổn định của tàu bị suy giảm.

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

Cơ sở lý thuyết của luận án dựa trên các nguyên lý thủy tĩnh học tàu thủy, động lực học vật rắn chuyển động trong chất lỏng và lý thuyết xử lý tín hiệu ngẫu nhiên:

  • Thuyết ổn định tàu thủy: Xem xét trạng thái cân bằng thông qua vị trí tương quan giữa trọng tâm $G$ và tâm nghiêng $M$. Khảo sát ba trạng thái cân bằng: bền ($GM > 0$), phiếm định ($GM = 0$) và không bền ($GM < 0$).
  • Ổn định tại góc nghiêng nhỏ: Mô men hồi phục được xác định theo biểu thức $M_{hp} = D \times GM \times \sin\theta$, trong đó chiều cao thế vững ban đầu $GM = KM - KG$. Trọng tâm tàu $KG$ được tổng hợp từ trọng lượng tàu không ($D_0, KG_0$) và các thành phần tải trọng riêng lẻ ($P_i, KG_i$).
  • Ổn định tại góc nghiêng lớn: Sử dụng đường cong cánh tay đòn ổn định tĩnh $G_0Z$ với quan hệ $G_0Z = KN - KG_0 \times \sin\theta$, trong đó giá trị $KN$ được trích xuất từ bảng đường cong hoành giao (Stability Cross Curves) theo lượng giãn nước $D$, và $KG_0$ là chiều cao trọng tâm sau khi đã hiệu chỉnh ảnh hưởng của mặt thoáng chất lỏng.
  • Tiêu chuẩn ổn định IS Code 2008 và tiêu chuẩn thời tiết: Đánh giá khả năng chống lật của tàu dưới tác động đồng thời của áp lực gió giật ($P = 540\text{ Pa}$) và sóng biển, dựa trên sự tương quan giữa diện tích hồi phục $b$ và diện tích nghiêng $a$.
  • Động lực học chuyển động lắc 6 bậc tự do: Mô hình hóa hệ phương trình vi phân chuyển động không gian theo các hệ trục tọa độ Euler, phân tách thành các chuyển động tịnh tiến (dọc - Surge, ngang - Sway, đứng - Heave) và chuyển động quay (chòng chành ngang - Roll, chúi dọc - Pitch, đảo mũi - Yaw).

Phương pháp nghiên cứu thực tế được triển khai bao gồm:

  • Phương pháp phân tích và tổng hợp: Hệ thống hóa hệ thống tiêu chuẩn kỹ thuật hàng hải quốc tế, cơ sở dữ liệu tai nạn và các công trình nghiên cứu tiền nhiệm.
  • Phương pháp mô hình hóa toán học: Mô hình hóa chuỗi thời gian góc nghiêng ngang $\theta(t)$ của tàu thành tín hiệu rời rạc $x(n) = r(n) + e(n)$, với $r(n)$ là thành phần dao động lắc riêng, $e(n)$ là tổng hợp giữa nhiễu ngẫu nhiên do sóng $w(n)$ và sai số đo đạc $d(n)$. Áp dụng thuật toán biến đổi Fourier rời rạc (DFT) và thuật toán biến đổi Fourier nhanh (FFT) để phân tích mật độ phổ công suất (PSD).
  • Phương pháp thuật toán dự đoán: Thiết lập kỹ thuật cửa sổ quan sát trượt có kích thước $V$ phút, độ dịch trượt $v$ giây chồng đè nhau. Áp dụng thuật toán trung bình hóa phổ năng lượng $Power_j = \text{AVG}(PSD_i)$ kết hợp mô hình hồi quy tuyến tính hoặc hàm nội suy đa thức $P(x)$ qua $q$ cửa sổ liên tiếp để dự báo tần số dao động riêng $F_{q+1}$ và chu kỳ lắc $T$.
  • Phương pháp chuyên gia: Tham vấn ý kiến chuyên môn từ các chuyên gia hàng hải và thuyền trưởng về các yếu tố rủi ro gây mất ổn định tàu hàng rời trên biển.
  • Phương pháp ứng dụng công nghệ thông tin: Lựa chọn môi trường lập trình để xây dựng phần mềm tự động hóa tính toán các bảng tra thủy tĩnh, vẽ đồ thị $G_0Z$ và kích hoạt cảnh báo ngưỡng an toàn.
  • Phương pháp thực nghiệm: Lắp đặt cụm cảm biến đo nghiêng và hệ thống thu nhận dữ liệu trực tiếp trên tàu để kiểm chứng độ tin cậy của thuật toán tính toán chu kỳ lắc và các thông số ổn định.

Nguồn tư liệu và dữ liệu: Hồ sơ kỹ thuật tàu (Bảng thủy tĩnh Hydrostatic Table, Bảng đường cong Cross Curves, Sơ đồ xếp hàng Cargo Stowage Plan, Bảng kiểm tra chu kỳ lắc Rolling Period Table); các số liệu thống kê tai nạn của Allianz; các tiêu chuẩn trong Bộ luật IS Code 2008 và SOLAS 74.

Nội dung chính theo từng chương

Chương 1: Cơ sở lý luận về ổn định tàu hàng rời theo thời gian thực

Chương 1 hệ thống hóa toàn diện các đặc trưng kết cấu của đội tàu hàng rời và nền tảng lý thuyết tính toán ổn định tàu. Tàu chuyên dụng chở hàng rời là loại tàu một boong, cấu trúc đáy đôi và mạn kép vững chắc, được thiết kế các két hông (top-side tanks) và két treo (hopper tanks) ở hai bên nhằm hạn chế diện tích mặt thoáng của hầm hàng và hỗ trợ điều chỉnh trọng tâm tàu. Tùy theo sức chở, tàu hàng rời được phân chia thành các nhóm kích thước tiêu chuẩn:

Phân cấp kích thước Phạm vi trọng tải toàn phần (DWT) Đặc điểm khai thác
Handysize 20.000 – 25.000 DWT (hoặc đến dưới 40.000 DWT) Phù hợp các cảng biển nhỏ, luồng lạch hạn chế
Handymax Từ 40.000 DWT Linh hoạt trong vận chuyển hàng rời thông dụng
Supramax Từ 50.000 DWT Sức chở trung bình, tính cơ động cao
Panamax 60.000 – 75.000 DWT Kích thước giới hạn đi qua kênh đào Panama cũ
Post-Panamax ("Baby capers") Dưới 125.000 DWT Kích thước vượt chuẩn Panamax cũ
Capesize Từ 125.000 DWT Gồm Newcastlemax (~185.000 DWT) và Setouchmax (~203.000 DWT)
VLOC (chuyển đổi từ VLCC) Trên 200.000 DWT Chuyên chở quặng sắt tải trọng siêu lớn

Về trang thiết bị làm hàng, tàu hàng rời được phân thành các nhóm: tàu có cẩu/băng tải tự làm hàng (Geared), tàu không có thiết bị làm hàng phụ thuộc hoàn toàn vào cảng (Gearless), tàu chở kết hợp quặng - hàng rời - dầu (OBO), tàu tự dỡ hàng trang bị gàu xúc/băng tải dọc mạn (Self-dischargers), và tàu xếp hàng rời dỡ đóng bao (BIBO - ví dụ tàu CHL INNOVATOR có công suất dỡ 300 tấn đường rời và đóng bao 50 kg trong một giờ).

Về mặt tiêu chuẩn kỹ thuật, chương tổng hợp các tiêu chí bắt buộc theo Bộ luật IS Code 2008:

  • Diện tích dưới đường cong cánh tay đòn ổn định tĩnh $G_0Z$ không được nhỏ hơn $0,055\text{ m-rad}$ tính đến góc nghiêng $30^\circ$, và không nhỏ hơn $0,090\text{ m-rad}$ tính đến góc nghiêng $40^\circ$ (hoặc góc ngập nước nếu góc này nhỏ hơn $40^\circ$).
  • Diện tích dưới đường cong $G_0Z$ nằm giữa khoảng $30^\circ$ và $40^\circ$ (hoặc góc ngập nước) không được nhỏ hơn $0,030\text{ m-rad}$.
  • Cánh tay đòn $G_0Z$ phải đạt độ lớn tối thiểu $0,20\text{ m}$ tại góc nghiêng $\ge 30^\circ$.
  • Góc nghiêng đạt giá trị cực đại của $G_0Z$ tối thiểu phải từ $25^\circ$ trở lên (khuyến nghị trên $30^\circ$).
  • Chiều cao thế vững ban đầu đã hiệu chỉnh mặt thoáng chất lỏng $G_0M$ không được nhỏ hơn $0,15\text{ m}$.

Đối với tiêu chuẩn thời tiết chịu tác động sóng gió, tay đòn nghiêng do gió liên tục $l_{w1}$ được tính theo công thức: $$l_{w1} = \frac{P \times A \times Z}{1000 \times g \times D}$$ với áp lực gió $P = 540\text{ Pa}$, gia tốc trọng trường $g = 9,81\text{ m/s}^2$, $A$ là diện tích hứng gió, $Z$ là cánh tay đòn hứng gió, và $D$ là lượng chiếm nước. Tay đòn gió giật xác định bằng $l_{w2} = 1,5 \times l_{w1}$. Góc nghiêng do gió tĩnh $\theta_0 \le 16^\circ$ hoặc không vượt quá 80% góc ngập mép boong. Biên độ lắc sóng $\theta_1$ được xác định theo: $$\theta_1 = 109 \times k \times X_1 \times X_2 \times \sqrt{r \times s}$$ trong đó hệ số $k$ phụ thuộc diện tích vây giảm lắc $A_k$ (tra theo tỷ số $100 \times A_k / (L_{wl} \times B)$); các hệ số $X_1, X_2$ tra theo tỷ số $B/d$ và hệ số béo thể tích $C_B$; hệ số $s$ phụ thuộc vào chu kỳ lắc ngang $T$.

Nội dung chương đi sâu phân tích ba yếu tố chính làm suy giảm ổn định của tàu hàng rời:

  1. Hiện tượng dịch chuyển khối hàng rời: Khi hàng bị xô dịch sang mạn một góc, trọng tâm hàng dời từ $g$ sang $g_2$, kéo theo trọng tâm tàu dời từ $G$ sang $G_2$. Các thành phần dịch chuyển ngang và đứng được tính bởi: $$GG_1 = \frac{P \times d_h}{D}; \quad G_1G_2 = \frac{P \times d_v}{D}$$ Nếu $G_x < 0$ (trọng tâm dời lên cao), chiều cao thế vững $GM$ sẽ suy giảm.
  2. Ảnh hưởng của két chất lỏng có mặt thoáng: Sự dịch chuyển chất lỏng trong két không đầy làm giảm chiều cao thế vững một lượng $GG_0$: $$GG_0 = \frac{\sum I_x \times \rho}{D}$$ với mô men quán tính mặt thoáng $I_x = (l \times b^3)/K$, hệ số hình dáng $K = 12$ đối với két hình chữ nhật, $K = 36$ với két tam giác vuông, và $K = 48$ với tam giác cân. Khi đó: $$G_0M = KM - KG - GG_0$$
  3. Hiện tượng hóa lỏng của hàng rời có tính nhão chảy: Độ ẩm vượt ngưỡng khiến hàng bị biến dạng sệt dưới tác động rung chấn của sóng biển, tạo ra mặt thoáng tự do lớn và dịch chuyển bất đối xứng sang một mạn. Quá trình này khiến chu kỳ lắc của tàu kéo dài bất thường và tàu bị nghiêng lệch cố định, dẫn tới lật chìm. Chương liệt kê các tai nạn nghiêm trọng:
  • Tàu Hui Long (20/05/2005) chở 5.185 tấn quặng fluorspar bị nghiêng $15^\circ$ rồi đắm tại Sungai Pakning, Indonesia.
  • Tàu M/v Asian Forest (18/07/2009) chở tinh quặng sắt lật tại vịnh Bengal, Ấn Độ.
  • Tàu M/v Jian Fu Star (27/10/2010) chở quặng niken đắm tại phía nam Đài Loan.
  • Tàu M/v Vinalines Queen (25/12/2011) của Việt Nam chở quặng niken bị đắm sau khi rời cảng Morowali, Indonesia.
  • Tàu Trans Summer (14/08/2013) chở 54.067 tấn quặng niken đắm tại Subaim, Indonesia.
  • Tàu Bulk Jupiter (02/01/2015) chở 46.400 tấn quặng bauxite từ Malaysia đi Trung Quốc bị nghiêng $45^\circ$ và chìm tại vùng biển Đông Nam Việt Nam.

Để khắc phục suy giảm ổn định, chương phân tích các giải pháp hiệu chỉnh:

  • Bơm thêm lượng nước ballast $P$ vào két: Biến thiên thế vững $\Delta GM$ được tính theo: $$\Delta GM = \frac{P}{D + P} \left( d_1 + \frac{\Delta d}{2} - GM - KG_{\text{chất lỏng}} \right)$$ Tùy thuộc vào giá trị $(d_1 + \Delta d/2 - GM)$ lớn hơn hay nhỏ hơn $KG_{\text{chất lỏng}}$, chiều cao thế vững sẽ tăng hoặc giảm.
  • Xả lượng nước ballast $P$ ra khỏi tàu: $$\Delta GM = \frac{P}{D - P} \left( d_1 - \frac{\Delta d}{2} - GM - KG_{\text{chất lỏng}} \right)$$
  • Nghiệp vụ đánh tẩy (san phẳng) mặt hàng: Quy định phân cấp góc nghỉ đối với hàng không gắn kết:
    • Góc nghỉ $\le 30^\circ$: Phải đánh tẩy toàn bộ bề mặt và chằng buộc cố định;
    • Góc nghỉ từ $30^\circ$ đến $35^\circ$: Chênh lệch chiều cao bề mặt đống hàng không vượt quá $1/10$ chiều rộng tàu hoặc không quá $0,15\text{ m}$;
    • Góc nghỉ $> 35^\circ$: Độ chênh lệch không vượt quá $1/10$ chiều rộng tàu hoặc không quá $0,20\text{ m}$.

Chương 2: Tính toán ổn định tàu hàng rời theo thời gian thực thông qua chu kỳ lắc ngang

Chương 2 thiết lập cơ sở toán học để xác định chu kỳ lắc ngang riêng và tính toán các chỉ số ổn định từ dữ liệu cảm biến đo góc nghiêng.

Chuyển động của tàu trên sóng được khảo sát trong hệ quy chiếu vật thể 6 bậc tự do thông qua hệ phương trình động lực học Euler: $$I_x \dot{\omega}_x - (I_y - I_z)\omega_y\omega_z = \sum M_x$$ $$I_y \dot{\omega}_y - (I_z - I_x)\omega_z\omega_x = \sum M_y$$ $$I_z \dot{\omega}_z - (I_x - I_y)\omega_x\omega_y = \sum M_z$$ trong đó $I_x, I_y, I_z$ là các mô men quán tính chính trung tâm của khối lượng tàu; $\omega_x, \omega_y, \omega_z$ là các thành phần vận tốc góc tương ứng với ba trục tọa độ liên kết $Gxyz$.

Tín hiệu đo độ nghiêng ngang thực tế $\theta(t)$ được số hóa thành chuỗi tín hiệu rời rạc: $$x(n) = r(n) + e(n) = r(n) + w(n) + d(n)$$ với $r(n)$ là dao động lắc ngang điều hòa của thân tàu, $w(n)$ là thành phần ngẫu nhiên do sóng kích động, và $d(n)$ là sai số thiết bị đo. Để bóc tách $r(n)$, tác giả triển khai quy trình xử lý tín hiệu:

  1. Làm mịn dữ liệu để triệt tiêu sai số đo $d(n)$.
  2. Phân chia chuỗi dữ liệu thành các cửa sổ quan sát có kích thước thời gian $V$ (phút), các cửa sổ dịch chuyển gối đầu nhau một khoảng trượt $v$ (giây) theo quy tắc $vx_i(k) = x(i \times v + k)$.
  3. Thực hiện phép biến đổi Fourier rời rạc (DFT) thông qua thuật toán biến đổi nhanh (FFT) trên từng đoạn: $$X_k = \sum_{n=0}^{N-1} x_n e^{-\frac{2\pi i}{N}kn}, \quad k = 0, \dots, N-1$$
  4. Xác định mật độ phổ năng lượng $PSD$. Do dao động của sóng $w(n)$ có tính ngẫu nhiên nên năng lượng phân tán rộng, trong khi dao động điều hòa của tàu $r(n)$ luôn tập trung năng lượng tại tần số riêng $f_r$. Tác giả tính toán phổ năng lượng trung bình đại diện cho cửa sổ thứ $j$: $$Power_j = \text{AVG}_{\forall i, f} (PSD_i^f)$$
  5. Áp dụng thuật toán dự đoán tần số dao động cho kỳ quan sát thứ $q+1$ dựa trên tập hợp phổ năng lượng qua $q$ cửa sổ liên tiếp bằng mô hình hồi quy tuyến tính: $$F_{q+1} = \text{predict}({ Power_j \mid j \in [1, q] })$$ hoặc thông qua hàm nội suy đa thức $P(x)$ với $F_{q+1} = P(q+1)$. Từ tần số dự báo, chu kỳ lắc ngang được xác định bằng $T = 1/F$.

Từ giá trị chu kỳ lắc ngang $T$ xác định được theo thời gian thực, hệ thống tiến hành tính toán chiều cao thế vững $G_0M$ theo công thức tổng quát trong Bộ luật IS Code: $$T = \frac{2 \times C \times B}{\sqrt{G_0M}} \implies G_0M = \left( \frac{2 \times C \times B}{T} \right)^2$$ trong đó $B$ là chiều rộng định hình của tàu, và $C$ là hệ số quán tính chuyển động lắc được chuẩn hóa theo công thức thực nghiệm của IMO: $$C = 0,373 + 0,023 \left(\frac{B}{d}\right) - 0,043 \left(\frac{L_{wl}}{100}\right)$$ với $d$ là mớn nước trung bình, $L_{wl}$ là chiều dài đường nước của tàu. Giá trị $G_0M$ tính được kết hợp cùng lượng giãn nước $D$ và bảng hoành giao $KN$ cho phép hệ thống tự động suy ra cao độ trọng tâm $KG_0 = KM - G_0M$, từ đó giải tích toàn bộ đường cong cánh tay đòn tĩnh $G_0Z(\theta) = KN(\theta) - KG_0 \times \sin\theta$ và kiểm tra tự động các tiêu chí diện tích diện tích dưới đường cong theo tiêu chuẩn IS Code 2008.

Chương 3: Thiết kế và xây dựng hệ thống thông báo ổn định tàu hàng rời theo thời gian thực

Chương 3 trình bày thiết kế kiến trúc tổng thể của hệ thống bao gồm hai phân hệ: phần cứng thu thập dữ liệu và phần mềm xử lý trung tâm:

  • Phân hệ phần cứng: Bao gồm các cảm biến đo góc nghiêng điện tử kỹ thuật số có độ chính xác cao được lắp đặt cố định tại trục dọc thân tàu, module chuyển đổi và truyền thông giao tiếp chuẩn công nghiệp để đưa dữ liệu góc nghiêng tức thời về máy tính tính toán trên buồng lái.
  • Phân hệ phần mềm: Được xây dựng dựa trên các module chức năng chính: module tiếp nhận và tiền xử lý tín hiệu thời gian thực; module phân tích phổ FFT và lọc nhiễu sóng biển; module thuật toán tính toán chu kỳ lắc riêng $T$ và chiều cao thế vững $G_0M$; module số hóa cơ sở dữ liệu hồ sơ tàu (bảng thủy tĩnh, bảng hoành giao $KN$); module giải tích đường cong $G_0Z$; và module giao diện đồ họa người dùng (GUI) hiển thị các thông số động lực học kèm hệ thống cảnh báo quang học - âm thanh khi các thông số vượt ngưỡng quy định của IS Code 2008.

Chương 4: Thực nghiệm hệ thống thông báo ổn định tàu hàng rời theo thời gian thực

Chương 4 mô tả quá trình triển khai thử nghiệm hệ thống trong môi trường vận hành thực tế trên tàu hàng rời:

  • Kiểm tra tính ổn định, độ nhạy và thời gian đáp ứng của cụm cảm biến đo nghiêng và đường truyền tín hiệu trong điều kiện rung động cơ học của buồng máy và sóng gió mặt boong.
  • Thu thập chuỗi số liệu dao động lắc ngang thực tế của tàu trong các điều kiện tải trọng khác nhau (chạy ballast và chở đầy tải).
  • Đánh giá độ hội tụ và tính chính xác của thuật toán lọc nhiễu phổ và thuật toán dự đoán tần số dao động thông qua việc so sánh kết quả tính toán $G_0M$ từ hệ thống với các giá trị thu được từ phương pháp tính toán thủy tĩnh truyền thống lập từ sơ đồ xếp hàng (Stowage Plan).
  • Hiệu chỉnh các tham số kích thước cửa sổ $V$ và độ dịch trượt $v$ nhằm tối ưu hóa độ chính xác và đảm bảo tốc độ phản hồi cảnh báo của hệ thống.

Kết quả và những đóng góp mới

Các đóng góp mới về mặt lý luận được luận án xác định:

  • Đề xuất phương pháp xác định chu kỳ dao động lắc ngang của tàu hàng rời từ tín hiệu của thiết bị đo góc nghiêng ngang theo thời gian thực dựa trên kỹ thuật phân tích phổ năng lượng và triệt tiêu thành phần kích động ngẫu nhiên của sóng biển.
  • Đề xuất thuật toán xác định chu kỳ lắc ngang kết hợp mô hình dự báo hồi quy/nội suy đa thức và thuật toán giải tích tự động các thông số ổn định ($G_0M$, đường cong $G_0Z$) theo thời gian thực.
  • Đề xuất mô hình kiến trúc và cấu trúc chức năng hoàn chỉnh cho hệ thống thông báo ổn định thời gian thực chuyên dụng trên tàu hàng rời.
  • Xây dựng phần mềm chuyên dụng có khả năng tính toán tự động, cập nhật liên tục đồ thị đường cong cánh tay đòn ổn định tĩnh $G_0Z$ và đưa ra các cảnh báo tự động khi các chỉ số ổn định suy giảm xuống dưới ngưỡng quy định của Bộ luật IS Code 2008.

Các đóng góp mới về mặt thực tiễn:

  • Chế tạo và tích hợp hoàn chỉnh hệ thống thông báo ổn định (kết hợp đồng bộ thiết bị phần cứng và phần mềm điều khiển), hoạt động như một công cụ hỗ trợ dẫn đường cho thuyền trưởng và sĩ quan hàng hải nhằm chủ động phòng ngừa các sự cố lật tàu do mất ổn định.
  • Cung cấp công cụ phục vụ công tác huấn luyện, đào tạo thực hành cho sĩ quan hàng hải và sinh viên chuyên ngành điều khiển tàu biển, nâng cao kỹ năng xử lý các tình huống suy giảm ổn định do dịch chuyển hoặc hóa lỏng hàng rời.

Kiến nghị và giải pháp đề xuất từ nghiên cứu:

  • Ứng dụng hệ thống giám sát thời gian thực trên các tàu chở hàng rời chuyên dụng, đặc biệt là các tàu thường xuyên vận chuyển các loại quặng có nguy cơ hóa lỏng cao như niken, bauxite, quặng sắt mịn.
  • Kết hợp dữ liệu cảnh báo thời gian thực của hệ thống với các quy trình điều chỉnh ballast đáy hoặc các phương án đánh tẩy bề mặt hầm hàng theo đúng góc nghỉ tiêu chuẩn khi phát hiện xu hướng suy giảm thế vững trong hành trình.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của nghiên cứu:

  • Phạm vi nghiên cứu của luận án giới hạn ở bài toán đánh giá ổn định nguyên vẹn trong điều kiện tàu nổi ở trạng thái cân bằng ban đầu, chưa xét đến các trường hợp tàu bị nghiêng ngang tĩnh trước khi hành trình.
  • Nghiên cứu chưa mở rộng khảo sát điều kiện ổn định khi tàu bị thủng vỏ hoặc ngập nước khoang cục bộ (ổn định tai nạn/hư hỏng).
  • Các thuật toán tính toán phụ thuộc vào độ chính xác của cơ sở dữ liệu thủy tĩnh số hóa ban đầu từ hồ sơ tàu do các xưởng đóng tàu cung cấp.

Hướng nghiên cứu phát triển tiếp theo:

  • Mở rộng thuật toán để đánh giá trạng thái ổn định của tàu trong điều kiện tồn tại góc nghiêng ngang tĩnh ban đầu kết hợp với hiện tượng mặt thoáng chất lỏng phức tạp.
  • Tích hợp thêm các bài toán đánh giá ổn định tai nạn theo quy chuẩn của SOLAS 74 Chương II-1 lên cùng nền tảng phần mềm thời gian thực.

Giá trị tham khảo

Luận án có giá trị tham khảo chuyên môn cho các nhóm đối tượng:

  • Sĩ quan quản lý và thuyền trưởng tàu hàng rời: Cung cấp phương pháp luận và công cụ trực quan để theo dõi liên tục trạng thái ổn định thực tế của tàu trên biển thay vì chỉ dựa vào các bảng tính toán tĩnh lập trước chuyến đi.
  • Giảng viên, nghiên cứu sinh và sinh viên chuyên ngành Điều khiển tàu biển, Kỹ thuật tàu thủy và Khoa học hàng hải: Tài liệu tham khảo về phương pháp mô hình hóa chuyển động lắc tàu, ứng dụng xử lý tín hiệu số (FFT, PSD) trong kỹ thuật hàng hải, và phân tích các cơ chế tai nạn do hóa lỏng hàng rời.
  • Các công ty khai thác tàu và cơ quan quản lý an toàn hàng hải: Tham khảo mô hình kiến trúc hệ thống để định hướng trang bị các giải pháp giám sát an toàn chủ động, hỗ trợ thẩm định và phòng ngừa rủi ro vận tải biển.

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

1. Vì sao việc kiểm tra ổn định tàu trước khi rời cảng không đủ đảm bảo an toàn cho tàu hàng rời trong suốt hành trình?
Trước khi rời cảng, tính ổn định chỉ được tính toán dựa trên trạng thái tĩnh của hàng hóa và chất dằn. Khi hành hải, các yếu tố ngoại cảnh và nội tại như sự xô dịch của hàng rời, sự hóa lỏng của quặng do rung chấn, sự hình thành mặt thoáng khi điều chỉnh nước dằn, hoặc khiếm khuyết vỏ tàu có thể làm trọng tâm tàu nâng cao hoặc dịch chuyển sang mạn, làm suy giảm nghiêm trọng chiều cao thế vững mà sĩ quan không thể phát hiện kịp thời nếu không có công cụ đo thời gian thực.

2. Tiêu chuẩn Bộ luật IS Code 2008 quy định các ngưỡng an toàn cụ thể nào đối với đường cong cánh tay đòn ổn định tĩnh $G_0Z$ và chiều cao thế vững $G_0M$?
Theo IS Code 2008, diện tích dưới đường cong $G_0Z$ phải đạt tối thiểu $0,055\text{ m-rad}$ đến góc $30^\circ$, tối thiểu $0,090\text{ m-rad}$ đến $40^\circ$ (hoặc góc ngập nước), và diện tích đoạn giữa $30^\circ$ và $40^\circ$ tối thiểu là $0,030\text{ m-rad}$. Độ lớn cánh tay đòn $G_0Z$ phải đạt ít nhất $0,20\text{ m}$ tại góc nghiêng $\ge 30^\circ$ và đạt cực đại tại góc $\ge 25^\circ$. Chiều cao thế vững ban đầu đã hiệu chỉnh mặt thoáng $G_0M$ không được nhỏ hơn $0,15\text{ m}$.

3. Cơ chế loại bỏ nhiễu sóng biển để trích xuất chu kỳ lắc riêng của tàu trong luận án được thực hiện như thế nào?
Nhiễu do sóng biển $w(n)$ có tính chất ngẫu nhiên, năng lượng phân tán rộng trên nhiều tần số. Ngược lại, dao động riêng của tàu $r(n)$ có tính điều hòa nên năng lượng tập trung tại tần số dao động chính. Bằng cách chia tín hiệu thành các cửa sổ trượt $V$ có độ dịch $v$ chồng đè, áp dụng biến đổi FFT và tính mật độ phổ năng lượng trung bình ($Power_j$), hệ thống duy trì được mức năng lượng tại tần số có xác suất xuất hiện cao nhất, từ đó triệt tiêu năng lượng của sóng và xác định chính xác tần số riêng qua thuật toán hồi quy hoặc nội suy.

4. Dựa vào công thức nào hệ thống có thể tính toán được chiều cao thế vững $G_0M$ từ chu kỳ lắc ngang đo được?
Hệ thống sử dụng công thức thực nghiệm được chuẩn hóa của IMO: $$G_0M = \left( \frac{2 \times C \times B}{T} \right)^2$$ với $T$ là chu kỳ lắc ngang thu được từ cảm biến và thuật toán xử lý tín hiệu; $B$ là chiều rộng định hình của tàu; $C$ là hệ số quán tính lắc được xác định dựa trên thông số mớn nước $d$ và chiều dài đường nước $L_{wl}$ theo công thức: $$C = 0,373 + 0,023(B/d) - 0,043(L_{wl}/100)$$

5. Hiện tượng hóa lỏng của hàng rời dẫn đến lật chìm tàu thông qua các giai đoạn cơ học nào?
Dưới tác động rung lắc liên tục của sóng, hàng có độ ẩm vượt ngưỡng bị hóa lỏng tạo ra mặt thoáng tự do, làm suy giảm chiều cao thế vững $GM$ và kéo dài chu kỳ lắc. Do tính nhớt cao, khi tàu nghiêng, khối hàng sệt chảy dồn sang mạn thấp nhưng không hồi lưu đều đặn khi tàu lắc ngược lại. Khối lượng hàng tích tụ làm dịch chuyển trọng tâm tàu sang mạn nghiêng, đẩy góc nghiêng tĩnh tăng dần (có thể vượt $30^\circ$ đến $45^\circ$). Khi gặp sóng gió bất lợi, tàu mất hoàn toàn mô men hồi phục và bị lật úp.

Kết luận

Luận án đã giải quyết bài toán cấp thiết trong kỹ thuật khai thác tàu hàng rời bằng việc xây dựng phương pháp luận và hệ thống hoàn chỉnh nhằm giám sát các thông số ổn định nguyên vẹn theo thời gian thực. Bằng sự kết hợp giữa mô hình động lực học chuyển động lắc, kỹ thuật phân tích phổ năng lượng FFT để lọc nhiễu sóng biển và thuật toán tính toán thế vững tự động, nghiên cứu cho phép chuyển đổi phương thức kiểm soát an toàn từ kiểm tra tĩnh định kỳ sang giám sát động liên tục. Kết quả thiết kế hệ thống phần cứng và phần mềm mang lại công cụ hỗ trợ cho sĩ quan hàng hải trong việc nhận diện sớm nguy cơ mất ổn định, đóng góp cơ sở khoa học và giải pháp kỹ thuật vào công tác đảm bảo an toàn sinh mạng và phương tiện trong vận tải biển.