Tổng quan nghiên cứu

Việt Nam thải ra hơn 70,7 triệu tấn rơm rạ khô mỗi năm sau các vụ thu hoạch lúa. Tuy nhiên, phần lớn lượng phụ phẩm này bị đốt trực tiếp trên đồng ruộng, phát thải ra môi trường từ 700 đến 4.100 mg khí methane ($CH_4$) và từ 19 đến 57 mg nitrous oxide ($N_2O$) trên mỗi kilôgam rơm rạ khô, đóng góp vào 10% tổng lượng phát thải khí nhà kính của ngành nông nghiệp. Thói quen xử lý này không chỉ gây ô nhiễm không khí nghiêm trọng do bụi mịn PM2.5 và PM10 mà còn làm thất thoát khoảng 45% lưu huỳnh, 40% nitơ và 33% kali trong đất. Trong bối cảnh thế giới thải ra hơn 300 triệu tấn rác thải nhựa mỗi năm với thời gian phân hủy kéo dài từ 100 đến 1.000 năm, việc chuyển đổi sinh khối rơm rạ thành vật liệu xanh có giá trị cao là một yêu cầu cấp bách.

Luận văn thạc sĩ chuyên ngành Kỹ thuật Hóa học mã số 8520301 của tác giả Nguyễn Hoài Nam, thực hiện tại Trường Đại học Bách khoa – ĐHQG TP.HCM dưới sự hướng dẫn của PGS. TS. Trần Tấn Việt (hoàn thành từ tháng 01/2024 đến tháng 05/2024 và bảo vệ vào tháng 07/2024), tập trung giải quyết bài toán này. Nghiên cứu đặt mục tiêu thiết kế và mô phỏng quy trình sản xuất tinh thể nanocellulose (Cellulose Nanocrystals - CNC) liên tục từ rơm rạ với lưu lượng nạp liệu 100 kg/h bằng phần mềm Aspen Plus. Kết quả đạt được hiệu suất thu hồi CNC đỉnh lên tới 90,32%, mở ra tiềm năng lớn để thay thế vật liệu polyme tổng hợp nhờ cơ tính vượt trội với độ bền kéo đạt 7.500–7.700 MPa và mô đun đàn hồi từ 130–250 GPa.

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 sự tích hợp giữa các mô hình động học hóa học và nhiệt động lực học quá trình:

  • Mô hình động học thủy phân axit giả đồng thể: Dựa trên nghiên cứu của Wang và các cộng sự, phản ứng chuyển hóa $\alpha$-cellulose thành tinh thể CNC và phản ứng phụ phân hủy glucose được mô tả qua hệ phương trình vi phân với năng lượng hoạt hóa $E_{a1} = 40,222$ kJ/mol (tạo CNC), $E_{a2} = 36,35$ kJ/mol (phân hủy CNC thành glucose) và $E_{a3} = 31,8$ kJ/mol (tạo glucose trực tiếp).
  • Mô hình động học tách Lignin và Silica: Áp dụng phương trình tốc độ Arrhenius theo Huang và các cộng sự để tính toán quá trình phân tách pha rắn trong môi trường kiềm nóng.
  • Gói nhiệt động NRTL (Non-Random Two-Liquid): Mô hình tính toán hệ số hoạt độ trong pha lỏng được lựa chọn để mô tả chính xác tương tác giữa các hợp chất phân cực mạnh ($H_2SO_4, H_2O, NaOH$) và các cấu tử phi phân cực hoặc chất rắn sinh khối ($Cellulose, Lignin, Hemicellulose, CNC$).

Cấu trúc lignocellulose trong rơm rạ chứa 20–40% cellulose, 15–24% hemicellulose và 10–18% lignin. Quá trình sản xuất CNC thực chất là loại bỏ các vùng vô định hình của sợi cellulose bằng axit mạnh để giữ lại các vùng tinh thể nano có đường kính từ 4–26 nm và chiều dài từ 40–260 nm.

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

  • Dữ liệu đầu vào và quy mô mô phỏng: Nguồn nguyên liệu rơm rạ đầu vào được chuẩn hóa ở công suất 100 kg/h với độ ẩm 5,58% và hàm lượng chất bốc đạt 65,8%. Dữ liệu thành phần sinh khối được phân tích chi tiết nhằm thiết lập cấu tử phi tiêu chuẩn trên Aspen Plus.
  • Công cụ và môi trường phân tích: Toàn bộ chuỗi 4 phân khu sản xuất được mô hình hóa trên Aspen Plus V. Lý do chọn phần mềm này là khả năng giải đồng thời các phương trình cân bằng vật chất và năng lượng phức tạp kết hợp kiểm soát điều kiện vận hành thiết bị phản ứng khuấy trộn liên tục (CSTR). Gói nhiệt động NRTL được lựa chọn thay cho UNIQUAC vì đáp ứng tốt hơn độ hội tụ của hệ điện ly và hỗn hợp đa cấu tử phức tạp.
  • Công cụ khảo sát độ nhạy (Sensity Analysis): Tiến hành quét các biến số vận hành gồm nồng độ $H_2SO_4$ (45 wt%, 60 wt%, 75 wt%), nhiệt độ phản ứng (30°C, 40°C, 50°C) và thời gian lưu phản ứng kéo dài từ 1 đến 5 giờ.
  • Thời gian và đối sánh: Nghiên cứu thực hiện từ tháng 01/2024 đến tháng 05/2024, sử dụng dữ liệu thực nghiệm đối chứng từ các công bố quốc tế uy tín nhằm hiệu chỉnh sai số mô hình dưới mức 5%.

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

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

Mô phỏng 4 phân khu công nghệ trên Aspen Plus đem lại các chỉ số kỹ thuật nổi bật:

  1. Hiệu quả xử lý sơ bộ và thủy phân kiềm: Tại cụm tiền xử lý, rơm rạ được rửa sạch và ngâm trong nước nóng 80°C với lưu lượng 200 kg/h trong 1 giờ để loại bỏ sáp và tạp chất vô cơ. Tiếp đó, quá trình xử lý bằng dung dịch NaOH 12 wt% ở 121°C trong 1 giờ đã loại bỏ hơn 85% lignin và phần lớn hemicellulose, chuyển hóa sinh khối thô thành dòng giàu cellulose.
  2. Hiệu suất tẩy trắng tạo $\alpha$-cellulose: Xử lý bằng $NaClO_2$ 5 wt% ở 75°C trong 90 phút với pH duy trì từ 3–5 bằng axit axetic 1 M đã loại bỏ triệt để lignin còn sót lại, tạo ra $\alpha$-cellulose có độ tinh khiết cao trước khi nạp vào cụm thủy phân axit.
  3. Ảnh hưởng của nồng độ axit đến hiệu suất CNC: Khi nồng độ axit $H_2SO_4$ tăng từ 45 wt% lên 60 wt% và 75 wt%, hiệu suất tạo CNC tăng trưởng rõ rệt. Tại nồng độ 75 wt%, quá trình thủy phân đạt đỉnh với hiệu suất thu hồi CNC đạt 90,32% ở 30°C sau 5 giờ phản ứng.
  4. Cân bằng vật chất tối ưu: Tỷ lệ rắn/lỏng (solid-to-liquor ratio) tối ưu được xác định ở mức 1:20 wt%, đảm bảo sự phân tán đồng đều của sợi cellulose trong pha lỏng và hạn chế tối đa sự nghẽn dòng trong thiết bị phản ứng CSTR.

Thảo luận kết quả

Dữ liệu mô phỏng từ Aspen Plus cho thấy nhiệt độ và nồng độ axit có tác động tương hỗ mật thiết đến độ tinh sạch và cấu trúc tinh thể của nanocellulose. Khi khảo sát trên đồ thị độ nhạy biến thiên theo thời gian từ 1 đến 5 giờ, ở nhiệt độ cao 50°C, tốc độ phản ứng ban đầu diễn ra rất nhanh nhưng phản ứng phụ phân hủy CNC thành đường đơn glucose ($k_2$) cũng tăng mạnh, làm giảm sản lượng CNC thu hồi sau 3 giờ. Ngược lại, ở nhiệt độ thấp 30°C kết hợp nồng độ axit 75 wt%, hằng số tốc độ hình thành CNC ($k_1$) chiếm ưu thế áp đảo, bảo toàn cấu trúc tinh thể hình que có kích thước nano đồng đều.

So sánh với các nghiên cứu trước đây như quy trình sản xuất CNC từ bột giấy Kraft của Gu và các cộng sự (hiệu suất khoảng 50%) hay quy trình thủy phân lá cọ dầu 1.000 kg/ngày của Qing và các cộng sự, quy trình mô phỏng rơm rạ trong luận văn này đạt hiệu suất vượt trội (90,32%). Các bảng đối chiếu dữ liệu giữa mô phỏng Aspen Plus và thực nghiệm của Thakur cùng Wang cho thấy độ lệch tuyệt đối chỉ dao động trong khoảng 2,1% đến 4,6%, chứng minh tính chính xác cao của mô hình tính toán nhiệt động NRTL.

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

Dựa trên các kết quả mô phỏng, tác giả đề xuất 4 nhóm giải pháp kỹ thuật cụ thể nhằm hiện thực hóa quy trình ở quy mô công nghiệp:

  1. Nâng quy mô thử nghiệm pilot lên 1.000 kg/h: Các kỹ sư thiết kế công nghệ và viện nghiên cứu cần triển khai hệ thống pilot từ quý 2/2025 đến quý 4/2025 với mục tiêu duy trì hiệu suất thu hồi CNC thực tế trên 88%, đồng thời kiểm soát chính xác thời gian lưu 5 giờ trong cụm phản ứng CSTR.
  2. Tích hợp hệ thống thu hồi và tái tuần hoàn $H_2SO_4$: Nhóm phát triển quy trình cần bổ sung tháp cô đặc axit chân không trong vòng 6 tháng tới nhằm tái sinh ít nhất 92% lượng axit sunfuric dư thừa từ công đoạn ly tâm tách rắn-lỏng (10.000 vòng/phút), giúp giảm đáng kể chi phí nguyên liệu hóa chất.
  3. Xây dựng cụm trung hòa tạo phụ phẩm thạch cao thương phẩm: Doanh nghiệp hợp tác sản xuất cần lắp đặt thiết bị phản ứng trung hòa bằng vôi tôi $Ca(OH)_2$ ở 25°C và 1 atm trong giai đoạn 2025–2026, đặt mục tiêu triệt tiêu 100% dòng thải axit ra môi trường và thu hồi thạch cao tinh khiết phục vụ ngành vật liệu xây dựng.
  4. Cài đặt hệ thống tự động hóa điều khiển thông số vận hành: Đội ngũ kỹ thuật nhà máy cần ứng dụng các cảm biến đo trực tuyến nồng độ và nhiệt độ trong vòng 12 tháng, khống chế biên độ sai số nhiệt độ dưới ±0,5°C và sai số nồng độ axit dưới 0,2% để ngăn chặn hiện tượng quá nhiệt gây cháy mẫu cellulose.

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

Luận văn là tài liệu tham khảo chuyên sâu và giá trị dành cho 4 nhóm đối tượng chính:

  1. Kỹ sư thiết kế quy trình hóa học và công nghệ vật liệu: Khai thác toàn bộ sơ đồ công nghệ (Process Flow Diagram), bảng cân bằng vật chất – năng lượng và thông số thiết lập mô hình NRTL trên Aspen Plus để ứng dụng thiết kế các nhà máy sản xuất polymer sinh học quy mô 100–1.000 kg/h.
  2. Học viên cao học và nhà nghiên cứu công nghệ nano sinh học: Sử dụng bộ dữ liệu động học chi tiết (năng lượng hoạt hóa $E_{a1} = 40,222$ kJ/mol, $E_{a2} = 36,35$ kJ/mol) và các phương trình phân tách lignin, silica để phát triển các đề tài mở rộng về biến tính bề mặt nanocellulose.
  3. Doanh nghiệp chế biến nông sản và sản xuất bao bì sinh học: Ứng dụng mô hình kinh tế tuần hoàn để chuyển đổi hàng triệu tấn rơm rạ phế thải thành màng nano composite tự phân hủy, phục vụ ngành bao bì thực phẩm và dược phẩm với độ bền kéo tương đương sợi carbon.
  4. Cơ quan quản lý môi trường và hoạch định chính sách kinh tế xanh: Tham khảo dữ liệu phát thải giảm thiểu (cắt giảm 700–4.100 mg $CH_4$/kg rơm rạ) để xây dựng chính sách khuyến khích tái chế sinh khối nông nghiệp bền vững tại các tỉnh đồng bằng.

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

Tại sao phần mềm Aspen Plus lại ưu tiên mô hình nhiệt động NRTL thay vì UNIQUAC cho mô phỏng này?
Mô hình NRTL (Non-Random Two-Liquid) tính toán chính xác hệ số hoạt độ trong pha lỏng có chứa hỗn hợp điện ly mạnh ($H_2SO_4, NaOH$) cùng các chất hữu cơ phân cực và không phân cực. Khi thử nghiệm mô hình UNIQUAC, phần mềm gặp lỗi phân kỳ do các cấu tử sinh khối không sẵn có thiếu hụt một số thông số tương tác năng lượng bề mặt, trong khi NRTL đảm bảo độ hội tụ tối ưu với sai số dưới 3%.

Nồng độ axit sunfuric và nhiệt độ nào đem lại hiệu suất tạo CNC cao nhất?
Kết quả khảo sát độ nhạy chỉ ra rằng điều kiện tối ưu đạt được ở nồng độ $H_2SO_4$ 75 wt%, nhiệt độ phản ứng 30°C và thời gian lưu 5 giờ. Ở chế độ này, hiệu suất chuyển hóa $\alpha$-cellulose thành tinh thể nanocellulose đạt đỉnh 90,32%, hạn chế tối đa phản ứng phụ phân hủy tinh thể thành glucose so với khi vận hành ở nhiệt độ 50°C.

Tỷ lệ rắn/lỏng (solid-to-liquor ratio) trong phản ứng thủy phân axit được thiết lập bao nhiêu?
Tỷ lệ rắn/lỏng được duy trì cố định ở mức 1:20 wt% (tương đương 1 g cellulose trong 20 ml dung dịch axit). Tỷ lệ này đảm bảo độ nhớt của bùn phản ứng trong thiết bị CSTR ở mức tối ưu, giúp axit khuếch tán đồng đều vào vùng vô định hình của sợi xenluloza mà không làm tiêu hao quá nhiều lượng hóa chất đầu vào.

Đặc tính cơ lý của nanocellulose từ rơm rạ có thể so sánh với những vật liệu nào?
Tinh thể nanocellulose thu được có độ bền kéo đạt 7.500–7.700 MPa và mô đun đàn hồi 130–250 GPa, tương đương với sợi carbon TORAYCA T1000G và vượt trội hơn thép không gỉ. Ngoài ra, vật liệu có độ bền nhiệt ổn định trong khoảng 120°C đến 220°C, lý tưởng làm pha gia cường cho màng composite kỹ thuật cao.

Mô hình mô phỏng xử lý giai đoạn tinh chế sản phẩm cuối cùng như thế nào?
Sau phản ứng thủy phân axit, hỗn hợp được rửa bằng nước cất qua hệ thống ly tâm tách pha lỏng ở 10.000 vòng/phút trong 15 phút để loại bỏ axit dư. Dòng huyền phù tiếp tục được thẩm tách đến pH trung tính, phân tán bằng sóng siêu âm và sấy đông khô (freeze-drying) nhằm thu hồi bột CNC với độ ẩm tiêu chuẩn dưới 2%.

Kết luận

  • Luận văn xây dựng thành công chuỗi mô phỏng liên tục 4 công đoạn: tiền xử lý, thủy phân kiềm (12 wt% NaOH, 121°C), tẩy trắng ($NaClO_2$ 5 wt%, 75°C) và thủy phân axit ($H_2SO_4$ 75 wt%, 30°C) với lưu lượng nạp 100 kg/h rơm rạ.
  • Xác lập điều kiện công nghệ tối ưu mang lại hiệu suất thu hồi tinh thể nanocellulose (CNC) đạt 90,32%, ăn khớp chặt chẽ với số liệu thực nghiệm (sai số dưới 5%).
  • Ứng dụng thành công gói nhiệt động NRTL và hệ phương trình động học Arrhenius với năng lượng hoạt hóa $E_{a1} = 40,222$ kJ/mol để dự báo chính xác diễn biến phản ứng trong thiết bị CSTR.
  • Cung cấp cơ sở khoa học và giải pháp công nghệ vững chắc để giải quyết triệt để bài toán ô nhiễm từ 70,7 triệu tấn rơm rạ phát thải hàng năm tại Việt Nam.
  • Đóng góp mô hình toán học hoàn chỉnh phục vụ nghiên cứu phát triển màng bao bì sinh học tự phân hủy và composite nano y sinh trong tương lai.

Quá trình chuyển đổi từ mô phỏng sang sản xuất thực tế cần được tiếp tục thử nghiệm trên các hệ thống pilot lớn hơn trong giai đoạn 2025–2026. Các kỹ sư và nhà nghiên cứu quan tâm có thể khai thác trực tiếp bộ thông số vận hành và tệp mô phỏng Aspen Plus trong luận văn để thúc đẩy các dự án công nghệ sinh học xanh phục vụ phát triển bền vững.