郭兆峰 班婷 馬艷 沈衛(wèi)強(qiáng)

摘要：為提高11JF-1.2型棉花秸稈飼料化發(fā)酵機(jī)的工作性能，保證輸出物料的質(zhì)量，對該機(jī)具工作參數(shù)進(jìn)行優(yōu)化。通過開展單因素試驗和正交試驗，并運用 Design-Expert 10.0軟件對試驗所得數(shù)據(jù)進(jìn)行響應(yīng)面優(yōu)化分析。試驗結(jié)果表明：各因素對評價指標(biāo)的貢獻(xiàn)率從大到小依次為主軸轉(zhuǎn)速、物料密度、加熱機(jī)溫度，通過軟件中的優(yōu)化求解功能得到最佳工作參數(shù)組合為主軸轉(zhuǎn)速5r/min、物料密度176.5kg/m3、加熱機(jī)溫度39.5℃。根據(jù)得到的最優(yōu)參數(shù)組合進(jìn)行試驗驗證，含水率方差平均值為1.009，發(fā)酵溫度平均值為36.9℃，優(yōu)化后機(jī)具性能提升明顯且滿足設(shè)計要求。

關(guān)鍵詞：棉花秸稈；飼料化；發(fā)酵機(jī)；參數(shù)優(yōu)化；正交試驗

中圖分類號：S226.5

文獻(xiàn)標(biāo)識碼：A

文章編號：20955553 （2023） 070111

07

Experiment and parameter optimization of 11JF-1.2 cotton straw

feeding-fermentation machine

Guo Zhaofeng1， 2， Ban Ting1， 2， Ma Yan1， 2， Shen Weiqiang1， 2

（1. Institute of Agricultural Machinery， Xinjiang Academy of Agricultural Science， Urumqi， 830091， China;

2. Key Laboratory of Intelligent Control Technology for Xinjiang Facility Agriculture， Urumqi， 830091， China）

Abstract： In order to improve the performance of 11JF-1.2 cotton straw feeding-fermentation machine and ensure the quality of output materials， the working parameters of the machine were optimized. By carrying out single-factor tests and orthogonal tests and using Design-Expert 10.0 software， response surface optimization analysis was carried out on the data obtained from the test. The test results show that the contribution rate of each factor to the evaluation index， from largest to smallest， was spindle speed， material density， and heater temperature. Through the optimization and solving function in the software， the optimal combination of working parameters was spindle speed at 5r/min， material density at 176.5kg/m3， and heater temperature at 39.5℃. According to the obtained optimal parameter combination， the average moisture content variance is 1.009， and the average fermentation temperature is 36.9℃. After optimization， the performance of the machine was significantly improved and met the design requirements.

Keywords： cotton straw； feeding； fermentation mechine； parameter optimization； orthogonal test

0 引言

農(nóng)作物秸稈是指在農(nóng)業(yè)生產(chǎn)過程中，稻谷、玉米、小麥等作物收獲后剩余的莖和少部分葉［1］。我國作為農(nóng)業(yè)大國，農(nóng)作物種植面積廣闊且豐富，因此收獲產(chǎn)生的農(nóng)作物秸稈數(shù)量大，種類較多以及分布范圍較廣。農(nóng)作物秸稈作為農(nóng)業(yè)生產(chǎn)之余的副產(chǎn)品，含有豐富的有機(jī)質(zhì)和微量元素，可以育肥土壤、增加土地生產(chǎn)力，也可作為畜牧、養(yǎng)殖業(yè)發(fā)展的有機(jī)飼料［27］。

我國是世界上最大的棉花生產(chǎn)大國之一，棉花生產(chǎn)基地主要位于地處歐亞大陸腹地的新疆。新疆氣候干燥、干旱少雨，適合大面積種植商品棉，近年來，新疆棉花種植面積和產(chǎn)量以及質(zhì)量都位于世界前列［8］。2021 年我國棉花種植總面積為3028.1khm2，新疆棉花種植面積達(dá)到2506.1khm2，產(chǎn)量約占全國產(chǎn)量的85％，約占全球棉花產(chǎn)量的20％。棉花摘除后，棉秸稈主要剩余莖、葉和棉桃殼，同一植株棉秸稈中棉莖占73%，棉桃殼占27%，有機(jī)物含量占比在91%以上，其中纖維素含量約40%左右［912］。棉秸稈作為一種豐富的可再生資源，應(yīng)用潛力巨大，新疆棉秸稈目前主要的處理方式為粉碎還田肥土和作為燃料使用，科學(xué)合理利用棉秸稈可推動農(nóng)業(yè)農(nóng)村可持續(xù)發(fā)展。

為了實現(xiàn)農(nóng)作物秸稈高效利用，秸稈的基質(zhì)化、飼料化和能源化等得到了諸多學(xué)者的關(guān)注和研究。吳慧等［13］將腐熟的棉花秸稈按基質(zhì)成分不同進(jìn)行分類配比，然后用在番茄幼苗上觀察長勢，結(jié)果表明棉花秸稈基質(zhì)對番茄幼苗的理化性質(zhì)以及其他方面長勢指標(biāo)優(yōu)于其他配比基質(zhì)。張艷萍等［14］利用腐熟的玉米秸稈培養(yǎng)篩選土壤微生物，獲得可降解秸稈纖維的真菌，可解決秸稈還田后普遍存在降解所需時間長的問題。章智鈞等［15］利用棉花秸稈等為原材料，研究不同配比基質(zhì)對草莓植株的影響，結(jié)果表明腐熟棉花秸稈可代替部分草炭作為草莓無土栽培基質(zhì)使用。 楊雪海等［16］使用同一品種不同時期的油菜秸稈作為試驗原料，通過營養(yǎng)成分檢驗、檢測、分析等方法分析油菜秸稈的飼用價值。Djibrillou等［17］研究了奶牛對未處理棉秸稈（UNTS）、尿素處理棉秸稈（UTS）、未處理棉秸稈+400g棉籽（UNTS+CS）的消化情況，通過方差分析顯示，尿素處理的棉秸稈可能是哺乳期馬拉地犬適宜的飼料，并發(fā)現(xiàn)添加棉籽的未處理秸稈也比單獨未處理秸稈有優(yōu)勢。

本研究針對新疆棉花種植面積廣、產(chǎn)量大、飼料化利用率低的特點和問題，研發(fā)了一種棉秸稈飼料化發(fā)酵設(shè)備，并對該設(shè)備工作參數(shù)進(jìn)行優(yōu)化，旨在充分利用新疆棉花秸稈資源，對深入研究棉花秸稈飼料化發(fā)酵技術(shù)有重要參考價值。

1 整機(jī)結(jié)構(gòu)與工作過程

1.1 整機(jī)結(jié)構(gòu)

11JF-1.2型棉花秸稈發(fā)酵設(shè)備配有上料系統(tǒng)、攪拌系統(tǒng)、加熱機(jī)送風(fēng)系統(tǒng)、出料提升系統(tǒng)以及電控系統(tǒng)等裝置，為棉秸稈飼料化加工發(fā)酵、工況參數(shù)優(yōu)化提供試驗平臺，設(shè)備簡圖如圖1所示，整機(jī)主要技術(shù)參數(shù)和攪拌系統(tǒng)主要參數(shù)分別如表1、表2所示。

1.2 工作過程

棉花秸稈發(fā)酵設(shè)備工作時，首先打開電源，通過電控表開啟風(fēng)機(jī)和加熱機(jī)，再啟動主軸攪拌系統(tǒng)，觀察設(shè)備正常運行1min后，通過上料箱向攪拌桶輸送粉碎后的棉花秸稈與配料混合物，此時風(fēng)機(jī)產(chǎn)生的氣流經(jīng)過加熱機(jī)形成熱風(fēng)，通過攪拌桶底部由負(fù)壓風(fēng)腔排出，攪拌桶內(nèi)混合物料溫度逐漸上升，待攪拌桶飽和后停止送料并閉合上蓋板，調(diào)節(jié)風(fēng)速和溫度并持續(xù)攪拌，當(dāng)達(dá)到要求發(fā)酵時間后打開出料口，物料通過出料口堆集、打包、封存發(fā)酵。攪拌系統(tǒng)結(jié)構(gòu)如圖2所示，主要由螺旋葉片、葉片支架、攪拌刀片、主軸以及圖中未表示出的攪拌桶組成。

2 試驗與分析

2.1 試驗方案

試驗地點位于新疆維吾爾自治區(qū)巴音郭楞蒙古自治州博湖縣進(jìn)行，試驗品種選擇采摘棉花后的新疆1號棉花秸稈，試驗前需要對棉花秸稈進(jìn)行前處理，即先將新疆1號棉花秸稈風(fēng)干或烘干至含水率為15%左右，使用秸稈粉碎機(jī)將棉花秸稈揉切成3～5cm小段，再加水、菌劑和麩皮充分?jǐn)嚢杌旌?，制成試驗所需物料。參考?biāo)準(zhǔn)GB/T 24675.6—2009《保護(hù)性耕作機(jī)械秸稈粉碎還田機(jī)》和NY/T 3020《農(nóng)作物秸稈綜合利用技術(shù)通則》、NY/T 1701《農(nóng)作物秸稈資源調(diào)查與評價技術(shù)規(guī)范》，設(shè)置主軸轉(zhuǎn)速、物料密度以及加熱機(jī)溫度為試驗因素，發(fā)酵均勻性Hf和發(fā)酵溫度Tf作為試驗指標(biāo)。

為保證喂入量恒定并減少試驗誤差，設(shè)計上料箱滿載重量為50±1kg，并在喂料過程中充分?jǐn)嚢枇铣匚锪匣旌衔?，使每次上料的密度差值在可接受的誤差范圍以內(nèi)。單因素試驗中，每組試驗只調(diào)整單個試驗因素參數(shù)，主軸轉(zhuǎn)速通過變頻調(diào)速器控制，物料密度是前處理階段的重點，需通過調(diào)整棉花秸稈小段、水、麩皮和菌劑的配比制成三種不同密度的試驗物料以備使用。改變加熱機(jī)溫度可以使發(fā)酵設(shè)備內(nèi)部溫度變化，物料攪拌過程中實時溫度可通過測溫槍測量。每組試驗結(jié)束后，必須嚴(yán)格清理發(fā)酵箱內(nèi)部的物料殘渣，避免物料混合造成密度誤差。對于發(fā)酵溫度，需要在試驗開始1h、1.5h和2h三個時刻使用測溫槍隨機(jī)測試三個位置的溫度并取平均值，每組試驗剛結(jié)束，需在出料口隨機(jī)5個位置取樣（每個位置取樣0.5kg），使用具備稱重功能的水分測定儀器（YLS16A鹵素水分測定儀）進(jìn)行含水率測定，由于試驗物料的復(fù)雜性，物料發(fā)酵均勻性無法實際測得，因此使用干基含水率方差來代替物料發(fā)酵的均勻性。含水率方差［18］計算公式如式（1）所示。

σ2=∑（X－μ）2N

（1）

式中：

X——每組樣本的含水率；

μ——樣本含水率的平均值；

N——樣本數(shù)量。

為了使試驗結(jié)果更加準(zhǔn)確、科學(xué)，每組單因素試驗需要重復(fù)進(jìn)行3次。

2.2 單因素試驗

圖3（a）為主軸轉(zhuǎn)速對發(fā)酵均勻性和發(fā)酵溫度的影響，設(shè)定加熱機(jī)溫度為39℃，物料密度為175kg/m3，主軸轉(zhuǎn)速3～7r/min范圍內(nèi)取5個值對秸稈發(fā)酵機(jī)發(fā)酵性能進(jìn)行評價。隨著主軸轉(zhuǎn)速逐漸升高，含水率方差呈下降趨勢，發(fā)酵溫度先上升后下降。含水率方差降低，說明棉秸稈物料發(fā)酵均勻性越好，發(fā)酵溫度隨主軸轉(zhuǎn)速增大呈二次曲線先上升后下降，歸因于主軸轉(zhuǎn)速越慢，物料混合均勻性差，測溫時局部出現(xiàn)極值溫度，隨著主軸轉(zhuǎn)速越快，物料混合均勻程度變高，發(fā)酵溫度更加趨于平穩(wěn)，不會出現(xiàn)極值溫度。圖3（b）為加熱機(jī)溫度對含水率方差和發(fā)酵溫度的影響，設(shè)定主軸轉(zhuǎn)速為5r/min，物料密度為175kg/m3，加熱機(jī)溫度在35℃～42℃之間取值。試驗結(jié)果表明，發(fā)酵溫度隨加熱機(jī)溫度增高呈正比例趨勢上升，當(dāng)加熱機(jī)溫度達(dá)到38℃后含水率方差趨于平穩(wěn)，歸因于溫度越高，物料失水率上升，含水率測定誤差減小，更加準(zhǔn)確。圖3（c）為物料密度對發(fā)酵性能的影響，當(dāng)物料密度為176kg/m3時含水率方差最低，物料均勻性最佳，發(fā)酵溫度隨密度增大呈逐漸下降趨勢。

2.3 正交試驗

為了深入研究主軸轉(zhuǎn)速、 加熱機(jī)溫度和物料密度三個試驗因素與秸稈發(fā)酵機(jī)評價指標(biāo)之間的交互作用，探究各因素影響規(guī)律并尋找較優(yōu)的工作參數(shù)組合，結(jié)合單因素試驗結(jié)果確定各因素水平和取值范圍，以含水率方差和發(fā)酵溫度為評價指標(biāo)進(jìn)行正交試驗。

依據(jù)單因素試驗結(jié)果，選取主軸轉(zhuǎn)速中間值為5r/min，物料密度中間值為175±0.5kg/m3， 加熱機(jī)溫度中間值為39℃，為提高試驗準(zhǔn)確性和科學(xué)性，對試驗因素進(jìn)行編碼如表3所示，試驗正交如表4所示。A、B、C為因素編碼值。

2.4 回歸模型建立與顯著性檢驗

運用Design-Expert 10軟件中Box-Behnken功能［19］設(shè)計試驗，并進(jìn)行多元回歸擬合分析，建立主軸轉(zhuǎn)速、物料密度和加熱機(jī)溫度三個因素與含水率方差、發(fā)酵溫度兩個評價指標(biāo)的響應(yīng)面正交回歸模型，再進(jìn)行方差分析，結(jié)果如表5所示。根據(jù)方差分析結(jié)果，含水率方差指標(biāo)模型為P=0.0039（P＜0.01），說明回歸模型達(dá)到了極顯著水平，擬合精度高，含水率方差和發(fā)酵溫度兩評價指標(biāo)失擬項P值分別為0.0648和0.0932（失擬項P＞0.05），表明失擬項不顯著。綜上，可以利用該響應(yīng)面正交回歸模型進(jìn)行參數(shù)優(yōu)化設(shè)計［20］。

在含水率方差指標(biāo)模型中，因素A、A2的P值均小于0.01，對回歸模型影響極顯著；因素B、AB、AC、BC、B2、C2的P值均大于0.05，對回歸模型影響不顯著；發(fā)酵溫度指標(biāo)模型中，因素A的P值小于0.01，對回歸模型影響極顯著；因素AB、AC、BC、B2、C2的P值均大于0.05，對回歸模型影響不顯著，由此可得模型回歸項優(yōu)化方程為

Y1=1.06－0.025A－0.016C+0.029A2

（2）

Y2=

36.68+0.34A+0.11B+0.15C－

0.27A2

（3）

貢獻(xiàn)率［21］Δj可以反映出單個參數(shù)對所建回歸模型的影響程度，Δj越大則影響程度就越大。貢獻(xiàn)率Δj可用式（4）、式（5）計算，即

δ=

0F≤1

1－1FF＞1

（4）

Δj=δj+12∑m

i=1

i≠jδij+δjj j=1，2，…，m

（5）

式中：

δ——考核值；

F——方差分析表中F值；

δj——一次項貢獻(xiàn)率，%；

δjj——二次項貢獻(xiàn)率，%；

δij——交互貢獻(xiàn)率，%。

根據(jù)上式計算得到各個因素貢獻(xiàn)率如表6所示，即回歸模型中每個因素對評價指標(biāo)的貢獻(xiàn)率從大到小順序為主軸轉(zhuǎn)速、物料密度、加熱機(jī)溫度。

圖4為各因素組合對含水率方差響應(yīng)曲面。

如圖4（a）所示，物料密度增高時，含水率方差呈先下降后上升趨勢，含水率方差與主軸轉(zhuǎn)速呈負(fù)相關(guān)關(guān)系。主軸轉(zhuǎn)速變化對響應(yīng)曲面影響更加明顯，歸因于當(dāng)主軸轉(zhuǎn)速增大時，單位時間內(nèi)物料被攪拌頻率增加，受熱均勻、混合程度高；如圖4（b）所示，主軸轉(zhuǎn)速和加熱機(jī)溫度增高時，含水率方差呈上升趨勢。加熱機(jī)溫度變化對含水率方差影響較小，主軸轉(zhuǎn)速對響應(yīng)曲面影響更加顯著；如圖4（c）所示，加熱機(jī)溫度與物料密度對含水率方差都呈現(xiàn)正相關(guān)關(guān)系，物料密度對含水率方差影響更明顯，歸因于當(dāng)主軸轉(zhuǎn)速不變、物料密度逐漸增大時，攪拌裝置單位時間內(nèi)需要處理的物料總量增加，會出現(xiàn)受熱不均勻、攪拌不到位的情況，影響發(fā)酵物料含水率變化。

以主軸轉(zhuǎn)速、物料密度以及加熱機(jī)溫度作為工作參數(shù)優(yōu)化對象，以物料含水率方差和發(fā)酵溫度作為優(yōu)化目標(biāo)，利用分析軟件中的優(yōu)化求解功能進(jìn)行較優(yōu)參數(shù)篩選，約束條件方程為

minY1（A，B，C）

Y2（37）

－1≤A≤1

－1≤B≤1

－1≤C≤1

（4）

得到棉花秸稈飼料化發(fā)酵設(shè)備最佳工作參數(shù)組合為：主軸轉(zhuǎn)速5.26r/min，物料密度176.54kg/m3，加熱機(jī)溫度為39.46℃，該參數(shù)組合下物料含水率方差和發(fā)酵溫度分別為1.005和37℃。

2.5 試驗驗證

為了方便操作控制并充分考慮各因素參數(shù)單位精度問題，對所得最佳參數(shù)進(jìn)行圓整，調(diào)整后的參數(shù)分別為：主軸轉(zhuǎn)速5r/min，物料密度176.5kg/m3，加熱機(jī)溫度為39.5℃。根據(jù)最優(yōu)參數(shù)進(jìn)行3次驗證試驗，試驗現(xiàn)場如圖5所示，統(tǒng)計計算含水率方差和發(fā)酵溫度并求平均值，如表7所示。

由表7可知，在最優(yōu)工作參數(shù)組合下含水率方差平均值為1.009，發(fā)酵溫度平均值為36.9℃，含水率方差接近1，說明該發(fā)酵設(shè)備輸出產(chǎn)物均勻性更加穩(wěn)定，優(yōu)化效果明顯，滿足機(jī)具設(shè)計要求［22］。

3 結(jié)論

1） 確定使用含水率方差代替物料均勻性作為試驗指標(biāo)，通過單因素試驗與正交試驗，得到了相關(guān)工作參數(shù)取值范圍為主軸轉(zhuǎn)速4～6r/min、物料密度175～177kg/m3、 加熱機(jī)溫度38℃～40℃，確定正交試驗因素為主軸轉(zhuǎn)速、物料密度、加熱機(jī)溫度。

2） 利用Design Expert軟件對試驗中獲取的數(shù)據(jù)進(jìn)行方差與響應(yīng)面分析，結(jié)果表明主軸轉(zhuǎn)速對含水率方差影響極顯著，各因素貢獻(xiàn)率順序從大到小為主軸轉(zhuǎn)速、物料密度、加熱機(jī)溫度。使用軟件中優(yōu)化求解功能對工作參數(shù)進(jìn)行優(yōu)化，得到較優(yōu)工作參數(shù)組合為主軸轉(zhuǎn)速5r/min，物料密度176.5kg/m3，加熱機(jī)溫度為39.5℃，該參數(shù)組合下平均含水率方差為1.009，平均發(fā)酵溫度為36.9℃。

參 考 文 獻(xiàn)

［1］ 那偉， 趙新穎， 祝延立， 等. 吉林省主要農(nóng)作物秸稈綜合利用現(xiàn)狀、問題及對策分析［J］. 黑龍江農(nóng)業(yè)科學(xué)， 2019（3）： 169-173.

Na Wei， Zhao Xinying， Zhu Yanli， et al. Current situation， problems and countermeasures for comprehensive utilization of main crop straw in Jilin Province ［J］. Heilongjiang Agricultural Sciences， 2019（3）： 169-173.

［2］ 張湘云， 劉利珍. 農(nóng)作物秸稈飼料化利用補(bǔ)貼政策的實施現(xiàn)狀及優(yōu)化建議［J］. 中國飼料， 2022（12）： 126-129.

Zhang Xiangyun， Liu Lizhen. Implementation status and optimization suggestions of subsidy policy for feed utilization of crop straw ［J］. China Feed， 2022（12）： 126-129.

［3］ Bao J C， Yu J H， Feng Z， et al. Situation of distribution and utilization of crop straw resources in seven western provinces， China ［J］. Chinese Journal of Applied Ecology， 2014， 25（1）： 181-187.

［4］ He K， Zhang J， Zeng Y. Rural households willingness to accept compensation for energy utilization of crop straw in China ［J］. Energy， 2018， 165： 562-571.

［5］ Wang B， Shen X， Chen S， et al. Distribution characteristics， resource utilization and popularizing demonstration of crop straw in southwest China： A comprehensive evaluation ［J］. Ecological Indicators， 2018， 93： 998-1004.

［6］ 田宜水， 趙立欣， 孫麗英， 等. 農(nóng)作物秸稈資源調(diào)查與評價方法研究［J］. 中國人口·資源與環(huán)境， 2011， 21（S1）： 583-586.

Tian Yishui， Zhao Lixin， Sun Liying， et al. Study on crop straw survey and evaluation method［J］. China Population， Resources and Environment， 2011， 21（S1）： 583-586.

［7］ 張永康， 李龍， 王娟， 等. 農(nóng)作物秸稈資源調(diào)查與評價［J］. 現(xiàn)代農(nóng)業(yè)科技， 2013（16）： 202-203.

Zhang Yongkang， Li Long， Wang Juan， et al. Investigation and evaluation of crop straw resources ［J］. Modern Agricultural Science and Technology， 2013（16）： 202-203.

［8］ 全國農(nóng)作物秸稈資源調(diào)查與評價報告［J］. 農(nóng)業(yè)工程技術(shù)（新能源產(chǎn)業(yè)）， 2011（2）： 2-5.

［9］ 王雪穎， 呂慧捷， 袁小偉， 等. 新疆棉花產(chǎn)業(yè)鏈發(fā)展現(xiàn)狀及展望［J］. 山西農(nóng)經(jīng)， 2022（9）： 1-4， 8.

［10］ Soni B， Barbary E， Mahmoud B. Chemical isolation and characterization of different cellulose nanofibers from cotton stalks［J］. Carbohydrate Polymers， 2015， 134： 581-589.

［11］ Al Afif R， Wendland M， Amon T， et al. Supercritical carbon dioxide enhanced pre-treatment of cotton stalks for methane production［J］. Energy， 2020， 194： 116903.

［12］ 趙樹琪， 李蔚， 戴寶生，等. 棉花秸稈綜合利用現(xiàn)狀分析［J］. 湖北農(nóng)業(yè)科學(xué)， 2017， 56（12）： 2201-2203.

Zhao Shuqi， Li Wei， Dai Baosheng， et al. The analysis of comprehensive utilization present situation of cotton stalk ［J］. Hubei Agricultural Sciences， 2017， 56（12）： 2201-2203.

［13］ 吳慧， 阿衣甫， 吳默涵. 不同配比棉花秸稈基質(zhì)對番茄幼苗生長的影響［J］.長江蔬菜， 2015（6）： 29-32.

Wu Hui， A Yifu， Wu Mohan. Effects of different mixed cotton straw substrates on growth of tomato seedlings ［J］. Journal of Changjiang Vegetables， 2015（6）： 29-32.

［14］ 張艷萍， 趙瑛， 張運暉. 玉米秸稈降解真菌的篩選鑒定［J］. 甘肅農(nóng)業(yè)科技， 2022， 53（5）： 55-59.

Zhang Yanping， Zhao Ying， Zhang Yunhui. Screening and identification of fungi for maize straw degradation ［J］. Gansu Agricultural Science and Technology， 2022， 53（5）： 55-59.

［15］ 章智鈞， 孫軍利， 張坤， 等. 利用棉稈混合基質(zhì)進(jìn)行草莓栽培試驗的研究［J］. 農(nóng)業(yè)科技與裝備， 2016（8）： 25-27， 30.

Zhang Zhijun， Sun Junli， Zhang Kun， et al. Research on experiment in cultivation of strawberry by mixed matrix of cotton stalks ［J］. Agricultural Science and Technology and Equipment， 2016（8）： 25-27， 30.

［16］ 楊雪海， 郭萬正， 黃少文， 等. 不同生長階段的油菜飼用營養(yǎng)價值評價研究［J］. 飼料工業(yè)， 2017， 38（21）： 19-22.

Yang Xuehai， Guo Wanzheng， Huang Shaowen， et al. Study on rape forage nutritive value evaluation in different growth stages ［J］. Feed Industry， 2017， 38（21）： 19-22.

［17］ Djibrillou O A， Pandey V S， Gouro S A， et al. Effect of urea-treated or untreated straw with cotton seed on performances of lactating Maradi （Red Sokoto） goats in Niger ［J］. Livestock Production Science， 1998， 55（2）： 117-125.

［18］ 王志剛， 胡勛鋒. 隨機(jī)變量數(shù)學(xué)期望與方差［J］. 海南大學(xué)學(xué)報（自然科學(xué)版）， 2009， 27（4）： 329-331， 335.

Wang Zhigang， Hu Xunfeng. Expected value and variance of random variable ［J］. Natural Science Journal of Hainan University， 2009， 27（4）： 329-331， 335.

［19］ 范玲， 唐佳蒙， 王鑫， 等. 正交試驗法與Box-Behnken響應(yīng)面法優(yōu)選白殘花—烏梅藥對提取工藝比較［J］. 中國藥業(yè)， 2022， 31（15）： 38-43.

Fan Ling， Tang Jiameng， Wang Xin， et al. Comparison of the orthogonal test and the BBD-RSM in optimizing the extraction process of Rosae Multiflorae Flos-Mume Fructus drug pair ［J］. China Pharmaceuticals， 2022， 31（15）： 38-43.

［20］ 馬艷， 班婷， 郭兆峰， 等. 基于超細(xì)粉碎機(jī)的棉花秸稈基質(zhì)對黃瓜穴盤育苗試驗研究［J］. 中國農(nóng)機(jī)化學(xué)報， 2020， 41（2）： 196-199， 205.

Ma Yan， Ban Ting， Guo Zhaofeng， et al. Experimental study on cotton straw substrate of cucumber seedling based on ultra-fine mill ［J］. Journal of Chinese Agricultural Mechanization， 2020， 41（2）： 196-199， 205.

［21］ 馬秋成， 盧安舸， 高連興， 等. 蓮子物料空氣動力學(xué)特性與殼仁分離裝置試驗［J］. 農(nóng)業(yè)工程學(xué)報， 2015， 31（6）： 297-303.

Ma Qiucheng， Lu Ange， Gao Lianxing， et al. Aerodynamic characteristics of lotus seed mixtures and test on pneumatic separating device for lotus seed kernel and contaminants ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2015， 31（6）： 297-303.

［22］ 郭兆峰， 班婷， 方躍， 等. 棉秸稈基質(zhì)化利用設(shè)備的研制及試驗［J］. 中國農(nóng)機(jī)化學(xué)報， 2020， 41（4）： 163-168.

Guo Zhaofeng， Ban Ting， Fang Yue， et al. Developing and test of cotton straw substrate utilization equipment ［J］. Journal of Chinese Agricultural Mechanization， 2020， 41（4）： 163-168.