Tổng quan nghiên cứu

Trong ngành công nghiệp hóa chất hiện đại, hơn 15% đến 20% các sự cố cháy nổ nhà xưởng nghiêm trọng trên quy mô toàn cầu xuất phát từ các phản ứng tỏa nhiệt mất kiểm soát. Hiện tượng này xảy ra do tốc độ tỏa nhiệt sinh ra theo hàm mũ nhiệt độ theo định luật Arrhenius, trong khi khả năng giải nhiệt của hệ thống làm mát bồn phản ứng chỉ tăng tuyến tính. Khi năng lượng giải phóng đột ngột vượt quá công suất làm mát, nguy cơ nổ nhiệt và phá hủy thiết bị là không thể tránh khỏi. Để giải quyết triệt để vấn đề này, việc xác định chính xác các thông số nhiệt động và động hóa học đóng vai trò then chốt trong việc đánh giá mức độ nguy hiểm và thiết kế quy trình vận hành an toàn.

Mục tiêu cụ thể của luận văn là xây dựng quy trình chuẩn xác định các thông số động học phản ứng tỏa nhiệt (bao gồm năng lượng hoạt hóa, thừa số trước hàm mũ, bậc phản ứng và mô hình động học) bằng kỹ thuật nhiệt lượng vi sai quét kết hợp phần mềm động học nhiệt chuyên dụng. Phạm vi nghiên cứu tập trung khảo sát phản ứng lưu hóa hệ cao su EPDM 3666 (có và không có lưu huỳnh) được thực hiện tại Trung tâm Nghiên cứu Vật liệu Polyme và Viện Khoa học Kỹ thuật Bảo hộ Lao động trong giai đoạn 2010 đến 2012. Ý nghĩa thực tiễn của công trình thể hiện qua việc tối ưu hóa biên dạng nhiệt độ gia công, giúp rút ngắn khoảng 18% đến 25% chu kỳ sản xuất, đồng thời cắt giảm hơn 90% nguy cơ xảy ra sự cố quá nhiệt trong quy trình chế biến vật liệu đàn hồi.

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 trên nền tảng vững chắc của hai hệ thống lý thuyết cốt lõi: thuyết động hóa học hình thức kinh điển phát triển từ các công trình của Arrhenius và Van 't Hoff từ những năm 1880, kết hợp với phương trình bán thực nghiệm Avrami mô tả động học các quá trình biến đổi pha trong chất rắn. Khung lý thuyết này làm rõ bản chất của tốc độ phản ứng, sự phụ thuộc của hằng số tốc độ vào nhiệt độ tuyệt đối và định luật tác dụng khối lượng của Guldberg và Waage.

Mô hình nghiên cứu triển khai theo 2 hướng tiếp cận: mô hình động học phi đẳng nhiệt và mô hình đẳng nhiệt. Trong phân tích động học nhiệt, 3 mô hình tự do được áp dụng gồm mô hình Friedman vi phân, mô hình Ozawa - Flynn - Wall (OFW) tích phân và mô hình phân tích theo tiêu chuẩn ASTM E698. Đi cùng với đó là mô hình cơ sở dựa trên việc giải hệ phương trình vi phân phi tuyến tính để nhận dạng cơ chế phản ứng đa giai đoạn.

Năm khái niệm chuyên ngành cốt lõi chi phối toàn bộ phân tích gồm:

  1. Năng lượng hoạt hóa ($E$): Hàng rào năng lượng tối thiểu mà phân tử phải vượt qua để phản ứng hóa học xảy ra.
  2. Thừa số trước hàm mũ ($\log A$): Đại lượng phản ánh tần số va chạm hiệu quả giữa các tiểu phân phản ứng.
  3. Biến thiên enthalpy ($\Delta H$): Lượng nhiệt giải phóng hoặc hấp thụ trong điều kiện áp suất không đổi.
  4. Nhiệt độ khởi phát ($T_{onset}$): Điểm nhiệt độ bắt đầu xuất hiện sự sai lệch của thông lượng nhiệt khỏi đường nền.
  5. Hiệu ứng tự xúc tác: Hiện tượng sản phẩm sinh ra đóng vai trò xúc tác làm gia tăng đột biến tốc độ phản ứng.

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

Nguồn dữ liệu thực nghiệm được thu thập từ 16 mẫu thử cao su EPDM 3666 có khối lượng tiêu chuẩn từ 5 mg đến 10 mg. Phương pháp chọn mẫu áp dụng là lấy mẫu ngẫu nhiên phân tầng từ mẻ cán luyện cao su trên máy cán hai trục tại Trung tâm Nghiên cứu Vật liệu Polyme, đảm bảo độ đồng nhất phụ gia và chất lưu hóa đạt trên 98,5%.

Lý do lựa chọn phương pháp phân tích: Nghiên cứu sử dụng kỹ thuật nhiệt lượng vi sai quét dòng nhiệt trên thiết bị NETZSCH DSC 204 F1 Phoenix. Thiết bị trang bị đầu dò với cặp nhiệt điện Alumel - Chromel có độ nhạy nhiệt cực cao, cho phép đo trực tiếp thông lượng nhiệt theo thời gian thực. Phương pháp quét nhiệt phi đẳng nhiệt được ưu tiên lựa chọn thay vì phương pháp đẳng nhiệt truyền thống vì giúp rút ngắn thời gian thí nghiệm từ hơn 10 giờ xuống dưới 3 giờ mỗi chuỗi đo, đồng thời loại bỏ sai số do phản ứng sớm trong giai đoạn gia nhiệt ban đầu.

Quá trình thực nghiệm tiến hành ở 4 tốc độ quét nhiệt tuyến tính gồm 5, 10, 15 và 20 K/phút trong môi trường khí trơ $N_2$ với lưu lượng ổn định 30 ml/phút, dải nhiệt độ khảo sát từ 30 °C đến 300 °C. Dữ liệu sau đó được trích xuất và xử lý hồi quy phi tuyến đa biến trên phần mềm NETZSCH Thermokinetics. Toàn bộ timeline nghiên cứu kéo dài 24 tháng, bao gồm 4 tháng chuẩn bị đơn phối liệu, 10 tháng phân tích DSC và 10 tháng mô phỏng tối ưu hóa.

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

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

Thực nghiệm đo quét nhiệt vi sai DSC và tính toán động học trên hệ cao su EPDM 3666 đã mang lại 4 phát hiện khoa học quan trọng:

Thứ nhất, phản ứng lưu hóa cao su EPDM thể hiện bản chất tỏa nhiệt mạnh với diện tích đỉnh nhiệt lượng phụ thuộc chặt chẽ vào thành phần phối liệu. Đối với mẫu cao su không sử dụng lưu huỳnh, biến thiên enthalpy $\Delta H$ đạt mức -15,4 J/g, trong khi mẫu sử dụng lưu huỳnh ghi nhận $\Delta H$ lên tới -42,8 J/g, tương đương mức tăng năng lượng tỏa ra hơn 177,9%.

Thứ hai, tốc độ gia nhiệt $\beta$ tác động rõ rệt đến vị trí các điểm đặc trưng của giản đồ DSC. Khi tăng tốc độ quét từ 5 K/phút lên 20 K/phút, nhiệt độ đỉnh cực đại ($T_{max}$) của hệ lưu hóa bằng lưu huỳnh dịch chuyển từ 182,5 °C lên 208,3 °C (tăng 25,8 °C), minh chứng cho độ trễ truyền nhiệt đặc trưng trong phân tích nhiệt chất rắn.

Thứ ba, việc xác định sơ bộ thông số động học bằng các mô hình tự do cho thấy năng lượng hoạt hóa trung bình của hệ có lưu huỳnh đạt từ 95,2 kJ/mol (theo mô hình Friedman) đến 101,4 kJ/mol (theo mô hình OFW). Ngược lại, hệ không lưu huỳnh đòi hỏi năng lượng hoạt hóa cao hơn đáng kể, dao động từ 128,4 kJ/mol đến 145,1 kJ/mol (cao hơn khoảng 34,9% đến 43,1%).

Thứ tư, phân tích theo mô hình cơ sở khẳng định phản ứng lưu hóa EPDM diễn ra theo cơ chế 2 giai đoạn nối tiếp có tính tự xúc tác với hệ số tương quan thống kê $R^2$ vượt mức 0,9985, bác bỏ giả định phản ứng đơn bậc truyền thống.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự suy giảm năng lượng hoạt hóa khi có mặt lưu huỳnh là do lưu huỳnh và các chất xúc tiến tạo thành các phức chất trung gian hoạt động, làm suy yếu các liên kết đôi trên mạch polyme EPDM, từ đó hạ thấp hàng rào năng lượng phản ứng. Điều này giúp phản ứng lưu hóa khởi phát ở nhiệt độ thấp hơn khoảng 20 °C đến 25 °C so với hệ không có lưu huỳnh, nhưng lại làm tăng mật độ dòng nhiệt tức thời, tiềm ẩn nguy cơ nổ nhiệt cục bộ nếu không giải nhiệt kịp thời.

So sánh với các công bố của Francis Stoessel về đánh giá an toàn nhiệt, việc áp dụng mô hình động học phi đẳng nhiệt kết hợp phần mềm NETZSCH Thermokinetics giúp giảm độ lệch chuẩn của thừa số trước hàm mũ $\log A$ từ 8,5% xuống dưới 1,2% so với việc chỉ tính toán dựa trên tiêu chuẩn ASTM E698 đơn điểm.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua 2 dạng công cụ chính: giản đồ DSC đa tốc độ quét biểu diễn sự phụ thuộc của thông lượng nhiệt vào nhiệt độ, và biểu đồ chuyển hóa nồng độ các cấu tử theo thời gian ở các mức nhiệt độ đẳng nhiệt từ 140 °C đến 200 °C. Cùng với đó, bảng phân loại mức độ nguy hiểm nhiệt dựa trên biến thiên enthalpy xác định hệ phản ứng EPDM thuộc nhóm nguy cơ nhiệt cấp độ 1 (thấp), hoàn toàn có thể kiểm soát an toàn nếu áp dụng đúng biên dạng gia nhiệt tối ưu.

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

Nhằm chuyển hóa các kết quả nghiên cứu lý thuyết vào thực tiễn sản xuất công nghiệp và kiểm soát an toàn, 4 giải pháp cụ thể được đề xuất:

  1. Áp dụng quy trình tối ưu hóa biên dạng nhiệt độ phi tuyến tính trong công nghệ lưu hóa cao su EPDM nhằm duy trì tốc độ sinh nhiệt dưới ngưỡng 1,5 W/kg, rút ngắn 20% tổng thời gian gia công khuôn ép trong vòng 6 tháng tới (do Đội ngũ Kỹ sư Công nghệ tại các nhà máy sản xuất cao su kỹ thuật thực hiện).
  2. Nâng cấp thiết bị và trang bị bổ sung module phần mềm NETZSCH Thermokinetics cho ít nhất 3 phòng thí nghiệm trọng điểm về an toàn hóa chất nhằm chuẩn hóa công tác xác định thông số động học trong lộ trình 12 tháng (do Viện Khoa học Kỹ thuật Bảo hộ Lao động chủ trì).
  3. Thiết lập bộ dữ liệu chuẩn quốc gia về ngưỡng nhiệt độ khởi phát ($T_{onset}$) và enthalpy tỏa nhiệt ($\Delta H$) cho hơn 50 hợp chất hữu cơ dễ phản ứng nhằm phục vụ công tác thẩm định an toàn nhà xưởng trước năm 2028 (do Cục Hóa chất phối hợp cùng các viện nghiên cứu chuyên ngành thực hiện).
  4. Tổ chức định kỳ 2 khóa đào tạo hàng năm về mô hình hóa động học nhiệt và đánh giá rủi ro nổ nhiệt cho trên 150 cán bộ quản lý an toàn và kỹ sư vận hành bồn phản ứng hóa chất (do các trường đại học chuyên ngành hóa học phối hợp với Bộ Công Thương triển khai).

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

Công trình nghiên cứu mang lại giá trị học thuật và ứng dụng sâu rộng cho 4 nhóm đối tượng chính:

  1. Kỹ sư công nghệ chế biến cao su và polyme: Tiếp thu phương pháp xác định chính xác năng lượng hoạt hóa và cơ chế lưu hóa 2 giai đoạn để tối ưu hóa thời gian lưu hóa đạt hiệu suất chuyển hóa trên 98%, giảm phế phẩm do hiện tượng quá nhiệt hay cháy xém vật liệu.
  2. Chuyên gia an toàn lao động và đánh giá rủi ro công nghiệp: Ứng dụng quy trình phân tích nhiệt DSC để sàng lọc nguy cơ nổ nhiệt cho các phản ứng tỏa nhiệt thể tích lớn trên 5.000 lít, từ đó thiết kế hệ thống giải nhiệt khẩn cấp đạt chuẩn an toàn quốc tế.
  3. Giảng viên và nhà nghiên cứu chuyên ngành Hóa lý - Hóa lý thuyết: Sử dụng khung so sánh thực nghiệm giữa 3 mô hình tự do (Friedman, OFW, ASTM E698) và mô hình cơ sở làm học liệu giảng dạy và tài liệu tham khảo cho các đề tài nghiên cứu động học chuyển pha.
  4. Học viên cao học và nghiên cứu sinh ngành Khoa học Vật liệu: Nắm vững kỹ thuật xử lý dữ liệu DSC phi đẳng nhiệt với độ tin cậy $R^2 > 0,998$, phục vụ trực tiếp cho việc thực hiện các luận văn liên quan đến phân hủy nhiệt, độ bền nhiệt và động học phản ứng pha rắn.

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

Kỹ thuật DSC dòng nhiệt có ưu điểm gì vượt trội so với phân tích nhiệt vi sai DTA truyền thống? Kỹ thuật DSC dòng nhiệt sử dụng cảm biến Alumel - Chromel tích hợp trong lò nhiệt chung, cho phép đo trực tiếp thông lượng nhiệt thực tế sinh ra hoặc mất đi theo thời gian. Nhờ cơ chế bù trừ và hiệu chuẩn phần cứng, DSC định lượng chính xác biến thiên enthalpy $\Delta H$ với độ sai số dưới 1,5%, vượt trội hoàn toàn so với DTA vốn chỉ ghi nhận sự chênh lệch nhiệt độ định tính.

Tại sao phương pháp động học phi đẳng nhiệt lại được ưu tiên hơn phương pháp đẳng nhiệt? Phương pháp phi đẳng nhiệt quét ở nhiều tốc độ gia nhiệt từ 5 K/phút đến 20 K/phút cho phép thu thập toàn bộ dữ liệu động học chỉ trong 3 giờ đến 4 giờ đo, tiết kiệm hơn 70% thời gian so với phương pháp đẳng nhiệt (thường mất trên 10 giờ). Hơn nữa, phương pháp này tránh được hiện tượng phản ứng bị mất mát tín hiệu trong giai đoạn gia nhiệt ban đầu.

Mô hình tự do Friedman và mô hình Ozawa - Flynn - Wall khác nhau như thế nào? Mô hình Friedman là phương pháp vi phân sử dụng logarit của tốc độ chuyển hóa vi phân tại từng mức độ chuyển hóa xác định mà không cần giả thiết hàm số, mang lại độ chính xác cao. Trong khi đó, mô hình OFW là phương pháp tích phân sử dụng công thức gần đúng tuyến tính của tốc độ quét theo nghịch đảo nhiệt độ $1/T$, dễ xử lý số liệu nhưng cần hiệu chỉnh theo ASTM E698 để giảm thiểu sai số từ 3% đến 5%.

Kết quả nghiên cứu trên cao su EPDM 3666 đóng góp gì cho an toàn sản xuất? Luận văn đã lượng hóa chính xác enthalpy tỏa nhiệt của hệ EPDM có lưu huỳnh là -42,8 J/g với năng lượng hoạt hóa 95,2 kJ/mol. Dữ liệu này giúp doanh nghiệp cài đặt chính xác ngưỡng nhiệt độ khống chế dưới 160 °C trong buồng cán ép, ngăn ngừa hoàn toàn nguy cơ nổ nhiệt cục bộ và tình trạng phân hủy chuỗi polyme do tích nhiệt.

Phần mềm NETZSCH Thermokinetics đóng vai trò gì trong việc dự đoán diễn biến phản ứng? Phần mềm giải hệ phương trình vi phân động học để mô phỏng chính xác sự biến đổi nồng độ chất phản ứng và sản phẩm theo bất kỳ chương trình nhiệt độ giả định nào. Điều này giúp doanh nghiệp tìm ra chế độ gia nhiệt tối ưu đạt độ chuyển hóa trên 98% mà không cần thực hiện hàng chục mẻ thử nghiệm thực tế tốn kém chi phí.

Kết luận

  • Xây dựng thành công quy trình toàn diện 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 quét nhiệt vi sai DSC dòng nhiệt kết hợp phần mềm mô phỏng chuyên dụng.
  • Lượng hóa chính xác các thông số nhiệt động của hệ lưu hóa cao su EPDM 3666 với biến thiên enthalpy dao động từ -15,4 J/g đến -42,8 J/g tùy thuộc vào hàm lượng lưu huỳnh.
  • Chứng minh phản ứng lưu hóa EPDM tuân theo cơ chế 2 giai đoạn tự xúc tác với hệ số tương quan $R^2 > 0,9985$, trong đó lưu huỳnh làm giảm năng lượng hoạt hóa từ khoảng 135 kJ/mol xuống còn xấp xỉ 98 kJ/mol.
  • Hoàn thiện mô hình tối ưu hóa biên dạng nhiệt độ gia công, giúp cắt giảm 18% đến 25% thời gian lưu hóa và loại bỏ triệt để nguy cơ quá nhiệt mất kiểm soát.
  • Cung cấp phương pháp luận khoa học và dữ liệu thực nghiệm tin cậy phục vụ công tác đánh giá an toàn nhiệt cho hơn 50 nhóm quy trình hóa chất công nghiệp tại Việt Nam.

Trong giai đoạn 2026 đến 2028, các nghiên cứu tiếp theo cần mở rộng áp dụng phương pháp này trên các hệ nhựa nhiệt rắn và vật liệu composite có enthalpy tỏa nhiệt cao trên -200 J/g. Các đơn vị sản xuất và cơ sở nghiên cứu hãy chủ động ứng dụng ngay kỹ thuật phân tích động học DSC vào dây chuyền chế biến nhằm tối đa hóa hiệu suất vận hành và đảm bảo an toàn tuyệt đối cho nhà xưởng.