Nghiên cứu phương pháp xác định các axit béo trong một số loại dầu mỡ động, thực vật Việt Nam bằng kỹ thuật sắc ký Võ Thị Việt Dung Trường Đại học Khoa học Tự nhiên Luận án TS. ngành: Hóa phân tích; Mã số: 62 44 29 01 Người hướng dẫn: PGS.TS. Nguyễn Xuân Trung, PGS. TS. Lưu Văn Bôi Năm bảo vệ: 2012 Abstract. Sử dụng phương pháp quy hoạch hóa thực nghiệm - xây dựng mô hình hồi quy bậc hai tâm trực giao để tìm điều kiện tối ưu cho quá trình xác định các axit béo trong dầu mỡ động thực vật. Kết quả nghiên cứu đã thu được tối ưu hơn nhiều so với các nghiên cứu phân tích axit béo trước đây. Đã sử dụng các kỹ thuật sắc ký hiện đại: sắc ký khí ion hóa ngọn lửa (GC/FID) và sắc ký lỏng siêu nhanh (UFLC) để phân tích đồng thời 37 metyl este axit béo trong hỗn hợp. Các nghiên cứu được tiến hành một cách có hệ thống từ khảo sát đơn lẻ đến mô hình hóa các yếu tố có ảnh hưởng tương hỗ. Quá trình xử lý thống kê và xác nhận giá trị sử dụng của phương pháp nghiên cứu được tiến hành đầy đủ, đáp ứng yêu cầu khi xây dựng phương pháp phân tích. Ứng dụng phương pháp nghiên cứu xác định các axit béo trong 8 mẫu dầu mỡ động thực vật thông dụng của Việt Nam gồm: mỡ cá basa, mỡ lợn, dầu mè, dầu đậu nành, dầu lạc, dầu dừa, dầu hạt cao su và dầu hạt jatropha. Phân tích đề xuất có khả năng xây dụng thành phương pháp tiêu chuẩn để xác định thành phần các axit béo trong dầu mỡ động thực vật, phục vụ nghiên cứu về dinh dưỡng, kiểm nghiệm thực phẩm, chế tạo nhiên liệu sinh học và các nghiên cứu khác. Keywords. Hóa phân tích; Axit béo; Kỹ thuật sắc ký; Dầu động vật; Dầu thực vật Content MỞ ĐẦU 1. Tính cấp thiết của luận án Dầu mỡ động, thực vật (DMĐTV) đóng vai trò to lớn trong việc cung cấp dinh dưỡng thiết yếu cho cơ thể sống dưới dạng thức ăn, thực phẩm chức năng, làm nguyên liệu chế biến thức ăn gia súc, nguyên liệu sản xuất nhiên liệu sinh học (biodiesel), … Dưới dạng thực phẩm thiết yếu, nếu sử dụng và tiêu thụ DMĐTV đúng cách, đúng hàm lượng chất béo sẽ cân bằng được dinh dưỡng và kiểm soát bệnh tật. Với các loại DMĐTV có hàm lượng dinh dưỡng thấp, dầu mỡ thải hoặc đã qua sử dụng, DMĐTV không ăn được như dầu hạt jatropha, hạt cao su, thầu dầu, tảo,… có thể sử dụng để làm nguyên liệu sản xuất biodiesel. Đây được xem là nguồn năng lượng thay thế đầy tiềm năng. Trên cả hai phương diện là thực phẩm dinh dưỡng và nhiên liệu sinh học thì tại Việt Nam hiện chưa có cơ sở dữ liệu đầy đủ về hàm lượng các axit béo trong dầu mỡ động, thực vật. Các quy trình phân tích đã công bố và áp dụng chủ yếu từ tài liệu tham khảo của nước ngoài. Các số liệu phân tích hàm lượng axit béo không đầy đủ, nên khó so sánh đánh giá toàn diện nguồn cung cấp axit béo cho nhu cầu làm thực phẩm hoặc cho nguyên liệu sản xuất nhiên liệu sinh học. Vì vậy đề tài: “Nghiên cứu phương pháp xác định các axit béo trong một số loại dầu mỡ động, thực vật Việt Nam bằng kỹ thuật sắc ký” sẽ giúp cho việc hoàn thiện quy trình phân tích axit béo, đáp ứng được yêu cầu phân tích cho các phòng thí nghiệm đầu ngành tại Việt Nam. 2. Mục tiêu và nội dung nghiên cứu Mục tiêu của đề tài là nghiên cứu tương đối toàn diện các vấn đề liên quan đến quá trình tách và xác định đồng thời các axit béo trong hỗn hợp, quá trình chuyển đổi DMĐTV thành metyl este axit béo (FAME) làm cơ sở khoa học để tách và xác định các axit béo trong 8 mẫu DMĐTV Việt Nam. Từ đó xây dựng phương pháp tiêu chuẩn xác định hàm lượng các axit béo trong DMĐTV ở nước ta. Do đó nội dung của luận án gồm: - Xây dựng mô hình thí nghiệm bậc hai tâm trực giao xét ảnh hưởng đồng thời các yếu tố đến khả năng tách và xác định các FAME bằng phương pháp sắc ký khí. - Dựa trên mô hình thí nghiệm ở trên tối ưu hóa và đánh giá ảnh hưởng các nhân tố khảo sát đến việc tách và xác định các FAME bằng phương pháp sắc ký khí detector ion hóa ngọn lửa (GC/FID) - Sử dụng hệ thống sắc ký lỏng siêu nhanh (UFLC) nghiên cứu tìm điều kiện tối ưu và đánh giá ảnh hưởng đồng thời các nhân tố khảo sát đến hiệu suất quá trình chuyển đổi este dầu mỡ. - Xác nhận giá trị sử dụng của phương pháp nghiên cứu. - Ứng dụng phương pháp nghiên cứu xác định các axit béo trong 8 mẫu DMĐTV Việt Nam. 3. Những điểm mới về khoa học của luận án 3.1. Lần đầu tiên ở Việt Nam sử dụng phương pháp quy hoạch hóa thực nghiệm theo mô hình thí nghiệm bậc hai tâm trực giao nghiên cứu ảnh hưởng đồng thời các yếu tố như tốc độ pha động, nhiệt độ ban đầu của cột tách, tốc độ tăng nhiệt, tỉ lệ bơm chia dòng đến khả năng tách và xác định các FAME bằng phương pháp sắc ký khí với detector ion hóa ngọn lửa. - Lần đầu tiên xây dựng mô hình thí nghiệm bậc hai tâm trực giao xét ảnh hưởng đồng thời các yếu tố đến hiệu suất quá trình chuyển đổi este DMĐTV như tỉ lệ mol metanol/dầu mỡ, hàm lượng chất xúc tác (axit hoặc kiềm), nhiệt độ, thời gian và môi trường phản ứng nhằm đạt hiêu suất chuyển đổi tối ưu ( 100%). 3.2. Dựa trên các điều kiện tối ưu hóa ở trên lần đầu tiên đã tách và xác định đồng thời 37 FAME trong hỗn hợp với độ lặp lại tốt, làm cơ sở để tiến hành xác định các axit béo trong DMĐTV. 3.3. Đã tiến hành nghiên cứu một cách có hệ thống các yếu tố ảnh hưởng đến việc tách, xác định các axit béo trong DMĐTV bằng phương pháp sắc ký khí. Kết quả phân tích các axit béo đã xử lý thống kê, xác nhận giá trị sử dụng của phương pháp, đáp ứng yêu cầu khi xây dựng phương pháp phân tích. 3.4. Lần đầu tiên đã tiến hành phân tích đồng thời các axit béo trong 8 mẫu DMĐTV Việt Nam: mỡ cá basa, mỡ lợn, dầu mè, dầu đậu nành, dầu lạc, dầu dừa, dầu hạt cao su và dầu hạt jatropha. 3.5. Phương pháp phân tích đề xuất có khả năng xây dựng thành phương pháp tiêu chuẩn xác định thành phần axit béo trong DMĐTV Việt Nam. 4. Bố cục của luận án Luận án gồm 151 trang với 42 bảng số liệu, 37 hình vẽ, 178 tài liệu tham khảo và 11 phần phụ lục được chia thành các phần như sau: mở đầu: 2 trang, chương 1. tổng quan tài liệu: 31 trang, chương 2. thực nghiệm và phương pháp nghiên cứu: 23 trang, chương 3. kết quả và thảo luận: 74 trang, kết luận và kiến nghị: 2 trang, tài liệu tham khảo: 18 trang. Chƣơng 1. TỔNG QUAN Chương này tổng hợp lý thuyết về chất béo và axit béo, các công trình nghiên cứu ở Việt Nam và nước ngoài về các phương pháp định lượng chất béo và axit béo, các phương pháp điều chế và các yếu tố ảnh hưởng đến quá trình điều chế FAME, phương pháp quy hoạch hóa thực nghiệm tìm điều kiện tối ưu tách và xác định các axit béo. Chƣơng 2. THỰC NGHIỆM VÀ PHƢƠNG PHÁP NGHIÊN CỨU 2.1. Nguyên liệu, hóa chất và thiết bị thí nghiệm Nguyên liệu DMĐTV được lựa chọn để nghiên cứu gồm: Mỡ cá basa An Giang do Công ty cổ phần xuất nhập khẩu thủy sản An Giang (Agifish) cung cấp, mỡ lợn từ các chợ Hà Nội, được nấu chảy thành dạng dầu, dầu đậu nành, dầu mè, dầu lạc, dầu dừa từ Hợp tác xã Nông nghiệp Quảng Ngãi, dầu được ép từ hạt đậu nành, hạt mè, hạt lạc và cơm dừa Quảng Ngãi, dầu hạt cao su và dầu hạt jatropha từ công ty Sản xuất – Xuất nhập khẩu Bình Dương, dầu được ép từ hạt cao su và hạt jatropha Đăk Nông. Chất chuẩn và nội chuẩn của Sigma (Mỹ), dung môi loại phân tích của Merck (Đức), các hóa chất khác: CH 3 OH công nghiệp 90% (Malaysia), KOH 82%, NaOH 96%, H 2 SO 4 96% (Trung Quốc),… Thiết bị sắc ky ́ lỏng siêu nhanh UFLC của Shimadzu (Nhật) với detector chỉ số khúc xạ RID–10A, máy sắc ký lỏng: LC–20AD, degasser: DGU–20A 3 , lò cột: CTO–20A, cột sắc ký: Cadenzal CD–C18 (250 mm × 4,6 mm ID × 3 μm) của Imtakt (Mỹ); thiết bị sắc ký khí Shimadzu GC–2010 (Kyoto, Nhật) với detector FID, cột mao quản phân cực mạnh SP TM - 2560 (100 m × 0,25 mm × 0,2 m) và cột mao quản ít phân cực Equity-5 (30 m × 0,25 mm × 0,25 m) của Supelco (Mỹ), máy khuấy từ gia nhiệt Velp 1200 rmp (Italia), lò vi sóng Sanyo 1200W (Nhật), bơm chân không IKA (Đức), máy rung siêu âm Telsonic 750W (Thụy Sỹ), 2.2. Phƣơng pháp nghiên cứu 2.2.1. Tách và xác định các FAME Sử dụng sắc ký khí với detector FID, cột mao quản phân cực mạnh SP TM -2560 (100 m × 0,25 mm × 0,25 m) có thành phần 100% poly(xyanoisopropyl siloxan), khí mang Heli, dung môi diclometan, nội chuẩn n-hexadecan. 2.2.2. Xác định hàm lƣợng axit béo tự do trong dầu mỡ Sử dụng phương pháp chuẩn độ với dung dịch NaOH 0,01M, chỉ thị phenolphatalein. 2.2.3. Chuyển đổi este dầu mỡ động thực vật Sử dụng phản ứng transeste hóa DMĐTV với metanol xúc tác kiềm và phản ứng este hóa DMĐTV với metanol xúc tác axit. 2.2.4. Xác định mức độ chuyển đổi este Sử dụng sắc ký lỏng siêu nhanh với cột sắc ký Cadenzal CD–C18 (250 mm × 4,6 mm × 3 μm), detector: RID–10A, nhiệt độ lò cột: 35 0 C, pha đô ̣ ng: Ace/ACN (7/3; v/v), tốc độ: 0,5 ml/phút, thể tích tiêm mẫu: 25µl. 2.2.5. Thu thập mẫu phân tích Sử dụng phương pháp lấy mẫu DMĐTV theo TCVN 2625:2005 (hay ISO 05555:2001). 2.2.6. Nghiên cứu tìm điều kiện tối ƣu cho quá trình tách và xác định các axit béo Sử dụng phương pháp mặt mục tiêu theo mô hình hồi quy bậc hai tâm trực giao với sự hỗ trợ của phần mềm Modde 5.0. 2.2.7. Đánh giá thống kê phƣơng pháp phân tích Thông qua việc xây dựng đường chuẩn, tính sai số, độ nhạy, độ lặp lại, độ đúng, hiệu suất thu hồi và so sánh kết quả phân tích với kết quả kiểm nghiệm Vilas (ISO 17025:2005). Chƣơng 3. KẾT QUẢ VÀ THẢO LUẬN 3.1. Nghiên cứu tối ƣu hóa điều kiện GC/FID xác định các FAME 3.1.1. Tối ƣu hóa điều kiện GC/FID theo phƣơng pháp đơn biến Trước hết chúng tôi khảo sát, so sánh khả năng tách các FAME (độ phân giải R) và độ lặp lại (độ lệch chuẩn tương đối RSD%) về tỉ lệ diện tích pic các FAME/IS để lựa chọn loại cột tách (cột mao quản phân cực mạnh SP TM -2560) và loại khí mang (khí He) cố định dùng cho GC/FID do trong quá trình phân tích việc thay đổi hai nhân tố này gặp nhiều khó khăn. Độ nhạy, độ ổn định (độ lặp lại) của hệ thống phân tích thay đổi đáng kể khi thay đổi hai yếu tố này. Sau đó chúng tôi tiếp tục khảo sát xác định khoảng biến thiên các điều kiện phân tích của GC/FID để tối ưu hóa các điều kiện này nhằm tách tốt các FAME. Kết quả khảo sát đơn biến tìm khoảng biến thiên các điều kiện phân tích GC/FID thu được như sau: - Tốc độ pha động (tốc độ khí mang) He: 18 – 25 cm/s - Tốc độ tăng nhiệt của cột tách: 2 – 4 0 C/phút - Nhiệt độ ban đầu của cột tách: 60 – 140 0 C - Tỉ lệ bơm chia dòng: 10:1 – 200:1 3.1.2. Tối ƣu hóa điều kiện GC/FID theo phƣơng pháp đa biến 3.1.2.1. Xây dựng mô hình hồi quy bậc hai tâm trực giao Sử dụng mô hình hồi quy bậc hai tâm trực giao để xây dựng phương trình hồi quy (PTHQ) và đánh giá ảnh hưởng đồng thời của 4 nhân tố khảo sát: tốc độ pha động, X 1 ; nhiệt độ ban đầu của cột tách, X 2 ; tốc độ tăng nhiệt của cột tách, X 3 ; tỉ lệ bơm chia dòng, X 4 đến khả năng tách (độ phân giải R) các FAME, y 1 và độ lệch chuẩn tương đối (RSD%) của tỉ lệ diện tích pic FAME/IS, y 2 . Kết quả thực nghiệm được trình bày trên bảng 3.2. Quan hệ giữa hàm mục tiêu và các nhân tố khảo sát được mô tả theo phương trình bậc hai: y = b 0 x 0 + b 1 x 1 + + b 12 x 1 x 2 + b 13 x 1 x 3 + + b 11 (x 1 2 – ) + b 22 (x 2 2 – ) +… Với tham số hồi quy  và cánh tay đòn sao d được tính sẵn cho mô hình hồi quy bậc hai tâm trực giao khi n = 4 là:  = 0,8000, d = 1,4142. Bảng 3.2. Kết quả thực nghiệm tìm điều kiện GC/FID tách và xác định các FAME theo mô hình hồi quy bậc hai tâm trực giao TT Biến mã hóa Biến số thực Kết quả x 1 x 2 x 3 x 4 X 1 X 2 X 3 X 4 y 1 y 2 1 -1 -1 -1 -1 18 60 2 10 0,5 9 1,76 2 1 -1 -1 -1 25 60 2 10 0,6 1,88 3 -1 1 -1 -1 18 140 2 10 1,0 2 1,75 4 1 1 -1 -1 25 140 2 10 0,2 1 1,79 5 -1 -1 1 -1 18 60 4 10 0 1,92 6 1 -1 1 -1 25 60 4 10 0 1,79 7 -1 1 1 -1 18 140 4 10 0,4 3 1,69 8 1 1 1 -1 25 140 4 10 0 1,75 9 -1 -1 -1 1 18 60 2 200 0,6 3,85 2 10 1 -1 -1 1 25 60 2 200 0,1 7 3,92 11 -1 1 -1 1 18 140 2 200 1,0 9 3,82 12 1 1 -1 1 25 140 2 200 0,2 4 3,88 13 -1 -1 1 1 18 60 4 200 0 3,84 14 1 -1 1 1 25 60 4 200 0 3,87 15 -1 1 1 1 18 140 4 200 0,4 7 3,84 16 1 1 1 1 25 140 4 200 0 3,92 17 - 1,4142 0 0 0 16,6 100 3 105 0,7 8 3,23 18 1,4142 0 0 0 26,4 100 3 105 0,1 5 3,18 19 0 - 1,4142 0 0 21,5 43,4 3 105 0,6 8 3,2 20 0 1,4142 0 0 21,5 156,6 3 105 0,9 9 3,23 21 0 0 - 1,4142 0 21,5 100 1,6 105 0,6 6 3,16 22 0 0 1,4142 0 21,5 100 4,4 105 0 3,19 23 0 0 0 - 1,4142 21,5 100 3 -29,3 - - 24 0 0 0 1,4142 21,5 100 3 239, 3 0,8 6 4,07 25 1 0 0 0 0 21,5 100 3 105 0,9 6 3,24 25 2 0 0 0 0 21,5 100 3 105 0,9 6 3,19 25 3 0 0 0 0 21,5 100 3 105 0,9 7 3,23 Ở đây, cần tiến hành N = 27 thí nghiệm để xác định các hệ số trong PTHQ. Thí nghiệm thứ 25 là thí nghiệm ở tâm được tiến hành 3 lần. Phương sai của các thí nghiệm được xác định theo thí nghiệm bổ sung ở tâm là 2 0(R) S = 7,43.10 -5 , 2 0(RSD) S = 7.10 -4 . Từ kết quả thí nghiệm, sử dụng phần mềm Modde 5.0 loại bỏ các hệ số hồi quy không phù hợp, tức là các hệ số có TN i t = i ib b /S < t p (f) hay TN i t < t (0,05) (2) = 4,3. Kết quả được chỉ ra trên bảng 3.3. Bảng 3.3. Kết quả phân tích ảnh hưởng của các nhân tố khảo sát đến khả năng tách các FAME (I) và độ lệch chuẩn tương đối của tỉ lệ diện tich pic FAME/IS (II) của GC/FID theo mô hình hồi quy bậc hai tâm trực giao y Hệ số hồi quy Phƣơng sai mẫu Chuẩn t Độ tin cậy (±) I II I II I II I II b 0 0,95 3,24 1,42.10 -16 8,01.10 - 19 8,0.10 7 3,6.10 9 0,03 0,06 x 1 -0,22 0,01 3,89.10 -13 0,30 3,7.10 5 0,018 0,01 0,03 x 2 0,12 -0,02 2,84.10 -10 0,17 7,0.10 3 0,049 0,01 0,03 x 3 -0,21 6,00.10 -4 6,10.10 -13 0,96 2,7.10 5 6,1.10 -4 0,01 0,03 x 4 0,01 1,04 0,12 1,45.10 - 16 0,024 8,6.10 7 0,01 0,03 x 1 .x 1 -0,24 -0,02 2,36.10 -11 0,29 4,9.10 4 0,037 0,02 0,04 x 2 .x 2 -0,05 -0,02 8,71.10 -5 0,41 5,7 0,031 0,02 0,04 x 3 .x 3 -0,31 -0,04 1,50.10 -12 0,08 2,5.10 5 0,14 0,02 0,04 x 4 .x 4 -0,05 -0,33 1,19.10 -5 1,15.10 -8 13 3,1.10 -3 0,02 0,05 x 1 .x 2 -0,11 0,01 3,17.10 -9 0,49 1,9.10 3 0,014 0,01 0,03 x 1 .x 3 0,10 -0,02 3,99.10 -9 0,26 1,6.10 3 0,039 0,01 0,03 x 1 .x 4 -0,01 0,01 0,30 0,49 0,013 0,014 0,01 0,03 x 2 .x 3 -0,01 3,10.10 -3 0,20 0,82 0,02 3,4.10 -3 0,01 0,03 x 2 .x 4 0,01 0,02 0,38 0,13 9,7.10 -3 0,056 0,01 0,03 x 3 .x 4 -0,01 1,90.10 -3 0,31 0,89 0,013 2,0.10 -3 0,01 0,03 N = 26 2 I R  0,9983 RSD I =0,025 P =0,95 2 II R  0,9984 RSD II =0,053 *Ghi chú: Những hệ số in nghiêng, đậm là hệ số không có ý nghĩa trong PTHQ. Những kết quả thu được ở bảng 3.3 cho thấy có ảnh hưởng bậc nhất tới khả năng tách các FAME y 1 là tốc độ pha động x 1 (tiêu cực), nhiệt độ ban đầu x 2 (tích cực) và tốc độ tăng nhiệt của cột tách x 3 (tiêu cực), còn tỉ lệ bơm chia dòng x 4 thì ảnh hưởng không đáng kể. Ngược lại tỉ lệ bơm chia dòng lại là nhân tố bậc nhất ảnh hưởng chủ yếu (tích cực) đến độ lệch chuẩn tương đối của tỉ lệ diện tích pic FAME/IS y 2 , ba nhân tố còn lại thì ít ảnh hưởng. Khả năng tách các FAME chịu ảnh hưởng tiêu cực ở bậc hai của cả bốn nhân tố khảo sát và ảnh hưởng đồng thời của tốc độ pha động và nhiệt độ ban đầu của cột tách, x 1 x 2 (tiêu cực), tốc độ pha động và tốc độ tăng nhiệt của cột tách, x 1 x 3 (tích cực). Trong khi đó độ lệch chuẩn tương đối của tỉ lệ diện tích pic FAME/IS chỉ chịu ảnh hưởng bậc hai của tỉ lệ bơm chia dòng, x 4 2 (tiêu cực), ảnh hưởng bậc hai và ảnh hưởng đồng thời của các nhân tố khác là không đáng kể. Từ kết quả thu được trên bảng 3.3, sau khi loại những nhân tố không phù hợp, có thể viết được PTHQ dạng mã hóa mô tả sự phụ thuộc của y 1 , y 2 vào các nhân tố x i như sau: y 1 = 0,95 – 0,22x 1 + 0,12x 2 – 0,21x 3 – 0,11x 1 x 2 + 0,10x 1 x 3 – 0,24(x 1 2 – 0,800) – 0,05(x 2 2 – 0,800) – 0,31(x 3 2 – 0,800) – 0,05(x 4 2 – 0,800) (3.1) y 2 = 3,24 + 1,04x 4 – 0,33(x 4 2 – 0,800) (3.2) Sử dụng chuẩn Fisher để kiểm tra sự phù hợp của mô hình: F mô hình I = I 2 phï hîp I 2 0 S S = 4 5 4 ,53.1 0 6 ,09 7,43.10    < F bảng = F 0,95 (15, 2) = 19,43. F mô hình II = II 2 phï hîp II 2 0 S S = 3 4 1,33 . 10 1,91 7.1 0    < F bảng = F 0,95 (23, 2) = 19,45. F mô hình < F bảng : mô hình thống kê mô tả đúng thực nghiệm. Do đó mô hình toán học mô tả quan hệ giữa khả năng tách các FAME và độ lệch chuẩn tương đối của tỉ lệ diện tích pic FAME/IS vào các biến thực như sau: Y 1 = -8,67 + 0,76X 1 + 0,03X 2 + 0,99X 3 + 0,01.10 -1 X 4 – 0,75.10 -3 X 1 X 2 + 0,03X 1 X 3 – 0,02X 1 2 – 0,33.10 -4 X 2 2 – 0,31X 3 2 – 0,05.10 -4 X 4 2 (3.3) Y 2 = 1,95 + 0,02X 4 – 3,65.10 -5 X 4 2 (3.4) Với o ii i i XX x    . Trong đó: X i : các giá trị dạng thực của nhân tố khảo sát, o i X : mức gốc (0) của X i , x i : các giá trị mã hóa của X i ,  i : khoảng biến thiên X i . Hai PTHQ ở trên phản ánh khá chính xác mô hình thực nghiệm, điều này được khẳng định qua các giá trị độ lệch chuẩn R 2 của mô hình (bảng 3.3). Mô hình có thể dự đoán khả năng tách các FAME và độ lệch chuẩn tương đối của tỉ lệ diện tích pic FAME/IS với độ chính xác cao, tương ứng là 99,83% và 99,84%, giữa giá trị thu được từ mô hình dự đoán và giá trị thực nghiệm chỉ có một sự sai khác nhỏ. 3.1.2.2. Xác định điều kiện tối ƣu của GC/FID Dựa trên mô hình hồi quy bậc hai tâm trực giao xây dựng được khi thực nghiệm khảo sát ảnh hưởng đồng thời của các nhân tố của GC/FID đến khả năng tách các FAME và độ lệch chuẩn tương đối của tỉ lệ diện tích pic FAME/IS để xác định điều kiện tối ưu. Bằng cách lấy đạo hàm của hàm hồi quy theo từng biến và cho bằng 0 để tìm ra cực trị của hàm (thông qua phần mềm Modde 5.0) tìm được điều kiện tối ưu để tách và xác định các FAME như sau: tốc độ khí mang He: 20,5 cm/s; tốc độ tăng nhiệt của cột tách: 2,4 0 C/phút; nhiệt độ ban đầu của cột tách: 140 0 C; tỉ lệ bơm chia dòng: 30:1. Ở điều kiện tối ưu trên, khả năng tách các FAME và độ lệch chuẩn tương đối của tỉ lệ diện tích pic FAME/IS có giá trị dự đoán theo mô hình tối ưu hóa lần lượt là 1,08 và 2,11%. Các giá trị trung bình tương ứng thu được khi tiến hành thực nghiệm 3 lần với chất chuẩn FAME 20–40–60 ppm là 1,08 và 2,12%. Giá trị các phương sai thực nghiệm theo điều kiện tối ưu thu được: I 2 tu S = 6,43.10 -5 , II 2 tu S = 5,33.10 -4 và của mô hình tương ứng là: I 2 0 S = 7,43.10 -5 , II 2 0 S = 7,00.10 -4 . Có F tính I = 7,43.10 - 5 /6,43.10 -5 = 1,16; F tính II = 7,00.10 -4 /5,33.10 -4 = 1,31 < F bảng (0,95; 2; 2) = 19. Như vậy độ chính xác của các số liệu theo mô hình và theo thực nghiệm ở điều kiện tối ưu là giống nhau. Kết quả này cho thấy mô hình dự đoán có độ chính xác cao, đảm bảo độ tin cậy. Đây là một công cụ hữu ích trong việc tiên đoán kết quả cũng như tối ưu hóa các điều kiện thực nghiệm. 3.2. Nghiên cứu chuyển đổi DMĐTV thành FAME 3.2.1. Kiểm tra hàm lƣợng axit béo tự do, lựa chọn quy trình chuyển đổi Hàm lượng axit béo tự do (%FFA) có ảnh hưởng đến quy trình chuyển đổi DMĐTV thành FAME. Các phản ứng chuyển đổi este dùng xúc tác kiềm cần có nguyên liệu đầu với hàm lượng FFA < 2,5%. Nếu hàm lượng FFA > 2,5% thì cần sử dụng nhiều kiềm hơn để trung hòa FFA, đồng thời sản phẩm tạo ra xà phòng và nước nhiều hơn. Xà phòng sinh ra làm tăng độ nhớt của hỗn hợp, tạo thành nhũ và làm cho việc tách glyxerin trở nên khó khăn. Phản ứng sinh ra nhiều nước có thể làm cho metyl este bị thủy phân tạo ra axit béo và metanol, do vậy tiêu tốn nhiều xúc tác hơn và hiệu suất phản ứng bị giảm đi đáng kể. Do đó với dầu mỡ có hàm lượng FFA > 2,5% cần chuyển đổi qua hai giai đoạn, giai đoạn đầu thực hiện phản ứng este hóa xúc tác axit để chuyển axit béo tự do thành FAME sau đó tiếp tục giai đoạn hai với phản ứng transeste hóa dầu mỡ sử dụng xúc tác kiềm. Vì vậy trước hết chúng tôi chuẩn độ kiểm tra hàm lượng axit béo tự do của nguyên liệu DMĐTV để lựa chọn quy trình chuyển đổi DMĐTV thành FAME thích hợp. Kết quả phân tích cho thấy hàm lượng FFA trong mỡ cá basa: 1,26%; mỡ lợn: 0,45%; dầu đậu nành: 0,21%; dầu mè: 1,68%; dầu lạc: 0,33% và dầu dừa: 0,21%. Các mẫu dầu mỡ này đều có hàm lượng FFA < 2,5%, hoàn toàn thích hợp tiến hành chuyển đổi thành FAME bằng phương pháp metyl este hóa xúc tác kiềm một giai đoạn. Riêng dầu hạt cao su và dầu hạt jatropha có hàm lượng FFA rất lớn, tương ứng 46,19% và 21,01%, phải tiến hành chuyển đổi qua hai giai đoạn, giai đoạn đầu este hóa với xúc tác axit, giai đoạn hai chuyển vị este với xúc tác kiềm. 3.2.2. Điều kiện chuyển đổi este theo phƣơng pháp đơn biến Qua khảo sát đơn biến các nhân tố ảnh hưởng đến hiệu suất phản ứng chuyển đổi este giữa mỡ cá basa và metanol sử dụng xúc tác kiềm và có mặt của dung môi chúng tôi thu được khoảng biến thiên các điều kiện phản ứng như sau: + Tỉ lệ mol metanol/mỡ cá (X 1 , mol/mol)): 4/1 – 7/1; + Nồng độ xúc tác KOH (X 2 , % m/m): 0,75 – 1,25%; + Nhiệt độ phản ứng (X 3 , 0 C): 30 – 60 0 C; + Thời gian phản ứng (X 4, phút): 10 – 40 phút.; + Hàm lượng dung môi axeton (X 5 , % m/m): 10 – 30%. Kết quả khảo sát đơn biến các nhân tố ảnh hưởng đến hiệu suất phản ứng este hóa xúc tác axit (giai đoạn 1) của quá trình chuyển đổi 2 giai đoạn điều chế FAME từ dầu hạt cao su thu được khoảng biến thiên các nhân tố như sau (tất cả các phản ứng đều thực hiện ở nhiệt độ sôi của hỗn hợp,  60 0 C): + Tỉ lệ mol metanol/dầu (X 1 , mol/mol): 3/1 – 9/1 + Hàm lượng xúc tác H 2 SO 4 (X 2 , % m/m): 0,5 – 2% + Thời gian phản ứng (X 3, giờ): 1 – 6 giờ + Hàm lượng dung môi axeton (X 4 , % m/m): 10 – 50% 3.2.3. Điều kiện chuyển đổi este theo phƣơng pháp đa biến * Xây dựng mô hình hồi quy bậc hai tâm trực giao Tương tự như quá trình thực nghiệm tối ưu hóa các điều kiện GC/FID để tách và xác định các FAME, với quá trình chuyển đổi DMĐTV, qua mô hình hồi quy bậc hai tâm trực giao, chúng tôi thu được các PTHQ dạng mã hóa sau khi loại bỏ các nhân tố ảnh hưởng không đáng kể như sau: - Hiệu suất phản ứng este hóa xúc tác kiềm điều chế FAME từ mỡ cá basa ( = 0,7698): y 1 = 94,06 + 5,88x 1 + 1,53x 2 + 2,18x 3 + 4,46x 4 + 0,99x 5 – 0,80x 1 x 2 – 0,77x 1 x 3 – 0,71x 1 x 4 – 4,76x 4 2 + 0,96x 5 2 (3.5) - Hiệu suất phản ứng este hóa xúc tác axit (giai đoạn 1) của quá trình chuyển đổi hai giai đoạn điều chế FAME từ dầu hạt cao su ( = 0,8000): y 2 = 105,36 + 9,24x 1 + 3,61x 2 + 3x 3 + 2,08x 4 – 1,18x 1 x 2 – 1,04x 1 x 4 – 1,47x 2 x 3 – 11,08x 1 2 – 3,2x 4 2 (3.6) Biến đổi o ii i i XX x    để đưa phương trình (3.5) và (3.6) về dạng thực tương ứng là : Y 1 = 38,41 + 7,18X 1 + 7,88X 2 + 0,30X 3 + 1,60X 4 – 0,28X 5 – 1,07X 1 X 2 – 0,03X 1 X 3 – 0,04X 1 X 4 – 0,03X 4 2 + 0,01X 5 2 (3.7) Y 2 = – 20,39 + 17,71X 1 + 4,43X 2 + 2,18X 3 + 0,68X 4 – 0,52X 1 X 2 – 0,02X 1 X 4 – 0,78X 2 X 3 – 1,23X 1 2 – 0,01X 4 2 (3.8) Kết quả phân tích cho thấy có ảnh hưởng tích cực nhất ở bậc 1 đến hiệu suất phản ứng chuyển đổi mỡ cá basa thành FAME là tỉ lệ mol metanol/mỡ cá (x 1 ), tiếp đến là thời gian phản ứng (x 4 ), nhiệt độ phản ứng (x 3 ), hàm lượng xúc tác (x 2 ) và hàm lượng dung môi (x 5 ) ít ảnh hưởng nhất. Ảnh hưởng tiêu cực nhất ở bậc hai đến hiệu suất phản ứng là nhân tố thời gian (x 4 2 ). Ảnh hưởng bậc hai của tỉ lệ mol metanol/mỡ cá (x 1 2 ), hàm lượng xúc tác (x 2 2 ), nhiệt độ phản ứng (x 3 2 ) và ảnh hưởng đồng thời của tỉ lệ mol metanol/mỡ cá và hàm lượng dung môi (x 1 x 5 ), hàm lượng xúc tác và nhiệt độ phản ứng (x 2 x 3 ), hàm lượng xúc tác và thời gian phản ứng (x 2 x 4 ), hàm lượng xúc tác và hàm lượng dung môi (x 2 x 5 ), nhiệt độ phản ứng và thời gian phản ứng (x 3 x 4 ), nhiệt độ phản ứng và hàm lượng dung môi (x 3 x 5 ), thời gian phản ứng và hàm lượng dung môi (x 4 x 5 ) đến hiệu suất phản ứng là không đáng kể. Sự tăng đồng thời của tỉ lệ mol metanol/mỡ cá với hàm lượng xúc tác (x 1 x 2 ), với nhiệt độ phản ứng (x 1 x 3 ) và với thời gian phản ứng (x 1 x 4 ) đều làm giảm hiệu suất phản ứng (ảnh hưởng tiêu cực). Kết quả này được thể hiện trên hình 3.14. Hình 3.14. Ảnh hưởng của các nhân tố khảo sát đến hiệu suất phản ứng chuyển đổi este xúc tác kiềm mỡ cá basa Đối với quá trình chuyển đổi hai giai đoạn, kết quả phân tích cho thấy cả bốn nhân tố khảo sát đều ảnh hưởng tích cực bậc 1 đến hiệu suất phản ứng este hóa xúc tác axit dầu hạt cao su trong đó tỉ lệ mol metanol/dầu (x 1 ) có ảnh hưởng tích cực nhiều nhất, tuy nhiên sự tăng đồng thời tỉ lệ metanol/dầu (x 1 x 1 ) lại làm giảm nhiều nhất hiệu suất phản ứng (ảnh hưởng tiêu cực). Các ảnh hưởng bậc hai và ảnh hưởng đồng thời của các nhân tố khảo sát cũng đều là ảnh hưởng tiêu cực. Sự thay đổi đồng thời hàm lượng xúc tác (x 2 x 2 ), thời gian phản ứng (x 3 x 3 ), tỉ lệ mol metanol/dầu và thời gian phản ứng (x 1 x 3 ), xúc tác và dung môi (x 2 x 4 ), thời gian phản ứng và dung môi (x 3 x 4 ) ít ảnh hưởng đến hiệu suất phản ứng. Ảnh hưởng của các nhân tố khảo sát đến hiệu suất phản ứng este hóa dầu hạt cao su được thể hiện trên hình 3.17. Hình 3.17. Sự biến thiên của hiệu suất phản ứng este hóa xúc tác axit dầu hạt cao su theo tỉ lệ mol metanol/dầu (3/1 – 9/1) và hàm lượng dung môi (10 – 50%) khi cố định hàm lượng xúc tác (1,25%) và thời gian phản ứng (3,5 giờ). Kết quả nghiên cứu cũng cho thấy độ tương quan giữa hiệu suất phản ứng chuyển đổi este mỡ cá basa thành FAME thu được từ thực nghiệm và tính toán từ PTHQ là 99,7%, của phản ứng este hóa xúc tác axit dầu hạt cao su là 98,7%. Như vậy giữa giá trị thu được từ mô hình dự đoán và giá trị thực nghiệm chỉ có một sự sai khác nhỏ. * Xác định điều kiện tối ƣu Dựa trên PTHQ xây dựng được, với sự hỗ trợ của phần mềm Modde 5.0, bằng phương pháp đạo hàm tìm được điều kiện tối ưu của quá trình chuyển đổi DMĐTV thành FAME từ mỡ cá basa: tỉ lệ mol metanol/mỡ cá: 6/1 (mol/mol); nồng độ xúc tác KOH: 0,98% (m/m); nhiệt độ phản ứng: 37,8 0 C; thời gian phản ứng: 23,5 phút; hàm lượng dung môi axeton: 30% (m/m). Ở điều kiện tối ưu này, hiệu suất chuyển đổi một giai đoạn các FAME dự đoán có giá trị là 99,9%, giá trị tương ứng thu được khi tiến hành thực nghiệm lặp lại ba lần với mẫu mỡ cá basa là 99,7%. Giá trị phương sai theo điều kiện tối ưu thu được là S tu 2 = 0,0233, theo mô hình là S 0 2 = 0,0351, F tn = 0,0351/0,0233 = 1,51 < F bảng (0,95, 2, 2) = 19, vậy độ chính xác của các số liệu theo mô hình và theo thực nghiệm ở điều kiện tối ưu là giống nhau. Điều kiện tối ưu của phản ứng este hóa xúc tác axit (giai đoạn một) của quá trình hai giai đoạn điều chế FAME từ dầu hạt cao su thu được như sau: tỉ lệ mol metanol/dầu: 7,5/1 (mol/mol); hàm lượng xúc tác H 2 SO 4 : 1,65% (m/m); thời gian phản ứng: 4 giờ; hàm lượng dung môi axeton: 35%; phản ứng thực hiện ở 60 o C. Ở điều kiện tối ưu này, hiệu suất phản ứng este hóa dự đoán có giá trị là 97,9%. Kết quả thực nghiệm thu được giá trị trung bình của 3 lần thí nghiệm ở điều kiện này là 98,3%. Giá trị phương sai theo điều kiện tối ưu thu được là S tu 2 = 0,034, của mô hình là S 0 2 = 0,350, F tn = 0,350/0,034 = 10,29 < F bảng (0,95, 2, 2) = 19, độ chính xác của các số liệu theo mô hình và theo thực nghiệm ở điều kiện tối ưu là giống nhau. Qua so sánh kết quả nghiên cứu với các kết quả chuyển đổi este DMĐTV theo hướng khảo sát đơn biến từng nhân tố với các loại xúc tác khác nhau cho thấy phản ứng chuyển đổi theo hướng tối ưu hóa đồng thời các nhân tố trong môi trường đồng thể đạt được hiệu suất cao hơn với lượng metanol cần dùng ít hơn, thời gian ngắn và nhiệt độ thấp hơn nhiều. Đó chính là ưu điểm của phương pháp tối ưu hóa chuyển đổi DMĐTV thành FAME này so với các nghiên cứu chuyển đổi trước đây. 3.3. Xác nhận giá trị sử dụng của phƣơng pháp phân tích 3.3.1. Xây dựng đƣờng chuẩn, xác định giới hạn phát hiện và giới hạn định lƣợng Để tiến hành xây dựng đường chuẩn dùng cho việc phân tích các FAME, một dãy chuẩn hỗn hợp 37 FAME có nồng độ từ 2–4–6 ppm đến 50–100–150 ppm được tiêm vào hệ GC/FID. Từ các kết quả phân tích, vẽ đồ thị sự phụ thuộc tỉ lệ diện tích pic FAME/IS vào nồng độ của các FAME thu được các đường chuẩn tương ứng. Dựa trên phương trình đường chuẩn: y = ax + b để tính LOD theo công thức: B LOD 3.S /b , trong đó S B : độ lệch chuẩn của mẫu trắng, nếu không làm thí nghiệm với mẫu trắng thì S B = S y (độ lệch chuẩn của phương trình hồi quy). LOQ = 3,3LOD. Kết quả phân tích được chỉ ra trên bảng 3.16. Từ kết quả phân tích cho thấy có sự tương quan tuyến tính chặt chẻ của tỉ lệ diện tích pic FAME/IS vào nồng độ các FAME (R  0,999) và hệ thống phân tích GC/FID có độ nhạy cao, phát hiện được các FAME đến 0,2 ppm và định lượng được đến 0,7 ppm. Bảng 3.16. Phương trình hồi quy, giới hạn phát hiện, giới hạn định lượng các FAME của phương pháp phân tích TT FAME Phƣơng trình hồi quy ( P =0,95) R LO D (pp m) LO Q (pp m) 1 C4:0 y = (0,2138  0,0334)x – (0,7961  1,6451) 0,9964 0,8 2,6 2 C6:0 y = (0,1066  0,0191)x – (0,3036  0,9336) 0,9953 1,6 5,2 3 C8:0 y = (0,2688  0,0347)x – (0,7646  1,7056) 0,9976 0,3 1,1 4 C10:0 y = (0,2776  0,0286)x – (0,5932  1,4071) 0,9984 0,7 2,5 5 C11:0 y = (0,2595  0,0220)x – (0,2424  0,5387) 0,9990 0,9 3,1 6 C12:0 y = (0,2646  0,0220)x – (0,4430  1,0873) 0,9990 1,0 3,3 7 C13:0 y = (0,2505  0,0216)x – (0,2174  0,5343) 0,9989 1,5 4,9 8 C14:0 y = (0,2599  0,0216)x – (0,4558  1,0707) 0,9990 1,1 3,5 9 C14:1 y = (0,2466  0,0159)x – (0,1023  0,3895) 0,9994 0,5 1,7 10 C15:0 y = (0,2512  0,0197)x – (0,1922  0,4856) 0,9991 0,9 3,1 11 C15:1 y = (0,2322  0,0146)x – (0,0498  0,3618) 0,9994 0,8 2,6 12 C16:0 y = (0,2581  0,0226)x – (0,5489  1,6721) 0,9989 2,0 6,5 [...]... jatropha 3.4.2 Xác định axit béo trong các mẫu DMĐTV Việt Nam Sau khi chuyển đổi hoàn toàn DMĐTV thành FAME, từ các điều kiện tối ưu tìm được của phương pháp GC/FID tách và xác định các FAME chúng tôi áp dụng để tách và xác định các axit béo trong 8 mẫu DMĐTV Việt Nam Kết quả phân tích thu được như sau: Trong mỡ cá basa có chứa đến 25 axit béo trong đó axit béo chiếm hàm lượng cao nhất là axit oleic (37%),... luận sau: 1 Nhằm xác định chính xác và chọn lọc đồng thời 37 axit béo, đã sử dụng mô hình hồi quy bậc hai tâm trực giao xây dựng phương trình hồi quy và đánh giá ảnh hưởng đồng thời các yếu tố đến khả năng tách và xác định các axit béo trong dầu mỡ động, thực vật Việt Nam 2 Trên cơ sở dữ kiện thực nghiệm đã tối ưu hóa điều kiện tách và xác định đồng thời 37 axit béo bằng phương pháp sắc ký khí detector... (25%) Hàm lượng các axit béo -3 là 7,6 mg/g chất béo, chiếm gần 1% và hàm lượng các axit béo -6 là 193 mg/g chất béo, chiếm 24% trong tổng số các axit béo Mỡ lợn chủ yếu chứa các axit béo mạch ngắn, không thấy sự có mặt của các axit béo mạch dài C22 và EPA, DHA như trong mỡ cá basa Trong dầu đậu nành phân tích được 15 axit béo, trong đó tổng hàm lượng các MUFA nhiều nhất (39%), tiếp đến là các PUFA (32%)... đến là dầu lạc, mỡ cá basa, mỡ lợn, dầu dừa ít có giá trị dinh dưỡng Dầu hạt cao su và dầu hạt jatropha có thể sử dụng làm nguyên liệu chính để sản xuất biodiesel KIẾN NGHỊ 1 Đề nghị sử dụng kết quả xây dựng quy trình phân tích và tiến tới phê duyệt thành phương pháp phân tích tiêu chuẩn xác định thành phần axit béo trong dầu mỡ động, thực vật 2 Áp dụng phương pháp xác định axit béo nghiên cứu để mở... là axit palmitic (28%), axit linoleic (14%), axit stearic (8%), axit myristic (5%), axit palmitoleic (1,4%), các axit béo khác đếu có hàm lượng thấp, dưới 1% Tổng các axit béo no (SFA) và axit béo không no một nối đôi (MUFA) gần xấp xỉ nhau, tương ứng là 43% và 38%, và gấp khoảng 2 lần các axit béo không no nhiều nối đôi (PUFA) (18%) Hàm lượng các axit béo -3 và -6 trong mỡ cá basa là 155 mg trong. .. là các SFA (29%) Dầu đậu nành chứa một lượng đáng kể các axit béo -3 (75 mg/g) và axit béo -6 (188 mg/g) Trong dầu mè có chứa 11 axit béo trong đó nhiều nhất là PUFA, 51%, chủ yếu là axit linoleic, tiếp đến là MUFA, 34%, chủ yếu là axit oleic và chỉ có 15% các SFA Dầu mè cũng chứa lượng đáng kể axit béo -6 (424 mg/g) và đặc biệt có chứa lượng bé thành phần của DHA (0,165 mg/g) mà các loại dầu thực. .. thực vật khác không có Trong dầu lạc có chứa 13 axit béo, nhiều nhất là các MUFA, 42%, chủ yếu là axit oleic, tiếp đến là các PUFA, 35%, chủ yếu là axit linoleic và chứa ít hơn các SFA, 22%, chủ yếu là axit palmitic Dầu lạc cũng chứa nhiều axit béo -6 (288 mg/g) nhưng lượng axit béo -3 thì rất ít (0,66 mg/g) Đối với dầu dừa thì khác với các dầu thực vật đã phân tích, thành phần chủ yếu lại là các. .. trò của DHA và ARA–Những lời khuyên đối với các bà mẹ mang thai và nuôi con bú”, Y học thực hành số 12, tr 92–95 26 Phạm Hùng Việt (2003), Cơ sở lý thuyết của phương pháp sắc ký khí, NXB Khoa học và Kỹ thuật Hà Nội, Hà Nội 27 Nguyễn Tường Vy, Trần Tử An, Đăng Đức Khanh, Trịnh Văn Lẩu (2005), Xác định các axit béo trong dầu màng hạt gấc bằng phương pháp sắc ký khí (GC/FID)”, Tạp chí kiểm nghiệm thuốc... độ chính xác của phương pháp nghiên cứu được chỉ ra trên bảng 3.28 Bảng 3.28 So sánh kết quả phân tích và kết quả kiểm nghiệm một số axit béo của mẫu mỡ cá basa và mẫu dầu mè Việt Nam Kết quả Kết quả kiểm Thông số phân tích nghiệm Vilas ttn (N=6) phân tích (mg/g) (mg/g) Mỡ cá basa Axit 227 238 2,570 palmitic Axit oleic 305 308 0,840 Axit linoleic 118 123 2,448 Dầu mè Axit 77 79 2,330 palmitic Axit oleic... phân tích Để xem xét tính chính xác của phương pháp nghiên cứu, chúng tôi tiến hành gửi kiểm nghiệm hàm lượng axit palmitic, axit oleic và axit linoleic trong một mẫu mỡ động vật (mỡ cá basa) và một mẫu dầu thực vật (dầu mè) đến phòng thí nghiệm được công nhận VILAS (theo ISO 17025:2005) (hệ thống công nhận phòng thí nghiệm của Việt Nam) , Viện Kiểm nghiệm An toàn Vệ sinh thực phẩm Quốc gia và so sánh . Nghiên cứu phương pháp xác định các axit béo trong một số loại dầu mỡ động, thực vật Việt Nam bằng kỹ thuật sắc ký Võ Thị Việt Dung . thành phương pháp phân tích tiêu chuẩn xác định thành phần axit béo trong dầu mỡ động, thực vật. 2. Áp dụng phương pháp xác định axit béo nghiên cứu để