能源甜菜酒精發(fā)酵的工藝優(yōu)化及數(shù)學(xué)模型.pdf

第 25 卷 第 6 期農(nóng) 業(yè) 工 程 學(xué) 報Vol.25No.6 2142009 年6 月Transactions of the CSAEJun. 2009 能源甜菜酒精發(fā)酵的工藝優(yōu)化及數(shù)學(xué)模型 史淑芝 1,2，程大友2，徐德昌2，馬鳳鳴1，魯兆新2， 董愛軍 2，代翠紅2，羅成飛2 （1東北農(nóng)業(yè)大學(xué)農(nóng)業(yè)部寒地作物生理生態(tài)重點開放試驗室，哈爾濱 150056； 2哈爾濱工業(yè)大學(xué)食品科學(xué)與工程學(xué)院，哈爾濱 150090） 摘要：利用菌種 ADY（Saccharomyces cerevisiae）對能源甜菜 NY0503 進行酒精發(fā)酵。應(yīng)用 Plackett- Burman 設(shè)計法從 底物濃度、料液比、加菌量、營養(yǎng)鹽添加量、加磷量、pH 值、轉(zhuǎn)速、發(fā)酵溫度和發(fā)酵時間 9 個因素中篩選出了加磷量、 發(fā)酵溫度和底物濃度為主要影響因素，并用響應(yīng)面分析法求出回歸方程。結(jié)果表明，經(jīng)優(yōu)化后最佳工藝條件為：底物濃 度 12%、料液比（質(zhì)量比）11、加菌量 15%、營養(yǎng)鹽添加量 0.5、加磷量 1.9%、pH 5.0、轉(zhuǎn)速 130 r/min、發(fā)酵溫度 31 和發(fā)酵時間 44 h，共進行 5 批次驗證試驗，結(jié)果為酒精轉(zhuǎn)化率的平均值為 95.48%，與模型的理論值 96.54%的差值僅占 理論值的 1.09%，進一步說明所建立的模型是切實可行的。 關(guān)鍵詞：酒精發(fā)酵，工藝優(yōu)化，數(shù)學(xué)模型，能源甜菜 doi：10.3969/j.issn.1002-6819.2009.06.040 中圖分類號：TS261.4文獻標識碼：A文章編號：1002-6819(2009)-6-0214-05 史淑芝，程大友，徐德昌，等.能源甜菜酒精發(fā)酵的工藝優(yōu)化及數(shù)學(xué)模型J. 農(nóng)業(yè)工程學(xué)報，2009，25(6)：214218. Shi Shuzhi, Cheng Dayou, Xu Dechang, et al. Technological optimization and mathematical model of alcohol fermentation of energy beetJ. Transactions of the CSAE, 2009,25(6)：214218.(in Chinese with English abstract) 0引言 隨著全世界能源的不斷擴展，可再生的生物質(zhì)能源 將是本世紀主要的能源物質(zhì)1- 4。從 20 世紀 70、80 年代 開始，世界上許多國家就開始尋求利用可再生的生物質(zhì) 進行能源開發(fā)5。而隨著我國經(jīng)濟的飛速發(fā)展，能源的消 耗越來越大，石油資源的短缺，長期供應(yīng)能力不足的問 題日顯突出6。 酒精是一種可再生能源,也是石化燃料最有希望的替 代產(chǎn)品7,8。這在世界上許多國家的人們都已達成共識。 但怎樣才能找到一種適宜的能夠產(chǎn)生酒精的可再生的生 物質(zhì)，又擺在人們面前。 甜菜作為糖料作物早已為人類所利用， 從 20 世紀 90 年代末開始，由于受世界食糖市場不穩(wěn)定和石油能源緊 張的影響，少數(shù)發(fā)達國家便開始嘗試選育能源甜菜的研 究9。由于能源甜菜是新興的可再生能源作物10，它具有 生物學(xué)產(chǎn)量高，含糖量適中等特點，是環(huán)保能源，又是 糖飼兼收的作物11。同時，利用能源甜菜發(fā)酵生產(chǎn)燃料 酒精還具有良好的環(huán)境效益和經(jīng)濟效益。 利用能源甜菜進行酒精發(fā)酵，主要是對篩選所得到 的其組分適于酒精發(fā)酵的甜菜塊根進行研究。本研究針 收稿日期：2008- 10- 28修訂日期：2009- 12- 30 基金項目： 農(nóng)業(yè)部寒地作物生理生態(tài)重點開放試驗室開放基金： 國家科技部 專項基金（NCSTE- 2006- JKZX- 022） 作者簡介：史淑芝（1967） ，女，博士研究生，副研究員，研究方向：植 物生理及農(nóng)產(chǎn)品深加工。 哈爾濱市南崗區(qū)海河路 202 號哈爾濱工業(yè)大學(xué)食 品學(xué)院 2649 信箱，150090。Email: 通訊作者：馬鳳鳴（1947） ，男，教授，博士生導(dǎo)師，從事作物生理研 究。哈爾濱市香坊區(qū)木材街 59 號東北農(nóng)業(yè)大學(xué)農(nóng)學(xué)院，150030。 Email: Fengming_ 對酒精發(fā)酵過程中存在的共性技術(shù)難題開展研究，研究 利用能源甜菜塊根發(fā)酵生成酒精的最佳工藝參數(shù)，建立 適宜的數(shù)學(xué)模型12，最終獲得較高的酒精轉(zhuǎn)化率，從而 為工業(yè)化生產(chǎn)和提高經(jīng)濟效益提供理論依據(jù)。 1材料與方法 1.1試驗材料 材料：能源甜菜 NY0503。 供試菌種：ADY（Saccharomyces cerevisiae） 。 試劑：AlCl36H2O，(NH4)2SO4，Mg SO4，KH2PO4， CaCl2，過磷酸鈣，鹽酸，蒸餾水。 儀器及用品：500 mL 容量瓶，100 mL 量筒，自制鋸 糊機， VENEMA 甜菜檢糖自動線， 天平， pH 計， DH5000 型恒溫培養(yǎng)箱，HZQ- F160A 型恒溫振蕩培養(yǎng)箱，EYELA N- 1000 旋轉(zhuǎn)蒸發(fā)儀，酒精計。 1.2試驗方法 1.2.1發(fā)酵底物制備 每個處理設(shè) 3 個平行試驗，按以下順序及方法進行。 1）制糊：取切削標準的 NY0503 甜菜塊根洗凈后用 自制的鋸糊機加工成糊狀，充分混勻。 2）蔗糖測定及配制底物濃度：用荷蘭進口的 VENEMA 甜菜檢糖自動線檢測樣品的蔗糖含量，具體做 法為：先配制 AlCl3溶液，使其含量為 85 g/L，待使用前 將其稀釋 47 倍達 1.81 g/L 后作為浸提液使用；浸糖時用 26 g 甜菜糊加入 177 mL 浸提液，充分攪拌后用甜菜檢糖 自動線檢測樣品的蔗糖含量，然后根據(jù)其蔗糖含量將樣 品分別配制成 12%和 15%的底物濃度。 3） 配制料液比： 取 100 g 甜菜糊放入 500 mL 三角瓶 內(nèi)， 用蒸餾水按試驗要求分別配成 10.5 和 11 的料液 第 6 期史淑芝等：能源甜菜酒精發(fā)酵的工藝優(yōu)化及數(shù)學(xué)模型215 比（質(zhì)量比） 。 4）殺菌：在每個三角瓶內(nèi)分別加入藥料比為 0.004% 的五氯苯酚鈉，混勻后將三角瓶放在 50的搖床內(nèi)以 180 r/min 的轉(zhuǎn)速殺菌 30 min13。 5）添加營養(yǎng)鹽：每種營養(yǎng)鹽的添加量見表 1。 表 1營養(yǎng)鹽的添加量 Table 1Addition amount of nutrition saltsgL- 1 水平(NH4)2SO4MgSO4KH2PO4CaCl2 0.50.552.52.5 111055 6）調(diào)整磷含量：將過磷酸鈣溶解，再加入到每個三 角瓶中使其終濃度達到試驗所要求的 1%或 2%水平。 7）調(diào)節(jié) pH 值：用酸度計對每個三角瓶內(nèi)的樣品用 濃鹽酸進行 pH 值調(diào)節(jié)，最終達到試驗方案的要求。 8）加菌：將活化好的菌種按試驗方案要求分別以 10%或 15%的菌量加到每個三角瓶中，然后混勻。 9）稱重：最后用天平將每個培養(yǎng)瓶進行稱量，并記 錄數(shù)據(jù)，以便進行發(fā)酵過程控制的 CO2失重測定。 1.2.2酒精發(fā)酵 首先，將已制備好的三角瓶放在 28的搖床培養(yǎng)箱 中，以 150 r/min 的轉(zhuǎn)速培養(yǎng) 4 h，并將瓶塞旋松，使部分 氧氣進入，以使酵母菌在有氧條件下快速繁殖。 然后，再將培養(yǎng)箱溫度分別調(diào)至 31和 35，轉(zhuǎn)速 分別調(diào)至 100 r/min 和 130 r/min，使培養(yǎng)瓶在近乎厭氧條 件下進行酒精發(fā)酵達到試驗所要求時間， 并且每隔 6 h 測 定一次 CO2失重量。 最后，以 CO2失重小于 0.1 g/h 時視為發(fā)酵終止，將 三角瓶取出，進行蒸餾14。 1.2.3酒精蒸餾及測定 將成熟的發(fā)酵醪全部倒入 500 mL 的圓底燒瓶中， 并 將旋轉(zhuǎn)蒸發(fā)儀的水浴鍋溫度調(diào)至 52，緩慢蒸發(fā)，以使 發(fā)酵所產(chǎn)生的酒精盡可能全部蒸出；最后將收集瓶中的 餾出液定容至 100 mL，混勻后倒入 100 mL 量筒中,用酒 精計 （標溫 20） 測酒精度（下緣為準） ，同時測定溫度， 然后換算成 20時的酒精度。每個處理做 3 次平行試 驗15。 1.2.4酒精轉(zhuǎn)化率計算方法 = mLg/mL g%53.8 100% 酒精轉(zhuǎn)化率 發(fā)酵醪體積() 發(fā)酵醪酒精含量(體積百分數(shù)) 相應(yīng)密度() 發(fā)酵底物質(zhì)量( ) 含糖率( ) 式中53.8100 g 蔗糖理論上轉(zhuǎn)化為酒精的質(zhì)量，g。 1.3試驗設(shè)計 應(yīng)用 SAS（Statistical Analysis System）V8.2 統(tǒng)計分 析軟件進行試驗設(shè)計和分析。 1.3.1Plackett- Burman（P- B）設(shè)計 Plackett- Burman（P- B）設(shè)計法是一種近飽和的兩水 平試驗設(shè)計方法。它能用最少試驗次數(shù)估計出因素的主 效應(yīng)，以從眾多的考察因素中快速有效地篩選出最為重 要的幾個因素，供進一步研究16。 應(yīng)用 Plackett- Burman（P- B）試驗設(shè)計法篩選重要因 素，并根據(jù)有關(guān)參考文獻和單因素試驗的結(jié)果9，設(shè)計 Plackett- Burman（P- B）試驗因素及水平，如表 2 所示。 表 2Plackett- Burman（P- B）設(shè)計的試驗因素與水平 Table 2Experimental factors and levels of Plackett- Burman(P- B) design 水平 因素 - 11 底物濃度/%1215 料液比1:0.51:1 加菌量/%1015 營養(yǎng)鹽0.51 加磷量/%12 pH 值4.55.5 轉(zhuǎn)速/rmin- 1100130 發(fā)酵溫度/3135 發(fā)酵時間/ h4460 注：1）營養(yǎng)鹽的添加為水平值，無單位；2）發(fā)酵時間在 4460 h 之間， CO2失重均小于 0.1 g/h，酒精轉(zhuǎn)化率稍有降低，但到 60 h 以后酒精轉(zhuǎn)化率 急劇下降，故設(shè)這兩個水平。 1.3.2Box- Behnken 設(shè)計及響應(yīng)面分析 根據(jù) Plackett- Burman（P- B）試驗設(shè)計法所篩選出的 3 個重要因素，利用 SAS V8.2 統(tǒng)計分析軟件進行試驗 Box- Behnken 設(shè)計，并進行響應(yīng)面分析，從而確定 3 個重 要因素的最佳水平。 2結(jié)果與分析 2.1根據(jù) PlackettBurman 設(shè)計法篩選重要因素 本文根據(jù)以前的研究和文獻報道17，選用試驗次數(shù) N=12的試驗設(shè)計，對底物濃度（X1）、料液比（X2）、 加菌量（X3）、營養(yǎng)鹽（X4）、加磷量（X5）、pH（X6）、 轉(zhuǎn)速（X7）、發(fā)酵溫度（X8）和發(fā)酵時間（X9）9個因素 進行考察，響應(yīng)值為酒精轉(zhuǎn)化率（Y）。發(fā)酵試驗設(shè)計及 試驗結(jié)果見表3。 表 3Plackett- Burman 試驗設(shè)計表及響應(yīng)值 Table 3Plackett- Burman design table and response values 序號X1X2X3X4X5X6X7X8X9Y/% 11- 11- 1- 1- 111176.15 211- 11- 1- 1- 11174.79 3- 111- 11- 1- 1- 1197.43 41- 111- 11- 1- 1- 182.28 511- 111- 1- 1- 1- 197.44 6111- 111- 11- 191.33 7- 1111- 111- 1190.15 8- 1- 1111- 111- 192.95 9- 1- 1- 1111- 11187.46 101- 1- 1- 1111- 1187.36 11- 11- 1- 1- 1111- 187.02 12- 1- 1- 1- 1- 1- 1- 1- 1- 189.10 注：底物濃度（X1）、料液比（X2）、加菌量（X3）、營養(yǎng)鹽（X4）、加磷 量（X5）、pH（X6）、轉(zhuǎn)速（X7）、發(fā)酵溫度（X8）和發(fā)酵時間（X9）9 個 因素進行考察，響應(yīng)值為酒精轉(zhuǎn)化率（Y）。 采用 SAS V8.2 進行各因素主效應(yīng)分析結(jié)果見表 4。 可見針對 ADY 菌種， 發(fā)酵起重要作用的 3 個因素依 216農(nóng)業(yè)工程學(xué)報2009 年 次為加磷量，底物濃度和發(fā)酵溫度，這 3 個因素均達到 了 95%的顯著水平；而加磷量的效應(yīng)值為正值，底物濃 度和發(fā)酵溫度的效應(yīng)值為負值，這在下一步水平設(shè)計時 需要充分考慮。 進而從表 4 的估計效應(yīng)和 t 值兩列中也可 以發(fā)現(xiàn)：發(fā)酵過程中加入磷量的多少對試驗結(jié)果影響極 大，另外發(fā)酵溫度和底物濃度對最終轉(zhuǎn)化率影響也較大。 表 4統(tǒng)計分析結(jié)果 Table 4Results of statistic analysis 因素 符號 因 素 估計 效應(yīng) 標準差t 值P 值排序 X1底物濃度- 5.79330.1812- 4.90480.03912 X2料液比3.810.18123.22560.08425 X3加菌量1.18670.18121.00470.42097 X4營養(yǎng)鹽- 0.55330.1812- 0.46850.68558 X5加磷量9.080.18127.68740.01651 X6pH 值- 0.37670.1812- 0.31890.78009 X7轉(zhuǎn)速1.44670.18121.22480.34536 X8發(fā)酵溫度- 5.67670.1812- 4.80600.04073 X9發(fā)酵時間- 4.46330.1812- 3.7780.06344 注：表中標準差是針對整個實驗的處理誤差，故為同一值。 2.2應(yīng)用響應(yīng)面分析法確定重要因素的最佳水平 2.2.1試驗因素水平設(shè)計 通過上述 Plackett- Burmen 設(shè)計得到的結(jié)果可知，各 因素水平可按表 4 中的效應(yīng)值選取適當?shù)乃?，并分?出：針對 ADY 菌種，加磷量，發(fā)酵溫度和底物濃度為影 響酒精發(fā)酵結(jié)果的 3 個主要因素，根據(jù)其效應(yīng)值的大小 和正負按排響應(yīng)面試驗 3 個因素和水平如表 5。 表 5響應(yīng)面分析試驗因素與水平 Table 5Experimental factors and levels of Response Surface Methodology(RSM) 水平 編號因素 - 101 X5加磷量/%11.82.5 X8發(fā)酵溫度/283133 X1底物濃度/%81215 2.2.2響應(yīng)面分析試驗設(shè)計及結(jié)果 由表6可見：本試驗共通過15次試驗來完成，其中12 次試驗為析因點分析試驗，并設(shè)置三次中心點重復(fù)試驗， 用來估計試驗誤差；以酒精的轉(zhuǎn)化率為響應(yīng)值，根據(jù)表6 中Box- Behnken設(shè)計的試驗結(jié)果， 運用SAS軟件對表6的結(jié) 果進行二次回歸分析，得到二次回歸方程Y= 96.14333 + 2.8075X5- 1.4375X8- 2.48X1- 11.96667X5X5- 0.565X5X8+ 0.29X5X1-5.72666X8X8-1.905X8X1- 9.481667X1X1，回歸 方程方差分析及模擬可信度分析結(jié)果見表7和表8。 從方差分析表7可以看出模型在 =0.01 水平上回歸 顯著，一次項、二次項對響應(yīng)值有顯著性影響，但交互 項影響不顯著，說明所選因素之間的交互作用不明顯， 采用該模型時可以不考慮因素間的交互作用；失擬項反 映的是試驗數(shù)據(jù)與模型不相符的情況， P=0.2805720.1， 說明失擬不顯著，因此證明模型選擇正確。 從表8可見：決定系數(shù)R2為98.64%，說明此模型可以 解釋98.64%試驗所得的酒精轉(zhuǎn)化率的變化，表明方程擬 合程度較好。而Y的變異系數(shù)CV則表示試驗的精確度， CV值越高，試驗的可靠性越低，本試驗中CV=2.043647， 數(shù)據(jù)較低，說明此試驗操作可信。綜上說明該回歸方程 為能源甜菜固態(tài)發(fā)酵產(chǎn)生酒精預(yù)測提供了一個合理的模 型。 表 6響應(yīng)面分析試驗設(shè)計及發(fā)酵結(jié)果 Table 6Design of RSM and fermentation results 序號X5X8X1響應(yīng)值 Y/% 1- 1- 1077.23 2- 11073.45 31- 1084.58 411078.54 50- 1- 182.45 60- 1180.26 701- 185.42 801175.61 9- 10- 174.44 1010- 178.87 11- 10169.94 1210175.53 1300096.52 1400097.08 1500094.83 表 7回歸方程的方差分析 Table 7Analysis of variance for regression equation 變異來源自由度平方和方差F值P值 X5163.0564563.0564522.646710.005065 X8116.5312516.531255.9371950.058897 X1149.203249.203217.671320.008459 X5X51528.7426528.7426189.89780.0001 X5X811.27691.27690.4585980.52834 X5X110.33640.33640.1208180.742303 X8X81121.0882121.088243.488790.001205 X8X1114.516114.51615.2134550.071244 X1X11331.9459331.9459119.21830.000112 模型91011.583112.398240.367780.000387 一次項3128.790942.930315.418410.00584 二次項3866.6631288.8877103.7540.0001 交互項316.12945.3764671.9309570.242783 誤差項513.921772.784353 失擬項311.17773.72592.7156050.280572 純誤差22.7440671.372033 所有項141025.505 表 8模型可信度分析 Table 8Fit statistics for the model 平均值決定系數(shù)R2校正后的R2模型誤差的平方根 Y的變異系數(shù)CV 81.6598.64%96.20%1.6686382.043647 注：表中的81.65為所有響應(yīng)值的平均值。 2.2.3響應(yīng)面分析及試驗因素的最佳水平確定 通過回歸方程繪制分析圖，研究所擬合的相應(yīng)曲面 的形狀， 響應(yīng)面立體分析圖和相應(yīng)等高線圖， 分別見圖1、 圖2和圖3。從圖1、2、3中及軟件分析可見，回歸方程存 第 6 期史淑芝等：能源甜菜酒精發(fā)酵的工藝優(yōu)化及數(shù)學(xué)模型217 在穩(wěn)定點，同時穩(wěn)定點為極大值。從圖1、圖2和圖3的等 高線圖中可以看出：其等高線圖近乎為圓形，說明任意 兩因素間交互作用不顯著；而從線條的疏密程度可以得 出：X5（加磷量）比X8（發(fā)酵溫度）和X1（底物濃度）對 響應(yīng)值所起的作用要大， 而X8和X1對響應(yīng)值所起的作用差 異不明顯。 圖1Y=f(X5, X8)響應(yīng)面立體分析圖和等高線圖 Fig.1Response surface and contours of adding phosphor and fermentation temperature 圖2Y=f(X5, X1)響應(yīng)面立體分析圖和等高線圖 Fig.2Response surface and contours of adding phosphor and substrate concentration 圖3Y=f(X8, X1)響應(yīng)面立體分析圖和等高線圖 Fig.3Response surface and contours of fermentation temperature and substrate concentration 通過嶺脊分析得到，Y 值的預(yù)測極大值為 96.54%， 而該預(yù)測值的標準誤為 0.9516，表明該試驗達到了顯著 水平；極大值所對應(yīng)的 3 個主要因素（X5，X8，X1）的編 碼值及對應(yīng)量見表 9。 表 9嶺脊分析因素與水平 Table 9Factors and levels for ridge analysis 因素編碼值Y/% X50.118521.88 X8- 0.1117730.66 X1- 0.1177411.53 2.2.4驗證試驗 根據(jù)實際情況分別將加磷量 （X5） ，發(fā)酵溫度 （X8） ， 底 物濃度（X1）的水平定為加磷量 1.9%,發(fā)酵溫度 31，底 物濃度 12%。在以上優(yōu)化條件下進行驗證試驗，共進行 5 批次 500 mL 搖瓶試驗，酒精轉(zhuǎn)化率平均結(jié)果為 95.48%， 該試驗值（95.48%）與模型的理論值（96.54%）的差值 僅占理論值的 1.09%， 可見該模型可以較好地反映出利用 ADY 菌種將能源甜菜 NY0503 塊根轉(zhuǎn)化成酒精的條件。 3結(jié)論 1）從底物濃度、料液比、加菌量、營養(yǎng)鹽添加量、 加磷量、pH、轉(zhuǎn)速、發(fā)酵溫度和發(fā)酵時間 9 個與發(fā)酵相 關(guān)的因素中篩選出加磷量、發(fā)酵溫度和底物濃度 3 個主 要影響因素。 2）利用響應(yīng)面分析法優(yōu)化后最佳工藝條件為：底物 218農(nóng)業(yè)工程學(xué)報2009 年 濃度 12%、料液比 11、加菌量 15%、營養(yǎng)鹽 0.5、加磷 量 1.9%、pH 5.0、轉(zhuǎn)速 130 r/min、發(fā)酵溫度 31和發(fā)酵 時間 44 h。 3）優(yōu)化條件下驗證試驗結(jié)果為：酒精轉(zhuǎn)化率的平均 值為 95.48%， 與模型的理論值 96.54%的差值僅占理論值 的 1.09%，進一步說明所建立的模型是切實可行的。 參考文獻 1路 明開發(fā)生物質(zhì)能源發(fā)展能源農(nóng)業(yè)N人民日報， 2005- 06- 20 2Berndes G, Hoogwijk M, Broek R. The contribution of biomass in the future global energy supply: a review of 17 studiesJ. Biomass and Bioenergy，2003，25(1)：128 3王海，盧旭東，張慧媛國內(nèi)外生物質(zhì)的開發(fā)與利用 J農(nóng)業(yè)工程學(xué)報，2006，22(增 1)：811 Wang Hai, Lu Xudong, Zhang Huiyuan. Exploitation and utilization of biomass at home and abroadJ. Transactions of the CSAE，2006，22(Supp 1)：811(in Chinese with English abstract) 4周中仁，吳文良生物質(zhì)能研究現(xiàn)狀及展望J農(nóng)業(yè)工程 學(xué)報，2005，21(12)：1215 Zhou Zhongren, Wu Wenliang. Status quo and prospects of biomass enernyJ. Transactions of the CSAE， 2005， 21(12)： 1215(in Chinese with English abstract) 5宋永芳 能源植物的開發(fā)與利用進展J 生物質(zhì)化學(xué)工程， 2006，40(6)：5153 Song Yongfang. Progress of exploitation and utilization on energy plantJ. Biomass Chemical Engineering，2006， 40(6)：5153(in Chinese with English abstract) 6孫永明，袁振宏，孫振鈞中國生物質(zhì)能源與生質(zhì)能利用 現(xiàn)狀與展望J可再生能源，2006，(2)：7882 Sun Yongming, Yuan Zhenhong, Sun Zhenjun. Status and future of bioenergy and biomass utilizationJ. Renewable Energy， 2006， (2)： 7882 (in Chinese with English abstract) 7Yan Lin, Shuzo Tanaka. Ethanol fermentation from biomass resources: current state and prospectsJ. Appl Microbiol Biotechnol，2006，69：627642 8孫振鈞中國生物質(zhì)產(chǎn)業(yè)及發(fā)展取向J農(nóng)業(yè)工程學(xué)報， 2004，20(5)：15 Sun Zhenjun .Biomass industry and its developmental trends in ChinaJ. Transactions of the CSAE， 2004， 20(5)： 15 (in Chinese with English abstract) 9BaiFW,AndersonWA,Moo- YoungM.Ethanol fermentation technologies from sugar and starch feedstocks J. BioechnologyAdvances，2008，26：89105 10 Hussy I, Hawkes F R, Dinsdale R, et al. Continuous fermentativehydrogenproductionfromsucroseand sugarbeetJ. Hydrogen Energy，2005，30：471483 11 史淑芝，程大友，馬鳳鳴，等生物質(zhì)能源作物能源甜 菜的開發(fā)利用J中國農(nóng)學(xué)通報，2007，(11)：416419 Shi Shuzhi, Cheng Dayou, Ma Fengming, et al. Exploitation and utilization of biomass energy crop- energy beetJ. Chinese Agricultural Science Bulletin，2007，(11)：416 419(in Chinese with English abstract) 12 Antonio De Le n- Rodr guez, Pilar Escalante- Minakata, Ana P. Barba de la Rosa, et al. Optimization of fermentation conditions for the production of the mescal from Agave salmiana using response surface methodologyJ. Chemical Engineering and Processing，2008，47：7682 13 章克昌酒精與蒸餾酒工藝學(xué)M北京：中國輕工業(yè)出 版社，2004 14 趙光鰲，金嶺南黃酒生產(chǎn)分析指南M北京：中國輕 工業(yè)出版社，1987：5181 15 郝林食品微生物試驗技術(shù)M北京：中國農(nóng)業(yè)出版 社，2006 16 虞悝響應(yīng)面法與多因子試驗J工業(yè)微生物，1989， 6：2130 17 歐宏宇，賈士儒SAS 軟件在微生物培養(yǎng)條件優(yōu)化中的應(yīng) 用J天津輕工業(yè)學(xué)院學(xué)報，2001，36：1417 Ou Hongyu, Jia Shiru. The application of SAS system in optimization of microbial culture conditionsJ. Journal of Tianjin Institute of Light Industry，2001，36：1417(in Chinese with English abstract) Technological optimization and mathematical model of alcohol fermentation of energy beet Shi Shuzhi1,2, Cheng Dayou2, Xu Dechang2, Ma Fengming1 , Lu Zhaoxin2, Dong Aijun2, Dai Cuihong2, Luo Chengfei2 (1. Key Laboratory of Cold Crop Physiology and Ecology, Min