芫菁科Meloidae昆虫分布广泛，是豆类作物常见害虫(李晓婷，2010；潘昭等，2011)，同时该科昆虫多具有重要药用价值，应用历史悠久(刘海波等，2019)，如大斑沟芫菁Hycleus phaleratus (Pallas, 1781)和眼斑沟芫菁Hycleus cichorii (Linnaeus, 1758)被《中华人民共和国药典》收录(任艳等，2018；国家药典委员会，2020)，其他有相关药用记载还包括绿芫菁Lytta caraganae (Pallas, 1781)、锯角豆芫菁Epicauta gorhami (Marseul, 1873)等(王雪梅，2007；王腾蛟等，2018)。芫菁科昆虫在新疆、甘肃、宁夏、河北、内蒙古等多地广泛分布(张润杨，2018；杨益春，2019)，而哈思山区是甘肃芫菁科昆虫聚集的典型代表区域(贺春贵等，2005)。该类昆虫分泌的斑蝥素及其衍生物具有抗肿瘤、增强免疫、杀虫等作用，医疗价值极高(Yan，1993；曾瑶波等，2016；颜道宇等，2017；李济森等，2019)。作为重要的药用昆虫资源，当前市场需求量日益增长，在人工养殖方面存在瓶颈，无法规模化，临床使用仍以野外收集为主，过度不合理采集使其野生数量日渐减少(陈仕江，2004；王雪梅等，2007；莫让瑜，2014)，而目前有关芫菁科昆虫分布特征的研究尚无详尽报道。

昆虫种群空间分布是昆虫种群的重要生态学属性之一(杨现明等，2020)，不仅对了解昆虫种群的爆发、扩散具有重要的意义，而且在抽样技术及监测、预报中也发挥了一定的指导作用(鲁武锋等，2018)。对芫菁科药用昆虫聚集典型代表区域进行系统调查，探索其空间分布规律和生活习性，对该类药用昆虫资源的保护和综合开发利用具有重要的作用和实际意义。因此，本文采用随机取样法对芫菁科药用昆虫的分布进行调查，利用聚集度指标法、Taylor幂函数法及Iwao回归分析法对昆虫的空间分布类型及抽样技术进行研究，以期为芫菁科药用昆虫资源的监测、保护和综合开发奠定基础，也为防治该科昆虫因数量过度造成农业危害提供参考。

1 材料与方法 1.1 研究地概况

哈思山区是以哈思山自然保护区为中心的农林牧交错带，地处甘肃省白银市靖远县境内，属陇中黄土高原农牧林区，地理坐标104°13′40″~104°35′00″E，36°58′20″~37°2′40″N，海拔2 400~3 017 m。年均气温6~7 ℃，年均降水量350~ 410 mm，降水极不均匀，多集中于7—9月(杨万军，2003；张启立等，2019)。全县现有耕地7.58万hm²，草原面积37.33万hm²，其中哈思山区草原主要为高寒草甸类、温性草原类和温性荒漠草原类(董四娟，2012)。该区域植被覆盖较好，但森林植被类型少，同时由于地处旱地农业向荒漠牧地过渡边界，荒漠化较严重，主要有石质荒漠、沙质荒漠和盐化荒漠(舒海元，1997；李连惠，2003)。

1.2 调查方法

本研究于2019年7月，当地芫菁科昆虫频繁活动时期，根据哈思山区概况及地形地貌特征，选择农田型、草地型、森林型和荒漠型4种具有代表性的生态类型样地各3块，每块样地设置5 m²(n=7) 正方形样方，样方间距约为500 m。部分采样生境如图 1所示，记录样方内芫菁科主要药用昆虫及其所分布植物的种类和数量。采用实验室自制昆虫捕捉装置进行样方内目标昆虫的捕捉，进行分类学鉴定。

芫菁科主要药用昆虫空间分布型测定取样方法：选取14个样方(农田型和草地型样地各7个样方)，采用随机抽样法取样，每个样方随机选取20株植株，记录植株上芫菁科昆虫数量，并进行空间分布型测定。

1.3 分析方法 1.3.1 空间分布型

采用聚集度指标法(王智，2007；贾彦霞等，2019)、Taylor幂函数法(Taylor，1961；杨紫涵等，2020)和Iwao回归分析法(Iwao，1968；刘延超等，2018)对芫菁科主要药用昆虫的空间分布型进行测定。

聚集度指标法：计算得到各样方昆虫的平均虫口密度(m)、样本方差(S²)、平均拥挤度m*。再按照聚集度指标的计算公式分别求出扩散系数(C)：C=S²/m、聚块性指数m*/m、负二项分布k值：k=m²/(S²-m)、丛生指数(I)：I=(S²-m)/m、聚集度指数(C_a)：C_a=1/k=(S²-m)/m²和L_a指标：L_a=m-m/S²+1。

Taylor幂函数法：根据Taylor(1961)提出的m与S²的回归关系公式：S²=lga+blgm，式中，a和b分别表示抽样因素和聚集特征。当lga=0、b=1时，种群为随机分布；当lga>0、b=1时，种群在一切密度下呈聚集分布，且不依赖于密度；当lga>0、b>1时，种群在一切密度下呈聚集分布，且依赖于密度聚集；当lga < 0、b < 1时，种群为均匀分布，且密度越高，分布越均匀。

Iwao回归分析法：m*-m回归模型为：m*=α+ βm，m*和m表示平均拥挤度与平均密度。式中，α为分布基本成分的分布性质(平均拥挤程度)，当α=0时，分布的基本成分为单个个体；当α>0时，个体间相互吸引，分布的基本成分为个体群；当α < 0时，个体间相互排斥。β为基本成分的空间分布型，当β=1时，种群为随机分布；当β＞1时，种群为聚集分布；当β < l时，种群为均匀分布。

1.3.2 聚集原因分析

应用Blackith(1958)提出的聚集均数λ对芫菁科主要药用昆虫的聚集因素进行分析。聚集均数公式为：λ=mγ/2k，式中，k为各样方昆虫负二项分布指数，m为平均值，γ为df=2k、P=0.05时的χ²值。当λ≥2时，聚集主要是昆虫的本身习性和环境条件综合影响的结果；当λ < 2时，聚集是由环境因素引起(Blackith，1958；阎雄飞等，2019)。

1.3.3 抽样技术

根据Iwao(1968)提出的理论抽样数计算公式：N=t²/D²[(α+1)/m+β-1]，计算芫菁科主要药用昆虫在不同虫口密度下的最适抽样数。式中，N为最适抽样数，D为允许误差，取t=1.96(95%可靠度)，m为平均虫口密度(张治科等，2018)。

采用Iwao(1968)提出的序贯抽样理论公式：为防治指标，t为一定置信度下的分布值，取t=1.96(置信概率为0.95)，α、β为m*-m回归模型参数。制定序贯抽样表，并以此进行防治决策(尹炯等，2018)。

1.4 数据统计分析

运用SPSS 20.0进行数据统计分析，运用Heat Mapper进行聚类热点图分析。

2 结果与分析 2.1 哈思山区芫菁科主要药用昆虫的分布 2.1.1 芫菁科主要药用昆虫种类及分布情况

对所有样地中的芫菁科药用昆虫进行鉴定，常见有西北斑芫菁Mylabris sibirica (Fischer von Waldheim, 1823)、丽斑芫菁M. speciosa (Pallas, 1781)、圆点斑芫菁M. aulica (Ménétriés, 1832)、暗头豆芫菁Epicauta obscurocephala (Reitter, 1873)、扁角豆芫菁E. impressicornis (Pic, 1913)和绿芫菁Lytta caraganae (Pallas, 1781)。调查发现，在不同类型样地中芫菁科主要药用昆虫数量及种类均存在一定差异(表 1)。该区域芫菁科昆虫分布于农田型样地和草地型样地，且农田型样地分布数量远多于草地型样地。西北斑芫菁、丽斑芫菁、圆点斑芫菁、暗头豆芫菁、扁角豆芫菁及绿芫菁在2种类型样地中均有分布，其中，西北斑芫菁在2种样地中最常见、数量最多，圆点斑芫菁次之。森林型和荒漠型样地中均未发现芫菁科昆虫分布。有研究报道，芫菁科昆虫啃食沙棘Hippophae rhamnoides、沙枣Elaeagnus angustifolia、锦鸡儿Caragana microphylla等的叶片(王文和，2015)，但本研究的沙棘林及其他森林型样地中均未发现芫菁科昆虫分布。

2.1.2 芫菁科主要药用昆虫在植物上的分布

调查发现, 芫菁科主要药用昆虫分布的植物种类较多(表 2)，共10科17种植株，集中于菊科Astera ceae、豆科Leguminosae和唇形科Lamiaceae植物。其中，西北斑芫菁分布的植物种类最丰富，共计11种，是区域内最常见的芫菁类昆虫。根据昆虫在不同植物上的分布情况，对6种芫菁类昆虫进行聚类分析，从分布热点图(图 2)可以看出，芫菁类昆虫可分为2类，西北斑芫菁和圆点斑芫菁为一类，绿芫菁、扁角豆芫菁、丽斑芫菁和暗头豆芫菁为另一类。这是由于西北斑芫菁和圆点斑芫菁对植物的选择相似，2种昆虫多分布在同一样地，均主要取食密花香薷Elsholtzia densa、紫苜蓿Medicago sativa、乳苣Mulgedium tataricum等植物。而其他4种芫菁分布相对零散，其中绿芫菁和扁角豆芫菁分布的植物种类相对较少，绿芫菁主要取食野豌豆Vicia sepium，该地区绿芫菁活动时间集中在豆类植物开花期间，是危害豆类作物生产最主要害虫；扁角豆芫菁主要分布在黄芪Astragalus membranaceus植株上，同时该地区大量种植紫苜蓿作为牧草，其在紫苜蓿地中也分布；丽斑芫菁主要分布在串铃草Phlomis mongolica、冰草Agropyron cristatum等上；暗头豆芫菁分布范围较广，在农田型和草地型样地中均有分布，主要分布在河北石头花Gypsophila tschiliensis和白花枝子花Dracocephalum heterophyllum的植株上，但总体数量较少。

2.2 芫菁科主要药用昆虫的空间分布型 2.2.1 聚集度指标测定结果

对各样方的调查数据进行聚集度指标计算(表 3)，各样方的m*>m、C>1、I>0、C_a>0、m*/m>1、L_a/m>1、k>0，各项指标表明芫菁科主要药用昆虫的空间分布型为聚集分布。

2.2.2 Taylor幂函数法则分析

将m和S²进行对数转换，分析二者之间线性关系(图 3：a)，芫菁科主要药用昆虫Taylor幂函数方程为：lgS²=0.204 2+1.103 8lgm(r=0.886)，式中lga=0.204 2>0、b=1.103 8>1，芫菁科主要药用昆虫种群在一切密度下呈聚集分布，且具密度依赖性，聚集强度随种群密度的升高而增加，与应用聚集度指标法测定的结果一致。

2.2.3 Iwao回归分析

将m和m*进行线性回归(图 3：b)，m*-m模型回归方程为：m*=0.642 5+ 1.069 7m(r=0.916)，由α=0.642 5>0、β=1.0697>1可知，芫菁科主要药用昆虫的空间分布型为聚集分布，种群分布的基本成分是个体群，个体间相互吸引，与采用聚集度指标法和Taylor幂法则分析法测定的种群空间分布型的结果一致。

2.2.4 聚集原因分析

根据昆虫的聚集原因分析(表 4)，各样方中芫菁科主要药用昆虫的λ>2，表明其聚集原因可能是由其本身的取食、产卵等习性及温度、日照等周围环境条件综合影响或其中之一引起。

2.3 抽样技术 2.3.1 最适抽样数

利用理论抽样数计算公式，取置信水平t=1.96(概率为0.95)，芫菁科主要药用昆虫理论抽样模型为：N= (6.310/m+ 0.268)/D²。当D为0.1、0.2和0.3时，m在1、2、4、6、8、10、15……50的条件下，计算得出不同密度条件下的最适抽样数(表 5)。在允许误差一致的情况下，随着平均虫口密度的增加，芫菁科主要药用昆虫的最适抽样数逐渐减少；在平均虫口密度相同的情况下，随着允许误差的增大，芫菁科主要药用昆虫的最适抽样数逐渐减少。

2.3.2 序贯抽样

根据序贯抽样理论公式，取置信水平t=1.96(概率为0.95)，得到芫菁科主要药用昆虫的序贯抽样模型为：T₁, T₀(n)=nm₀±1.96 。本研究中芫菁科主要药用昆虫防治指标参考鲁武锋等(2018)的研究结果，取m₀=10，得序贯抽样模型为：T₁, T₀=10n± 9.480 。当前尚无相关文献报道芫菁科昆虫的防治指标，本研究设定防治指标(m₀)为2、4、6、8、10、12头/株，建立不同防治指标下芫菁科主要药用昆虫的序贯抽样模型(表 6)。取抽样株数为10、20、30、40、50、60、70、80、90、100时，建立序贯抽样表(表 7)。在样地中调查时，若昆虫累计数量超过T₁，需对样地进行防治；若昆虫累计数量低于T₀，可定为不需防治样地；若介于T₁和T₀之间，则应继续调查，直至抽样调查株数达到理论抽样模型中抽样数。

3 讨论

在甘肃哈思山区不同类型样地中，调查发现的芫菁科药用昆虫共3属6种，分布于农田型和草地型样地中，西北斑芫菁、丽斑芫菁、圆点斑芫菁、暗头豆芫菁、扁角豆芫菁及绿芫菁在2种类型样地中均有分布。在调查过程中发现，该区域草地型样地中适宜芫菁科昆虫取食的植物较农田型样地少，多分布在禁牧区和保护区草地，整体分布数量远少于农田型样地。此外，有研究发现芫菁科昆虫与蝗虫、蜜蜂存在着明显的相关性(胡子渊，2020)，其幼虫主要取食蝗虫卵或蜜蜂幼虫，而农田中密集分布的作物及其他植物通常吸引了大量的蝗虫和蜜蜂。这可能是造成该区域芫菁科主要药用昆虫在农田型和草地型样地中分布数量差异的原因。

昆虫的空间分布型是由昆虫种群的生物学、生态学特性和所处生境条件所决定的空间分布方式(王小武等，2017；鲁武锋等，2018)。本研究通过聚集度指标法、Taylor幂函数法及Iwao回归法分析哈思山区芫菁科药用昆虫的分布特征，结果均表明芫菁科药用昆虫的空间分布型为聚集分布，分布的基本成分是个体群，与李亚林(2010)研究中描述芫菁科昆虫具有群聚性相符合。芫菁科昆虫为复变态昆虫，雌虫产卵量大，平均产卵达500余粒(刘洋洋，2016)。幼虫为肉食性，主要取食蝗虫卵，部分取食蜂类幼虫。成虫为植食性，喜食植物的花和嫩叶，具有集群危害的特性，调查过程中甚至发现不同种类的芫菁科昆虫会聚集在一起取食，如西北斑芫菁与圆点斑芫菁多分布在同一样地，这可能属于营养群集或“婚期”群集(杨玉霞，2007)。芫菁科药用昆虫的聚集原因究竟是本身的取食、产卵、交配等习性及温度、日照等周围环境条件综合影响，还是其中之一所引起，仍需进一步研究。在空间分布型研究的基础上，本研究对芫菁科昆虫的抽样技术进行分析；在同一允许误差内，随着平均虫口密度的增加，最适抽样数逐渐减少。不同虫口密度下的最适理论抽样数随成虫数量的增大而减小。因而在调查时，通过预估植株上芫菁科药用昆虫的虫口密度，确定抽样的数量，当虫口密度较低时，应适当增加长植株的抽样株数，反之，则可适当减少。同时还确定了芫菁科昆虫的理论抽样模型，建立了不同防治指标下芫菁科昆虫的序贯抽样模型及一定允许误差和不同虫口密度下的序贯抽样结果表。

昆虫是动物药的重要来源之一，是中药的重要组成，占有独特地位。然而，与药用植物相比，有关药用昆虫资源方面的研究相对薄弱(韦桂宁等，2019；张星贤等，2019)。由于其动态的生活史，在时间和空间上的不确定导致研究的不易性，严重制约药用昆虫资源的开发利用。芫菁科药用昆虫种类众多，形成的中药或民间药种类丰富，如中药品种斑蝥、青娘子、葛上亭长、地胆等(黄文华，2001)，但是现代药学研究集中在大斑沟芫菁和眼斑沟芫菁分泌的斑蝥素及其衍生物方面(刘亚楠，2013)。一方面可能因过度利用造成2种芫菁野生资源一定程度的破坏，另一方面却使得同类资源大量浪费。因此对于药用昆虫而言，针对其同科同属或者近缘种属大类资源进行系统研究，既能够为医药产业提供丰富的资源，也不易造成资源的过度使用和破坏。

本研究建立了芫菁科主要药用昆虫农田型样地理论抽样模型及序贯抽样模型，分析了该地区芫菁科主要药用昆虫分布特征和生活习性，研究结果为该类药用昆虫资源的监测提供参考，一方面有利于其可持续应用，一方面可为农业生产中该类昆虫数量过度造成的危害提供防治理论依据。