腐食酪螨在平菇种植场的空间分布型及抽样技术分析

张萌 吕玉华 李子玲 刘斌 赖茜茹

【摘   要】 对广西南宁平菇种植场危害平菇的腐食酪螨的空间分布情况进行调研，对不同的抽样技术进行对比，以Iwao 回归线性方程、Taylor 回归线性方程以及聚集指标来确定腐食酪螨的空间分布型为聚集分布。抽样技术的对比可以得出“对角线”和“平行线式”是菇场中较为理想的抽样方式。

【关键词】 腐食酪螨;空间分布;抽样技术

腐食酪螨（Tyrophagus putrescentiae）是一种在世界各地广泛分布且个体微小的螨虫，是我国南方地区的优势粉螨种类，属于粉螨科（Acaridae）、食酪螨属（Tyrophagus）[1， 2]，常生长在蛋白质和脂肪含量高的储藏物之中，食用储藏粮食，会使得储藏物的品质下降[3， 4]，并可能造成严重经济损失[5]。腐食酪螨也是食用菌的主要害虫之一[6]，喜爱食用食用菌的孢子和子实体，通过其边缘锯齿状的螯肢剪碎食用菌子实体上的菌肉，食用子实体[7]，当螨的数量过多时，会严重影响子实体的产量和品质。近来年，随着产业结构调整和扶贫攻坚工作的开展，广西平菇种植规模迅速扩大，腐食酪螨的危害也愈发地显著。为更有效地做好防控工作，减少其对平菇生产的影响，需对腐食酪螨的生物学、生态学特性进行全面的研究，而空间分布型是其中相当重要的一个方面。关于腐食酪螨空间分布型的研究目前已有一些报道，如湛孝东等对腐食酪螨自然种群消长动态及空间分布型研究表明，腐食酪螨在仓储面粉中的空间格局为聚集分布[8];陶莉等认为，在仓储小麦中腐食酪螨是以个体群为基本成份呈聚集分布，且密度越高，聚集度越大[9]。平菇种植一般在菇棚、菇房等人工设施内进行，腐食酪螨在其间所处的环境与上述研究的仓库环境有极大的不同，但有关腐食酪螨在食用菌种植场特别是平菇种植场的空间分布型的研究尚未见有报道。为确定腐食酪螨在平菇种植场的最佳田间取样方式，同时为制定安全、有效的防控方案提供参考，本文对腐食酪螨在平菇种植场的空间分布型及抽样技术进行了初步探究。

1  研究方法

1.1  調查方法

调查于2019年3-4月腐食酪螨冬后恢复活动并繁殖增长的时段进行，调查场地位于广西南宁市三塘镇宇辉平菇种植场。该平菇种植场为钢架结构，覆盖双层遮阳网，面积约180m2。种植场内菌包成行排列于地面，每行堆叠5层菌包，行距1m。样点设于菌包行间地面，隔2行菌包设1样点行，每样点行隔1米取1个样点，每样点行设4个样点，共设21行样点行84个样点。每次调查时在各样点放1张100cm2正方形粘虫板，有胶面朝上，24h后记录粘虫板粘附的螨数，共调查3次，每次间隔约15d。

1.2  测定方法

1.2.1 聚集度指标  以每平方厘米粘虫板粘附的腐食酪螨数量作为样点的平均数m〔m=M/（N*100），其中M为腐食酪螨的总数，N为样本总数〕。采用以下5种常用聚集度指标测定腐食酪螨的空间分布型：扩散系数（C= S2/m）、负二项分布值（K=m2/（S2- m）、Cassie（1962）指标（CA=（ S2-m）/ m2）、丛生指标（I=S2/m-1）、平均拥挤度（m*=m+m/K）、聚快性指标（m*/m）。

各聚集度判断指标如下：C=1，为随机分布，C>1，为聚集分布，C<1，为均匀分布;CA=0，为随机分布，CA>0，为聚集分布，CA<0，为均匀分布;I=0，为随机分布，I<0为均匀分布，I>0，为聚集分布;K>0且K<8为聚集分布，K≥8，为随机分布，K<0为均匀分布;m*/m=1，为随机分布，m*/m<1，为均匀分布，m*/m>1，为聚集分布。

1.2.2 Taylor幂函数法则  采用Taylor发现的样本的方差（S2）与抽样平均数（m）之间的函数关系，关系式为：lgS2=lga+blgm，其中判断标准为：当lga=0，b=1时，为随机分布;当lga>0，b=1时，为聚集分布，聚集度不因种群密度的改变而变化;lga>0，b>1时，为聚集分布，聚集度随种群的密度增大而增大;当lga<0，b<1时，分布均匀。

1.2.3 Iwao 回归分析法  Lloyd（1969）方法中用m*表示平均拥挤度指标，Iwao发展了Lloyd的方法，确定了m*与m的关系为m*=α+βm，α为分布基本成分的平均拥挤度：当α=0时，基本成分是个体;当α>0时，基本成分是个体群;当α<0时，基本成分是个体。β为基本成分的空间分布型：当β=1时，为随机分布;当β>1时，为聚集分布;当β<1时，为均匀分布。

1.2.4  理论抽样数  根据理论抽样数公式（Iwao）：n=（t2/D2）[（ α+1）/m+β-1]，可以确定不同密度虫口的最佳理论抽样数，公式中n为最佳抽样数，t为一定概率下的置信水平（当P=0.95时，t=1.96），m为平均螨虫密度，D为允许误差值，一般取值为0.1～0.3，α、β为Iwao回归分析法中m*=α+βm的值[10]。

1.3  数据处理方法

试验数据均由Excel2010及SPSS25.0处理分析。

2  结果与分析

2.1  腐食酪螨在平菇种植场的田间分布型

根据3次调查数据计算的方差（S2）、扩散指数（C）、聚集度指数（CA）、丛生指标（I）、平均拥挤度（m*）和聚块性指数（m*/m）、负二项分布指数（K）等聚集度指标，详见表2-1。

由表2-1可知各次平调查各项聚集度指标均满足I>0、C>1、K<1、m*/m> 1、CA > 0，说明其腐食酪螨空间分布型属于聚集分布。此外，由于K值与聚集程度有关，K值越小说明聚集程度越高，因此还可以判断腐食酪螨属于高度的聚集分布。

2.2  Iwao的m*-m的线性回归指数分析

对表2-1的平均拥挤度m*和平均数m进行线性回归分析，结果如图2-1所示。Iwao线性回归分析表明，平均拥挤度（m*）与平均数m之间的线性关系式为m*= 3.396m + 8.9681（R2 = 0.9959，P<0.01），从线性回归式可以得出：α=3.396>0，表示菇场中腐食酪螨个体之间是存在相互吸引的，而β=8.9681>1，则表示菇场的腐食酪螨为聚集分布，这个结果和聚集指标的检验结果相一致。

2.3  Taylor的lgS2-lgm回归分析

对表2-1的方差S2和平均数m进行线性回归分析，结果如图2-2所示。Taylor线性回归分析表明，方差S2与平均数m之间的线性关系式为lgS2= 1.5137lgm + 0.9888（R2 = 0.9977，P<0.01），从线性回归式可以得出：lgα=0.9888>0且β=1.5137>1，表示菇场的腐食酪螨为聚集分布，这个结果和聚集指标的检验结果相一致。

2.4   腐食酪螨理论抽样数

将m*-m回归方程（Iwao）中α和β带入最适理论抽样数公式中，则平菇种植场中腐食酪螨最适理论抽样式为：N=t2（9.9681/m+2.396）/D2 ，取t=1.96，允许误差D为0.1、0.2和0.3，根据此求出腐食酪螨在不同虫口密度下的最适抽样数量，如表2-2。结果表明，在相同允许误差的条件下，随着腐食酪螨虫口密度的增大，所需抽样数逐渐减少;且在相同虫口密度下，允许误差越大，所需抽样数越少。

2.5  腐食酪螨不同调查取样方法比较

使用对角线式、平行线式、棋盘式、Z字形式和五点取样式这五种方式在菇场实况图上取样，在这五种抽样方法中，对角线式和平行线式这两种方法最好，效果最佳，按照这两种抽样方法进行抽样所得到的平均数在P=0.01水平上，t测验无显著差异，他们的平均相对误差率为9.48%、14.57%，因此“对角线式”和“平行线式”取样方法作为平菇种植场中腐食酪螨的较为合理的抽样方法。

3  小结与讨论

本研究调查时间集中在3-4月，为菇场腐食酪螨危害严重期，在此期间对腐食酪螨的空间分布的研究，使用了扩散指数（C）、聚集度指数（CA）、丛生指标（I）、平均拥挤度（m*）和聚块性指数（m*/m）等聚集指标，测定分析腐食酪螨在菇场的空间分布型为聚集分布，这与以往文献所报道的密度对螨类种群分布型的影响最为明显相一致[11， 12]。对Iwao 回归分析及Taylor 回归分析的判定结果也是相一致，都说明腐食酪螨在菇场中呈聚集分布，且随着密度增加其聚集程度也增大，符合负二项分布，这与之前文献所报道在仓储小麦[13]和面粉中[14]的腐食酪螨空间分布型相一致。

通过五种抽样方法的比较可知，“对角线式”和“平行线式”取样方法可以作为在平菇种植场调查腐食酪螨的有效取样方式，为后期进行有效抽样统计提供科学与理论依据。

参考文献：

[1] 于 晓， 范青海. 腐食酪螨的发生与防治[J]. 福建农业科技， 2002（06）：49-50.

[2] Hubert J.， Munzbergova Z.， Kucerova Z.， 等. Comparison of com munities of stored product mites in grain mass and grain residues     in the Czech Republic[J]. EXP APPL ACAROL， 2006，39（2）： 149- 158.

[3] 沈兆鹏. 中国储粮螨类种类及其危害[J]. 武汉食品工业学院学报， 1996（01）：44-47.

[4] 李隆术， 轩静渊， 范青海. 四川省食品螨类名录[J]. 西南农业 大学学报， 1992（01）：27-38.

[5] Erban Tomas， Rybanska Dagmar， Hubert Jan. Population Growthof the Generalist Mite Tyrophagus putrescentiae （Acari： Acaridida） Following Adaptation to High- or Low-Fat and High- orLow-Protein Diets and the Effect of Dietary Switch[J]. ENVIRONENTOMOL， 2015，44（6）：1599-1604.

[6] 何 嘉， 張 陶， 李正跃， 等. 我国食用菌害虫研究现状[J]. 中国食用菌， 2005（01）：21-24.

[7] 陆云华. 江西食用菌螨类分类检索[J]. 宜春师专学报， 1997（05）：84-86.

[8] 湛孝东， 郭 伟， 柴 强， 等. 仓储面粉中腐食酪螨自然种群消长动态及空间分布型研究[J]. 中国血吸虫病防治杂志， 2017，29       （05）：587-591.

[9] 陶 莉， 李朝品. 腐食酪螨种群消长及空间分布型研究[J]. 南京医科大学学报（自然科学版）， 2006（10）：944-947.

[10] 洪 波， 张 锋， 李英梅， 等. 食芽象甲成虫在陕北枣园的空间分布格局[J]. 植物保护， 2017，43（06）：113-117.

[11] 谢寿安， 吕淑杰， 袁 锋， 等. 我国昆虫生态学的研究与展望 [J]. 西北林学院学报， 2002（01）：51-54.

[12] 王炳英， 张飞萍. 竹小爪螨空间分布型的研究[J]. 西南林学院学报， 2003（03）：41-43.

[13] 陶 莉， 李朝品. 腐食酪螨种群消长及空间分布型研究[J]. 南京医科大学学报（自然科学版）， 2006（10）：944-947.

[14] 湛孝东， 郭 伟， 柴 强， 等. 仓储面粉中腐食酪螨自然种群消长动态及空间分布型研究[J]. 中国血吸虫病防治杂志， 2017，       29（05）：587-591.

（编辑：李丹）