基于D-optimal法优化香菇菌种培养基配方的研究

摘要：为了筛选和优化香菇原种及栽培种的培养基质配方，采用D-optimal设计方法，以麦粒和木屑不同配比为原料优化香菇原种培养基质，以木屑、玉米芯、麸皮不同配比为原料优化香菇栽培种培养基质，以香菇品种L808作为供试菌种，分别以其菌丝萌发期、菌丝长速、满袋期为评价指标，通过对各评价指标的测量，建立了各配比基质与香菇培养基质配方响应值之间的回归模型，从而科学的优化出香菇原种及栽培种栽培基质的配方。试验结果表明，香菇原种栽培基质最优配方为50%麦粒+50%木屑；香菇栽培种栽培基质最优配方为37.69%玉米芯+23.33%麸皮+38.98%木屑。在以上2个配方的栽培基质接种后，香菇菌丝的生长旺盛，萌发期短、满袋期短，且理化性质较优，说明优化得到的栽培基质配方具有较高的可行性，该设计方法也在优化培养料配比上是科学并且可行的。

关键词：D-optimal法；香菇；原种；栽培种；培养基质；配方

中图分类号：S646.1 文献标志码：A 文章编号：2097-2172（2024）08-0724-10

doi：10.3969/j.issn.2097-2172.2024.08.006

Optimization of Lentinula edodes Spawn Medium Formula Based on

D-optimal Method

REN Aiming 1， Bao Yzheng 2， Han Aimin 2， LI Tong 1， LIU Mingjun 1， WANG Xiaowei 1，

YANGJianjie 1， Yang Qin 1， YANG Renlu 3， FU Aifang 3

（1. Institute of Vegetables， Gansu Academy of Agricultural Sciences， Lanzhou Gansu 730070， China; 2. Agricultural Technology Extension Centre of Tanchang County， Tanchang Gansu 742500， China; 3. Tanchang County Xingtan

Fungi Industry Co.， Ltd.， Tanchang Gansu 742500， China）

Abstract： In order to screen and optimize the culture medium formula of Lentinula edodes original and cultivated species， the D-optimal design method was used to optimize the original culture medium of Lentinula edodes with different ratios of wheat grain and sawdust as raw materials， and the cultivation medium of Lentinula edodes with different ratios of sawdust， corncob and bran as raw materials. The Lentinula edodes variety L808 was used as the test strain， and the mycelial germination period， mycelial growth rate and full bag period were used as the evaluation indexes， respectively. Through the measurement of the indexes， the regression model between the ratio matrix and the response value of the Lentinula edodes culture medium formula was established， so as to scientifically optimize the formula of the original and cultivated species of Lentinula edodes. The results showed that the optimal formula of Lentinula edodes original culture was： 50% wheat grain + 50% sawdust; the optimal formula of Lentinula edodes cultivation was： 37.69% corncob+23.33% bran+38.98% sawdust. After inoculating the cultivation substrates of the above two formulas， the mycelium of Lentinula edodes grew vigorously， with a short germination period and a short full bag period， and the physicochemical properties were also superior， indicating that the optimized cultivation substrate formula is highly feasible. This design method is scientific and feasible in optimizing the ratio of cultivation materials.

Key words： D-optimal design; Lentinula edodes; Original species; Cultivated species; Culture medium; Formula

香菇（Lentinula edodes）是目前国内种植规模最大、产量最高的食用菌［1 ］。菌种是香菇栽培中的关键，香菇菌种的优劣通常取决于纯度高、活力强的香菇母种以及原种和栽培种培养料的科学合理选择［2 ］。近年来，香菇菌种的培养以及菌棒的栽培主要原料为各种阔叶树木的木屑，由于全国乡村振兴战略的开展，香菇栽培面积的扩大，导致栽培原材料木屑缺乏严重［3 ］，为了缓解木屑资源短缺的问题，培养料中可添加玉米芯、棉籽壳、秸秆等代替木屑培养香菇原种、栽培种及栽培香菇［4 - 6 ］，这也是解决目前面临资源短缺困境的有效途径。同时，使用部分麦粒、玉米芯和麸皮代替木屑培养香菇菌种能有效利用多种农业废弃物，缓解对木屑资源的采用，以及减少栽培种植成本，促进农业的可持续发展［3 ］。但培养基原料的不同，也会影响原种、栽培种的菌丝的生长速度及活力［7 ］，并影响后期香菇的营养成分、品质、产量［8 ］。因此进行科学的筛选并优化菌种的培养料配方并培养出优质的菌种是香菇栽培过程中的关键途径。杨建杰等［9 ］根据对多年文献的汇总，选择出多个优良的配方，再结合当地地理优势筛选并优化出香菇的栽培种配方。夏敏等［10 ］通过有限的配方均匀设计筛选并优化出了玉米秸秆代料栽培香菇的优良配方。周韬［11 ］采用正交试验设计优化出了香菇液体菌种的优良配方。但是目前采用D-optimal设计优化原种及栽培种配方却鲜有研究，更多的使用在蔬菜基质的优化试验和药品处方的优化方面［12 - 13 ］。

本试验将通过D-optimal设计对香菇的原种及栽培种的培养料配方进行科学的优化，即用麦粒代替木屑设计不同的配比对原种培养基质配方进行优化，用玉米芯和麸皮代替部分木屑设计不同配比的栽培种培养料配方进行优化，从而获得一个优良的原种培养料配方和一个栽培种培养料配方，最终为有效地缓解目前木屑资源短缺的问题提供科学参考。

1 材料与方法

1.1 材料与设备

供试菌种为香菇L808试管种及麦粒种，由甘肃省农业科学院蔬菜研究所食用菌研究室提供。供试木屑购买于陕西省汉中市宁强县，供试麦粒、玉米芯和麸皮均购买于甘肃省白银市景泰县。各供试原料的理化性质如表1所示。

供试设备有QY-8型电脑程序控制装袋机（启源机械有限公司）、免锅炉节能环保灭菌锅（诸城市良工机械有限公司）、ST20型pH仪（奥豪斯仪器常州有限公司）。

1.2 试验方法

1.2.1 混料配方设计 根据前期研究及相关文 献［14 - 16 ］，本试验将原种所用的麦粒用量设置为0～50.0%，木屑用量设置为50.0%～100.0%；将栽培种所用的木屑的用量设置为20.0%～100.0%，玉米芯的用量设置为0～40.0%，麸皮的用量设置为0～40.0%。根据Design-Expert 8.0.6软件中的Mixture-Optimal进行设计试验，原种栽培料共设计出13个混料组合配方（表2），栽培种栽培料共设计出16个混料组合配方（表3）。

1.2.2 制作方法 分别按照表2、表3中不同混料组合配方将原料混合均匀，使用装料机将原料装入17 cm×33 cm的菌袋，每个配方制作15袋，完成装袋后转移至灭菌锅121 ℃灭菌3 h，灭菌完成后转移至已消毒的接种室中冷却备用。随机从灭菌后每个配方中的菌袋中取出6袋用于测量培养基质的容重、总孔隙度、通气孔隙度、持水孔隙度、pH、含水量、含氮量；不同混料组合配方的原种培养基质中均使用试管种接种，栽培种培养基质中均使用麦粒种接种。接种后转移至养菌室中培养，并观察记录菌丝长速、萌发期、满袋期。每个处理均重复6次。

1.2.3 指标测量 采用环刀法并参考倪琳［17 ］的方法测量培养料的容重、总孔隙度、通气孔隙度、持水孔隙度；参考王涛等［18 ］方法并稍作修改，测量培养料pH。即称取5 g培养料，放入50 mL去离子水中充分震荡，静置30 min后取上清液使用pH仪测量pH。含水量用烘干恒重法测量，含氮量用凯式法测量。菌丝萌发期、满袋期通过感官观察进行记录，菌丝长速测量参考常婷婷等［19 ］的方法采用划线法测量。

1.2.4 验证实验 验证试验的制作方法同1.2.2，测量方法同1.2.3。

1.3 数据整理与分析

数据采用Microsoft Excel 2016软件进行计算、整理；使用SPSS 26.0统计软件单因素ANOVA检测对数据进行方差分析并使用Duncan法做差异显著性分析；使用Design Expert 8.0.6软件分析模型、进行拟合分析、建立回归方程模型，得出等高线图和响应曲面图，利用Optimization功能，进行多目标的优化分析。

2 结果与分析

2.1 不同混料组合配方栽培基质的理化性质

2.1.1 原种混料组合配方基质的理化性质 由表4可知，不同原种混料组合配方栽培基质的容重以配方3、8最高，均为0.310 g/cm3，均显著高于其余配方（P < 0.05）；配方1、4、6、7、9、12的容重次之，均为0.290 g/cm3，均与配方2、5、11差异显著（P < 0.05），均与配方10、13差异不显著（P > 0.05）；配方5的容重最小，仅为0.260 g/cm3，与配方2、11均差异不显著（P > 0.05），与其余配方均差异显著（P < 0.05）。总孔隙度以配方10、12、11的较高，分别为75.01%、74.23%、71.91%，与其他配方均差异显著（P < 0.05）。通气孔隙度以配方8最高，为22.33%，与配方13差异不显著（P > 0.05），与其他配方均差异显著（P < 0.05）；配方13次之，为18.78%，与配方3、6、9、10、12均差异不显著（P > 0.05），与配方1、2、4、5、7、11均差异显著（P < 0.05）；配方5的通气孔隙度最小，为9.36%，与配方1、2、4、7均差异不显著（P > 0.05），与其余配方均差异显著（P < 0.05）。持水孔隙度以配方10的最高，为60.43%，显著高于其他配方（P < 0.05）；配方11次之，为57.72%，与配方12差异不显著（P > 0.05），与其余配方均差异显著（P < 0.05）；配方8最小，为31.44%，与其余配方均差异显著（P < 0.05）。pH以配方2的最高，为7.71，均显著高于其他配方（P < 0.05）。含水量以配方1的最高，为496.8 g/kg，除与配方12的差异显著（P < 0.05）外，与其余配方均无显著性差异（P > 0.05）。含氮量以配方13、9、8的较高，分别为13.6、13.3、13.1 g/kg，均与配方3差异不显著（P > 0.05），均与其余配方差异显著（P < 0.05）；配方5含氮量最低，为3.6 g/kg，与配方2、10、11差异不显著（P > 0.05），与其余配方均差异显著（P < 0.05）。

2.1.2 栽培种混料组合配方栽培基质的理化性质

由表5可知，不同栽培种混料组合配方栽培基质的容重以配方10最大，为0.32 g/cm3，与其他配方均呈现显著性差异（P < 0.05）；配方4、11的容重次之，均为0.30 g/cm3，均与配方14、16差异不显著（P > 0.05），均与其余配方差异显著（P < 0.05）；配方2、3容重最小，均为0.25 g/cm3，均与其余配方差异显著（P < 0.05）。总孔隙度以配方11最大，为79.02%，与配方12差异不显著（P > 0.05），与其余配方均差异显著（P < 0.05）；配方12次之，为77.47%，与配方3、5、15均差异显著（P < 0.05），与其余配方均差异不显著（P > 0.05）。通气孔隙度以配方11最大，为24.19%，与配方1、2、5、6、7、8、10、12、13均差异不显著（P > 0.05），与其余配方均差异显著（P < 0.05）；配方2次之，为23.69%，与配方3、4、9、14、15均差异显著 （P < 0.05），与其余配方均差异不显著（P > 0.05）。持水孔隙度以配方4最大，为60.05%，均显著高于其他配方（P < 0.05）；配方9次之，为55.68%，与配方5、6均差异显著（P < 0.05），与其余配方均差异不显著（P > 0.05）；配方5的持水孔隙度最小，为48.46%，与配方1、2、3、6、7、8均差异不显著（P > 0.05），与其余配方均差异显著（P < 0.05）。pH以配方9最高，为7.53，与配方1、4、8、11均差异不显著（P > 0.05），与其余配方均差异显著（P < 0.05）；配方2最低，为6.71，与配方7差异不显著（P > 0.05），与其余配方均差异显著（P < 0.05）。含水量以配方7最高，为475.6 g/kg，与配方5差异不显著（P > 0.05），与其余配方均差异显著（P < 0.05）；配方5次之，为464.3 g/kg，与配方1、4、10、14、16均差异不显著（P > 0.05），与其余配方均差异显著（P < 0.05）；配方2最小，为395.7 g/kg，与配方3差异不显著（P &gt; 0.05），与其余配方均差异显著（P < 0.05）。含氮量以配方11最高，为15.0 g/kg，与配方2、16、10 均差异不显著（P > 0.05），与其余配方均差异显著（P < 0.05）；配方2、16、10较高，分别为14.3、14.3、13.0 g/kg，均与配方11差异不显著（P > 0.05），但均与其余配方差异显著（P < 0.05）；配方7最低，为3.5 g/kg，与配方3、5、6均差异不显著，与其余配方均差异显著（P < 0.05）。

2.2 预测模型及回归方程的建立

2.2.1 原种配方的预测模型及回归方程的建立 原种混料各配方基质香菇生长情况见表6，采用Design-Expert软件分别对香菇萌发期、菌丝长速、满袋期各指标进行回归拟合分析。根据Design-Expert建议原种的萌发期、菌丝长速、满袋期选用线性分析模型，得到的预测方程如表7所示。

根据表7可知，对于香菇原种的萌发期指标而言，预测方程的P < 0.01，表示选用的模型高度显著；回归系数R2=0.986 3，说明香菇原种萌发期的变异中有98.63%是由变量（A麦粒、B木屑）引起的，即表明该模型与试验数据有98.63%的符合度；校正相关系数R2Adj=0.985 1、失拟项（0.522 3） > 0.05不显著，说明该模型与实际拟合较好，具有较高的可信度。回归方程的自变量系数K，反应各自变量对方程的贡献程度，由原种萌发期方程可知，KB（6.00） > KA（2.46），表明2种培养料对萌发期的贡献程度为木屑优于麦粒（B > A）。

对于香菇原种的菌丝长速指标而言（表7），预测方程的P < 0.01，表示选用的模型高度显著；回归系数R2=0.968 8，说明香菇原种菌丝长速的变异中有96.88%是由变量（A麦粒、B木屑）引起的，即表明该模型与试验数据有96.88%的符合度；校正相关系数R2Adj=0.965 9，失拟项（0.051 5） > 0.05不显著，说明该模型与实际拟合较好，具有较高的可信度。回归方程的自变量系数K，反应各自变量对方程的贡献程度，由原种菌丝长速方程可知，KA（0.67） > KB（0.21），表明2种培养料对菌丝长速的贡献程度为麦粒优于木屑（A > B）。说明就菌丝长速而言麦粒的贡献程度高于木屑，即可以使用麦粒代替木屑作为培养香菇原种的培养料。

对于香菇原种的满袋期指标而言（表7），预测方程的P < 0.01，表示选用的模型高度显著；回归系数R2=0.960 9，说明香菇原种满袋期的变异中有96.09%是由变量（A麦粒、B木屑）引起的，即表明该模型与试验数据有96.09%的符合度；校正相关系数R2Adj=0.957 4，失拟项（0.421 3） > 0.05不显著，说明该模型与实际拟合较好，具有较高的可信度。回归方程的自变量系数K，反应各自变量对方程的贡献程度，由原种满袋期方程可知，KB（53.87） > KA（35.26），表明2种培养料对满袋期的贡献程度为木屑优于麦粒（B > A）。

2.2.2 栽培种配方的预测模型及回归方程的建立

栽培种混料各配方基质的香菇生长情况见表8，采用Design-Expert软件分别对香菇萌发期、菌丝长速、满袋期进行回归拟合分析。根据Design-Expert建议原种的萌发期、菌丝长速、满袋期选用线性分析模型，得到的预测方程如表9所示。

根据表9可知，对于香菇栽培种的萌发期指标而言，预测方程的P < 0.01，表示选用的模型高度显著；回归系数R2=0.951 5，说明香菇栽培种萌发期的变异中有95.15%是由变量（A玉米芯、B麸皮、C木屑）引起的，即表明该模型与试验数据有95.15%的符合度；校正相关系数R2Adj=0.927 3、失拟项（0.496 8） > 0.05不显著，说明该模型与实际拟合较好，具有较高的可信度。回归方程的自变量系数K，反应各自变量对方程的贡献程度，由栽培种萌发期方程可知，KB（5.89） > KC（5.84） > KA（2.94），表明3种培养料对萌发期的贡献程度由高到低依次为麸皮、木屑、玉米芯（B > C > A）。二次项系数KAB（-4.08） > KAC（-4.68） > KBC（-5.50），说明玉米芯（A）和麸皮（B）搭配对栽培种萌发期的贡献最大。

对于香菇栽培种的菌丝长速指标而言，预测方程的P < 0.01，表示选用的模型高度显著；回归系数R2=0.988 7，说明香菇栽培种菌丝长速的变异中有98.87%是由变量（A玉米芯、B麸皮、C木屑）引起的，即表明该模型与试验数据有98.87%的符合度；校正相关系数R2Adj=0.971 7、失拟项（0.869 4） > 0.05不显著，说明该模型与实际拟合较好，具有较高的可信度。回归方程的自变量系数K，反应各自变量对方程的贡献程度，由栽培种菌丝长速方程可知，KB（2.73） > KA（0.85） > KC（0.18），表明3种培养料对萌发期的贡献程度由大到小依次为麸皮、玉米芯、木屑（B > A > C）。二次项系数KAC（-0.52） > KBC（-4.29） > KAB（-4.84），说明玉米芯和木屑搭配对栽培种菌丝长速的贡献最大（表9）。

对于香菇栽培种的满袋期指标而言，预测方程的P < 0.01，表示选用的模型高度显著；回归系数R2=0.991 2，说明香菇栽培种满袋期的变异中有99.12%是由变量（A玉米芯、B麸皮、C木屑）引起的，即表明该模型与试验数据有99.12%的符合度；校正相关系数R2Adj=0.977 9、失拟项（0.366 5） > 0.05不显著，说明该模型与实际拟合较好，具有较高的可信度。回归方程的自变量系数K，反应各自变量对方程的贡献程度，由栽培种满袋期方程可知，KB（256.00） > KA（158.13） > KC（54.91），表明3种培养料对萌发期的贡献程度由大到小依次为麸皮、玉米芯、木屑（B > A > C）。二次项系数KAC（-260.19） > KBC（-457.81） > KAB（-663.17），说明玉米芯和木屑搭配对栽培种满袋期的贡献最大（表9）。

2.3 不同原料配比对香菇原种和栽培种的影响

2.3.1 不同原料配比对香菇原种的影响 原种混合料为麦粒和木屑2种，其不同配比与原种萌发期、菌丝长速、满袋期的影响曲线见图1、图2、图3。根据图1可知，随着麦粒含量的增加，随之木屑含量越小，菌丝的萌发期逐渐变短，配方3、8、13的试验基质均为麦粒50%和木屑50%，其萌发期分别为2.50、2.50、2.17 d，均短于其他配方（表6）；由图2可知，随着麦粒含量的增加，原种菌丝长速越快，配方3、8、13的菌丝长速分别为0.64 、0.67 、0.63 cm/d，均高于其他配方（表6）。由图3可知随着麦粒含量的增加，香菇原种的满袋期越短，配方3、8、13的满袋期分别为35.75 、35.58 、36.80 d，均短于其他配方（表6）。根据建立的回归方程及响应图最终预测认为香菇原种栽培基质的最佳优化配方为50%麦粒+50%木屑。

2.3.2 不同原料配比对香菇栽培种的影响 栽培种混合原料为玉米芯、麸皮、木屑，其不同配比与香菇栽培种菌丝萌发期、菌丝长速、满袋期的影响曲线见图4、图5、图6。由图4可知不同原料配比与栽培种萌发期的响应面图为曲面，故三者之间存在交互作用。配方12的萌发期最短，为2.67 d；而配方5、7的萌发期最长，为5.75 d（表8）。曲面波动较小，但受木屑的影响较大，木屑添加量越少，栽培种的萌发期就越短。由图5可知不同原料配比与栽培种菌丝长速的响应面图为曲面，故三者之间存在交互作用。配方2号的菌丝长速最快为0.57 cm/d，而配方5的栽培种的菌丝长速最慢，为0.17 cm/d（表8）。受玉米芯的影响曲面波动较大，玉米芯添加量越多，栽培种的菌丝长速越快。由图6可知不同原料配比与栽培种满袋的响应面图为曲面，故三者之间存在交互作用。配方4、8的满袋期最短，均为37.67 d；而配方5的栽满袋期最长，为55.17 d（表8）。受木屑的影响曲面波动较大，木屑添加量越少，栽培种的满袋期越短。根据建立的回归方程及响应面图最终预测认为香菇栽培种栽培基质的最佳优化配方为37.69%玉米芯+23.33%麸皮+38.98%木屑。

2.4 验证实验

为了让香菇原种和栽培种建立模型中的预测值，更加拟合实际，故再次通过优化的配方进行验证，验证实验结果表明，香菇原种的萌发期、菌丝长速、满袋期的实际值与预测值均差异较小，误差率分别为1.60%、1.49%、0.74%；香菇栽培种的萌发期、菌丝长速、满袋期的实际值与预测值差异较小，误差率分别为0.67%、2.00%、4.64%。

3 讨论与结论

优质的培养料及科学的配方是菌种生产中的关键，培养料的多元化既能改善菌丝生长所需的营养物质，又可以缓解原料单一造成的资源缺乏问题［20 - 22 ］。香菇作为木腐菌，其菌丝生长过程中将分解木材中的木质素、半纤维素、纤维素以及有机氮提供其生长［23 ］，为了解决木材资源缺乏的问题，大量的研究学者将使用木质素及纤维素含量较高的玉米芯、秸秆等代替了木材［24 ］。但很多研究人员在香菇种植过程中均凭借多年的经验来添加代替原料，未能科学客观的指出相关的原料配比，香菇培养料种类的不同，其中所含有的营养物质也会随之改变，进而导致菌种质量的优劣有所差别［25 - 27 ］。杨建杰等［28 ］采用单纯形格子试验设计方法，以木屑、棉籽壳、玉米芯和麦草为原料优化出平菇的优质栽培配方，改试验设计方法为最先出现的、最基本的一种混料回归设计，但是最后试验结果将需人工计算［29 ］。而本试验中所采用的D-optimal最优混料设计试验则与其原理类似，优点在于除了数据的测量与录入以外，其他计算均可在计算机软件中完成，大大节省了工作效率以及提高了人工计算带来的误差，从而更科学的得出试验结果［30 - 31 ］。

沈霞［32 ］采用D-optimal最优混料设计法，以芦苇、五节芒、斑茅和芒萁为原料，代替全部木屑，经过不同的配比，优化出灰树花的栽培基质配方，得出最优的配方为45%五节芒+30%芒萁+5%麦麸+18%玉米粉+2%石膏，但对不同配比的基质理化性质未作研究。本试验通过测定不同配比培养料的理化性质，更加详细深入地研究了香菇原种及栽培种的培养基质配方，经分析得出香菇原种和栽培种栽培基质的最优配方分别为50%麦粒+50%木屑和37.69%玉米芯+23.33%麸皮+ 38.98%木屑，在以上2个配方的栽培基质接种后，香菇菌丝的生长旺盛，萌发期短、满袋期短，且理化性质优，说明优化得到的栽培基质配方具有较高的可行性，该设计方法也在优化培养料配比上是科学并且可行的。由于液体菌种的发展迅速，目前国内外学者通过D-optimal最优混料设计法对菌种培养料配方优化鲜有研究，关于适宜西北地区香菇菌棒混料培养基的筛选和优化也有待进一步研究。

参考文献：

［1］ 汪文琦，冯继虎，段 辉，等. 康县香菇产业现状及发展建议［J］. 寒旱农业科学，2024，3（3）：208-211.

［1］ 戈永杰，尚晓冬，谭 琦. 通过限制培养基中氮元素含量促进香菇分泌木质素过氧化物酶［J］. 食用菌学报，2022，29（6）：25-34.

［2］ 卓 英，谭 琦，陈明杰，等. 香菇主要栽培菌株遗传多样性的AFLP分析［J］. 菌物学报，2006（2）：203-210.

［3］ 黄 毅，郑永德，张长本. 菌林矛盾成为制约香菇、黑木耳产业发展的瓶颈［J］. 食药用菌，2018，26（6）：333-337.

［4］ 李 军，蔡 婧，常 堃，等. 玉米秆玉米芯和玉米粉反季节栽培香菇试验［J］. 食用菌，2017，39（6）：37-39.

［5］ 覃中凤，鄂肖勍，何志强，等. 香菇栽培料配方优选试验［J］. 现代农业科技，2021，787（5）：56-57；62.

［6］ 王洪奇，陈 顺，于广峰，等. 玉米秸秆代料栽培香菇菌种的筛选［J］. 微生物学杂志，2016，36（6）：73-75.

［7］ ANDRADE M， SILVA J， MINHONI M， et al. Mycelial growth of two Lentinula edodes strains in culture media prepared with sawdust extracts from seven eucalyptus species and three eucalyptus clones［J］. Acta Scientiarum Agronomy， 2008， 30（3）： 333-337.

［8］ XU X F， LI A Z， CHENG S M， et al. Reappraisal of phylogenetic status and genetic diversity analysis of Asian population of Lentinula edodes［J］. Progress in Natural Science： Materials International， 2006（3）： 274-280.

［9］ 杨建杰，张桂香，杨 琴，等. 玉米芯栽培香菇的配方优化及胞外酶活性研究［J］. 中国食用菌，2021，40（4）：26-31.

［10］ 夏 敏，余明玉，杜瑞卿，等. 配方均匀设计对玉米秸秆代料栽培香菇的配方优化［J］. 南京农业大学学报，2011，34（4）：138-142.

［11］ 周 韬. 香菇液体菌种培养基优化及其菌种质量评价体系的建立［D］. 武汉：华中农业大学，2017.

［12］ 王攀磊，潘艳华，郭玉蓉，等. 基于混料均匀设计的无土栽培基质研究与应用［J］. 云南农业大学学报（自然科学），2018，33（5）：797-803.

［13］ 王永洁，孙铭忆，邓莉莉，等. 基于AHP-CRITIC的D-最优混料设计优化复方黄芪乳膏处方［J］. 中国实验方剂学杂志，2020，26（3）：110-116.

［14］ 张桂香，杨建杰，刘明军，等. 不同出菇温度下香菇各潮次菇产品的品质变化［J］. 甘肃农业科技，2020（10）：28-31.

［15］ 曲 扬. 不同菌株及栽培模式对香菇CO2释放量的影响［D］. 上海：上海海洋大学，2020.

［16］ 张清凤，刘元剑，杨静康，等. 苹果枝木屑作栽培料生产香菇效果试验［J］. 上海蔬菜，2022（5）：73-74.

［17］ 倪 琳. 不同比例腐解秸秆及生物质炭代替泥炭用作番茄育苗基质研究［D］. 长春：吉林农业大学，2021.

［18］ 王 涛，雷锦桂，黄语燕，等. D-最优混料设计优化草菇菇渣复合基质的黄瓜栽培配方［J］. 土壤，2022，

54（4）：723-732.

［19］ 常婷婷，赵 妍，杨焕玲，等. 香菇复配氨基酸培养基优化及其对香菇菌丝生长生理的影响［J］. 微生物学通报，2022，49（2）：545-555.

［20］ DELMAS J， MAMOUN M. Study of some growth and fruiting factors of Pleurotus cornucopiae［J］. Comptes Rendus Des Seances De Lacademie Dagriculture De France，1980， 66（3）： 294-301.

［21］ 段应策，胡姿仪，杨 帆，等. pH和缓冲作用对香菇菌丝生长的影响［J］. 中国农业科学，2020，53（22）：4683-4690.

［22］ 赵林超. 氮源浓度及乙酸钠对香菇生长的影响研究［D］. 南京：南京农业大学，2019.

［23］ 张金霞，蔡为明，黄晨阳，等. 中国食用菌栽培学［M］. 北京：中国农业出版社，2020.

［24］ 于 惠. 玉米秸秆预处理及用作食用菌栽培基质的研究［D］. 武汉：华中农业大学，2018.

［25］ SILVA E， MACHUCA A， MILAGRES A. Effect of cereal brans on Lentinula edodes growth and enzyme activities during cultivation on forestry waste［J］. Letters in Applied Microbiology， 2005， 40（4）： 283-288.

［26］ 季泽洋. 不同培养料对香菇菌丝生长、产量和营养成分的影响［D］. 合肥：安徽农业大学，2018.

［27］ POONAM B， RANJANA P， TEJO P. Enhanced antioxidant properties as a function of selenium uptake by edible mushrooms cultivated on selenium-accumulated waste post-harvest wheat and paddy residues［J］. International Journal of Recycling of Organic Waste in Agriculture， 2014， 3（4）： 127-132.

［28］ 杨建杰，张桂香，杨 琴，等. 玉米芯栽培平菇培养料配方优化试验［J］. 甘肃农业科技，2017，504（12）：31-34.

［29］ 熊子康，刘力惟，宁建辉，等. 混料切片试验设计［J］. 系统科学与数学，2020，40（2）：262-274.

［30］ 罗嘉成，张崇岐. 混料试验渐近D-最优设计的聚类算法［J］. 数理统计与管理，2022，41（3）：402-412.

［31］ 王浩宇，张崇岐. 基于特殊混料模型D-最优设计搜索的交换点式门限接受算法［J］. 系统科学与数学，2020，40（2）：210-224.

［32］ 沈 霞. 灰树花菌草栽培基质配方优化［J］. 农学学报，2021，11（8）：64-69.