牛紅志，孔曉英，李連華，孫永明，袁振宏，王瑤，周賢友

（1 中國科學(xué)院廣州能源研究所，中國科學(xué)院可再生能源重點(diǎn)實驗室，廣東 廣州 510640；2 中國科學(xué)院大學(xué)，北京 100049）

引 言

稻殼是稻谷加工過程中的主要副產(chǎn)物之一，是一種量大、價廉的可再生資源。截至2013年，我國稻谷年產(chǎn)量連續(xù)6年突破1.9×108t[1]，而稻殼約占稻谷重量的17%～20%[2-3]，即每年的稻殼產(chǎn)量超過3.2×107t。稻殼中含纖維素28.6%～36.10%、半纖維素16.48%～28.6%、木質(zhì)素14.30%～35.90%、粗蛋白2.5%～3.1%、粗脂肪0.7%～1.0%，二氧化硅等礦物質(zhì)含量占3%～22%[4-7]。稻殼中木質(zhì)纖維素含量高導(dǎo)致不適宜直接作為飼料，而礦物質(zhì)含量高又導(dǎo)致其直燃困難且引起爐內(nèi)結(jié)渣，從而在一定程度上限制了稻殼資源的規(guī)模化利用。

厭氧發(fā)酵制備生物燃?xì)馐寝r(nóng)業(yè)廢棄物規(guī)模化處理和資源化利用的有效途徑[8-12]。有機(jī)廢棄物經(jīng)過厭氧發(fā)酵處理，不僅可以獲得戶用、發(fā)電或者車用氣體燃料[13]，還可制備液態(tài)或固態(tài)有機(jī)肥料[11]或作為飼料添加劑[7]，沼渣也可作為生產(chǎn)二氧化硅等材料的原料[14]。目前稻殼厭氧發(fā)酵研究主要在厭氧發(fā)酵工藝及產(chǎn)氣潛力方面[15-17]，還沒有系統(tǒng)地進(jìn)行物質(zhì)與能量分布的研究。物質(zhì)流分析是在特定的系統(tǒng)范圍內(nèi)對特定的物質(zhì)進(jìn)行工業(yè)代謝研究的有效手段，展示系統(tǒng)內(nèi)各個過程中不同物質(zhì)的流動情況及其相互關(guān)系，評估各個過程對系統(tǒng)產(chǎn)生的影響[18-19]。

本研究以稻殼作為原料，采用批式中溫厭氧發(fā)酵工藝研究稻殼制備生物燃?xì)膺^程中的產(chǎn)氣特性，并結(jié)合物質(zhì)流分析方法研究厭氧發(fā)酵過程中的C、N 元素流動狀況以及物質(zhì)和能量的轉(zhuǎn)化利用效率，為稻殼資源的統(tǒng)籌管理和能源化利用提供參考。

1 實驗材料與方法

1.1 實驗原料

實驗原料為稻殼，取自武漢，經(jīng)粉碎過篩（20目，粒徑＜0.85 mm）備用；接種物為本實驗室長期馴化的中溫厭氧發(fā)酵菌種。

1.2 實驗裝置及操作

實驗裝置如圖1所示，2.5 L 玻璃反應(yīng)器置于35℃±1℃的水浴鍋中，側(cè)壁上有2 個出口，上出口為取氣樣口，下出口為取液樣口，上部連接集氣瓶，集氣瓶后連通集液瓶，通過排飽和鹽水法測定生物燃?xì)猱a(chǎn)量。接種物加入量為1800 ml，按發(fā)酵總TS 濃度的5%加入稻殼（稻殼與接種物的總質(zhì)量比為1:22），并添加2.5%（以稻殼總質(zhì)量計）的NH4HCO3作為緩沖劑和補(bǔ)充氮源。充入高純N2排出反應(yīng)器頂部的空氣。實驗期間每天手動攪拌2 次。當(dāng)日產(chǎn)氣率持續(xù)7 d 低于1 ml·(gVSRH)-1·d-1后停止實驗，共運(yùn)行39 d。發(fā)酵剩余物稱量后過濾，濾渣于105℃烘干作為沼渣（biogas residue）分析，濾液作為沼液（biogas slurry）分析。

圖1 厭氧發(fā)酵實驗裝置Fig.1 Experimental sketch of anaerobic fermentation

1.3 數(shù)據(jù)測試與分析方法

總固體（TS）含量和揮發(fā)性固體（VS）含量分別采用 105℃烘干和 550℃煅燒法測定；熱值（calorific value,CV）由量熱儀C2000（德國IKA@公司）測定；固體中的C、H、N 和S 元素含量采用Vario EL cube 元素分析儀（德國Elementar 公司）測定；纖維素（CL）和半纖維素（HCL）測定參考美國可再生能源實驗室標(biāo)準(zhǔn)測定方法（LAP,NREL）[17,20]。

液體中的總碳（TC）、總無機(jī)碳（TIC）和總氮（TNb）采用Vario TOC 分析儀（德國Elementar公司）測定，樣品經(jīng)12000 r·min-1離心10 min，取上清液作為分析樣品。

沼氣中的CH4、CO2、N2和H2等氣體含量用島津GC2014 型高效氣相色譜測定，TCD 檢測器，Porapak Q 色譜柱，載氣為Ar，柱箱和檢測器溫度分別為50℃和120℃，測樣時間5 min。

耗電量用優(yōu)利德（UNI-T）UT230C 功率插座測定，測量范圍0.1～2200 W，測量精度1 級。

物質(zhì)流分析方法參見“Practical handbook of material flow analysis”[18]。物質(zhì)流分析系統(tǒng)模型應(yīng)用Stan2.5 軟件（TU Vienna）建立，物質(zhì)流分析結(jié)果應(yīng)用e!Sankey 軟件（ifu Hamburg）繪圖。

2 結(jié)果與討論

2.1 發(fā)酵前后反應(yīng)理化性質(zhì)變化

稻殼和沼渣的理化性質(zhì)見表1。稻殼原料C 含量為46.26%，熱值為19.20 kJ·g-1，C/N 為25.84，處于適宜厭氧發(fā)酵的C/N（20～30）范圍[21-22]。厭氧發(fā)酵前后C 元素的相對含量由46.26%下降到42.93%，而N 元素的相對含量由1.78%上升到1.86%，S 元素的相對含量由0.15%上升到0.27%。這是因為厭氧發(fā)酵過程中C 元素降解并轉(zhuǎn)移到生物燃?xì)庵?，而N 和S 仍主要儲存在發(fā)酵剩余物中。

接種物和沼液的理化性質(zhì)見表1和表2。接種物中TC、TNb 含量分別為1055 mg·L-1和710 mg·L-1，沼液中TC、TNb 含量分別為1525 mg·L-1和 1120 mg·L-1，而二者總有機(jī)碳（TOC,TOC=TC-TIC,305 mg·L-1和352 mg·L-1）和TC/TNb（1.49 和1.36）差異不大。

稻殼厭氧發(fā)酵前后組分的去除率分析見表3。TS 和VS 去除率分別為46.63%和57.39%；C 元素的去除率為32.77%，不考慮接種物和沼液含C 量時去除率為35.73%；N 元素去除率為10.17%，不考慮接種物和沼液中含N 量時去除率為42.88%；S元素去除率為負(fù)值，主要原因是僅考察了原料和沼渣中的S 含量，接種物中的S 可引起沼渣的S 含量增加；纖維素和半纖維素的去除率分別為42.00%和54.65%，沼渣中纖維素相對含量21.29%高于半纖維素的12.64%，說明在稻殼厭氧發(fā)酵過程中半纖維素相對更容易被降解轉(zhuǎn)化為生物燃?xì)狻?/p>

2.2 發(fā)酵過程的產(chǎn)氣性能變化[23]

稻殼厭氧發(fā)酵制備生物燃?xì)膺^程中CH4體積分?jǐn)?shù)及CH4/CO2體積比隨發(fā)酵過程的變化如圖2所示。CH4體積分?jǐn)?shù)在第5 天達(dá)到65.27%，而后稍有下降并穩(wěn)定在58%～60%之間。在發(fā)酵過程中，CH4/CO2體積比由0.69 迅速升高到2.03，之后下降并穩(wěn)定在1.60～1.70 之間。這主要與厭氧發(fā)酵所處的階段有關(guān)。在發(fā)酵初期的1～3 d，原料迅速水解酸化，此過程產(chǎn)生的氣體中CO2含量高，故CH4/CO2偏低；在發(fā)酵中后期，由于產(chǎn)甲烷菌的生長繁殖，氣體中CH4含量升高，CO2含量降低，故CH4/CO2值升高并保持穩(wěn)定。

厭氧發(fā)酵39 d 生物燃?xì)饫鄯e產(chǎn)氣量（TVbiogas）為24.29 L，生物燃?xì)夂虲H4產(chǎn)氣率隨發(fā)酵時間的變化及結(jié)果見圖3和表4。稻殼厭氧發(fā)酵制備生物燃?xì)獾漠a(chǎn)氣效果明顯優(yōu)于Contreras 等[24]的研究結(jié)果：生物燃?xì)猱a(chǎn)氣率由44 ml·(gVSRH)-1提高到297.41 ml·(gVSRH)-1，達(dá)到80%產(chǎn)氣率由21 d 縮短到13 d，CH4產(chǎn)率由19 ml·(gVSRH)-1提高到164.40 ml·(gVSRH)-1，CH4含量由45%提高到55.28%。這主要是因為稻殼原料的理化性質(zhì)不同，Contreras等選用的稻殼原料的VS 為77.8%，C 含量為37.7%，C/N 為99:1，不是適宜的厭氧發(fā)酵C/N比例。生物燃?xì)庵蠧H4含量與Nijaguna[25]的研究結(jié)果一致。

表1 接種物、調(diào)節(jié)劑、稻殼和沼渣的理化性質(zhì)Table 1 Physicochemical parameters of inoculum,additive,rice hulls and biogas residue

表2 接種物和沼液的理化性質(zhì)Table 2 Physicochemical parameters of inoculum and biogas slurry

表3 厭氧發(fā)酵過程不同物質(zhì)的去除率Table 3 Removal rates of involved substances in anaerobic fermentation process

圖2 CH4 體積分?jǐn)?shù)和CH4/CO2Fig.2 CH4 volume percentage and CH4/CO2 rate

圖3 生物燃?xì)夂图淄榈漠a(chǎn)氣率變化Fig.3 Gas production changes of biogas and CH4

理論產(chǎn)甲烷潛力（theoretical methane potential,TMP）是在標(biāo)準(zhǔn)狀態(tài)下原料完全降解達(dá)到的最大產(chǎn)甲烷率，可由Buswell 方程[26-27]計算得出。根據(jù)元素分析結(jié)果[表2，O%=VS%-(N%+C%+H%+S%)]推出稻殼厭氧發(fā)酵過程的近似化學(xué)反應(yīng)方程式如下：

計算得稻殼理論產(chǎn)CH4率為527.53 ml·(gVS)-1，甲烷與二氧化碳體積比為1.26。實驗中總CH4與總CO2體積比（1.25）與理論值接近，實際產(chǎn)CH4率僅占理論的31.16%。這是由于部分難降解的有機(jī)物難以被微生物降解利用，厭氧發(fā)酵過程中微生物生長會消耗一部分有機(jī)物。

2.3 物質(zhì)流分析系統(tǒng)的建立

為了闡明厭氧發(fā)酵過程中的物質(zhì)和能量留存和利用狀況，以稻殼單獨(dú)厭氧發(fā)酵制備生物燃?xì)獾膶嶒炦^程為模型建立了一個開放型的物質(zhì)流分析系統(tǒng)（MFA system），如圖4所示。本研究僅考慮原料的粉碎和厭氧發(fā)酵過程，并假定粉碎和加熱需要耗電并最終主要以熱量形式散失到周圍環(huán)境中。系統(tǒng)空間邊界包括原料粉碎和厭氧發(fā)酵兩個環(huán)節(jié)和10 條物質(zhì)流，系統(tǒng)時間邊界為一個發(fā)酵周期。出入系統(tǒng)的物質(zhì)：稻殼原料、接種物和調(diào)節(jié)劑作為系統(tǒng)的輸入流物質(zhì)，沼氣和發(fā)酵剩余物作為系統(tǒng)的輸出流物質(zhì)；出入系統(tǒng)的能量：化學(xué)能以物質(zhì)為載體隨物質(zhì)流動，電能是輸入流，熱能是輸出流。系統(tǒng)內(nèi)部的粉碎原料物質(zhì)流既是原料粉碎環(huán)節(jié)的輸出流又是厭氧發(fā)酵環(huán)節(jié)的輸入流。

圖4 稻殼厭氧發(fā)酵的物質(zhì)流分析系統(tǒng)模型Fig.4 System model for MFA of fermentation process (cy:system cycle,similarly hereinafter)

表4 稻殼原料制備生物燃?xì)獾漠a(chǎn)氣效果Table 4 Biogas and CH4 production results from rice hulls

2.4 發(fā)酵過程的物質(zhì)流動分析

物質(zhì)流動分析的對象是出入系統(tǒng)的實體物質(zhì)，分析的層次包含宏觀物質(zhì)層面和微觀元素層面，微觀元素分析建立在宏觀物質(zhì)分析的基礎(chǔ)上。稻殼單獨(dú)厭氧發(fā)酵系統(tǒng)的物質(zhì)流動分析數(shù)據(jù)見表5，宏觀物質(zhì)分析包含原料、接種物、調(diào)節(jié)劑、粉碎原料、生物燃?xì)夂褪Ｓ辔锏龋⒂^元素分析包含C、N 元素。生物燃?xì)獾拿芏龋é裝iogas）按氣體成分CH455%和CO245%推算為1.28 g·L-1。數(shù)據(jù)輸入建立的系統(tǒng)模型中，采用Stan 軟件，應(yīng)用IAL-IMPL2013 模 型對實驗數(shù)據(jù)進(jìn)行95%置信優(yōu)化 得到結(jié)果，桑基圖呈現(xiàn)物質(zhì)流動分析結(jié)果如圖5～圖7所示。

圖5 宏觀物質(zhì)的物質(zhì)流動分析結(jié)果Fig.5 MFA results of goods (per system cycle)

圖6 微觀C 元素的物質(zhì)流分析結(jié)果 Fig.6 MFA results of C (per system cycle)

圖7 微觀N 元素的物質(zhì)流分析結(jié)果Fig.7 MFA results of N (per system cycle)

表5 稻殼厭氧發(fā)酵系統(tǒng)物質(zhì)流動分析數(shù)據(jù)Table 5 Materials data for MFA of anaerobic fermentation process

表6 稻殼厭氧發(fā)酵系統(tǒng)能量流動分析數(shù)據(jù)Table 6 Energy data for MFA of anaerobic fermentation process

物質(zhì)流動分析結(jié)果顯示：① 稻殼原料有30.0%的物質(zhì)降解轉(zhuǎn)化為生物燃?xì)?，其余物質(zhì)進(jìn)入剩余物 中；② 原料中C 元素有30.8%流入生物燃?xì)?，其余流入剩余物，為沼?.4%、沼渣62.9%；③ 原料中N 元素主要流入剩余物中，其在剩余物中分布為沼液63.2%、沼渣36.8%。

綜上所述，稻殼厭氧發(fā)酵制備生物燃?xì)膺^程中原料及其中的C 元素轉(zhuǎn)化為生物燃?xì)獾睦寐蕛H在30%左右，N 元素主要流入剩余物中，因此發(fā)酵剩余物中含有超過原料65%的C 資源和90%的N 資源。合理管理和利用沼渣、沼液資源，既能提高資源利用效率，又能避免造成環(huán)境的污染。

2.5 發(fā)酵過程的能量流動分析

能量流動分析的對象是出入系統(tǒng)的能量，包括物質(zhì)能量和非物質(zhì)能量。物質(zhì)能量是實體物質(zhì)本身攜帶的化學(xué)能，以熱值作為依據(jù)；非物質(zhì)能量是沒有特定物質(zhì)載體的能量，如電能或熱能。因此能量的流動分析必須以物質(zhì)的流動分析為基礎(chǔ)。厭氧發(fā)酵系統(tǒng)的物質(zhì)能量對象有原料、接種物。粉碎原料、沼氣和剩余物是物質(zhì)能量，根據(jù)其單位熱值計算；調(diào)節(jié)劑因其用量極少，能量不計。電耗和散熱是非物質(zhì)能量，由于儀器工作過程功率不穩(wěn)定，用測電儀跟蹤記錄設(shè)備的耗電量，散熱能量由電耗能量確定。稻殼單獨(dú)厭氧發(fā)酵系統(tǒng)的能量流動分析數(shù)據(jù)見表6，分析同物質(zhì)流動分析。能量流分析結(jié)果如圖8和圖9所示。

圖8 能量流動分析結(jié)果Fig.8 MFA results of energy (per system cycle)

圖9 能量流動分析結(jié)果（僅化學(xué)能）Fig.9 MFA results of chemical energy (per system cycle)

能量流動分析結(jié)果顯示：① 原料中的能量有33.7%進(jìn)入生物燃?xì)庵?，其余主要流入剩余物中；?厭氧發(fā)酵實驗過程中耗電量遠(yuǎn)遠(yuǎn)大于原料的能量，主要是散熱耗能，占系統(tǒng)總能量的97.3%。

因此，稻殼厭氧發(fā)酵制備生物燃?xì)膺^程中必須注意剩余物能量的回收利用，尤其是含有原料約60%能量的沼渣，降低發(fā)酵過程的耗能關(guān)鍵在于減少系統(tǒng)的散熱，加強(qiáng)發(fā)酵過程的控溫節(jié)能措施。

3 結(jié) 論

（1）在稻殼厭氧發(fā)酵過程中，累積產(chǎn)氣率和產(chǎn)CH4分 別 為 297.41 ml·(gVSRH)-1和 164.40 ml·(gVSRH)-1，生物燃?xì)庵蠧H4平均含量為55.28%，CH4產(chǎn)率占理論的31.16%。

（2）稻殼制備生物燃?xì)獾奈镔|(zhì)轉(zhuǎn)化效率約為30.0%，發(fā)酵剩余物中含有超過原料65%的C 元素和90%的N 元素。合理管理和利用沼渣、沼液資源，是提高資源利用效率和減少環(huán)境污染的共同要求。

（3）稻殼制備生物燃?xì)獾哪芰哭D(zhuǎn)化效率約為33.7%，散熱導(dǎo)致厭氧發(fā)酵過程能量浪費(fèi)嚴(yán)重。稻殼厭氧發(fā)酵制備生物燃?xì)膺^程中，應(yīng)加強(qiáng)沼渣中能量的回收利用，降低系統(tǒng)的散熱耗能，達(dá)到提高原料能量轉(zhuǎn)化效率和節(jié)能減排的雙重目的。

在資源、能源和環(huán)境問題不斷加劇的背景下，農(nóng)業(yè)加工廢棄物的能源化和生態(tài)化利用顯得尤為重要。應(yīng)用物質(zhì)流分析方法研究了稻殼資源能源化利用過程的物質(zhì)和能量轉(zhuǎn)化效率，結(jié)果顯示厭氧發(fā)酵制備生物燃?xì)鈫我毁Y源轉(zhuǎn)化利用途徑效率不高，發(fā)酵剩余物還含有超過原料60%的資源、能源，若肆意排放會變成嚴(yán)重的污染源。因此，稻殼資源制備生物燃?xì)饽茉椿猛緩巾毰c農(nóng)業(yè)、園林和材料制造等行業(yè)鏈接起來，建立一條資源轉(zhuǎn)化和能源利用效率高、環(huán)境污染少的綠色循環(huán)生態(tài)鏈。

[1]National data(National Bureau of Statistics of China).http://data.stats.gov.cn/index

[2]Zheng G J,Zhou Y J,Zhang J,Cheng K K,Zhao X B,Zhang T,et al.Pretreatment of rice hulls for cellulase production by solid substrate fermentation [J].Journal of Wood Chemistry and Technology,2007,27:65-71

[3]Lim J S,Manan Z A,Alwi S R W,Hashim H.A review on utilisation of biomass from rice industry as a source of renewable energy [J].Renewable & Sustainable Energy Reviews,2012,16:3084-3094

[4]Vadiveloo J,Nurfariza B,Fadel J G.Nutritional improvement of rice husks [J].Animal Feed Science and Technology,2009,151:299-305

[5]Thomsen S T,Kadar Z,Schmidt J E.Compositional analysis and projected biofuel potentials from common West African agricultural residues [J].Biomass & Bioenergy,2014,63:210-217

[6]Imamoglu E,Dalay M C,Sukan F V.Regional differences in rice hulls supply for bioethanol production [J].Applied Biochemistry and Biotechnology,2013,171:2065-2074

[7]Oduguwa O O,Edema M O,Ayeni A O.Physico-chemical and microbiological analyses of fermented corn cob,rice bran and cowpea husk for use in composite rabbit feed [J].Bioresour.Technol.,2008,99:1816-1820

[8]Demirbas A,Ozturk T.Anaerobic digestion of agricultural solid residues [J].International Journal of Green Energy,2004,1:483-494

[9]Li Lianhua(李連華),Ma Longlong(馬隆龍),Yuan Zhenhong(袁振宏),Liu Xiaofeng(劉曉風(fēng)).Study on anaerobic digestion of straw stalk [J].Journal of Agro-Environment Science(農(nóng)業(yè)環(huán)境科學(xué)學(xué)報),2007,26(1):335-338

[10]Comino E,Rosso M,Riggio V.Investigation of increasing organic loading rate in the co-digestion of energy crops and cow manure mix [J].Bioresource Technology,2010,101:3013-3019

[11]Grundmann P,Ehlers M H,Uckert G.Responses of agricultural bioenergy sectors in Brandenburg (Germany) to climate,economic and legal changes:an application of Holling’s adaptive cycle [J].Energ.Policy,2012,48:118-129

[12]Benbelkacem H,Bayard R,Abdelhay A,Zhang Y,Gourdon R.Effect of leachate injection modes on municipal solid waste degradation in anaerobic bioreactor [J].Bioresource Technology,2010,101:5206-5212

[13]Zhen Feng(甄峰),Li Dong(李東),Sun Yongming(孫永明),Li Zhibing(李志兵),Yuan Zhenhong(袁振宏).High value application and purification technology of biogas [J].Environment Science & Technology(環(huán)境科學(xué)與技術(shù)),2012,35(11):103-108

[14]Rohatgi K,Prasad S V,Rohatgi P K.Release of silica-rich particles from rice husk by microbial fermentation [J].Journal of Materials Science Letters,1987,6:829-831

[15]Contreras L M,Schelle H,Sebrango C R,Pereda I.Methane potential and biodegradability of rice straw,rice husk and rice residues from the drying process [J].Water Science and Technology,2012,65:1142-1149

[16]Okeh O C,Onwosi C O,Odibo F J C.Biogas production from rice husks generated from various rice mills in Ebonyi State,Nigeria [J].Renewable Energy,2014,62:204-208

[17]Li L H,Li D,Sun Y M,Ma L L,Yuan Z H,Kong X Y.Effect of temperature and solid concentration on anaerobic digestion of rice straw in South China [J].International Journal of Hydrogen Energy,2010,35:7261-7266

[18]Brunner P H,Rechberger H.Practical Handbook for Material Flow Analysis[M].Boca Raton,Fla.London:CRC,2004

[19]Rechberger H,Cencic O,Fruhwirth R.Uncertainty in material flow analysis [J].J.Ind.Ecol.,2014,18:159-160

[20]Yu Qiang,Zhuang Xinshu,Lv Shuangliang.Liquid hot water pretreatment of sugarcane bagasse and its comparison with chemical pretreatment methdes for the sugar recovery and structural hanges [J].Bioenergy Research,2013,129:592-598

[21]Kitani O,Hall C W.Biomass Handbook[M].New York:Gordon and Breach Science Publishers,1989.https://www.ecn.nl/phyllis2/ Browse/ Standard/ECN-Phyllis#rice husk

[22]Weiland P.Bases of the Methane Fermentation — Biology and Substrate//Biogas als Regenerative Energie — Stand und Perspektiven [M].Tagung Hanover,2001:19-33

[23]Mustonen S,Raiko R,Luukkanen J.Bioenergy consumption and biogas potential in Cambodian households [J].Sustainability,2013,5(5):1875-1892

[24]Contreras L M,Schelle H,Sebrango C R,et al.Methane potential and biodegradability of rice straw,rice husk and rice residues from the drying process [J].Water Science and Technology,2012,65(6):1142-1149

[25]Nijaguna B.Biogas Technology[M].New Delhi,India:New Age International (P) Ltd.,2002

[26]Buswell A M,Muller H F.Mechanics of methane fermentation [J].Journal of Industrial Engineering Chemistry,1955,44(3):550-559

[27]Li Dong(李東),Ye Jingqing(葉景清),Zhen Feng(甄峰).Effects of anaerobic co-digestion of different proportions between rice straw and chicken manure on biogas yield rate [J].Transactions of Chinese Society of Agricultural Engineering(農(nóng)業(yè)機(jī)械學(xué)報),2013,29(2):232-238