面包蝦中副溶血性弧菌在pH影響下生長(zhǎng)預(yù)測(cè)模型的建立

賀 群，黃 和，*，劉 亞，蔣志紅，勵(lì)建榮

（1.廣東省水產(chǎn)品加工與安全重點(diǎn)實(shí)驗(yàn)室，廣東普通高等學(xué)校水產(chǎn)品深加工重點(diǎn)實(shí)驗(yàn)室，廣東海洋大學(xué)食品科技學(xué)院，廣東湛江 524088；2.渤海大學(xué)，遼寧省食品安全重點(diǎn)實(shí)驗(yàn)室，遼寧錦州 121013）

面包蝦中副溶血性弧菌在pH影響下生長(zhǎng)預(yù)測(cè)模型的建立

賀 群1，黃 和1，*，劉 亞1，蔣志紅1，勵(lì)建榮2

（1.廣東省水產(chǎn)品加工與安全重點(diǎn)實(shí)驗(yàn)室，廣東普通高等學(xué)校水產(chǎn)品深加工重點(diǎn)實(shí)驗(yàn)室，廣東海洋大學(xué)食品科技學(xué)院，廣東湛江 524088；2.渤海大學(xué)，遼寧省食品安全重點(diǎn)實(shí)驗(yàn)室，遼寧錦州 121013）

為探討pH對(duì)面包蝦生產(chǎn)過程中副溶血性弧菌的生長(zhǎng)規(guī)律的影響，測(cè)定了不同pH（5.5、6.0、6.5、7.0、7.5、8.0、8.5）下副溶血性弧菌的生長(zhǎng)曲線，運(yùn)用Origin 8.0軟件和DMFit軟件，分別采用修正Gompertz模型、修正Logistic模型、Baranyi模型對(duì)生長(zhǎng)曲線進(jìn)行擬合，選擇決定系數(shù)最高的方程建立初級(jí)模型，并求出最大比生長(zhǎng)速率。應(yīng)用修正Arrhenius模型擬合最大比生長(zhǎng)速率和pH之間的曲線建立二級(jí)模型。結(jié)果表明，通過模型檢驗(yàn)比較，修正Compertz模型的擬合度較好（決定系數(shù)=0.99），進(jìn)一步擬合所得二級(jí)模型擬合度較好（決定系數(shù)=0.98），偏差因子為0.99，準(zhǔn)確因子1.0，均在可接受的范圍內(nèi)，表明模型能較好的預(yù)測(cè)面包蝦中副溶血性弧菌在pH為5.5～8.5環(huán)境下的生長(zhǎng)動(dòng)態(tài)。

面包蝦，副溶血性弧菌，生長(zhǎng)模型

副溶血性弧菌（Vibrio parahaemolyticus）是革蘭氏陰性兼性厭氧菌，具有嗜鹽性，最適生長(zhǎng)pH是7～8，對(duì)酸敏感[1]。副溶血性弧菌主要中毒食品是海產(chǎn)品和鹽腌制品，感染臨床癥狀主要是腹瀉、胃寒發(fā)熱、惡心、嘔吐等典型胃腸炎癥狀。來自國(guó)家食源性疾病監(jiān)測(cè)網(wǎng)的數(shù)據(jù)顯示，自1998年以來，副溶血性弧菌引起的中毒事件已經(jīng)高居我國(guó)微生物中毒的首位[2]。姬華[3]、吳英嬌[4]、田金玲[5]、盧曉鳳等[6]研究了對(duì)蝦、蛤肉和模擬蟹肉中副溶血性弧菌在不同培養(yǎng)條件下的變化規(guī)律，發(fā)現(xiàn)Gompertz方程擬合度較高。Zwietering等[7]經(jīng)過比較發(fā)現(xiàn)修正Gompertz方程更能充分描述細(xì)菌生長(zhǎng)曲線。Miles[8]、Gooch[9]、Soon等[10]研究了不同培養(yǎng)條件下致病性和非致病性副溶血性弧菌在肉湯和牡蠣中的生長(zhǎng)情況，發(fā)現(xiàn)初級(jí)模型都能較好的符合修正Gompertz方程。

副溶血性弧菌的感染大部分是因?yàn)槭秤昧宋词焱富蛘咴谏a(chǎn)加工流通過程中再次污染的海產(chǎn)品所導(dǎo)致的，尤其是蝦[11]。我國(guó)的對(duì)蝦產(chǎn)品生產(chǎn)量居世界首位，也在世界蝦產(chǎn)品交易中占較大份額，在過去，我國(guó)對(duì)蝦產(chǎn)品的出口主要是以冷凍蝦為主，但進(jìn)入21世紀(jì)后，蝦制品出口量有了較大幅度增長(zhǎng)，并很快占據(jù)了主要地位[12]。在2005年初美國(guó)對(duì)我國(guó)初級(jí)冷凍蝦類產(chǎn)品開始征收高額反傾銷稅后，我國(guó)冷凍蝦出口嚴(yán)重受阻，也直接導(dǎo)致了我國(guó)蝦類產(chǎn)品出口的降低，面包蝦作為一種鮮蝦的深加工產(chǎn)品，占到了我國(guó)蝦類產(chǎn)品出口增長(zhǎng)的一大半，帶動(dòng)了我國(guó)整個(gè)蝦類產(chǎn)品的出口。面包蝦在生產(chǎn)過程中除去低溫幾乎沒有任何的抑菌工序，因此，本文通過研究不同pH影響下的面包蝦中副溶血性弧菌的生長(zhǎng)動(dòng)態(tài)，確定最適初級(jí)模型，建立生長(zhǎng)動(dòng)力學(xué)參數(shù)最大比生長(zhǎng)速率和遲滯期與溫度之間的二級(jí)模型，探討面包蝦中副溶血性弧菌生長(zhǎng)規(guī)律與pH環(huán)境之間的關(guān)系，對(duì)完善面包蝦的加工技術(shù)具有理論指導(dǎo)意義。

1 材料與方法

1.1 材料與儀器

蝴蝶形面包蝦 湛江國(guó)聯(lián)水產(chǎn)股份有限公司提供；副溶血性弧菌（ATCC17802） 廣東省微生物研究所，菌株在含3%NaCl和20%甘油的腦心浸液肉湯（BHI）中于-80℃保存，菌株接種于2.5%NaCl胰蛋白胨大豆肉湯（TSB）中，37℃培養(yǎng)箱中培養(yǎng)16～20h，活化3代備用；腦心浸液肉湯（BHI）、TCBS瓊脂、3%氯化鈉胰蛋白胨大豆瓊脂（2.5%NaCl TSA）、胰蛋白胨大豆肉湯（TSB）、3%氯化鈉堿性蛋白胨水、營(yíng)養(yǎng)瓊脂（NA） 均購(gòu)自青島海博生物技術(shù)有限公司；NaCl、Na2HPO4·2H2O 、NaH2PO4·2H2O 分 析 純 ，廣 東 省 光華科技股份有限公司。

LS-B50L型立式壓力蒸汽滅菌器 上海華線醫(yī)用核子儀器有限公司；XW-80A型旋渦混合器 海門其林貝爾儀器制造有限公司；拍擊式均質(zhì)器 西班牙IUL公司；SPX-250B-Z型生化培養(yǎng)箱（5～50℃，± 1.5℃） 上海博訊實(shí)業(yè)有限公司醫(yī)療設(shè)備廠；Sartorius PB-10型pH計(jì) 德國(guó)賽多利斯公司；電子天平 日本島津公司；1mL、200μL移液槍 德國(guó)Eppendorf公司；BCD-277型冰箱 博西華家用電器有限公司。

1.2 菌懸液的制備

將副溶血性弧菌劃線接種于TCBS平板，挑取典型綠色斗笠狀菌落，接種至含有9mL 2.5%NaCl TSB的試管中，37℃、120r/min培養(yǎng)12h，利用無菌的2.5% NaCl TSB培養(yǎng)基的調(diào)節(jié)菌液OD600為0.5[3]，用無菌3% NaCl堿性蛋白胨水梯度稀釋，選擇10-4稀釋度作為添加的菌懸液。

1.3 實(shí)驗(yàn)用無菌面包蝦的制備[3，13]

將面包蝦分裝于高壓平口袋中，約40g/包，真空包裝，面包蝦在400MPa、28℃下處理20min。超高壓處理后準(zhǔn)確稱取10g面包蝦于無菌均質(zhì)袋中，加入90mL無菌3%氯化鈉堿性蛋白胨水，均質(zhì)機(jī)拍打90s，吸取100μL面包蝦勻漿液，采用傾注平板法倒3%NaCl NA和2.5%NaCl TSA平板，37℃培養(yǎng)48h，不長(zhǎng)菌代表超高壓處理殺死了所有的雜菌。將無菌的面包蝦置于-18℃冰箱中備用。

1.4 實(shí)驗(yàn)方法

1.4.1 副溶血性弧菌的接種、培養(yǎng)及計(jì)數(shù)[3]無菌稱取10g無菌面包蝦置于直徑為90mm無菌培養(yǎng)皿中，用 不同pH的磷 酸緩 沖 液 浸 泡1～2h后 撈 出 瀝 干[5]，每份蝦仁用無菌槍頭均勻注射接種100μL菌懸液，接種后使面包蝦的初始菌數(shù)達(dá)到102～103CFU/g。

接種完畢后，置于30℃生化培養(yǎng)箱中培養(yǎng)，在選定的時(shí)間取出樣品于無菌均質(zhì)袋中，加入90mL無菌3%氯化鈉堿性蛋白胨水，均質(zhì)機(jī)拍打90s，吸取1mL面包蝦勻漿液注入含有9mL無菌3%氯化鈉堿性蛋白胨水試管中，渦旋混勻后進(jìn)行逐級(jí)稀釋，選擇3個(gè)連續(xù)的適宜稀釋度，吸取100μL采用傾注法倒2.5% NaCl TSA平板，37℃培養(yǎng)24h后計(jì)數(shù)[3，14]。

1.4.2 生 長(zhǎng) 模 型 的 建 立 一 級(jí) 模 型 的 擬 合 采 用Origin 8.0 軟 件（美 國(guó) Microcal公 司）和 DMFit軟 件（http：//modelling.combase.cc/DMFit.aspx），分別用修正Gompertz模型、修正Logistic模型、Baranyi模型[15-17]擬合面包蝦在不同溫度下實(shí)際測(cè)得的副溶血性弧菌濃度隨時(shí)間變化的生長(zhǎng)曲線，再做對(duì)比選出最適模型，3種模型的方程分別為：

式（1）、式（2）中，logNt是在時(shí)間t時(shí)菌落計(jì)數(shù)的對(duì)數(shù)值（lgcfu/mL），A是初始菌量，C是初始菌量和最大菌量之間的差值（lgcfu/mL），M是生長(zhǎng)速率達(dá)到最大值時(shí)的時(shí)間，B是在時(shí)間M時(shí)的相對(duì)最大生長(zhǎng)速率（h-1）。式（3）中，y0和ymax分別是細(xì)菌的初始菌量和最大菌量；μm是最大比生長(zhǎng)速率[lg（cfu/mL）/h-1]；tlag即遲滯期（h）。修正Gompertz模型、修正Logistic模型生長(zhǎng)參數(shù)參照姬華[3]來計(jì)算。

Davey[18]將 自 然 對(duì) 數(shù) 方 程 與Arrhenius關(guān) 系 式 聯(lián)用，用來描述微生物生長(zhǎng)速率與環(huán)境因子之間的關(guān)系，其方程為：

式（4）中，k是生長(zhǎng)速率，即GR，是回歸系數(shù)，根據(jù)初級(jí)最適模型求得GR，將GR和對(duì)應(yīng)的pH代入方程式，擬合出，然后代入原方程，求得二級(jí)模型。二級(jí)模型采用Origin 8.0軟件擬合。

1.4.3 模型的數(shù)學(xué)檢驗(yàn) 一級(jí)模型的擬合優(yōu)度采用決定系數(shù)R2評(píng)價(jià)。二級(jí)模型的殘差平方和（RSS）、平方根誤差（RMSE）、偏差因子（Bf）和準(zhǔn)確性因子（Af）來評(píng)定[16]，表達(dá)式如下：

式中，μobs和μpred分別最大比生長(zhǎng)速率的觀測(cè)值和預(yù)測(cè)值，n為實(shí)驗(yàn)次數(shù)。

2 結(jié)果與分析

2.1 副溶血性弧菌一級(jí)生長(zhǎng)模型的建立

修正Gompertz模型、修正Logistic模型對(duì)副溶血性弧菌在pH為5.5、6、6.5、7、7.5、8、8.5的生長(zhǎng)曲線進(jìn)行擬合，由Origin 8.0軟件計(jì)算出擬合方程的各個(gè)參數(shù)，如表1所示。

表1 修正Gompertz模型、修正Logistic模型參數(shù)及統(tǒng)計(jì)分析Table 1 Statistic analysis of growth parameters of modified Gompertz model and modified Logistic model

使用Baranyi模型對(duì)副溶血性弧菌在pH為5.5、6.0、6.5、7.0、7.5、8.0、8.5的生長(zhǎng)曲線進(jìn)行擬合，DMFit軟件同時(shí)還計(jì)算出了模型的生長(zhǎng)速率μm以及擬合曲線的標(biāo)準(zhǔn)差s和決定系數(shù)R2，所得的參數(shù)及統(tǒng)計(jì)分析如表2所示。

表2 Baranyi模型參數(shù)及統(tǒng)計(jì)分析Table 2 Statistic analysis of growth parameters of Baranyi&Robert model

實(shí)驗(yàn)結(jié)果表明，修正Gompertz模型、修正Logistic模型、Baranyi模型擬合，均能得到較高的決定系數(shù)，分別為0.99、0.98、0.99，選用修正Gompertz模型得出的 參 數(shù) 進(jìn) 一 步 建 立 二 級(jí) 模 型 ，以lg（Nt/N0）對(duì) 時(shí) 間t作圖，副溶血性弧菌的一級(jí)模型擬合曲線如圖1所示。pH為5.5～8.5之間修正Gompertz模型的最大比生長(zhǎng)速 率 μm[（lg（cfu/mL）h）分 別 為 0.5164、0.5244、0.5363、0.5510、0.5728、0.6008、0.6182，隨著pH的升高，微生物的最大比生長(zhǎng)速率也隨之升高。

2.2 副溶血性弧菌二級(jí)生長(zhǎng)模型的建立

將不同pH下修正Gompertz模型擬合得到的最大生長(zhǎng)速率GR代入式（4）擬合，即得到最大生長(zhǎng)速率-pH的二級(jí)模型，擬合曲線見圖2。

得 到 的 擬 合 方 程 為 ：lnGR=-0.053-0.081×pH+ 0.010×（pH）2式（9）

從圖2中可以看出，pH5.5～8.5范圍內(nèi)，副溶血性弧菌的比生長(zhǎng)速率隨著pH的升高而增加，擬合方程的決定系數(shù)R2達(dá)到了0.98。

比生長(zhǎng)速率的預(yù)測(cè)值與實(shí)測(cè)值的殘差平方和（RSS）、平 方 根 誤 差（RMSE）、偏 差 因 子（Bf）、準(zhǔn) 確 性因子（Af）見表3。

表3 比生長(zhǎng)速率的預(yù)測(cè)值與實(shí)測(cè)值的統(tǒng)計(jì)學(xué)分析Table 3 Statistic analysis of μmbetween predicted value and observed value

從實(shí)驗(yàn)結(jié)果可以看出，二級(jí)模型具有較高的R2和較低的RSS值，表明Ratkowsky平方根模型對(duì)生長(zhǎng)速率的擬合度較好，通過數(shù)學(xué)參數(shù)檢驗(yàn)，比較μm的預(yù)測(cè)值與實(shí)測(cè)值，可以計(jì)算出偏差因子和準(zhǔn)確因子，從上表可以看出，RMSE值較小，預(yù)測(cè)值和實(shí)測(cè)值的偏差很小，偏差因子Bf接近于1；準(zhǔn)確因子Af也接近于1，表明模型的擬合效果好，準(zhǔn)確度高，能很好地預(yù)測(cè)面包蝦中副溶血性弧菌的生長(zhǎng)。

3 結(jié)論與討論

本實(shí)驗(yàn)研究了面包蝦中副溶血性弧菌在不同pH（5.5～8.5）影響下的生長(zhǎng)動(dòng)態(tài)，選擇擬合度最好的修正Gompertz方程（R2=0.99）來建立初級(jí)模型，求得副溶血性弧菌的生長(zhǎng)動(dòng)力學(xué)參數(shù)最大比生長(zhǎng)速率。采用自然對(duì)數(shù)方程與Arrhenius聯(lián)用方程來擬合最大比生長(zhǎng)速率與環(huán)境因子pH之間的規(guī)律，擬合度高（R2= 0.99）。經(jīng)過數(shù)學(xué)驗(yàn)證，模型具有較低的RSS值和RMSE值，偏差因子和準(zhǔn)確因子分別為0.99和1.0，均在可接受的范圍內(nèi)，表明模型能較好的描述不同pH影響下面包蝦中副溶血性弧菌的生長(zhǎng)規(guī)律。

王璐華[19]采用濁度法建立了3%氯化鈉堿性蛋白水中中副溶血性弧菌在不同pH影響下的預(yù)測(cè)模型，結(jié)果表明pH在5.5～10.0之間時(shí)，副溶血性弧菌的最大比生長(zhǎng)速率是隨著pH的升高而先增大后減小，在pH為8.5時(shí)達(dá)到最大。李濤[20]研究了TSB中副溶血性弧菌在不同pH影響下的生長(zhǎng)情況，在pH為8時(shí)達(dá)到最大。這些都與本文結(jié)果類似，但考慮到在調(diào)整pH的方法和所選擇的培養(yǎng)基上存在的區(qū)別以及菌株之間本身存在的特性差異，在具體的生長(zhǎng)動(dòng)力學(xué)參數(shù)上可能存在差距。

面包蝦作為一種鮮蝦的深加工產(chǎn)品，其也占據(jù)了蝦制品出口中一定的份額。然而在實(shí)際的生產(chǎn)過程中，除了依靠低溫抑制微生物的生長(zhǎng)，通過調(diào)整pH、鹽度等能影響微生物生長(zhǎng)的因素的方法也是完善面包蝦加工技術(shù)的又一途徑。副溶血性弧菌作為一種重要的食源性致病菌，在蝦體本身尤其是夏季更是檢出率頗高，因此，在面包蝦的加工環(huán)節(jié)，通過有效的手段降低其繁殖速率和數(shù)量，也是保證面包蝦安全生產(chǎn)的必要措施。

圖1 pH為5.5、6、6.5、7、7.5、8、8.5條件下面包蝦中副溶血性弧菌的一級(jí)生長(zhǎng)模型Fig.1 The primary growth model for Vibrio parahaemolyticus in Bread Shrimp at 5.5，6.0，6.5，7.0，7.5，8.0，8.5 of pH

圖2 lnGR和pH二級(jí)模型擬合曲線圖Fig.2 Fitted curve of growth rate and pH of the secondary model

[1]俞鶯，寧喜斌.對(duì)蝦副溶血性弧菌的風(fēng)險(xiǎn)評(píng)估[J].現(xiàn)代食品科技，2006，22（3）：184-196.

[2]劉秀梅.食源性疾病監(jiān)控技術(shù)的研究[J].中國(guó)食品衛(wèi)生雜志，2004，16（1）：3-9.

[3]姬華.對(duì)蝦中食源性弧菌預(yù)測(cè)模型建立及風(fēng)險(xiǎn)評(píng)估[D].無錫：江南大學(xué)，2012.

[4]吳英嬌.雜色蛤中副溶血性弧菌的定量風(fēng)險(xiǎn)評(píng)估[D].福州：福建農(nóng)林大學(xué)，2011.

[5]盧曉鳳.蛤肉中副溶血性弧菌風(fēng)險(xiǎn)評(píng)估體系初探[D].泰安：山東農(nóng)業(yè)大學(xué)，2007.

[6]田金玲.模擬蟹肉中副溶血性弧菌生長(zhǎng)預(yù)測(cè)模型的建立及應(yīng)用[D].湛江：廣東海洋大學(xué)，2010.

[7]Zwietering M H，Jongenburger L，Rombouts F M，et al. Modeling of the bacterial growth curve[J].Appl Environ Microbiol，1990，56（6）：1875-1881.

[8]Miles D W，Ross T，Olley J，et al.Development and evaluation of a predictive model for the effect of temperature and water activity on the growth rate of Vibrio parahaemolyticus[J].Int J Food Microbiol，1997，16，38（2-3）：133-142.

[9]Gooch J A，DePaola A，Bowers J，et al.Growth and survival of Vibrio parahaemolyticus in postharvest American oysters[J].J Food Prot，2002，65（6）：970-974.

[10]Soon K S，Min K J，Jung Y J，et al.A model of the effect of temperature on the growth of pathogenic and nonpathogenic Vibrio parahaemolyticus isolated from oysters in Korea[J].FoodMicrobiol，2008，25（5）：635-641.

[11]Alam M J，Miyoshi S I，&Shinoda S.Studies on pathogenic Vibro parahaemolyticus during a warm weather season in the Seto Inland Sea，Japan[J].Environmental Microbiology，2003，5（8）：706-710.

[12]陳文漢，寧凌.中國(guó)對(duì)蝦產(chǎn)業(yè)現(xiàn)狀及相關(guān)研究綜述[J].漁業(yè)經(jīng)濟(jì)研究，2008（2）：21-27.

[13]殷允旭.蝦仁的超高壓殺菌工藝研究[D].安徽：合肥工業(yè)大學(xué)，2007.

[14]中華人民共和國(guó)國(guó)家標(biāo)準(zhǔn)，食品微生物學(xué)檢驗(yàn)副溶血性弧菌檢驗(yàn)[S].GB/T 4789.7-2008.

[15]Gibson A M，Bratchell N，Roberts T A.The effect of sodium chloride and temperature on the rate and extent of growth of Clostridium botulinum type A in pasteurized pork slurry[J]. Application of Bacteriol，1987，62（6）：479-490.

[16]Gibson A M，Bratchell N，Roberts T A.Predicting microbial growth：growth reponses of salmonellae in a laboratory medium as affected by pH，sodium chloride and storage temperature.[J]. International Journal of Food Microbiology，1988 ，6（2）：155-178.

[17]Baranyi J，Roberts T A.Mathematics of predictive food microbiology[J].Int J Food Microbiol，1995，26（2）：199-218.

[18]Davey K R.A terminology for models in predictive microbiology[J].Food Microbiology，1992，9（4）：353-354.

[19]王璐華.水產(chǎn)品副溶血性弧菌預(yù)警模型的研究[D].上海：上海海洋大學(xué)，2008.

[20]李濤.副溶血性弧菌生長(zhǎng)特性及雙因素預(yù)警模型的研究[D].上海：上海海洋大學(xué)，2010.

Establishment of growth predictive model of Vibrio parahaemolyticus in bread shrimp under different pH condition

HE Qun1，HUANG He1，*，LIU Ya1，JIANG Zhi-hong1，LI Jian-rong2
（1.Guangdong Provincial Key Laboratory of Aquatic Product Processing and Safety，Key Laboratory of Advanced Processing of Aquatic Products of Guangdong Higher Education Institution，College of Food Science and Technology，Guangdong Ocean University，Zhanjiang 524088，China；2.Food Safety Key laboratory of Liaoning Province，Bohai University，Jinzhou 121013，China）

To study the impact of pH to the growth law of Vibrio parahaemolyticus in the process of bread shimp，growth curves under different pH （ 5.5 ， 6.0 ， 6.5 ， 7.0 ， 7.5 ， 8.0 and 8.5 ， respectively ） were determined.The modified Compertz，modified Logistic and Baranyi models，respectively，were applied to fit the growth curves by Origin 8.0 software and DMFit software，the maximum specific growth rates was calculated by the primary model with the highest coefficient of determination.The modified Arrhenius model was used to fit the curves between the maximum specific growth rates and pH for the establishment of the secongary model.Results showed that the modified Compertz model（R2=0.99） was more accurate and useful in fitting V.parahaemolyticus growth curves than the modified Logistic and Baranyi&Robert models，respectively.And the secondary model had a high goodness-of-fit（R2=0.98 ）.The values of bias factor and accuracy factor were also acceptable statistically that were close to 1.0.The predictive models were highly reliable.These could effectively predict the growth of Vibrio parahaemolyticus under above pH condition.

bread shrimp；Vibrio parahaemolyticus；growth model

TS201.3

A

1002-0306（2014）20-0107-05

10.13386/j.issn1002-0306.2014.20.014

2014-03-17

賀群（1987-），女，在讀碩士研究生，研究方向：食品質(zhì)量與安全。

* 通訊作者：黃和（1962-），男，教授，研究方向：食品質(zhì)量與安全。

“十二五”國(guó)家科技支撐計(jì)劃項(xiàng)目（2012BAD29B06）；廣東省省部產(chǎn)學(xué)研結(jié)合項(xiàng)目（2011B090400154）；湛江市科技計(jì)劃項(xiàng)目（湛科[2010]174號(hào)）。