食用大豆油儲藏過程中品質(zhì)變化的預(yù)測

唐瑞麗，袁先雯，馮燕玲，韓洪玲，董　月，高瑀瓏*，袁　建，汪海峰，鞠興榮

（南京財(cái)經(jīng)大學(xué)食品科學(xué)與工程學(xué)院，江蘇省現(xiàn)代糧食流通與安全協(xié)同創(chuàng)新中心，江蘇高校糧油質(zhì)量安全控制及深加工重點(diǎn)實(shí)驗(yàn)室，江蘇 南京 210023）

食用大豆油儲藏過程中品質(zhì)變化的預(yù)測

唐瑞麗，袁先雯，馮燕玲，韓洪玲，董月，高瑀瓏*，袁建，汪海峰，鞠興榮

（南京財(cái)經(jīng)大學(xué)食品科學(xué)與工程學(xué)院，江蘇省現(xiàn)代糧食流通與安全協(xié)同創(chuàng)新中心，江蘇高校糧油質(zhì)量安全控制及深加工重點(diǎn)實(shí)驗(yàn)室，江蘇 南京 210023）

以一級大豆油為原料，儲藏溫度設(shè)定為40 ℃，研究儲藏時(shí)間、抗氧化劑用量、水分含量對其品質(zhì)的影響。在單因素試驗(yàn)的基礎(chǔ)上，以儲藏時(shí)間、抗氧化劑用量、水分含量3 個(gè)因素為自變量，大豆油品質(zhì)指標(biāo)（過氧化值和氧化穩(wěn)定指數(shù)）為響應(yīng)值，采用響應(yīng)面Box-Behnken試驗(yàn)設(shè)計(jì)法，建立了大豆油品質(zhì)隨儲藏時(shí)間、抗氧化劑用量、水分含量變化的二次多項(xiàng)式預(yù)測模型。結(jié)果表明：儲藏時(shí)間、抗氧化劑用量、水分含量對大豆油的過氧化值和氧化穩(wěn)定性指數(shù)的影響顯著；大豆油品質(zhì)指標(biāo)預(yù)測模型的方差分析表明，模型皆極顯著（P＜0.001），實(shí)驗(yàn)誤差小，可以此模型來預(yù)測大豆油過氧化值和氧化穩(wěn)定指數(shù)隨儲藏時(shí)間、抗氧化劑用量、水分含量的變化，具有較高的擬合優(yōu)度。

大豆油；過氧化值；氧化穩(wěn)定指數(shù)；響應(yīng)面法；預(yù)測模型

大豆油原料豐富，營養(yǎng)價(jià)值高，是我國居民主要的食用油。大豆油含有較高的不飽和脂肪酸［1］，在儲藏、運(yùn)輸、消費(fèi)過程中，容易受外界條件諸如溫度、水分、氧氣等因素的影響，發(fā)生氧化酸敗，降低其品質(zhì)［2-3］，酸敗的植物油不僅口感差，營養(yǎng)價(jià)值低，甚至可能引起疾病，還可能致癌，給人體的健康帶來極大危害［4］。

為了延緩大豆油變質(zhì)，常用方法有低溫儲藏、避光、隔氧、添加抗氧化劑等。超市桶裝油一般都添加了叔丁基對苯二酚（tert-butylhydroquinone，TBHQ）。合成抗氧化劑的安全性一直受到質(zhì)疑，越來越多的天然抗氧化劑得到重視［5-6］，我國國家標(biāo)準(zhǔn)規(guī)定TBHQ在食用油中的添加量不超過0.2 g/kg。植物甾醇不僅具有抗氧化能力［7］，是一種新型的天然抗氧化劑，而且可以減少人體中的膽固醇水平［8］，被美國食品藥品監(jiān)督管理局認(rèn)可，可以作為食品抗氧化劑使用［9-10］。氧化穩(wěn)定指數(shù)（oxidative stability index，OSI）值可反映油脂氧化穩(wěn)定性，過氧化值（peroxide value，POV）是衡量油脂氧化初期的氧化程度的重要指標(biāo)［11］，為了提高食用大豆油的抗氧化性，降低TBHQ的在大豆油中添加量，本研究以添加天然植物甾醇和TBHQ為復(fù)合抗氧化劑，以儲藏時(shí)間、抗氧化劑用量、水分含量3 個(gè)因素為自變量，建立了大豆油儲藏過程中品質(zhì)指標(biāo)POV和OSI值的預(yù)測模型，以期為檢測食用油品質(zhì)安全提供重要理論依據(jù)。

我國夏季溫度較高，油品極易氧化變質(zhì)，因此，研究高溫儲藏條件下（40 ℃）大豆油品質(zhì)變化的規(guī)律非常必要。響應(yīng)面分析法是一種函數(shù)估算工具［12］，能準(zhǔn)確地用于研究各因素、因素之間與響應(yīng)值之間的相互關(guān)系［13-14］，目前已經(jīng)在工藝優(yōu)化、模型預(yù)測等方面得到越來越多的應(yīng)用［15］。應(yīng)用響應(yīng)面法預(yù)測研究大豆油品質(zhì)變化，為食用油在儲藏過程中的安全預(yù)測提供理論依據(jù)和技術(shù)支持，對研究食用油脂儲藏具有重要理論價(jià)值和實(shí)踐指導(dǎo)意義。

1　材料與方法

1.1材料與試劑

一級大豆油（未加任何抗氧化劑，澄清透明、水分含量0.05%、酸值0.08 mg KOH/g、POV1.18 mmol/kg）中儲糧鎮(zhèn)江糧油有限公司；TBHQ 廣東廣業(yè)清怡食品科技有限公司；植物甾醇（含量為95%，其中β-谷甾醇45%、豆甾醇20%、菜甾醇30%（菜籽甾醇和菜油甾醇）） 西安開來生物工程有限公司。

冰乙酸、碘化鉀、異辛烷、硫代硫酸鈉標(biāo)準(zhǔn)溶液、淀粉指示劑、無水硫酸鈉、甲醇（均為分析純） 南京丁貝生物技術(shù)有限公司；卡爾·費(fèi)休試劑 天津科密歐化學(xué)試劑有限公司；實(shí)驗(yàn)用水為蒸餾水。

1.2儀器與設(shè)備

TP-214型電子天平 北京丹佛儀器有限公司；GNP-9160型隔水式恒溫培養(yǎng)箱 上海三發(fā)科學(xué)儀器有限公司；HZQ-F160型全溫振蕩培養(yǎng)箱 太倉市實(shí)驗(yàn)設(shè)備廠；OSI-24型氧化穩(wěn)定性測定儀 美國Ominion公司；BIC-250型人工氣候箱 上海博訊實(shí)業(yè)有限公司醫(yī)療設(shè)備廠；101-3AS型電熱鼓風(fēng)干燥箱 上海蘇進(jìn)儀器設(shè)備廠；SB25-12D型超聲波清洗器 寧波新芝生物科技股份有限公司；電子萬用爐 北京永光明醫(yī)療儀器有限公司；卡爾·費(fèi)休水分測定儀 瑞士萬通中國；DS-1型高速組織搗碎機(jī) 上海標(biāo)本模型廠。

1.3方法

1.3.1一級大豆油樣品干燥處理

油樣干燥：無水硫酸鈉使用量為2%，加入一定量的一級大豆油中，充分混合，靜置1 d，取上層油樣，測定水分含量約為0，即為充分干燥的油樣。

1.3.2一級大豆油樣品抗氧化劑TBHQ和植物甾醇的添加

于棕色具塞瓶，分別加入一定量干燥的油樣，按照單因素試驗(yàn)和響應(yīng)面試驗(yàn)設(shè)計(jì)，抗氧化劑TBHQ和植物甾醇的添加均按照1∶1的比例加入大豆油樣品中，最高抗氧化劑用量為0.4 g/kg（即TBHQ與植物甾醇各為0.2 g/kg）。

1.3.3大豆油樣品的加水

準(zhǔn)確稱取1.3.1節(jié)所制備的添加了抗氧化劑的大豆油樣，置于高速組織搗碎機(jī)中，按照單因素試驗(yàn)和響應(yīng)面試驗(yàn)設(shè)計(jì)，加入適量的蒸餾水，均質(zhì)5 min，確保水和油脂充分均質(zhì)，快速裝入具塞棕色試劑瓶中。

1.3.4一級大豆油樣品的儲藏實(shí)驗(yàn)

將添加過抗氧化劑和水大豆油樣品，按照單因素試驗(yàn)和響應(yīng)面試驗(yàn)設(shè)計(jì)給定條件進(jìn)行恒溫40 ℃儲藏實(shí)驗(yàn)，盛油容器采用具塞棕色瓶500 mL，儲油量350 g，旋緊瓶塞，閉口儲藏。取樣方法：按照指標(biāo)測定方法要求取油樣，用一次性吸管快速取樣，然后旋緊瓶塞，混合均勻，放回培養(yǎng)箱。所有實(shí)驗(yàn)數(shù)據(jù)均重復(fù)測定3 次。

1.3.5指標(biāo)測定

1.3.5.1POV的測定

參照GB/T 5538—2005《動植物油脂：過氧化值測定》的方法來測定樣品的POV。

1.3.5.2OSI值的測定

參照GB/T 21121—2007《動植物油脂：氧化穩(wěn)定性的測定》的方法來測定樣品的OSI值。準(zhǔn)確稱取（5±0.2）g油樣于玻璃導(dǎo)管中，壓縮空氣壓力30 psi，OSI內(nèi)部調(diào)節(jié)器提供空氣壓力5.5 psi，測定溫度110 ℃，測試水量50 mL，平行測定3 次，取平均值。

1.3.5.3水分含量的測定

參照GB/T 26626—2011《動植物油脂水分含量測定：卡爾費(fèi)休法》（無吡啶）的方法測定樣品的原始含水量。

1.3.6單因素試驗(yàn)設(shè)計(jì)

以儲藏時(shí)間、抗氧化劑用量、水分含量3 個(gè)因素進(jìn)行單因素試驗(yàn)。大豆油水分含量0.05%、抗氧化劑含量0.2 g/kg（TBHQ和植物甾醇各0.1 g/kg）時(shí)，研究儲藏時(shí)間1、2、3、4、5、6 個(gè)月對大豆油POV和OSI值的影響；水分含量0.05%、儲藏時(shí)間6 個(gè)月時(shí)，研究抗氧化劑0、0.1、0.2、0.25、0.3、0.4 g/kg對POV和OSI值的影響；抗氧化劑含量0.2 g/kg（TBHQ和植物甾醇各0.1 g/kg）、儲藏時(shí)間6 個(gè)月時(shí)，研究水分含量0%、0.05%、0.1%、0.15%、0.20%、0.30%對POV和OSI值的影響。

1.3.7響應(yīng)面試驗(yàn)設(shè)計(jì)

在單因素試驗(yàn)基礎(chǔ)上，本研究采用Box-Behnken的設(shè)計(jì)原理［16-18］，進(jìn)行三因素三水平響應(yīng)面模型預(yù)測試驗(yàn)，以儲藏時(shí)間、抗氧化劑用量、水分含量3 個(gè)因素作為自變量，分別以X1、X2、X3表示，并以+1、0、-1分別代表自變量的高、中、低水平，按方程xi=（Xi-X0）/ΔX 對自變量進(jìn)行編碼。其中，xi為自變量的編碼值，Xi為自變量的真實(shí)值，X0為試驗(yàn)中心點(diǎn)處自變量的真實(shí)值，ΔX為自變量的變化步長，試驗(yàn)因素水平和編碼見表1。

表1　Box-Behnken試驗(yàn)因素與水平Table1　Codes and levels of factors used in Box-Behnken design

以大豆油過氧化值POV（Y1）和大豆油氧化穩(wěn)定性指數(shù)OSI值（Y2）作為響應(yīng)值。設(shè)POV（Y1）和OSI值（Y2）的預(yù)測模型由最小二乘法擬合的二次多項(xiàng)方程。

1.4數(shù)據(jù)處理

所有實(shí)驗(yàn)數(shù)據(jù)均重復(fù)測定3 次，利用Design-Expert（Version 6.0.5）軟件來設(shè)計(jì)試驗(yàn)，并對實(shí)驗(yàn)數(shù)據(jù)進(jìn)行模型擬合及方差分析；利用SPSS（Version 10.0）對數(shù)據(jù)進(jìn)行相關(guān)性分析。

2　結(jié)果與分析

2.1單因素試驗(yàn)結(jié)果

2.1.1儲藏時(shí)間對大豆油POV和OSI值的影響

圖1　儲藏時(shí)間對大豆油POV和OSI值的影響Fig.1　Effect of storage time on POV and OSI in soybean oil

由圖1可知，40 ℃儲藏，大豆油的POV隨著儲藏時(shí)間的延長顯著增加（P＜0.05），OSI值隨著儲藏時(shí)間的延長而顯著降低（P＜0.05）；儲藏6 個(gè)月時(shí)變化趨于平緩（P＞0.05），考慮到進(jìn)一步研究40 ℃儲藏大豆油品質(zhì)變化的預(yù)測模型，儲藏時(shí)間最長期限選擇6 個(gè)月。

2.1.2抗氧化劑用量對大豆油POV和OSI值的影響

圖2　抗氧化劑用量對大豆油POV和OSI值的影響Fig.2　Effects of the content of antioxidants on POV and OSI in soybean oil

由圖2可知，40 ℃儲藏6 個(gè)月，隨著抗氧化劑含量的增加，大豆油的POV明顯降低（P＜0.05），OSI值顯著增大（P＜0.05）；抗氧化劑含量為0.4 g/kg時(shí)，POV趨于平緩（P＞0.05）。綜合考慮抗氧化劑用量的作用效果和國標(biāo)規(guī)定的TBHQ在油脂中的最大使用添加量為0.2 g/kg，本研究選擇最高抗氧化劑是0.4 g/kg（TBHQ與植物甾醇各為0.2 g/kg）。吳時(shí)敏等［7］研究了植物甾醇在儲藏過程中對菜籽高級烹調(diào)油的抗氧化作用，發(fā)現(xiàn)植物甾醇添加量不大于0.10%（1.0 g/kg）時(shí)，隨著添加量的增大，抗氧化效果增強(qiáng)，與本研究的結(jié)果一致。

2.1.3水分含量對大豆油POV和OSI值的影響

圖3　水分含量對大豆油POV和OSI值的影響Fig.3　The effects of the moisture contents on POV and OSI in soybean

由圖3可知，40 ℃儲藏6 個(gè)月，隨著水分含量的增加，大豆油的POV明顯增大（P＜0.05），OSI值顯著降低（P＜0.05）；水分含量為0.3%時(shí)，大豆油POV變化不顯著（P＞0.05）?？紤]到大豆油容納的水分能力有限，保證大豆油與添加水分充分達(dá)到均一狀態(tài)，選擇最大水分含量為0.3%。

2.2預(yù)測模型的建立及其顯著性檢驗(yàn)

根據(jù)Design-Expert軟件所安排的試驗(yàn)組合，大豆油儲藏溫度為恒溫40 ℃，每組試驗(yàn)做3 個(gè)平行。試驗(yàn)結(jié)果見表2。

基于表2試驗(yàn)結(jié)果，運(yùn)用Design-Expert軟件進(jìn)行數(shù)據(jù)分析，得到大豆油POV和OSI值回歸模型的方差分析、因素的顯著性檢驗(yàn)結(jié)果，分別見表3～6。

表2　響應(yīng)面試驗(yàn)設(shè)計(jì)及其結(jié)果Table2　Box-Behnken design arrangement and corresponding experimental results

表3　POV回歸模型方差分析Table3　Analysis of variance for the regression equation for peroxide value

表4　OSI值回歸模型方差分析Table4　Analysis of variance for the regression equation for OSI

表5　POV回歸方程系數(shù)顯著性檢驗(yàn)Table5　Significance test for the regression coefficient for peroxide value

表6　OSI值回歸方程系數(shù)顯著性檢驗(yàn)Table6　Significance test for the regression coefficient for oxidative stability index

表5所對應(yīng)的POV預(yù)測模型方程為（1）：

由表3方差分析可知：POV預(yù)測模型方程的F= 36.47＞F0.01（9，4）=14.66（P＜0.000 1），表明方程（1）極顯著；失擬項(xiàng)F=2.26＜F0.05（9，3）= 8.81，失擬項(xiàng)P=0.223 6 ＞ 0.05，不顯著。方程（1）一次項(xiàng)（x1、x2、x3）和二次項(xiàng)對結(jié)果影響極顯著（P＜0.01）。方程的校正決定系數(shù)為0.952 3，表明該方程能解釋95.23%響應(yīng)值的變化，僅有總變異的4.77%不能用此方程來解釋；相關(guān)系數(shù)R2為0.979 1，表明模型有很好的擬合度，試驗(yàn)誤差小，該預(yù)測模型是合適的。

表6所對應(yīng)的OSI預(yù)測模型方程為（2）：

由表6方差分析可知，模型方程F值為19.42，P=0.000 4，表明該方程極顯著；失擬項(xiàng)F值為1.37＜F0.05（9，3）=8.81，失擬項(xiàng)P = 0.371 8＞0.05，不顯著。方程（2）一次項(xiàng)和二次項(xiàng)對OSI值的影響極顯著（P＜0.01）。方程的校正決定系數(shù)為0.912 0，表明該方程能解釋91.20%響應(yīng)值的變化；相關(guān)系數(shù)R2為0.961 5，表明模型擬合度好，試驗(yàn)誤差小。

2.3預(yù)測模型的驗(yàn)證

為了進(jìn)一步驗(yàn)證大豆油POV和OSI值預(yù)測模型方程（1）和（2）的合適性和有效性，進(jìn)行了6 組驗(yàn)證實(shí)驗(yàn)，結(jié)果見表7。

表7　POV和OSI值預(yù)測模型的驗(yàn)證結(jié)果Table7　Verified results from the prediction models for peroxide value and oxidative stability index

利用SPSS對表7試驗(yàn)數(shù)據(jù)進(jìn)行相關(guān)性分析，大豆油的POV實(shí)測值和預(yù)測值的相關(guān)系數(shù)為0.998 2，OSI的實(shí)測值和預(yù)測值的相關(guān)系數(shù)為0.992 4，表明模型方程（1）和（2）是合適有效的，能夠反映大豆油POV和OSI值與儲藏時(shí)間、抗氧化劑用量和水分含量的關(guān)系，可以用于40 ℃儲藏條件時(shí)大豆油POV與OSI值的理論預(yù)測，能較好地預(yù)測POV和OSI值隨研究因素的變化。崔楠等［19］在響應(yīng)面法研究中也通過對最優(yōu)條件進(jìn)行驗(yàn)證實(shí)驗(yàn)來證實(shí)數(shù)據(jù)參數(shù)的準(zhǔn)確性。于明等［20］在響應(yīng)面法研究中對優(yōu)化條件進(jìn)行5 組實(shí)驗(yàn)來驗(yàn)證所得數(shù)學(xué)模型的可靠性。本研究采用6 組不同數(shù)據(jù)，每組3 個(gè)平行，進(jìn)行驗(yàn)證實(shí)驗(yàn)，證實(shí)了預(yù)測模型是可靠的。

2.4大豆油POV的響應(yīng)面圖形分析

儲藏時(shí)間（X1）、抗氧化劑用量（X2）、水分含量（X3）及其交互作用對POV影響的響應(yīng)面見圖4。響應(yīng)曲面的陡峭程度可以反映該因素對響應(yīng)值影響的大小，曲面越陡，影響越大，曲面越平坦，影響越?。?1-22］。

圖4　任兩因素及其交互作用對POV影響的響應(yīng)面圖Fig.4　Response surface plots showing the effects of two experimental parameters and their interaction on peroxide value

由圖4可以看出，隨著儲藏時(shí)間的延長和水分含量的增加，大豆油的POV增大，隨著抗氧化劑含量的增加，大豆油的POV降低，該圖直觀分析與方差分析的結(jié)果相一致。

2.5大豆油OSI值的響應(yīng)面分析

儲藏時(shí)間（X1）、抗氧化劑用量（X2）、水分含量（X3）及其交互作用的OSI值響應(yīng)面圖如圖5所示。

圖5　任兩因素及其交互作用對OSI值影響的響應(yīng)面圖Fig.5　Response surface plots showing the effects of two experimental parameters and their interactive effects on oxidative stability index

由圖5響應(yīng)面曲線可以看出，隨著儲藏時(shí)間延長、水分含量的增加，OSI值顯著降低，大豆油穩(wěn)定性降低。隨抗氧化劑用量增大，OSI值增大，大豆油儲藏穩(wěn)定性在增加。該圖直觀分析與方差分析的結(jié)果相一致。

圖4和圖5結(jié)果均表明高溫條件下，儲藏時(shí)間曲面最陡，儲藏時(shí)間對響應(yīng)值的影響最大，可見，大豆油品在較高溫度條件下不能長期儲存。袁建等［23］在模擬油罐儲藏大豆油氧化穩(wěn)定性研究中指出對大豆油POV和OSI值的影響因素儲藏溫度＞儲藏時(shí)間＞含水量。王茜茜［24］在研究中指出各因素對菜籽油儲藏品質(zhì)影響的順序是：儲藏溫度＞儲藏時(shí)間＞含水量。這與本研究結(jié)果中對大豆油POV和OSI響應(yīng)值影響儲藏時(shí)間＞含水量結(jié)果相一致。為保障儲藏期大豆油品質(zhì)和貨架期，大豆油應(yīng)在較低溫度和較低水分含量條件下儲存，并添加適量的抗氧化劑。

3　結(jié) 論

在單因素試驗(yàn)的基礎(chǔ)上，采用Box-Behnken原理進(jìn)行試驗(yàn)設(shè)計(jì)，分別建立了大豆油POV和OSI值的預(yù)測模型方程，模型的試驗(yàn)數(shù)據(jù)擬合程度好，試驗(yàn)誤差小，可以分別用模型方程（1）和（2）分析和預(yù)測POV和OSI值隨儲藏時(shí)間（X1）、抗氧化劑用量（X2）和水分含量（X3）的變化；模型方程（1）和（2）的系數(shù)顯著性檢驗(yàn)表明，大豆油在40 ℃儲藏的過程中X1、X2和X3顯著地影響POV和OSI值；并對模型方程進(jìn)行驗(yàn)證，結(jié)果表明，模型方程（1）和（2）合適有效；響應(yīng)面直觀分析與方差分析的結(jié)果是一致的。POV是衡量油脂氧化初期的氧化程度的重要指標(biāo)，OSI值是反映油脂氧化穩(wěn)定性的指標(biāo)，對不同儲藏時(shí)期的大豆油的POV和OSI值進(jìn)行理論預(yù)測，這不僅為食用油品質(zhì)檢測提供理論依據(jù)，而且對保障食用油品質(zhì)安全具有重要意義。

［1］ GERTZ C， KLOSTERMANN S， KOCHHAR S P. Testing and comparing oxidative stability of vegetable oils and fats at frying temperature［J］. European Journal of Lipid Science and Technology，2000， 102（8/9）: 543-551. DOI:10.1002/1438-9312（200009）102:8/9＜543::AID-EJLT543＞3.0.CO；2-V.

［2］ LERCKER G， RODRIGUEZ M T. Cholesterol oxidation mechanism［M］//GUARDIOLA F， DUTTA P C， CODONY R， et al. Cholesterol and phytosterol oxidation products: analysis， occurrence and biological effects. AOCS Press， Champaign， 2002: 1-26.

［3］ FARHOOSH R， EINAFSHAR S， SHARAYEI P. The effect of commercial refining steps on the rancidity measures of soybean and canola oils［J］. Food Chemistry， 2009， 115（3）: 933-938. DOI:10.1016/ j.foodchem.2009.01.035.

［4］ LERCKER G， RODRIGUEZ M T. Cholesterol oxidation mechanism［M］//GUARDIOLA F， DUTTA P C， CODONY R， et al. Cholesterol and phytosterol oxidation products: analysis， occurrence and biological effects. Champaign: AOCS Press， 2002: 1-26.

［5］ HOU Dexing. Potential mechanism of cancer chemoprevention by anthocyanin［J］. Current Molecular Medicines， 2003， 3（2）: 149-159. DOI:10.2174/1566524033361555.

［6］ ZHANG Ying， YANG Lei， ZU Yuangang， et al. Oxidative stability of sunflower oil supplemented with carnosic acid compared with synthetic antioxidants during accelerated storage［J］. Food Chemistry，2010， 118（3）: 656-662. DOI:10.1016/j.foodchem.2009.05.038.

［7］ 吳時(shí)敏， 吳謀成， 馬莉. 植物甾醇在菜籽高級烹調(diào)油中的抗氧化作用（Ⅰ）: 常溫下抗氧化作用的研究［J］. 中國油脂， 2003， 28（4）: 52-54. DOI:10.3321/j.issn:1003-7969.2003.04.017.

［8］ DEMONTY I， RAS R T， van der KNAAP H C， et al. Continuous dose-response relationship of the LDL-cholesterol-lowering effect of phytosterol intake［J］. The Journal of Nutrition， 2009， 139（2）: 205-207. DOI:10.3945/jn.108.095125.

［9］ CHOUDHARY S P， TRAN L S. Phytosterols: perspectives in human nutrition and clinical therapy［J］. Current Medicinal Chemistry， 2011，18（29）: 4557-4567. DOI:10.2174/092986711797287593.

［10］ CARR T P， ASH M M， BROWN A W. Cholesterol-lowering phytosterols: factors affecting their use and efficacy［J］. Nutrition and Dietary Supplements， 2010， 2: 59-72. DOI:10.2147/NDS.S10974.

［11］ 闞建全. 食品化學(xué)［M］. 2版. 北京: 中國農(nóng)業(yè)大學(xué)出版社， 2002: 112-134.

［12］ RESCONI V C， KEENAN D F， GOUGH S， et al. Response surface methodology analysis of rice starch and fructo-oligosaccharides as substitutes for phosphate and dextrose in whole muscle cooked hams［J］. LWT-Food Science and Technology， 2015， 64（2）: 946-958. DOI:10.1016/j.lwt.2015.06.053.

［13］ MOHAMMADI R， MOHAMMADIFAR M A， MORTAZAVIAN A M， et al. Extraction optimization of pepsin-soluble collagen from eggshell membrane by response surface methodology（RSM）［J］. Food Chemistry， 2016， 190: 186-193. DOI:10.1016/ j.foodchem.2015.05.073.

［14］ NAHR F K， MOKARRAM R R， HEJAZI M A， et al. Optimization of the nanocellulose based cryoprotective medium to enhance the viability of freeze dried Lactobacillus plantarum using response surface methodology［J］. LWT-Food Science and Technology， 2015，64（1）: 326-332. DOI:10.1016/j.lwt.2015.06.004.

［15］ KOGKAKI E A， NATSKOULIS P I， PANAGOU E Z. Modeling the effect of natamycin， pine-resin and environmental factors on the growth and OTA production by Aspergillus carbonarius using response surface methodology［J］. Food Research International， 2016，79: 19-28. DOI:10.1016/j.foodres.2015.11.023.

［16］ GAO Yulong. Efficacy of high hydrostatic pressure and mild heat to reduce Geobacillus stearothermophilus as 1.1923 spores in model food systems［J］. Journal of Food Safety， 2010， 30（1）: 124-140. DOI:10.1111/j.1745-4565.2009.00194.x.

［17］ GAO Yulong， QIU Weifen， WU Ding， et al. Assessment of Clostridium perfringens spore response to high hydrostatic pressure and heat with nisin［J］. Applied Biochemistry and Biotechnology，2011， 164（7）: 1083-1095. DOI:10.1007/s12010-011-9196-0.

［18］ GAO Yulong， JU Xingrong. Inactivation of Bacillus coagulans spores subjected to combinations of high-pressure processing and nisin［J］. Transactions of the ASABE， 2011， 54（1）: 385-392. DOI:10.13031/2013.36235

［19］ 崔楠， 陶曉赟， 李娟， 等. 響應(yīng)面法優(yōu)化黃粉蟲蛋白制備ACE抑制肽的條件［J］. 食品科學(xué)， 2014， 35（15）: 156-160. DOI:10.7506/spkx1002-6630-201415032.

［20］ 于明， 朱文學(xué)， 郭菡， 等. 響應(yīng)面分析法優(yōu)化超聲提取槐豆膠工藝［J］.食品科學(xué)， 2013， 34（2）: 114-118.

［21］ 牟建樓， 王頡， 孫劍鋒. 響應(yīng)面法優(yōu)化纖維素酶提取蘋果渣中水溶性膳食纖維［J］. 食品科學(xué)， 2012， 33（8）: 95-98.

［22］ 王若蘭， 牛會敏， 鄧立陽， 等. 花生殼水溶性膳食纖維超聲提取工藝響應(yīng)面優(yōu)化［J］. 河南工業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2013， 34（1）: 6-11.

［23］ 袁建， 何海艷， 何榮， 等. 模擬油罐儲藏大豆油氧化穩(wěn)定性研究［J］. 中國糧油學(xué)報(bào)， 2013， 28（3）: 92-98. DOI:10.3969/ j.issn.1003-0174.2013.03.018.

［24］ 王茜茜. 菜籽油儲藏穩(wěn)定性研究［D］. 南京: 南京財(cái)經(jīng)大學(xué)， 2012.

Prediction of the Quality Changes of Soybean Oil during Storage

TANG Ruili， YUAN Xianwen， FENG Yanling， HAN Hongling， DONG Yue， GAO Yulong*， YUAN Jian， WANG Haifeng， JU Xingrong
（Collaborative Innovation Center for Modern Grain Circulation and Safety， Key Laboratory of Grains and Oils Quality Control and Processing， College of Food Science and Engineering， Nanjing University of Finance and Economics， Nanjing 210023， China）

In this investigation， the effects of storage time， the amount of antioxidants and moisture content on the quality of soybean oil （commercially available， grade 1） were studied during storage at 40 ℃. These three factors were investigated by one-factor-at-a-time method and then they were taken as independent variables for further optimization by response surface methodology using Box-Behnken design. The quality parameters peroxide value and oxidative stability index were selected as response variables. A quadratic polynomial predictive model for each response was developed. The results showed that all the three factors had a significant impact on peroxide value and oxidative stability index in soybean oil. Analysis of variance demonstrated that the two models were highly statistically significant （P ＜ 0.001） and had small error. They had high goodness of fit at the same time.

soybean oil； peroxide value （POV）； oxidative stability index （OSI）； response surface methodology； prediction model

10.7506/spkx1002-6630-201618041

TS221

A

1002-6630（2016）18-0256-06

唐瑞麗， 袁先雯， 馮燕玲， 等. 食用大豆油儲藏過程中品質(zhì)變化的預(yù)測［J］. 食品科學(xué)， 2016， 37（18）: 256-261. DOI:10.7506/spkx1002-6630-201618041. http://www.spkx.net.cn

TANG Ruili， YUAN Xianwen， FENG Yanling， et al. Prediction of the quality changes of soybean oil during storage［J］. Food Science， 2016， 37（18）: 256-261. （in Chinese with English abstract） DOI:10.7506/spkx1002-6630-201618041. http://www.spkx.net.cn

2015-12-31

糧食公益性行業(yè)科研專項(xiàng)（201413007-05）；江蘇省自然科學(xué)基金項(xiàng)目（BK20131435）；高校自然科學(xué)研究項(xiàng)目（12KJB550003）；國家自然科學(xué)基金面上項(xiàng)目（31371864）；江蘇高校優(yōu)勢學(xué)科建設(shè)工程資助項(xiàng)目（PADA）

唐瑞麗（1989—），女，碩士研究生，研究方向?yàn)橛椭瑑Σ嘏c生物技術(shù)。E-mail：tangruilishipin@qq.com

高瑀瓏（1974—），男，教授，博士，研究方向?yàn)樯锛夹g(shù)。E-mail：yulonggao19762001@163.com