Tổng quan nghiên cứu

Cơm dừa là nguồn nguyên liệu nông sản chiến lược tại vùng Đồng bằng sông Cửu Long, chiếm khoảng 22% tổng khối lượng trái dừa khô bên cạnh 25% nước dừa, 35% vỏ dừa và 12% gáo dừa. Hai chỉ tiêu cốt lõi quyết định giá trị kinh tế và chất lượng thương phẩm của cơm dừa là hàm lượng lipid tổng và hàm lượng ẩm. Trong thực tế, dải biến thiên của hàm ẩm dao động từ 0,39% đến 82,21% và hàm lượng lipid dao động từ 2,8867% đến 68,4363%. Tuy nhiên, việc thu mua cơm dừa phục vụ sản xuất dầu dừa, sữa dừa hay cơm dừa nạo sấy tại các cơ sở chế biến hiện nay chủ yếu vẫn dựa trên kinh nghiệm cảm quan, dễ dẫn đến tranh chấp thương mại và sai lệch chất lượng.

Các phương pháp phân tích phòng thí nghiệm theo tiêu chuẩn như trích ly Soxhlet hay sấy đến khối lượng không đổi thường đòi hỏi thời gian kéo dài từ 4 đến 8 giờ, tiêu tốn lượng lớn dung môi diethyl ether độc hại và chi phí vận hành cao. Do đó, nghiên cứu này được thực hiện nhằm xây dựng phương pháp đánh giá nhanh, không phá hủy mẫu thông qua thiết bị quang phổ cận hồng ngoại cầm tay SCiO trong dải bước sóng 740 – 1070 nm kết hợp thuật toán toán hóa học đa biến. Đề tài được triển khai trên 345 mẫu cơm dừa thu thập tại 4 tỉnh Bến Tre, Cần Thơ, Trà Vinh và Vĩnh Long trong giai đoạn từ tháng 7/2016 đến tháng 6/2017. Nghiên cứu mang lại ý nghĩa thực tiễn to lớn khi rút ngắn thời gian phân tích xuống dưới 1 phút cho mỗi mẫu, đạt hệ số tương quan suy đoán vượt 0,93 và cắt giảm tối đa chi phí hóa chất kiểm nghiệm.

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 lý thuyết quang phổ cận hồng ngoại (NIRS), kỹ thuật toán hóa học Chemometrics và các mô hình hồi quy đa biến tuyến tính:

  • Lý thuyết quang phổ cận hồng ngoại: Vùng quang phổ cận hồng ngoại trải dài từ bước sóng 750 – 2500 nm, ghi nhận sự hấp thụ bức xạ tương ứng với các dao động âm bội và dao động kết hợp của các liên kết phân tử C-H, O-H và N-H. Trong dải bước sóng hẹp 740 – 1070 nm, các đỉnh hấp thụ đặc trưng của nước và liên kết béo trong mô sinh học thể hiện rõ nét, cho phép trích xuất thông tin định lượng mà không cần phá hủy cấu trúc mẫu.
  • Mô hình hồi quy bình phương tối thiểu từng phần (PLS): Thuật toán PLS đóng vai trò cầu nối toán học giữa ma trận dữ liệu phổ đa biến và ma trận nồng độ thuộc tính phân tích thực nghiệm, giúp giải quyết triệt để vấn đề chồng chéo tín hiệu quang phổ.
  • Phương pháp tiền xử lý dữ liệu phổ: Ứng dụng biến thiên đã chuẩn hóa (SNV) nhằm loại bỏ hiện tượng tán xạ ánh sáng và hiệu ứng kích thước hạt, kết hợp cùng phép biến đổi đạo hàm bậc một Savitzky-Golay để triệt tiêu độ trôi đường nền, tăng cường độ phân giải của các vân hấp thụ.
  • Phân tích nhân tố khám phá (EFA): Được sử dụng để tối ưu hóa số lượng biến tiềm ẩn, cân bằng giữa khả năng giải thích phương sai tích lũy và việc kiểm soát sai số dự đoán.

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

Nghiên cứu sử dụng tổng cộng 345 mẫu cơm dừa từ các trái dừa trưởng thành từ 12 tháng tuổi trở lên tại Bến Tre, Cần Thơ, Trà Vinh và Vĩnh Long. Toàn bộ mẫu được chia ngẫu nhiên thành hai tập dữ liệu độc lập: 300 mẫu chuẩn dùng để huấn luyện mô hình hiệu chuẩn và 45 mẫu kiểm chứng dùng để đánh giá độ chuẩn xác ngoại suy.

Phương pháp chọn mẫu bao trùm các mức độ sấy ẩm (sấy ở 50 độ C trong các khoảng 1, 2, 3, 4 giờ đến khối lượng không đổi) và các mức độ ép tách béo khác nhau nhằm tạo ra dải nồng độ liên tục, phong phú. Phương pháp phân tích đối chứng bao gồm: xác định độ ẩm theo TCVN 8949:2011 (ISO 665:2000) bằng cân sấy ẩm hồng ngoại Scaltec và định lượng lipid tổng theo TCVN 8948:2011 (ISO 659:2009) bằng phương pháp trích ly Soxhlet với dung môi diethyl ether.

Phổ cận hồng ngoại được thu nhận bằng thiết bị SCiO (Consumer Physics – Mỹ) trên 331 điểm bước sóng từ 740 đến 1070 nm, thực hiện quét 15 lần trên mỗi mẫu để lấy giá trị phổ trung bình. Toàn bộ ma trận dữ liệu được xử lý trên phần mềm Unscrambler X phiên bản 10.5. Việc lựa chọn phương pháp PLS kết hợp tiền xử lý quang học được xem là giải pháp tối ưu vì vừa bảo toàn được 99% lượng thông tin phổ hữu ích, vừa loại trừ các tác nhân gây nhiễu vật lý từ bề mặt thô nhám của cơm dừa nạo.

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ô hình hóa đa biến đã đem lại các kết quả định lượng cụ thể:

  • Xây dựng dải nồng độ mẫu chuẩn toàn diện: Tập 300 mẫu chuẩn sở hữu hàm lượng lipid trung bình 31,4173% ± 16,3171% (dao động từ 2,8867% đến 68,4363%) và hàm ẩm trung bình 34,2338% ± 24,3137% (dao động từ 0,39% đến 82,21%). Tập 45 mẫu kiểm chứng có giá trị nằm hoàn toàn trong phạm vi đường chuẩn với hàm lượng lipid từ 5,7078% đến 66,0093% và độ ẩm từ 1,68% đến 75,30%.
  • Tối ưu hóa mô hình định lượng lipid tổng: Mô hình PLS đạt hiệu quả cao nhất ở 11 factors tối ưu, cho hệ số tương quan hiệu chuẩn r = 0,9367 (R² = 0,8775), sai số căn bậc hai trung bình hiệu chuẩn RMSEC = 5,7020%, sai số kiểm chứng chéo RMSECV = 6,1698% và sai số dự đoán trên tập kiểm chứng RMSEP = 4,6531%. Ở 36 factors, phương sai tích lũy của dải phổ giải thích được tới 99,9987% thông tin.
  • Tối ưu hóa mô hình định lượng hàm ẩm: Mô hình suy đoán độ ẩm đạt độ chính xác vượt trội với 7 factors tối ưu, hệ số tương quan r = 0,9613, sai số RMSEC = 6,6839%, RMSECV = 7,0562% và sai số dự đoán RMSEP = 5,6975%. Số lượng 10 factors giải thích được 92,5644% biến thiên của hàm ẩm và 99,4136% thông tin quang phổ.
  • Khẳng định độ tin cậy qua kiểm định thống kê: Kiểm định t-Student và chuẩn Fisher so sánh giữa kết quả dự đoán của thiết bị NIR cầm tay và phương pháp phân tích tiêu chuẩn trên 45 mẫu thử nghiệm độc lập cho thấy sự tương đồng cao, không có sai khác mang ý nghĩa thống kê (p > 0,05).

Thảo luận kết quả

Mối tương quan tuyến tính chặt chẽ giữa dữ liệu quang phổ và hai thành phần chất béo, độ ẩm xuất phát từ cơ chế hấp thụ năng lượng quang học đặc thù của các liên kết phân tử. Tại vùng bước sóng gần 970 nm, các liên kết O-H của phân tử nước hấp thụ mạnh mẽ, trong khi dải sóng 920 – 930 nm ghi nhận dao động âm bội bậc ba của liên kết C-H trong các chuỗi acid béo no và không no của cơm dừa.

Dữ liệu thực nghiệm được trực quan hóa xuất sắc thông qua các biểu đồ phân tán (scatter plot) giữa giá trị đo thực tế và giá trị suy đoán, phản ánh độ chụm cao xung quanh đường hồi quy tuyến tính chuẩn. Biểu đồ diễn biến sai số RMSECV theo từng factor giúp xác lập chính xác điểm dừng tối ưu: việc chọn dưới 10 factors ở chỉ tiêu lipid gây ra hiện tượng thiếu khớp (underfitting) do mất dữ liệu phổ, trong khi chọn trên 12 factors khiến RMSECV tăng trở lại từ 6,1698% lên 6,6035% do mô hình bị quá khớp (overfitting) bởi các tín hiệu nhiễu.

Kết quả này hoàn toàn tương thích với các nghiên cứu quốc tế tiêu biểu như nghiên cứu của J. Wold trên thịt bò (r = 0,98) hay F. Pedreschi trên khoai tây chiên (r = 0,99). Điều này chứng minh thiết bị NIR cầm tay dạng nhỏ gọn kết hợp Chemometrics hoàn toàn đủ năng lực thay thế các phép đo hóa ướt truyền thống trong kiểm soát chất lượng nông sản.

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

Nhằm chuyển hóa các kết quả nghiên cứu thành giá trị kinh tế thực tiễn cho ngành dừa, các khuyến nghị chiến lược được đề xuất như sau:

  • Đầu tư trang bị máy quang phổ NIR cầm tay tại hệ thống trạm thu mua: Các doanh nghiệp sản xuất dầu dừa, cơm dừa nạo sấy tại Đồng bằng sông Cửu Long cần triển khai trang bị thiết bị NIR cầm tay tại 100% điểm cân thu mua trong nửa đầu năm 2027. Mục tiêu rút ngắn thời gian kiểm nghiệm nguyên liệu từ 6 giờ xuống dưới 1 phút/mẫu, giảm thiểu 85% chi phí dung môi hóa chất kiểm tra.
  • Xây dựng quy trình tái hiệu chuẩn mô hình hồi quy định kỳ: Bộ phận Đảm bảo chất lượng (QA/QC) của các nhà máy thực phẩm cần định kỳ 6 tháng một lần thực hiện cập nhật ma trận phổ với tối thiểu 50 mẫu dừa theo từng vụ mùa thu hoạch, duy trì sai số dự đoán RMSEP luôn ổn định dưới ngưỡng 5,0%.
  • Chuẩn hóa hướng dẫn thao tác xử lý mẫu tại hiện trường: Ban kỹ thuật sản xuất cần ban hành quy chuẩn đồng nhất mẫu cơm dừa (nghiền nhỏ đồng đều, kiểm soát nhiệt độ mẫu đo trong khoảng 20 – 25 độ C) áp dụng từ quý 3 năm 2027 nhằm nâng cao độ lặp lại của phép quét phổ lên trên 98,5%.
  • Hoàn thiện khung định giá thu mua nông sản theo hàm lượng lipid thực tế: Hiệp hội Dừa Việt Nam cần phối hợp với các cơ quan quản lý ban hành bảng biểu giá thu mua linh hoạt dựa trên hàm lượng lipid định lượng bằng thiết bị đo nhanh, thay thế hoàn toàn phương pháp phân loại cảm quan trong năm 2028, giúp gia tăng giá trị chuỗi cung ứng thêm 12 – 15%.

Đố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 thiết thực cho nhiều nhóm đối tượng chuyên môn:

  • Doanh nghiệp sản xuất và chế biến dừa: Các nhà máy ép dầu dừa nguyên chất, chế biến sữa dừa và cơm dừa nạo sấy xuất khẩu có thể ứng dụng trực tiếp thiết bị NIR cầm tay để nghiệm thu nhanh nguyên liệu đầu vào, giảm thiểu thất thoát chất lượng và tiết kiệm hàng trăm triệu đồng chi phí phân tích hàng năm.
  • Chuyên viên QA/QC và kỹ sư công nghệ thực phẩm: Tài liệu cung cấp quy trình chi tiết về kỹ thuật xử lý dữ liệu Chemometrics trên Unscrambler X 10.5, phương pháp biến đổi phổ SNV và đạo hàm Savitzky-Golay, tạo tiền đề xây dựng mô hình kiểm tra nhanh cho các nông sản khác như hạt điều, cà phê, ca cao.
  • Giảng viên, nghiên cứu sinh và học viên cao học: Luận văn là tài liệu tham khảo chuẩn mực về phương pháp kết hợp giữa hóa phân tích thực nghiệm và toán tin ứng dụng, cung cấp dữ liệu đối sánh chi tiết theo các tiêu chuẩn TCVN 8948:2011 và TCVN 8949:2011.
  • Cơ quan kiểm nghiệm và quản lý chất lượng nông nghiệp: Cung cấp cơ sở khoa học vững chắc để các trung tâm kỹ thuật tiêu chuẩn đo lường chất lượng xem xét, ban hành các quy trình kiểm tra nhanh tại hiện trường phục vụ công tác giám sát chuỗi cung ứng thực phẩm an toàn.

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

Phương pháp quang phổ cận hồng ngoại NIR cầm tay có độ chính xác tương đương phương pháp Soxhlet truyền thống không? Phương pháp NIR kết hợp mô hình PLS đạt độ đúng và độ lặp lại tương đương với phương pháp Soxhlet tiêu chuẩn. Kiểm định Fisher và t-Student trên 45 mẫu kiểm chứng độc lập cho thấy sự tương đồng cao với hệ số tương quan r = 0,9367 và sai số dự đoán RMSEP chỉ 4,6531%, hoàn toàn đáp ứng các tiêu chuẩn kiểm nghiệm công nghiệp.

Tại sao cần phải thực hiện tiền xử lý phổ bằng đạo hàm bậc một và chuẩn hóa biến thiên SNV? Do mẫu cơm dừa nạo có bề mặt gồ ghề và kích thước hạt không đồng nhất, ánh sáng phản xạ bị tán xạ mạnh. Phép biến đổi SNV loại bỏ hiện tượng sai lệch do độ tán xạ ánh sáng, còn đạo hàm bậc một Savitzky-Golay triệt tiêu trôi đường nền, giúp làm nổi bật các đỉnh hấp thụ đặc trưng của lipid và nước trong dải bước sóng 740 – 1070 nm.

Hiện tượng quá khớp (overfitting) và thiếu khớp (underfitting) trong mô hình PLS được kiểm soát như thế nào? Nghiên cứu kiểm soát hiện tượng này thông qua giá trị sai số kiểm chứng chéo RMSECV. Ở chỉ tiêu lipid, việc sử dụng dưới 10 factors làm mất thông tin hữu ích (underfitting), trong khi vượt quá 12 factors khiến RMSECV tăng từ 6,1698% lên 6,6035% do mô hình hấp thụ cả tín hiệu nhiễu (overfitting). Mức 11 factors được chứng minh là điểm tối ưu cân bằng hoàn hảo.

Dải đo của mô hình có bao quát được mọi loại cơm dừa trên thị trường hiện nay không? Mô hình bao phủ dải nồng độ rất rộng: hàm lượng lipid từ 2,8867% đến 68,4363% và độ ẩm từ 0,39% đến 82,21%. Nhờ đó, phương pháp đáp ứng tốt mọi trạng thái mẫu thực tế, từ cơm dừa tươi có độ ẩm cao (trên 80%), cơm dừa trưởng thành (36 – 41% béo) cho đến copra sấy khô phục vụ ép dầu (trên 60% béo).

Thời gian và chi phí cho mỗi lần đo mẫu bằng thiết bị NIR cầm tay SCiO là bao nhiêu? Mỗi mẫu cơm dừa chỉ mất từ 30 đến 45 giây để hoàn tất 15 lần quét phổ trung bình và hiển thị ngay kết quả độ ẩm, chất béo trên màn hình. Thiết bị vận hành không sử dụng hóa chất hay dung môi hữu cơ độc hại, giúp giảm trên 80% chi phí vận hành so với phương pháp trích ly hóa học mất từ 4 đến 6 giờ.

Kết luận

  • Xây dựng thành công quy trình phân tích nhanh đồng thời hàm lượng lipid tổng và hàm lượng ẩm trong cơm dừa bằng thiết bị quang phổ cận hồng ngoại NIR cầm tay SCiO trong dải bước sóng 740 – 1070 nm.
  • Xác lập mô hình hồi quy PLS tối ưu với 11 factors cho hàm lượng lipid (r = 0,9367, RMSEP = 4,6531%) và 7 factors cho hàm lượng ẩm (r = 0,9613, RMSEP = 5,6975%).
  • Ứng dụng xuất sắc kỹ thuật tiền xử lý phổ kết hợp đạo hàm bậc một Savitzky-Golay và chuẩn hóa SNV trên bộ dữ liệu 300 mẫu chuẩn và 45 mẫu kiểm chứng.
  • Chứng minh phương pháp đo nhanh đạt độ tin cậy tương đương các phương pháp tiêu chuẩn quốc gia TCVN 8948:2011 và TCVN 8949:2011 nhưng rút ngắn thời gian phân tích từ vài giờ xuống dưới 1 phút.
  • Loại bỏ hoàn toàn việc sử dụng dung môi hữu cơ độc hại diethyl ether, đóng góp giải pháp công nghệ xanh thân thiện với môi trường và an toàn cho người thao tác.
  • Lộ trình phát triển tiếp theo trong giai đoạn 2027 – 2028: Mở rộng mô hình hiệu chuẩn cho các chỉ tiêu acid béo tự do (FFA), protein và phát triển thuật toán phân tích tích hợp trên nền tảng điện toán đám mây.
  • Hãy chủ động tiếp cận và ứng dụng công nghệ quang phổ NIR cầm tay ngay hôm nay để tối ưu hóa quy trình kiểm soát chất lượng, nâng cao hiệu quả kinh tế và nâng tầm vị thế nông sản dừa Việt Nam trên trường quốc tế!