Tổng quan nghiên cứu

Trong ngành công nghiệp hóa chất hiện đại, các sự cố cháy nổ bắt nguồn từ phản ứng tỏa nhiệt mất kiểm soát chiếm tới hơn 60% tổng số các tai nạn nghiêm trọng tại nhà xưởng và kho chứa. Khi tốc độ sinh nhiệt của phản ứng hóa học tăng theo hàm số mũ theo định luật Arrhenius trong khi khả năng tản nhiệt của hệ thống làm mát chỉ tăng tuyến tính theo nhiệt độ, hiện tượng tích tụ nhiệt đột ngột sẽ dẫn đến nổ nhiệt cục bộ. Luận văn thạc sĩ khoa học chuyên ngành Hóa lý thuyết và Hóa lý tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội đã tập trung giải quyết bài toán cốt lõi: xác định chính xác các thông số nhiệt động và động học của phản ứng tỏa nhiệt nhằm dự báo an toàn và tối ưu hóa quy trình công nghệ.

Mục tiêu cụ thể của công trình là đo đạc các đặc trưng nhiệt và xây dựng mô hình động học hoàn chỉnh cho phản ứng lưu hóa cao su Etylen Propylen Dien Monome (EPDM 3666) có phối trộn 0,5 phần khối lượng lưu huỳnh và hệ không sử dụng lưu huỳnh. Nghiên cứu được thực hiện tại Phòng Thí nghiệm An toàn Hóa chất thuộc Viện Nghiên cứu Khoa học Kỹ thuật Bảo hộ Lao động trong giai đoạn từ năm 2010 đến năm 2012. Phạm vi khảo sát nhiệt độ trải rộng từ 30°C đến 310°C dưới các tốc độ gia nhiệt biến thiên từ 5 K/phút đến 20 K/phút. Kết quả nghiên cứu cung cấp cơ sở dữ liệu quan trọng giúp kiểm soát quá trình truyền nhiệt, nâng cao độ bền vật liệu và giảm thiểu rủi ro quá nhiệt trong công nghệ sản xuất polyme với khả năng phòng ngừa sự cố đạt trên 95%.

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 được xây dựng dựa trên nền tảng vững chắc của nhiệt động lực học và động hóa học hình thức:

  1. Phương trình Arrhenius và Năng lượng hoạt hóa: Xác lập mối quan hệ phụ thuộc giữa hằng số tốc độ phản ứng với nhiệt độ tuyệt đối thông qua biểu thức $k = A \exp(-E/RT)$, trong đó năng lượng hoạt hóa đại diện cho rào cản năng lượng tối thiểu để chuyển hóa chất phản ứng thành sản phẩm.
  2. Phương trình biến đổi pha bán thực nghiệm Avrami: Mô tả động học quá trình khâu mạch và kết tinh pha rắn theo hàm thời gian $x = 1 - \exp[-(kt)^n]$, tạo tiền đề cho việc phân tích động học bất đẳng nhiệt.
  3. Định luật Hess về nhiệt phản ứng: Xác định sự thay đổi entanpy phản ứng ($\Delta H < 0$) thông qua tích phân diện tích đường cong nhiệt lượng quét vi sai so với đường nền.
  4. Hệ thống các khái niệm then chốt: Năng lượng hoạt hóa ($E$), thừa số trước hàm mũ ($\log A$), bậc phản ứng ($n$), nhiệt độ khởi phát tỏa nhiệt ($T_{onset}$), và biến thiên entanpy lưu hóa ($\Delta H$).

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

Nguồn dữ liệu thực nghiệm được thu thập trực tiếp thông qua hệ thống đo nhiệt lượng vi sai quét dòng nhiệt tiên tiến DSC 204 F1 Phoenix của hãng NETZSCH (Đức). Các mẫu cao su EPDM 3666 chưa lưu hóa được cán luyện và phối trộn tại Trung tâm Nghiên cứu Vật liệu Polyme thuộc Trường Đại học Bách Khoa Hà Nội theo đơn phối liệu chuẩn gồm 100 phần khối lượng cao su nền cùng các chất hoạt hóa kẽm oxit (5 phần khối lượng), axit stearic (1 phần khối lượng) và hệ xúc tiến chuyên dụng.

  • Quy mô và phương pháp chọn mẫu: Mẫu đo có khối lượng tối ưu từ 3 mg đến 5 mg, được cắt nhỏ, dàn phẳng đáy chén nhôm đường kính 6 mm và dập nắp kín nhằm loại trừ hiện tượng bay hơi của phụ gia hữu cơ có điểm sôi thấp. Khối lượng mẫu được kiểm soát chính xác đến ±0,1 mg trên cân phân tích Shimadzu AUY 220. Phương pháp lấy mẫu đại diện ngẫu nhiên từ các tấm cao su cán dày 3 mm đảm bảo tính đồng nhất cấu trúc.
  • Quy trình phân tích: Phép đo quét nhiệt được tiến hành trong môi trường khí trơ Nitơ có độ tinh khiết 99,999% với lưu lượng dòng 20 ml/phút trong buồng phản ứng 250 ml. Tốc độ quét nhiệt bất đẳng nhiệt được thiết lập lần lượt tại 4 mức: 5, 10, 15 và 20 K/phút. Việc lựa chọn phương pháp phân tích nhiệt bất đẳng nhiệt cho phép rút ngắn hơn 80% thời gian thực nghiệm so với phương pháp đẳng nhiệt truyền thống, đồng thời ghi nhận toàn diện các biến đổi hóa lý theo dải nhiệt độ liên tục.
  • Xử lý số liệu: Dữ liệu dạng ASCII từ phần mềm Proteus Analysis được trích xuất sang phần mềm NETZSCH Thermokinetics để phân tích theo mô hình tự do (Friedman, Ozawa-Flynn-Wall, ASTM E698) và mô hình cơ sở phi tuyến đa bước lựa chọn từ 77 mô hình động học lý thuyết.

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

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

Quá trình phân tích thực nghiệm và mô phỏng động học nhiệt đã làm sáng tỏ nhiều đặc tính quan trọng của hệ phản ứng:

  • Đặc trưng nhiệt động học của quá trình lưu hóa: Giản đồ DSC ghi nhận phản ứng lưu hóa cao su EPDM bằng lưu huỳnh diễn ra ở vùng nhiệt độ thấp hơn đáng kể, khởi phát từ khoảng 140°C và đạt đỉnh cực đại ở khoảng 165°C đến 185°C tùy theo tốc độ quét nhiệt. Hiệu ứng tỏa nhiệt đo được có biến thiên entanpy $\Delta H$ đạt xấp xỉ 42,6 J/g. Ngược lại, mẫu không chứa lưu huỳnh chỉ bắt đầu phản ứng ở nhiệt độ trên 190°C và đạt cực đại tại 235°C đến 255°C với nhiệt lượng tỏa ra chỉ đạt khoảng 21,8 J/g, thấp hơn gần 49% so với hệ có lưu huỳnh.
  • Giá trị năng lượng hoạt hóa: Phân tích mô hình tự do Ozawa-Flynn-Wall (OFW) và Friedman chỉ ra rằng hệ lưu hóa có lưu huỳnh có năng lượng hoạt hóa $E$ dao động từ 112 kJ/mol đến 125 kJ/mol. Trong khi đó, hệ không sử dụng lưu huỳnh đòi hỏi rào cản năng lượng cao hơn rõ rệt, đạt mức 158 kJ/mol đến 172 kJ/mol, tức cao hơn khoảng 37% đến 40%.
  • Cơ chế phản ứng theo mô hình cơ sở: Hệ thống hồi quy phi tuyến tính đa bước đã xác định phản ứng lưu hóa EPDM diễn ra theo cơ chế phản ứng hai giai đoạn nối tiếp. Giai đoạn thứ nhất tuân theo quy luật phản ứng bậc $n$ với $n_1 \approx 1,05$, tiếp nối bởi giai đoạn thứ hai với $n_2 \approx 0,82$. Hệ số tương quan hồi quy đạt mức $r = 0,9992$, khẳng định mức độ trùng khớp rất cao giữa mô hình toán học và thực nghiệm.
  • Tối ưu hóa thời gian chu trình công nghệ: Từ các thông số động học thu được, phần mềm đã mô phỏng việc tối ưu hóa biên dạng nhiệt độ, giúp rút ngắn thời gian đóng rắn hoàn toàn từ 40 phút xuống còn 18 phút ở chế độ gia nhiệt đa bậc, tiết kiệm tới 55% thời gian gia công mà vẫn đảm bảo độ chuyển hóa đạt trên 98%.

Thảo luận kết quả

Sự khác biệt rõ rệt về mặt nhiệt động học giữa hai hệ mẫu bắt nguồn từ cơ chế hóa học của liên kết ngang. Khi có mặt 0,5 phần khối lượng lưu huỳnh kết hợp cùng hệ xúc tiến dithiocarbamate và thiazole, các gốc tự do polysunfua được hình thành nhanh chóng ở nhiệt độ thấp, tấn công vào các vị trí nối đôi không no của phân tử dien trong mạch EPDM. Quá trình này giải phóng năng lượng liên kết lớn hơn và diễn ra nhanh hơn so với cơ chế tự khâu mạch nhiệt thuần túy của mẫu không có lưu huỳnh vốn đòi hỏi bẻ gãy các liên kết cacbon kém hoạt động hơn.

Khi so sánh với các công bố quốc tế, kết quả này hoàn toàn tương đồng với nghiên cứu của tác giả J. Opfermann trên hệ cao su EPDM lưu hóa bằng peroxit có năng lượng hoạt hóa giai đoạn một là 180 kJ/mol, cũng như tương thích với các phân tích nhiệt của A. Germain về phân hủy thuốc trừ sâu photpho hữu cơ với nhiệt tỏa ra từ 730 J/g đến 1100 J/g. Trên đồ thị biểu diễn dòng nhiệt theo nhiệt độ, các đỉnh phản ứng có xu hướng dịch chuyển tịnh tiến về phía nhiệt độ cao hơn khi tốc độ quét tăng từ 5 K/phút lên 20 K/phút do quán tính truyền nhiệt của vật liệu mẫu. Đồ thị chuyển hóa theo thời gian thể hiện rõ đường cong dạng chữ S đặc trưng cho các phản ứng khâu mạch tự xúc tác cục bộ, chứng minh tính tin cậy tuyệt đối của phương pháp vi sai quét nhiệt DSC.

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

Dựa trên kết quả đo đạc chính xác và mô hình động học đã thiết lập, luận văn đưa ra 4 nhóm giải pháp thực tiễn:

  1. Chuẩn hóa quy trình kiểm soát an toàn nhiệt đầu vào: Áp dụng phương pháp quét nhiệt DSC 3 tốc độ (5, 10, 15 K/phút) để sàng lọc nhiệt lượng tỏa ra của 100% các lô nguyên liệu phối trộn trước khi đưa vào dây chuyền ép khuôn lớn. Chủ thể thực hiện là Bộ phận Quản lý Chất lượng (QC) và Phòng Thí nghiệm R&D tại các nhà máy sản xuất cao su kỹ thuật, hoàn thành trong vòng 3 tháng đầu triển khai, hướng tới mục tiêu duy trì tỷ lệ sai hỏng nhiệt dưới 0,2%.
  2. Thiết kế hệ thống giải nhiệt công nghiệp chủ động: Nâng cấp công suất bộ trao đổi nhiệt tại các bồn phản ứng và máy ép thủy lực sao cho công suất làm mát danh định luôn lớn hơn tốc độ tỏa nhiệt cực đại của phản ứng ít nhất 1,5 lần. Bộ phận Kỹ thuật Cơ điện và An toàn Lao động chịu trách nhiệm thực thi trong thời hạn 6 tháng, đảm bảo triệt tiêu hoàn toàn nguy cơ mất kiểm soát nhiệt và đạt chỉ tiêu 0 vụ sự cố cháy nổ.
  3. Chuyển đổi biên dạng gia nhiệt sang chế độ tối ưu hóa đa bước: Thay thế quy trình gia nhiệt đẳng nhiệt 160°C cố định bằng chương trình nhiệt 2 giai đoạn (gia nhiệt nhanh lên 120°C sau đó nâng có kiểm soát lên 155°C theo mô hình Thermokinetics). Kỹ sư công nghệ chế biến thực hiện áp dụng trong vòng 2 tháng, nhằm mục tiêu cắt giảm 35% thời gian lưu hóa khuôn và giảm 18% chi phí tiêu hao năng lượng điện.
  4. Xây dựng cơ sở dữ liệu số hóa về thông số động học hóa chất: Phối hợp giữa Viện Nghiên cứu Khoa học Kỹ thuật Bảo hộ Lao động và các trường đại học khối kỹ thuật để thành lập ngân hàng dữ liệu an toàn nhiệt cho hơn 50 loại vật liệu polyme và hóa chất công nghiệp phổ biến theo tiêu chuẩn ISO/IEC 17025 trong lộ trình 12 tháng.

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

Nội dung và kết quả nghiên cứu mang lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng trọng tâm:

  1. Kỹ sư công nghệ polyme và gia công cao su: Sử dụng dữ liệu thực nghiệm về đơn phối liệu EPDM 3666 và các thông số $E$, $\log A$, $n$ để thiết lập chính xác chế độ lưu hóa ép khuôn cho các sản phẩm gioăng kính ô tô, ống dẫn nhiệt động cơ và màng chống thấm công trình.
  2. Chuyên viên an toàn hóa chất và bảo hộ lao động: Ứng dụng quy trình đánh giá mối nguy nhiệt bằng thiết bị DSC để tính toán các kịch bản nổ nhiệt, từ đó xây dựng phương án ứng phó sự cố khẩn cấp và thiết kế van xả áp an toàn cho bồn công nghệ.
  3. Cán bộ nghiên cứu, giảng viên và học viên sau đại học: Khai thác các phương pháp toán học xử lý dữ liệu động học nhiệt phi tuyến (mô hình OFW, Friedman, ASTM E698) làm tài liệu giảng dạy và nghiên cứu chuyên sâu trong chuyên ngành Hóa lý và Khoa học Vật liệu.
  4. Ban giám đốc và nhà quản lý sản xuất tại các doanh nghiệp hóa chất: Ứng dụng mô hình tối ưu hóa biên dạng nhiệt để rút ngắn thời gian chiếm dụng thiết bị, tăng năng suất chuyền ép lên 25% và giảm thiểu tiêu hao năng lượng vận hành.

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

Tại sao kỹ thuật DSC kết hợp phần mềm động học nhiệt lại vượt trội hơn phương pháp thử nghiệm truyền thống?
Kỹ thuật DSC sử dụng lượng mẫu cực nhỏ từ 3 đến 5 mg, giúp loại bỏ nguy cơ cháy nổ trong phòng thí nghiệm. Việc kết hợp phần mềm NETZSCH Thermokinetics cho phép mô phỏng chính xác hàng trăm kịch bản nhiệt độ phức tạp mà không cần thực hiện nhiều mẻ đo tốn kém, tiết kiệm hơn 70% thời gian nghiên cứu.

Sự có mặt của lưu huỳnh làm thay đổi động học phản ứng lưu hóa EPDM như thế nào?
Phối trộn 0,5 phần khối lượng lưu huỳnh giúp giảm năng lượng hoạt hóa từ khoảng 165 kJ/mol xuống xấp xỉ 118 kJ/mol. Sự sụt giảm rào cản năng lượng này giúp phản ứng khởi phát sớm hơn gần 50°C và làm tăng tốc độ tạo liên kết không gian 3 chiều lên hơn 40%.

Mô hình động học tự do khác biệt gì so với mô hình cơ sở trong phân tích nhiệt?
Mô hình tự do như Friedman và OFW ước tính nhanh năng lượng hoạt hóa theo từng mức độ chuyển hóa mà không cần giả định trước cơ chế. Ngược lại, mô hình cơ sở giải hệ phương trình vi phân cho 77 dạng phản ứng nối tiếp hay song song, xác định chính xác bậc phản ứng và thừa số tần số.

Làm thế nào để ứng dụng thông số động học vào việc tối ưu hóa năng lượng sản xuất?
Dựa trên phương trình động học đã xác lập, kỹ sư có thể lập trình đường cong gia nhiệt phi tuyến trên máy ép tự động. Quy trình này giúp phản ứng đạt độ chuyển hóa 98% trong thời gian ngắn nhất mà không cần duy trì nhiệt độ cao liên tục, giúp tiết kiệm 18% điện năng.

Phương pháp này có thể áp dụng để đánh giá nguy cơ cháy nổ cho các hóa chất khác không?
Phương pháp này hoàn toàn có thể áp dụng hiệu quả cho hầu hết các hợp chất hữu cơ dễ phân hủy nhiệt như peroxit, thuốc nổ hữu cơ, phân bón và nhựa tổng hợp. Việc xác định nhiệt độ bắt đầu tỏa nhiệt $T_{onset}$ và biến thiên entanpy $\Delta H$ cho phép phân loại mức độ nguy hiểm nhiệt theo các tiêu chuẩn quốc tế.

Kết luận

Nghiên cứu xác định các thông số động học của phản ứng tỏa nhiệt bằng kỹ thuật nhiệt lượng vi sai quét DSC đã mang lại những đóng góp khoa học và thực tiễn vững chắc:

  • Xác định thành công toàn bộ hệ thống thông số nhiệt động của phản ứng lưu hóa cao su EPDM 3666 với biến thiên entanpy $\Delta H$ đạt 42,6 J/g đối với hệ có lưu huỳnh và 21,8 J/g đối với hệ không lưu huỳnh.
  • Khẳng định tính ưu việt của việc kết hợp mô hình động học tự do (OFW, Friedman) và mô hình cơ sở phi tuyến đa bước, xác định chính xác năng lượng hoạt hóa của hệ có lưu huỳnh đạt 112 đến 125 kJ/mol.
  • Chứng minh phản ứng khâu mạch cao su EPDM diễn ra theo cơ chế 2 giai đoạn nối tiếp với hệ số tương quan tuyến tính rất cao đạt $r = 0,9992$.
  • Thiết lập thành công biên dạng nhiệt độ tối ưu hóa, hỗ trợ rút ngắn hơn 50% thời gian chu trình lưu hóa trong sản xuất công nghiệp.
  • Hoàn thiện quy trình chuẩn hóa ứng dụng thiết bị DSC 204 F1 Phoenix phục vụ công tác đánh giá mối nguy hiểm nổ nhiệt tại Việt Nam.

Trong giai đoạn 12 tháng tới, định hướng tiếp theo là mở rộng phương pháp nghiên cứu động học nhiệt sang các hệ vật liệu compozit polyme nền cao su thiên nhiên và nhựa nhiệt rắn, đồng thời tiến hành thử nghiệm kiểm chứng trên các thiết bị phản ứng quy mô pilot 50 lít. Các doanh nghiệp và viện nghiên cứu quan tâm có thể liên hệ hợp tác để chuyển giao quy trình phân tích và tối ưu hóa nhiệt động lực học tiên tiến này.