全 文 ：烟蚜在烤烟田分布动态的地统计学分析∗
夏鹏亮１，２　 刘映红２　 樊　 俊１　 谭　 军１∗∗
（ １湖北省烟草公司恩施州公司， 湖北恩施 ４４５０００； ２西南大学植物保护学院， 重庆 ４００７１５）
摘　 要　 烟蚜属半翅目蚜科，是烟田的重要迁飞性害虫，成蚜和若蚜聚集于烟叶幼嫩部位，吸
食植物汁液，滋生霉污病并传播烟草病毒病害，造成产量和质量损失．本文在湖北恩施地区选
择烤烟田对烟叶从移栽至中部烟叶开始采收阶段的田间有翅烟蚜和无翅烟蚜的空间分布和
发生动态进行调查，并通过地统计学方法分析其半变异函数特性，模拟其田间迁移模式．结果
表明： 恩施烟区有翅烟蚜种群动态为双峰曲线，烟叶移栽后第 ３ 周和烟叶打顶后第 ２ 周为发
生高峰，呈现随机、聚集、随机、聚集、随机 ５步过程；无翅烟蚜种群动态为单峰曲线，烟叶打顶
前为发生高峰，呈现随机、聚集、随机 ３步过程；人为因素和寄主特性对其种群密度影响较大；
空间分布模拟插值图可以清晰反映出田间烟蚜的发生动态，结合 Ｐｅａｒｓｏｎ 相关分析发现，有翅
蚜迁入高峰期，烟蚜种群数量小并高度集中，构成了田间烟蚜的虫口基数，是烟蚜防治的关键
时期．
关键词　 空间动态； 半变异函数； Ｋｒｉｇｉｎｇ插值； 相关分析
文章编号　 １００１－９３３２（２０１５）０２－０５４８－０７　 中图分类号　 Ｓ４３１．３　 文献标识码　 Ａ
Ｇｅｏｓｔａｔｉｓｔｉｃａｌ ａｎａｌｙｓｉｓ ｏｎ ｄｉｓｔｒｉｂｕｔｉｏｎ ｄｙｎａｍｉｃｓ ｏｆ Ｍｙｚｕｓ ｐｅｒｓｉｃａｅ（Ｓｕｌｚｅｒ） ｉｎ ｆｌｕｅ⁃ｃｕｒｅｄ ｔｏｂａｃ⁃
ｃｏ ｆｉｅｌｄ． ＸＩＡ Ｐｅｎｇ⁃ｌｉａｎｇ１，２， ＬＩＵ Ｙｉｎｇ⁃ｈｏｎｇ２， ＦＡＮ Ｊｕｎ１， ＴＡＮ Ｊｕｎ１ （１Ｅｎｓｈｉ Ｂｒａｎｃｈ Ｃｏｍｐａｎｙ ｏｆ Ｈｕ⁃
ｂｅｉ Ｐｒｏｖｉｎｃｉａｌ Ｔｏｂａｃｃｏ Ｃｏｒｐｏｒａｔｉｏｎ， Ｅｎｓｈｉ ４４５０００， Ｈｕｂｅｉ， Ｃｈｉｎａ； ２Ｃｏｌｌｅｇｅ ｏｆ Ｐｌａｎｔ Ｐｒｏｔｅｃｔｉｏｎ，
Ｓｏｕｔｈｗｅｓｔ Ｕｎｉｖｅｒｓｉｔｙ， Ｃｈｏｎｇｑｉｎｇ ４００７１５， Ｃｈｉｎａ） ． ⁃Ｃｈｉｎ． Ｊ． Ａｐｐｌ． Ｅｃｏｌ．， ２０１５， ２６（２）： ５４８－５５４．
Ａｂｓｔｒａｃｔ： Ｍｙｚｕｓ ｐｅｒｓｉｃａｅ ｂｅｌｏｎｇｉｎｇ ｔｏ Ａｐｈｉｄｉｄａｅ， Ｈｅｍｉｐｔｅｒａ， ｉｓ ａｎ ｉｍｐｏｒｔａｎｔ ｍｉｇｒａｔｏｒｙ ｐｅｓｔ ｉｎ
ｔｏｂａｃｃｏ ｆｉｅｌｄ． Ａｓ ｎｙｍｐｈ ａｎｄ ａｄｕｌｔ， ｉｔ ｓｕｃｋｓ ｔｈｅ ｊｕｉｃｅ， ｂｒｅｅｄｓ ｔｈｅ ｍｉｌｄｅｗ ｓｔａｉｎｓ ｄｉｓｅａｓｅ， ｓｐｒｅａｄｓ ｔｏ⁃
ｂａｃｃｏ ｖｉｒｕｓ ｄｉｓｅａｓｅｓ ａｎｄ ｃａｕｓｅｓ ｈｕｇｅ ｌｏｓｓｅｓ ｔｏ ｔｈｅ ｙｉｅｌｄ ａｎｄ ｑｕａｌｉｔｙ． Ｔｈｅ ｄｉｓｔｒｉｂｕｔｉｏｎ ｐａｔｔｅｒｎ ａｎｄ ｄｙ⁃
ｎａｍｉｃｓ ｏｆ ｗｉｎｇｅｄ ａｎｄ ｗｉｎｇｌｅｓｓ ａｐｈｉｄｓ ｉｎ ｔｈｅ ｆｉｅｌｄ ｗｅｒｅ ｉｎｖｅｓｔｉｇａｔｅｄ ｆｒｏｍ ｔｈｅ ｔｒａｎｓｐｌａｎｔｉｎｇ ｏｆ ｔｏｂａｃｃｏ
ｔｏ ｔｈｅ ｈａｒｖｅｓｔｉｎｇ ｓｔａｇｅ ｏｆ ｍｉｄ⁃ｐｌａｃｅ ｔｏｂａｃｃｏ ｌｅａｖｅｓ ｉｎ Ｅｎｓｈｉ， Ｈｕｂｅｉ． Ｔｈｅ ｓｅｍｉｖａｒｉａｂｌｅ ｆｕｎｃｔｉｏｎ
ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｗｅｒｅ ａｎａｌｙｚｅｄ ｂｙ ｇｅｏｓｔａｔｉｓｔｉｃａｌ ｍｅｔｈｏｄ， ａｎｄ ｔｈｅ ｆｉｅｌｄ ｍｉｇｒａｔｉｏｎ ｐａｔｔｅｒｎ ｗｅｒｅ ｓｉｍｕｌａ⁃
ｔｅｄ． Ｔｈｅ ｒｅｓｕｌｔｓ ｓｈｏｗｅｄ ｔｈａｔ ｔｈｅ ｐｏｐｕｌａｔｉｏｎ ｄｙｎａｍｉｃｓ ｏｆ ｗｉｎｇｅｄ ａｐｈｉｄｓ ｉｎ Ｅｎｓｈｉ ｗｅｒｅ ｏｆ ｂｉｍｏｄａｌ
ｃｕｒｖｅ， ｗｉｔｈ ｔｗｏ ｐｅａｋｓ ａｔ ３ ｗｅｅｋｓ ａｆｔｅｒ ｔｒａｎｓｐｌａｎｔｉｎｇ ａｎｄ ２ ｗｅｅｋｓ ａｆｔｅｒ ｍｕｌｔｉ⁃ｔｏｐｐｉｎｇ ｏｆ ｔｏｂａｃｃｏ ｌｅａｖ⁃
ｅｓ， ａｎｄ ｔｈｅｒｅ ｗｅｒｅ ｆｉｖｅ⁃ｓｔｅｐ ｐｒｏｃｅｓｓ ｓｕｃｈ ａｓ ｒａｎｄｏｍ， ａｇｇｒｅｇａｔｉｏｎ， ｒａｎｄｏｍ， ａｇｇｒｅｇａｔｉｏｎ ａｎｄ ｒａｎ⁃
ｄｏｍ． Ｔｈｅ ｐｏｐｕｌａｔｉｏｎ ｄｙｎａｍｉｃｓ ｏｆ ｗｉｎｇｌｅｓｓ ｐｅａｃｈ ａｐｈｉｄｓ ｗｅｒｅ ｏｆ ｓｉｎｇｌｅ⁃ｐｅａｋ ｃｕｒｖｅ， ｇｅｔｔｉｎｇ ｉｔｓ ｐｅａｋ
ｂｅｆｏｒｅ ｍｕｌｔｉ⁃ｔｏｐｐｉｎｇ， ａｎｄ ｈａｄ ｒａｎｄｏｍ， ａｇｇｒｅｇａｔｉｏｎ， ｒａｎｄｏｍ ｔｈｒｅｅ⁃ｓｔｅｐ ｐｒｏｃｅｓｓ． Ｈｕｍａｎ ｆａｃｔｏｒｓ ａｎｄ
ｔｈｅ ｈｏｓｔｓ ｈａｄ ｃｏｎｓｉｄｅｒａｂｌｅ ｅｆｆｅｃｔｓ ｏｎ ｔｈｅ ｐｏｐｕｌａｔｉｏｎ ｄｅｎｓｉｔｙ． Ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ｓｉｍｕｌａｔｉｏｎ⁃ｉｎｔｅｒｐｏ⁃
ｌａｔｉｏｎ⁃ｆｉｇｕｒｅ ｃｏｕｌｄ ｃｌｅａｒｌｙ ｒｅｆｌｅｃｔ ｔｈｅ ｄｙｎａｍｉｃｓ ｏｆ ｔｏｂａｃｃｏ ａｐｈｉｄｓ． Ｃｏｍｂｉｎｅｄ ｗｉｔｈ ｔｈｅ Ｐｅａｒｓｏｎ ｃｏｒｒｅ⁃
ｌａｔｉｏｎ ａｎａｌｙｓｉｓ， ｗｅ ｆｏｕｎｄ ｔｈａｔ ｔｈｅ ｐｏｐｕｌａｔｉｏｎ ｄｅｎｓｉｔｙ ｗａｓ ｌｏｗ ａｎｄ ｈｉｇｈｌｙ ｃｏｎｃｅｎｔｒａｔｅｄ ａｓ ｗｉｎｇｅｄ ｔｙｐｅ
ｉｎ ｔｈｅ ｉｍｍｉｇｒａｔｉｏｎ ｐｅｒｉｏｄ， ｗｈｉｃｈ ｗａｓ ｔｈｅ ｋｅｙ ｐｅｒｉｏｄ ｆｏｒ ｔｈｅ ｍａｎａｇｅｍｅｎｔ ｏｆ ｐｅａｃｈ ａｐｈｉｄｓ．
Ｋｅｙ ｗｏｒｄｓ： ｓｐａｔｉａｌ ｄｙｎａｍｉｃ； ｓｅｍｉｖａｒｉａｂｌｅ ｆｕｎｃｔｉｏｎ； Ｋｒｉｇｉｎｇ ｉｎｔｅｒｐｏｌａｔｉｏｎ； ｃｏｒｒｅｌａｔｉｏｎ ａｎａｌｙｓｉｓ．
∗湖北省烟草公司重点项目（０２７Ｙ２０１３⁃００５）和中国烟草总公司科
技重点项目（１１０２０１２０２０１４）资助．
∗∗通讯作者． Ｅ⁃ｍａｉｌ： １１３９５６００５４＠ ｑｑ．ｃｏｍ
２０１４⁃０５⁃３０收稿，２０１４⁃１１⁃２６接受．
　 　 烟蚜（Ｍｙｚｕｓ ｐｅｒｓｉｃａｅ）属半翅目蚜科，是烟田重
要的迁飞性害虫，广泛分布于全国各大烟区，是烤烟
上最为重要的害虫之一．春季烟叶移栽后迁入烟田，
成蚜和若蚜聚集于烟叶幼嫩部位，吸食植物汁液，滋
生霉污病并传播烟草病毒病害，造成产量和质量损
失［１］ ．昆虫种群的空间分布包括水平空间分布与垂
直空间分布两个方面， 目前对水平空间分布的研究
较多［２－３］ ．黄顶成等［４］认为，烟蚜种群时间动态可分
为发展、剧增、大发生和衰退 ４ 个阶段．李厥鲁［５］将
应 用 生 态 学 报　 ２０１５年 ２月　 第 ２６卷　 第 ２期　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　
Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ， Ｆｅｂ． ２０１５， ２６（２）： ５４８－５５４
我国不同生态烟区烟蚜种群数量消长型划分为单峰
型和双峰型两种类型．邓建华等［６］利用经典统计学
方法研究表明，烟蚜分布型为聚集分布，聚集度对密
度具有依赖性．经典统计学方法将调查对象作为独
立变量处理［７－８］，而昆虫种群间总是存在多样化的
内在联系．地统计学是在地质分析和统计分析相结
合的基础上形成的分析空间相关变量的理论和方
法，承认空间相关性的存在，基于半方差函数，研究
空间分布中既有随机性又有结构性的自然现象，可
最大限度利用稀疏的或无规律的空间数据［９－１４］ ．本
文通过地统计学的方法模拟烟蚜有翅、无翅蚜的田
间分布格局和发生动态，以期为该害虫的发生规律、
抽样检测和综合防治提供科学依据．
１　 研究地区与研究方法
１􀆰 １　 研究区概况
调查田设于湖北省恩施州咸丰县丁寨乡双拱桥
村（２９°３７′２５″ Ｎ，１０９°２７′３６″ Ｅ，海拔 ８１５ ｍ），烟叶品
种为烤烟型云烟 ８７，垄沟东西向排列．烟叶生长全过
程不施用任何农药，施肥、灌溉等按常规管理，调查
田四周均为同品种烟叶［１５－１７］ ．
１􀆰 ２　 调查方法
将 ３０ ｍ×３６ ｍ 的调查田（５０ 株×３０ 行，其中株
距 ０．６ ｍ，行距 １．２ ｍ），按栅格式分成 １００个样方，每
个样方 １５株（５株×３ 行）；在每个样方中央处设 ２０
ｃｍ×１０ ｃｍ黄色粘虫板一张（高度与当时最大烟叶
等高），共计 １００ 张，每周更换黄色粘虫板，用以检
测有翅蚜的发生动态，以每个样点 １５株烟叶上无翅
蚜的平均值代表每株无翅蚜带蚜量．调查时间为
２０１２年烟叶移栽至中部烟叶开始采收，每周 １ 次，
共 １７次，分别对应移栽后第 ０ ～ １６ 周（期间烟叶共
打顶 ３次） ［１８］ ．
１􀆰 ３　 半变异函数分析
半变异函数是指分割距离为 ｈ 的两位点 ｘ 和
ｘ＋ｈ处的调查值 Ｚ（ｘ）和 Ｚ（ ｘ＋ｈ）之间的变异，用其
增量［Ｚ（ｘ）－Ｚ（ｘ＋ｈ）］平方的数学期望表示．其计算
公式为［１９］：
γ（ｈ） ＝ １
２Ｎ（ｈ）∑
Ｎ（ｈ）
ｉ ＝ １
［Ｚ（ｘｉ） － Ｚ（ｘｉ ＋ ｈ）］ ２
式中：γ（ｈ）为半方差函数；Ｎ（ｈ）为间隔距离为 ｈ 的
样本对的数目；ｈ 是一个有方向特性和长度特性的
向量；Ｚ（ｘｉ）指在点 ｘｉ处的调查值；而 ｘｉ指样本在二
维空间中存在的位置参数，是向量点［２０－２３］ ．常见的
半方差函数模型有常数模型、球形模型、指数模型、
高斯模型等，半方差函数模型有 ３个最重要的参数：
块金常数、基台值和变程，块金常数是指变异函数曲
线被延伸到间隔距离为零时的截距，它反映了区域
化变量内部随机性的可能程度；基台值是指达到平
衡时的变异函数值，它反映了变量变化幅度的大小；
变程是当变异函数的值达到平衡时的间隔距离，它
反映了区域化变量影响范围的大小．Ｃ ／ （Ｃ０＋Ｃ）的大
小判定系统内变量的空间相关程度，＜０．２５ 时，其空
间相关性较弱，结构性变异较小，随机性变异较大；
介于 ０． ２５ ～ ０􀆰 ７５，则有中等程度的空间相关性；
＞０􀆰 ７５，则具有较强的空间相关性，其结构性变异较
大，随机性变异较小［２４－２７］ ．
１􀆰 ４　 空间数据插值
Ｋｒｉｇｉｎｇ插值是一种用来估计观测样点间内插
值的地统计学方法，在已知某变量的半变异函数模
拟模型条件下，可以利用样点观测值对研究区域上
未取样点的区域化变量值进行最小误差估值，任意
一点的估计值是通过该点影响范围内的 ｎ个有效样
本值 Ｚ（ｘｉ）的线性组合得到，即：
Ｚ∗Ｖ ＝∑
ｎ
ｉ ＝ １
λ ｉＺ（ｘｉ）
式中：λ ｉ是与样点观察值 Ｚ（ｘｉ）有关的加权系数，用
来表示各样点值 Ｚ（ｘｉ）对估计值 Ｚ∗Ｖ 的贡献［２８
－３０］ ．
本文基础数据统计通过 Ｅｘｃｅｌ ２００３ 完成，半变
异函数分析和克里格插值模拟通过 ＤＰＳ ７．０５ 软件
完成．
２　 结果与分析
２􀆰 １　 烟蚜田间分布动态基本特征
田间有翅蚜的发生为双峰曲线，烟叶移栽后第
３周和烟叶打顶后第 ２ 周为发生高峰．移栽后第 ２
周，第一次发现有翅蚜，密度仅为每黄板 ０．４头；第 ３
周有翅蚜上升为每黄板 ５．９ 头，为烟叶生长季第一
个迁移高峰；第 ４ 周有翅蚜数量下降为每黄板 ３．２
头；之后一个月，数量基本稳定，田间一直存在一定
数量的有翅蚜．
烟叶第一次打顶后第 ８天和第 １５天，有翅蚜上
升为每黄板 １０．３ 和 ２９．４ 头，为第二次高峰；此后叶
片颜色变深、厚度增加、蜡质层加厚，不再适宜烟蚜
生长，有翅蚜大量迁移，数量迅速下降至每黄板 １ 头
以下．
田间无翅蚜的发生为单峰曲线，第 ８ 周为发生
高峰．无翅蚜发生晚于有翅蚜 １周，平均每株 ０．８头，
此后数量急剧增加，中间出现过多次降雨均未明显
９４５２期　 　 　 　 　 　 　 　 　 　 　 　 　 　 夏鹏亮等： 烟蚜在烤烟田分布动态的地统计学分析　 　 　 　 　 　
图 １　 有翅蚜和无翅蚜的种群动态曲线
Ｆｉｇ．１　 Ｐｏｐｕｌａｔｉｏｎ ｄｙｎａｍｉｃｓ ｏｆ ｗｉｎｇｅｄ ａｎｄ ｗｉｎｇｌｅｓｓ ａｐｈｉｄ．
Ｙ： 有翅蚜 Ｗｉｎｇｅｄ ａｐｈｉｄ； Ｗ： 无翅蚜 Ｗｉｎｇｌｅｓｓ ａｐｈｉｄ． 下同 Ｔｈｅ ｓａｍｅ
ｂｅｌｏｗ．
影响其增长速度，于第 ８周达到高峰，平均每株虫量
为 ２８６．１头．田间烟叶第 １ 次打顶前，无翅烟蚜主要
分布于顶部幼嫩叶片，打顶后烟蚜数量急剧下降，第
９周和第 １０ 周虫量分别为每株 １０３．２ 和 １２８．３ 头，
种群恢复速度明显低于打顶前．烟叶第 ２ 次打顶，即
第 １１周和第 １２周，平均每株虫量为 １８．３和 ８．２头，
此时烟叶逐渐不适宜烟蚜取食，种群已经失去自然
恢复能力；第 ３次打顶后，种群平均数量低于每株 １
头．有翅、无翅烟蚜田间动态曲线见图 １．
２􀆰 ２　 烟蚜田间分布的空间结构特点和模型参数
通过半变异函数拟合，得到有翅蚜和无翅蚜不
同时期的最佳空间分布模型．结果表明，烟蚜密度较
大时（有翅蚜密度大于每板３．６头，无翅蚜密度大于
每株 ４．６头），多为球形分布；密度较小时（有翅蚜密
度小于每板 ３．３ 头，无翅蚜密度小于每株 ０．８ 头），
多为随机分布．总体上，无翅蚜呈现随机、聚集、再随
机 ３步过程，有翅蚜呈现随机、聚集、随机、聚集、再
随机 ５步过程．烟苗移栽后第 ２ 周有翅蚜开始迁移
至烟田，迁移过程中，有翅蚜随机降落，呈随机分布．
第 ３周为迁移第 １次高峰，大量有翅蚜迁移，降落后
选择适宜生存场所，在宜居地点形成聚集区域，成聚
集分布．无翅蚜在发生初期，受有翅蚜影响，表现为
随机分布，１周后即成为聚集分布，新生代无翅蚜繁
殖能力极强，迅速形成聚集种群，并长期维持聚集状
态．第 ４周后，有翅蚜主要来源于无翅蚜的种群分
化，有翅蚜数量较少且随机产生，呈随机分布状态．
烟叶第 １次打顶后，烟叶物理化学特性不再适宜烟
蚜取食，聚集分布的无翅蚜分化产生大量有翅蚜，呈
聚集分布状态；第 １３周以后，经过 ３次打顶，田间无
翅蚜数量急剧下降，无翅蚜随机分布于残留烟杈和
腋芽上，由此产生的有翅蚜亦呈随机分布状态．总体
而言，块金值、基台值随烟蚜种群密度的增大而增
大，变异幅度随着密度增大而增加；变程为 １２ ～ ５６
ｍ，在此范围内，烟蚜数量之间存在相关关系，其变
程随虫口密度增大而增大，表明空间依赖距离随密
度增大而增大；随机程度与密度呈负相关，说明密
度增大，由空间自相关引起的空间相关性增大，由
随机性引起的空间相关性减小；决定系数均较大
（≥０．７５），说明半变异函数拟合效果良好．其结构特
点和模型参数见表 １和表 ２．
表 １　 有翅蚜田间分布参数
Ｔａｂｌｅ １　 Ｄｉｓｔｒｉｂｕｔｉｏｎ ｐａｒａｍｅｔｅｒｓ ｏｆ ｗｉｎｇｅｄ ａｐｈｉｄ
周
Ｗｅｅｋ
密度
Ｄｅｎｓｉｔｙ
（ ｉｎｄ·ｐｌａｔｅ－１）
块金值
Ｎｕｇｇｅｔ
基台值
Ｓｔｉｌｌ
随机程度
Ｒａｎｄｏｍ ｄｅｇｒｅｅ
变程
Ｒａｎｇｅ
（ｍ）
决定系数
Ｒ２
模型
Ｍｏｄｅｌ
分布
Ｄｉｓｔｒｉｂｕｔｉｏｎ
０ ０
１ ０
２ ０．３５ － － － － － 随机 Ｒａｎｄｏｍ 随机 Ｒａｎｄｏｍ
３ ５．８６ ７．３８ ３６．３９ ０．２０ １２．９６ ０．９７ 球形 Ｓｐｈｅｒｅ 聚集 Ｃｌｕｍｐｅｄ
４ ３．１５ － － － － － 随机 Ｒａｎｄｏｍ 随机 Ｒａｎｄｏｍ
５ ２．２２ － － － － － 随机 Ｒａｎｄｏｍ 随机 Ｒａｎｄｏｍ
６ １．２８ － － － － － 随机 Ｒａｎｄｏｍ 随机 Ｒａｎｄｏｍ
７ ３．２５ － － － － － 随机 Ｒａｎｄｏｍ 随机 Ｒａｎｄｏｍ
８ ２．１１ ３．１２ ３．８５ ０．８１ ２５．３６ ０．７５ 球形 Ｓｐｈｅｒｅ 聚集 Ｃｌｕｍｐｅｄ
９ ３．６１ ５．７６ ９．２９ ０．６２ ２９．３１ ０．８６ 球形 Ｓｐｈｅｒｅ 聚集 Ｃｌｕｍｐｅｄ
１０ １０．３１ １５．３６ ３７．４６ ０．４１ ２５．３６ ０．８６ 球形 Ｓｐｈｅｒｅ 聚集 Ｃｌｕｍｐｅｄ
１１ ２９．３９ ５１．３２ １７６．９７ ０．２９ ５６．３１ ０．９５ 球形 Ｓｐｈｅｒｅ 聚集 Ｃｌｕｍｐｅｄ
１２ ５．６０ １１．２６ ２６．１９ ０．４３ ４４．８６ ０．７９ 球形 Ｓｐｈｅｒｅ 聚集 Ｃｌｕｍｐｅｄ
１３ １．８６ ３．４１ ６．２０ ０．５５ １２．３９ ０．７５ 高斯 Ｇａｕｓｓ 聚集 Ｃｌｕｍｐｅｄ
１４ ０．１２ － － － － － 随机 Ｒａｎｄｏｍ 随机 Ｒａｎｄｏｍ
１５ ０．０３ － － － － － 随机 Ｒａｎｄｏｍ 随机 Ｒａｎｄｏｍ
１６ ０．０２ － － － － － 随机 Ｒａｎｄｏｍ 随机 Ｒａｎｄｏｍ
０５５ 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 应　 用　 生　 态　 学　 报　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 ２６卷
表 ２　 无翅蚜田间分布参数
Ｔａｂｌｅ ２　 Ｄｉｓｔｒｉｂｕｔｉｏｎ ｐａｒａｍｅｔｅｒｓ ｏｆ ｗｉｎｇｌｅｓｓ ａｐｈｉｄ
周
Ｗｅｅｋ
密度
Ｄｅｎｓｉｔｙ
（ ｉｎｄ·ｐａｔｅ－１）
块金值
Ｎｕｇｇｅｔ
基台值
Ｓｔｉｌｌ
随机程度
Ｒａｎｄｏｍ ｄｅｇｒｅｅ
变程
Ｒａｎｇｅ
（ｍ）
决定系数
Ｒ２
模型
Ｍｏｄｅｌ
分布
Ｄｉｓｔｒｉｂｕｔｉｏｎ
０ ０
１ ０
２ ０
３ ０．７５ － － － － － 随机 Ｒａｎｄｏｍ 随机 Ｒａｎｄｏｍ
４ ４．５６ ３．４１ ３．８８ ０．８８ １３．２３ ０．８７ 高斯 Ｇａｕｓｓ 聚集 Ｃｌｕｍｐｅｄ
５ １７．５３ １９．２０ ２５．９５ ０．７４ １２．８５ ０．９３ 高斯 Ｇａｕｓｓ 聚集 Ｃｌｕｍｐｅｄ
６ ４１．３６ １５６．３５ ２９５．００ ０．５３ １８．５３ ０．７８ 球形 Ｓｐｈｅｒｅ 聚集 Ｃｌｕｍｐｅｄ
７ １０２．２２ ７５４．８８ ２６０３．０３ ０．２９ １５．３６ ０．８５ 球形 Ｓｐｈｅｒｅ 聚集 Ｃｌｕｍｐｅｄ
８ ２８６．１４ ２６０９．７１ １２４２７．１９ ０．２１ ３１．３６ ０．９４ 球形 Ｓｐｈｅｒｅ 聚集 Ｃｌｕｍｐｅｄ
９ １０３．１６ １１８９．３６ ３０２５．３６ ０．３９ ３４．９９ ０．７７ 球形 Ｓｐｈｅｒｅ 聚集 Ｃｌｕｍｐｅｄ
１０ １２８．３３ １５３２．２１ ３６０９．８８ ０．４２ ２９．３３ ０．８８ 球形 Ｓｐｈｅｒｅ 聚集 Ｃｌｕｍｐｅｄ
１１ １８．３３ ５５．３９ ９３．８８ ０．５９ ４４．３９ ０．７９ 球形 Ｓｐｈｅｒｅ 聚集 Ｃｌｕｍｐｅｄ
１２ ８．１５ ９．８８ １２．８３ ０．７７ ３９．５４ ０．８１ 球形 Ｓｐｈｅｒｅ 聚集 Ｃｌｕｍｐｅｄ
１３ ０．８４ － － － － － 随机 Ｒａｎｄｏｍ 随机 Ｒａｎｄｏｍ
１４ ０．５５ － － － － － 随机 Ｒａｎｄｏｍ 随机 Ｒａｎｄｏｍ
１５ ０．３２ － － － － － 随机 Ｒａｎｄｏｍ 随机 Ｒａｎｄｏｍ
１６ ０．１８ － － － － － 随机 Ｒａｎｄｏｍ 随机 Ｒａｎｄｏｍ
２􀆰 ３　 烟蚜田间分布动态模拟
利用克里格插值法，根据最佳模型方程及参数，
模拟有翅蚜和无翅蚜在不同时期的田间分布格局，
如图 ２所示．
图 ２　 不同时间烟蚜空间格局的 Ｋｒｉｇｉｎｇ插值分布
Ｆｉｇ．２　 Ｋｒｉｇｉｎｇ ｄｉｓｔｒｉｂｕｔｉｏｎｓ ｏｆ ｓｐａｔｉａｌ ｐａｔｔｅｒｎ ｏｆ ｐｅａｃｈ ａｐｈｉｄ ａｔ ｄｉｆｆｅｒｅｎｔ ｓｔａｇｅｓ．
数字代表调查时间 ０～１６周 Ｔｈｅ ｎｕｍｂｅｒ ｓｔｏｏｄ ｂｙ ｔｈｅ ｉｎｖｅｓｔｉｇａｔｉｎｇ ｔｉｍｅ ａｓ ｔｈｅ ０－１６ｔｈ ｗｅｅｋ．
１５５２期　 　 　 　 　 　 　 　 　 　 　 　 　 　 夏鹏亮等： 烟蚜在烤烟田分布动态的地统计学分析　 　 　 　 　 　
表 ３　 有翅蚜和无翅蚜不同调查时期空间结构的相关性
Ｔａｂｌｅ ３　 Ｃｏｒｒｅｌａｔｉｏｎ ｏｆ ｓｐａｔｉａｌ ｓｔｒｕｃｔｕｒｅ ｂｅｔｗｅｅｎ ｗｉｎｇｅｄ ａｐｈｉｄ ａｎｄ ｗｉｎｇｌｅｓｓ ａｐｈｉｄ ｉｎ ｄｉｆｆｅｒｅｎｔ ｓｔａｇｅｓ
ｘ２ ｘ３ ｘ４ ｘ５ ｘ６ ｘ７ ｘ８ ｘ９ ｘ１０ ｘ１１ ｘ１２ ｘ１３ ｘ１４ ｘ１５ ｘ１６ ｘ１７ ｘ１８ ｘ１９ ｘ２０ ｘ２１ ｘ２２ ｘ２３ ｘ２４ ｘ２５ ｘ２６ ｘ２７ ｘ２８ ｘ２９
ｘ１ ∗∗ ∗∗ ∗∗ ∗∗ ∗∗ ∗∗ ∗ ∗∗
ｘ２ ∗∗ ∗∗ ∗∗ ∗ ∗∗ ∗∗ ∗
ｘ３ ∗∗ ∗∗ ∗∗ ∗∗ ∗ ∗∗ ∗
ｘ４ ∗∗ ∗∗ ∗∗ ∗
ｘ５ ∗ ∗∗ ∗∗ ∗∗
ｘ６ ∗∗ ∗∗ ∗∗ ∗∗ ∗∗ ∗∗ ∗
ｘ７ ∗∗ ∗∗ ∗∗ ∗∗ ∗∗ ∗∗
ｘ８ ∗
ｘ９ ∗∗ ∗∗ ∗∗ ∗∗ ∗∗
ｘ１０ ∗∗ ∗∗ ∗∗
ｘ１１ ∗∗ ∗∗ ∗
ｘ１２ ∗∗ ∗ ∗∗
ｘ１３ ∗∗ ∗
ｘ１４ ∗ ∗∗ ∗ ∗∗
ｘ１５ ∗ ∗∗
ｘ１６ ∗∗ ∗
ｘ１７ ∗∗ ∗ ∗
ｘ１８ ∗∗
ｘ１９
ｘ２０
ｘ２１
ｘ２２ ∗
ｘ２３
ｘ２４
ｘ２５
ｘ２６
ｘ２７ ∗∗
ｘ２８
∗Ｐ＜０．０５； ∗∗ Ｐ＜０．０１． ｘ１ ～ ｘ１４为无翅蚜第 ３～１６周共 １４次调查结果，ｘ１５ ～ ｘ２９为有翅蚜第 ２～ １６ 周共 １５ 次调查结果 Ｔｈｅ ｘ１ －ｘ１４ ｍｅａｎｔ ｔｈｅ １４
ｓｕｒｖｅｙ ｒｅｓｕｌｔｓ ｏｆ ｗｉｎｇｌｅｓｓ ａｐｈｉｄ ｉｎ ｔｈｅ ３－１６ｔｈ ｗｅｅｋ， ａｎｄ ｘ１５－ｘ２９ ｍｅａｎｔ ｔｈｅ １５ ｓｕｒｖｅｙ ｒｅｓｕｌｔｓ ｏｆ ｗｉｎｇｅｄ ａｐｈｉｄ ｉｎ ｔｈｅ ２－１６ｔｈ ｗｅｅｋ．
　 　 有翅蚜在烟苗移栽后第 ３周出现第一次数量高
峰，在迁移过程中有翅蚜为寻找宜居场所，受各种生
态环境因素影响，主要分布在田块西半部；受其影
响，在第 ４周无翅蚜也主要分布于田块西半部，并逐
步向东半部扩散，至第 ８周扩散至整个烟田，并长期
维持．第 ８周后的有翅蚜主要来源于无翅蚜的种群
随机分化，也分散于整个田块．
２􀆰 ４　 不同时期烟蚜田间分布的相关性
用 Ｐｅａｒｓｏｎ相关法分析不同时期烟蚜田间分布
的相关性（表 ３）．结果表明，有翅蚜第 ２～３周与无翅
蚜第 ３～７周空间分布显著相关，表明越冬代有翅蚜
影响了无翅蚜的前期田间分布；无翅蚜在第 ３ ～ ７ 周
５次调查的空间分布均显著相关；第 ８～１３周 ６次调
查空间分布也显著相关．结果与图 ２ 空间模拟图完
全吻合，前 ５次调查烟蚜主要分布于西半部，后 ６ 次
调查无翅烟蚜扩散至整个田块．最后两次调查时，烟
叶采收过半，叶片肥厚不适合烟蚜生长，仅少量存在
于新生腋芽上，加之烟叶采收期，频繁人为清除，烟
蚜种群结构极不稳定，与其他时期空间结构无显著
相关性．
３　 讨　 　 论
烟蚜的空间格局参数，受海拔、气候、地理位置
影响较大，不同的调查时间和地点，其空间参数有所
不同．有翅蚜负责田间迁移，保障种群的延续，无翅
蚜通过大量繁殖，保障种群基数的扩大．有翅、无翅
蚜配合，保障烟蚜在多变的自然、人文条件下繁衍不
息．整个烟叶生育期，无翅蚜分布曲线为单峰型，不
同于贵州等烟区的发生特点．恩施烟区烟叶旺长期
的气候、寄主环境适合烟蚜生长，无翅蚜数量没有出
现低谷．无翅蚜的发生受人为因素和寄主因素影响
较大，打顶过程带走了大量烟蚜，同时打顶后烟叶理
化性质不适宜烟蚜取食，是烟蚜迁出烟田的主要原
因．空间模拟插值图可以清晰反映出田间昆虫的发
生动态，与 Ｐｅａｒｓｏｎ相关分析结合，可以分析空间结
构的变化规律．
２５５ 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 应　 用　 生　 态　 学　 报　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 ２６卷
黄顶成等［４］认为，烟蚜种群时间动态可分为发
展、剧增、大发生和衰退 ４ 个阶段，与本文无翅烟蚜
的发生特征相吻合．李厥鲁［５］将我国不同生态烟区
烟蚜种群数量消长型划分为单峰型和双峰型两种类
型．本研究发现，恩施烟区属于单峰型．邓建华等［６］
利用经典统计学方法研究表明，烟蚜分布型为聚集
分布，聚集度对密度具有依赖性，本文通过地统计学
的方法，进一步证明了聚集度对密度的依赖性．
有翅蚜迁入高峰，种群数量小并高度集中于部
分地区，构成了田间烟蚜的虫口基数，是烟蚜防治的
关键时期．有翅蚜迁出后，进入菜园和桃园继续为
害，成为翌年烟田烟蚜的基本来源，烟叶采收后，对
周边菜园和桃园进行统一灭蚜，有助于减少翌年烟
蚜种群基数．地统计学方法承认空间相关性的存在，
补充了经典统计学方法中所有样本为独立变量的不
足，可有效用于昆虫种群分布规律研究．
参考文献
［１］　 Ｚｈｕ Ｘ⁃Ｃ （朱贤朝）， Ｗａｎｇ Ｙ⁃Ｔ （王彦亭）， Ｗａｎｇ Ｚ⁃Ｆ
（王智发）． Ｔｈｅ Ｍａｎｕａｌ ｏｆ Ｃｈｉｎｅｓｅ Ｔｏｂａｃｃｏ Ｄｉｓｅａｓｅｓ
ａｎｄ Ｉｎｓｅｃｔｓ． Ｂｅｉｊｉｎｇ： Ｃｈｉｎｅｓｅ Ａｇｒｉｃｕｌｔｕｒｅ Ｐｒｅｓｓ， ２００２：
１０４－１０５ （ｉｎ Ｃｈｉｎｅｓｅ）
［２］　 Ｚｈａｎｇ Ｙ⁃Ｃ （张运慈）， Ｃｈｅｎ Ｘ （陈 　 新）， Ｈｅ Ｚ⁃Ｌ
（贺钟麟）， ｅｔ ａｌ． Ｔｈｅ ｒｅｓｅａｒｃｈ ｏｆ ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ
ａｎｄ ｓａｍｐｌｉｎｇ ｔｅｃｈｎｏｌｏｇｙ ｏｆ Ｍｙｚｕｓ ｐｅｒｃｉｃａｅ． Ｃｈｉｎａ Ｔｏｂａｃ⁃
ｃｏ （中国烟草）， １９８８， ４（１）： １－３ （ｉｎ Ｃｈｉｎｅｓｅ）
［３］　 Ｃｈｅｎ Ｙ⁃Ｎ （陈永年）， Ｗｅｎ Ｌ⁃Ｚ （文礼章）， Ｐａｎ Ｔ （潘
桐）． Ｔｈｅ ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｐｅａｃｈ ａｐｈｉｄ Ｍｙｚｕｓ ｐｅｒｓｉ⁃
ｃａｅ ａｎｄ ｉｔｓ ｓａｍｐｌｉｎｇ ｔｅｃｈｎｏｌｏｇｙ． Ｊｏｕｒｎａｌ ｏｆ Ｈｅｎａｎ Ａｇｒｉ⁃
ｃｕｌｔｕｒａｌ Ｕｎｉｖｅｒｓｉｔｙ （河南农业大学学报）， １９９０， ２４
（４）： ４３６－４４３ （ｉｎ Ｃｈｉｎｅｓｅ）
［４］　 Ｈｕａｎｇ Ｄ⁃Ｃ （黄顶成）， Ｚｈａｎｇ Ｒ⁃Ｚ （张润志）， Ｗａｎｇ Ｂ
（王 　 斌）． Ｔｈｅ ｐｏｐｕｌａｔｉｏｎ ｄｙｎａｍｉｃｓ ａｎｄ ｐｒｅｖｅｎｔｉｏｎ
ｃｏｕｎｔｅｒｍｅａｓｕｒｅｓ ｏｆ Ｍｙｚｕｓ ｐｅｒｃｉｃａｅ ｉｎ ｐｒｏｔｅｃｔｉｖｅ ａｒｅａ．
Ｐｌａｎｔ Ｐｒｏｔｅｃｔｉｏｎ （植物保护）， ２００９， ３５（６）： １０８－１１３
（ｉｎ Ｃｈｉｎｅｓｅ）
［５］ 　 Ｌｉ Ｊ⁃Ｌ （李厥鲁）． Ｔｈｅ ｏｃｃｕｒｒｅｎｃｅ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ａｎｄ
ｔｈｅ ｐｒｅｖｅｎｔｉｏｎ ｄｉｖｉｓｉｏｎ ｏｆ ｐｅａｃｈ ａｐｈｉｄ ｉｎ ｄｉｆｆｅｒｅｎｔ ｅｃｏ⁃
ｌｏｇｉｃａｌ ｒｅｇｉｏｎ． Ｃｈｉｎａ Ｔｏｂａｃｃｏ （中国烟草）， １９９４（４）：
１－８ （ｉｎ Ｃｈｉｎｅｓｅ）
［６］　 Ｄｅｎｇ Ｊ⁃Ｈ （邓建华）， Ｌｉ Ｔ⁃Ｆ （李天飞）， Ｗｕ Ｘ⁃Ｆ （吴
兴富）． Ｔｈｅ ｖｅｒｔｉｃａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ａｎｄ ｄｉｓｔｒｉｂｕｔｉｏｎ ｔｙｐｅ ｏｆ
ｐｅａｃｈ ａｐｈｉｄ ｏｎ ｔｏｂａｃｃｏ ｐｌａｎｔｓ． Ｚｏｏｌｏｇｉｃａｌ Ｒｅｓｅａｒｃｈ （动
物学研究）， １９９８， １９（６）： ４４６－４５２ （ｉｎ Ｃｈｉｎｅｓｅ）
［７］　 Ｓｈａｎｇ Ｓ⁃Ｈ （商胜华）， Ｃｈｅｎ Ｑ⁃Ｙ （陈庆园）， Ｘｕ Ｍ⁃Ｌ
（徐卯林）， ｅｔ ａｌ． Ｔｈｅ ｐｒｅｌｉｍｉｎａｒｙ ｓｔｕｄｙ ｏｆ ｔｈｅ ｏｃｃｕｒ⁃
ｒｅｎｃｅ ｒｅｇｕｌａｒｉｔｙ ａｎｄ ｐｒｅｄｉｃｔｉｏｎ ｍｏｄｅｌ ｏｆ Ｍｙｚｕｓ ｐｅｒｃｉｃａｅ
ｉｎ Ｇｕｉｚｈｏｕ Ｐｒｏｖｉｎｃｅ． Ｐｌａｎｔ Ｐｒｏｔｅｃｔｉｏｎ （植物保护），
２０１０， ３６（５）： ８６－９１ （ｉｎ Ｃｈｉｎｅｓｅ）
［８］　 Ｓｈｉ Ｘ⁃Ｌ （施晓丽）， Ｂｉ Ｓ⁃Ｄ （毕守东）， Ｇｅｎｇ Ｊ⁃Ｇ （耿
继光）， ｅｔ ａｌ． Ｒｅｌａｔｉｏｎｓｈｉｐｓ ｂｅｔｗｅｅｎ ｍａｉｎ ｉｎｓｅｃｔ ｐｅｓｔｓ
ａｎｄ ｔｈｅｉｒ ｐｒｅｄａｔｏｒｙ ｎａｔｕｒａｌ ｅｎｅｍｉｅｓ ｉｎ ‘５１８’ ｎｅｃｔａｒｉｎｅ
ｏｒｃｈａｒｄ． Ａｃｔａ Ｅｃｏｌｏｇｉｃａ Ｓｉｎｉｃａ （生态学报）， ２０１１， ３１
（１５）： ４３７２－４３８４ （ｉｎ Ｃｈｉｎｅｓｅ）
［９］　 Ｌｉ Ｙ⁃Ｃ （李友常）， Ｘｉａ Ｎ⁃Ｂ （夏乃斌）， Ｔｕ Ｑ⁃Ｈ （屠泉
洪）， ｅｔ ａｌ． Ａ ｇｅｏｓｔａｔｉｓｔｉｃａｌ ａｎａｌｙｓｉｓ ｏｆ ｓｐａｔｉａｌ ｐａｔｔｅｒｎｓ ｏｆ
Ａｎｏｐｉｏｐｈｏｒａ ｇｌａｂｒｉｐｅｎｎｉｓ ｉｎ ｐｏｐｌａｒｓ． Ａｃｔａ Ｅｃｏｌｏｇｉｃａ Ｓｉｎｉ⁃
ｃａ （生态学报）， １９９７， １７（４）： ３９３－４０１ （ｉｎ Ｃｈｉｎｅｓｅ）
［１０］　 Ｚｏｎｇ Ｓ⁃Ｘ （宗世祥）， Ｌｕｏ Ｙ⁃Ｑ （骆有庆）， Ｘｕ Ｚ⁃Ｃ （许
志春）， ｅｔ ａｌ． Ｇｅｏｓｔａｔｉｓｔｉｃａｌ ａｎａｌｙｓｉｓ ｏｎ ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕ⁃
ｔｉｏｎ ｏｆ Ｈｏｌｃｏｃｅｒｕｓ ｈｉｐｐｏｐｈａｅｃｏｌｕｓ ｅｇｇｓ ａｎｄ ｌａｒｖａ． Ａｃｔａ
Ｅｃｏｌｏｇｉｃａ Ｓｉｎｉｃａ （生态学报）， ２００５， ２５（４）： ８３１－８３６
（ｉｎ Ｃｈｉｎｅｓｅ）
［１１］　 Ｗａｎｇ Ｔ （王　 涛）， Ｗｅｎ Ｊ⁃Ｂ （温俊宝）， Ｌｕｏ Ｙ⁃Ｑ （骆
有庆）， ｅｔ ａｌ． Ｇｅｏｓｔａｔｉｓｔｉｃａｌ ｓｔｕｄｙ ｏｎ ｔｈｅ ｓｐａｔｉａｌ ｐａｔｔｅｒｎ
ｏｆ Ａｎｏｐｌｏｐｈｏｒａ ｇｌａｂｒｉｐｅｎｎｉｓ ｉｎ ｔｈｒｅｅ ｔｙｐｅｓ ｏｆ ｓｔａｎｄｓ． Ａｃｔａ
Ｅｃｏｌｏｇｉｃａ Ｓｉｎｉｃａ （生态学报）， ２００６， ２６（９）： ３０４１－
３０４８ （ｉｎ Ｃｈｉｎｅｓｅ）
［１２］　 Ｈｕａｎｇ Ｓ⁃Ｓ （黄寿山）， Ｈｕ Ｈ⁃Ｊ （胡慧建）， Ｌｉａｎｇ Ｇ⁃Ｗ
（梁广文）． Ｔｈｅ ｓｐａｔｉａｌ ｐａｔｔｅｒｎ ｏｆ ｄｉｐａｕｓｉｎｇ ｌａｒｖａ ｏｆ
ｓｔｒｉｐｅｄ ｒｉｃｅ ｂｏｒｅｒ， Ｃｈｉｌｏ ｓｕｐｐｅｒｓｓａｌｉｓ Ｗａｌｋｅｒ， ｗｉｔｈ ｇｅｏ⁃
ｇｒａｐｈｉｃｓ． Ａｃｔａ Ｅｃｏｌｏｇｉｃａ Ｓｉｎｉｃａ （生态学报）， １９９９， １９
（２）： ２５０－２５３ （ｉｎ Ｃｈｉｎｅｓｅ）
［１３］　 Ｓｈｉ Ｇ⁃Ｓ （石根生）， Ｌｉ Ｄ⁃Ｍ （李典谟）． Ｇｅｏｓｔａｔｉｓｔｉｃｓ
ａｎａｌｙｓｉｓ ｏｆ ｓｐａｔｉａｌ ｐａｔｔｅｒｎ ｏｆ Ｄｅｎｄｒｏｌｉｍｕｓ ｐｕｎｃｔａｔｕｓ． Ｃｈｉ⁃
ｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用生态学报），
１９９７， ８（６）： ６１２－６１６ （ｉｎ Ｃｈｉｎｅｓｅ）
［１４］　 Ｙｕａｎ Ｚ⁃Ｍ （袁哲明）， Ｆｕ Ｗ （付　 威）， Ｌｉ Ｆ⁃Ｙ （李方
一）． Ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ｐａｔｔｅｒｎ ｏｆ Ｃｈｉｌｏ ｓｕｐｐｒｅｓｓａｌｉｓ ａｎ⁃
ａｌｙｚｅｄ ｂｙ ｃｌａｓｓｉｃａｌ ｍｅｔｈｏｄ ａｎｄ ｇｅｏｓｔａｔｉｓｔｉｃｓ． Ｃｈｉｎｅｓｅ
Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用生态学报）， ２００４， １５
（４）： ６１０－６１４ （ｉｎ Ｃｈｉｎｅｓｅ）
［１５］　 Ｙａｎ Ｘ⁃Ｈ （闫香慧）， Ｈｕａｎｇ Ｙ （黄　 燕）． Ｇｅｏｓｔａｔｉｓｔｉ⁃
ｃａｌ ａｎａｌｙｓｉｓ ｏｎ ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｍａｃｒｏｐｔｅｒｏｕｓ ｗｈｉｔｅ⁃
ｂａｃｋｅｄ ｐｌａｎｔｈｏｐｐｅｒ， Ｓｏｇａｔｅｌｌａ ｆｕｒｃｉｆｅｒａ （ Ｈｏｒｖáｔｈ ）．
Ｊｏｕｒｎａｌ ｏｆ Ｃｈｉｎａ Ｗｅｓｔ Ｎｏｒｍａｌ Ｕｎｉｖｅｒｓｉｔｙ （Ｎａｔｕｒａｌ Ｓｃｉ⁃
ｅｎｃｅ） （西华师范大学学报：自然科学版）， ２０１２， ３３
（１）： ５４－６０ （ｉｎ Ｃｈｉｎｅｓｅ）
［１６］　 Ｚｏｕ Ｙ⁃Ｄ （邹运鼎）， Ｂｉ Ｓ⁃Ｄ （毕守东）， Ｚｈｏｕ Ｘ⁃Ｚ （周
夏芝）， ｅｔ ａｌ． Ｇｅｏｓｔａｔｉｓｔｉｃａｌ ａｎａｌｙｓｉｓ ｏｎ ｓｐａｔｉａｌ ｐａｔｔｅｒｎｓ
ｏｆ Ｅｒｙｔｈｒｏｎｅｕ ｒａｓｕｄｒａ ａｎｄ Ｅｒｉｇｏｎｉｄｉｕｍ ｇｒａｍｉｎｉｃｏｌａ． Ｃｈｉ⁃
ｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用生态学报），
２００２， １３（１２）： １６４５－１６４８ （ｉｎ Ｃｈｉｎｅｓｅ）
［１７］　 Ｙａｎ Ｘ⁃Ｈ （闫香慧）， Ｚｈａｏ Ｚ⁃Ｍ （赵志模）， Ｌｉｕ Ｈ （刘
怀）． Ｇｅｏｓｔａｔｉｓｔｉｃａｌ ａｎａｌｙｓｉｓ ｏｎ ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｔｈｅ
ｎｙｍｐｈ ｏｆ ｗｈｉｔｅｂａｃｋｅｄ ｐｌａｎｔｈｏｐｐｅｒ， Ｓｏｇａｔｅｌｌａ ｆｕｒｃｉｆｅｒａ
（Ｈｏｒｖáｔｈ）． Ｓｃｉｅｎｔｉａ Ａｇｒｉｃｕｌｔｕｒａ Ｓｉｎｉｃａ （中国农业科
学）， ２０１０， ４３（３）： ４９７－５０６ （ｉｎ Ｃｈｉｎｅｓｅ）
［１８］ 　 Ｚｏｕ Ｙ⁃Ｄ （邹运鼎）， Ｂｉ Ｓ⁃Ｄ （毕守东）， Ｗａｎｇ Ｘ⁃Ｓ
（王祥胜）， ｅｔ ａｌ． Ｇｅｏｓｔａｔｉｓｃａｌ ａｎａｌｙｓｉｓ ｏｎ ｓｐａｔｉａｌ ｐａｔ⁃
ｔｅｒｎｓ ｏｆ Ｍａｃｒｏｓｉｐｈｕｍ ｇｒａｎａｒｉｕｍ ａｎｄ Ａｐｈｉｄｉｕｓ． Ｃｈｉｎｅｓｅ
Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用生态学报）， ２００１， １２
（６）： ８８７－８９１ （ｉｎ Ｃｈｉｎｅｓｅ）
［１９］　 Ｗｉｌｌｉａｍｓ Ｌ， Ｓｃｈｏｔｚｋｏ ＤＪ， ＭｃＣａｆｆｒｅｙ ＪＰ． Ｇｅｏｓｔａｔｉｓｔｉｃａｌ
ｄｅｓｃｒｉｐｔｉｏｎ ｏｆ ｔｈｅ ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ Ｌｉｍｏｎｉｕｓ ｃａｌｉｆｏｒ⁃
ｎｉｃｕｓ （Ｃｏｌｅｏｐｔｅｒａ： Ｅｌａｔｅｒｉｄａｅ） ｗｉｒｅｗｏｒｍｓ ｉｎ ｔｈｅ ｎｏｒｔｈ⁃
ｗｅｓｔｅｒｎ Ｕｎｉｔｅｄ Ｓｔａｔｅｓ， ｗｉｔｈ ｃｏｍｍｅｎｔｓ ｏｎ ｓａｍｐｌｉｎｇ．
Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｅｎｔｏｍｏｌｏｇｙ， １９９２， ２１： ９８３－９９５
［２０］　 Ｄａｖｉｄ ＦＮ， Ｍｏｏｒｅ ＰＧ． Ｎｏｔｅｓ ｏｎ ｃｏｎｔａｇｉｏｕｓ ｄｉｓｔｒｉｂｕｔｉｏｎｓ
３５５２期　 　 　 　 　 　 　 　 　 　 　 　 　 　 夏鹏亮等： 烟蚜在烤烟田分布动态的地统计学分析　 　 　 　 　 　
ｉｎ ｐｌａｎｔ ｐｏｐｕｌａｔｉｏｎｓ． Ａｎｎａｌｓ ｏｆ Ｂｏｔａｎｙ， １９５４， １８： ４７－
５３
［２１］　 Ａｒｂｏｕｓ ＡＧ， Ｋｅｒｒｉｃｈ ＪＥ． Ａｃｃｉｄｅｎｔ ｓｔａｔｉｓｔｉｃｓ ａｎｄ ｔｈｅ ｃｏｎ⁃
ｃｅｐｔ ｏｆ ａｃｃｉｄｅｎｔ⁃ｐｒｏｎｅｎｅｓｓ． Ｂｉｏｍｅｔｒｉｃｓ， １９５１， ７： ３４０－
４３２
［２２］　 Ｗｅｂｓｔｅｒ Ｒ． Ｅｓｔｉｍａｔｉｎｇ ｔｈｅ ｓｐａｔｉａｌ ｓｃａｌｅｓ ｏｆ ｒｅｇｉｏｎａｌｉｚｅｄ
ｖａｒｉａｂｌｅｓ ｂｙ ｎｅｓｔｅｄ ｓａｍｐｌｉｎｇ， ｈｉｅｒａｒｃｈｉｃａｌ ａｎａｌｙｓｉｓ ｏｆ
ｖａｒｉａｎｃｅ ａｎｄ ｒｅｓｉｄｕａｌ ｍａｘｉｍｕｍ ｌｉｋｅｌｉｈｏｏｄ． Ｃｏｍｐｕｔｅｒｓ
ａｎｄ Ｇｅｏｓｃｉｅｎｃｅｓ， ２００６， ３２： １３２０－１３３３
［２３］　 Ｄｕｔｔａ Ｓ， Ｂｈａｔｔａｃｈａｒｙａ ＢＫ， Ｒａｊａｋ ＤＲ， ｅｔ ａｌ． Ｍｏｄｅｌｌｉｎｇ
ｒｅｇｉｏｎａｌ ｌｅｖｅｌ ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ａｐｈｉｄ （Ｌｉｐａｐｈｉｓ ｅｒｙ⁃
ｓｉｍｉ） ｇｒｏｗｔｈ ｉｎ Ｉｎｄｉａｎ ｍｕｓｔａｒｄ ｕｓｉｎｇ ｓａｔｅｌｌｉｔｅ⁃ｂａｓｅｄ ｒｅ⁃
ｍｏｔｅ ｓｅｎｓｉｎｇ ｄａｔａ． Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｐｅｓｔ Ｍａｎａｇｅ⁃
ｍｅｎｔ， ２００８， ５４： ５１－６２
［２４］　 Ｊｉａｎｇ Ｙ⁃Ｊ （蒋勇军）， Ｙｕａｎ Ｄ⁃Ｘ （袁道先）， Ｘｉｅ Ｓ⁃Ｙ
（谢世友）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｓｏｉｌ ｏｒｇａｎｉｃ ｍａｔ⁃
ｔｅｒ ｃｏｎｔｅｎｔ ｉｎ ａ ｔｙｐｉｃａｌ ｋａｒｓｔ ｗａｔｅｒｓｈｅｄ： Ｃａｓｅ ｓｔｕｄｙ ｏｆ
Ｘｉａｏｊｉａｎｇ ｗａｔｅｒｓｈｅｄ， Ｙｕｎｎａｎ Ｐｒｏｖｉｎｃｅ． Ａｃｔａ Ｅｃｏｌｏｇｉｃａ
Ｓｉｎｉｃａ （生态学报）， ２００７， ２７（５）： ２０４０－ ２０４７ （ ｉｎ
Ｃｈｉｎｅｓｅ）
［２５］　 Ｌｕ Ｐ （路　 鹏）， Ｐｅｎｇ Ｐ⁃Ｑ （彭佩钦）， Ｓｏｎｇ Ｂ⁃Ｌ （宋
变兰）， ｅｔ ａｌ． Ｇｅｏｓｔａｔｉｓｔｉｃａｌ ａｎｄ ＧＩＳ ａｎａｌｙｓｅｓ ｏｎ ｓｏｉｌ ｔｏ⁃
ｔａｌ Ｐ ｉｎ ｔｈｅ ｔｙｐｉｃａｌ ａｒｅａ ｏｆ Ｄｏｎｇｔｉｎｇ Ｌａｋｅ Ｐｌａｉｎ． Ｓｃｉｅｎｔｉａ
Ａｇｒｉｃｕｌｔｕｒａ Ｓｉｎｉｃａ （中国农业科学）， ２００５， ３８（６）：
１２０４－１２１２ （ｉｎ Ｃｈｉｎｅｓｅ）
［２６］　 Ｙｕａｎ Ｚ⁃Ｍ （袁哲明）， Ｌｉ Ｆ⁃Ｙ （李方一）， Ｈｕ Ｘ⁃Ｙ （胡
湘粤）， ｅｔ ａｌ． Ｇｅｏｓｔａｔｉｓｔｉｃａｌ ａｎａｌｙｓｉｓ ｏｎ ｔｅｍｐｏｒａｌ ｐａｔ⁃
ｔｅｒｎｓ ｏｆ Ｃｈｉｌｏ ｓｕｐｐｒｅｓｓａｌｉｓ ｐｏｐｕｌａｔｉｏｎ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ
Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用生态学报）， ２００６， １７ （ ４）：
６７３－６７７ （ｉｎ Ｃｈｉｎｅｓｅ）
［２７］　 Ｆｕｅｎｔｅｓ Ｍ， Ｃｈａｕｄｈｕｒｉ Ａ， Ｈｏｌｌａｎｄ Ｄ． Ｂａｙｅｓｉａｎ ｅｎｔｒｏｐｙ
ｆｏｒ ｓｐａｔｉａｌ ｓａｍｐｌｉｎｇ ｄｅｓｉｇｎ ｏｆ ｅｎｖｉｒｏｎｍｅｎｔａｌ ｄａｔａ． Ｅｎｖｉ⁃
ｒｏｎｍｅｎｔａｌ ａｎｄ Ｅｃｏｌｏｇｉｃａｌ Ｓｔａｔｉｓｔｉｃｓ， ２００７， １４： ３２３－３４０
［２８］ 　 Ｚｈｏｕ Ｑ （周 　 强）， Ｚｈａｎｇ Ｒ⁃Ｊ （张润杰）， Ｇｕ Ｄ⁃Ｘ
（古德祥）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｓｔｒｕｃｔｕｒｅ ａｎｄ ｐａｔｔｅｒｎ ｏｆ Ｎｉ⁃
ｌａｐａｒｖａ ｔａｌｕｇｅｎｓ ｐｏｐｕｌａｔｉｏｎ ｉｎ ｌａｒｇｅ⁃ｓｃａｌｅ． Ｃｈｉｎｅｓｅ Ｊｏｕｒ⁃
ｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用生态学报）， ２００１， １２
（２）： ２４９－２５２ （ｉｎ Ｃｈｉｎｅｓｅ）
［２９］　 Ｔａｎｇ Ｑ⁃Ｙ （唐启义）， Ｆｅｎｇ Ｇ⁃Ｍ （冯明光）． Ｐｒａｃｔｉｃａｌ
Ｓｔａｔｉｓｔｉｃａｌ Ａｎａｌｙｓｉｓ ａｎｄ ＤＰＳ Ｄａｔａ Ｐｒｏｃｅｓｓｉｎｇ Ｓｙｓｔｅｍ．
Ｂｅｉｊｉｎｇ： Ｓｃｉｅｎｃｅ Ｐｒｅｓｓ， ２００２： ２５０－２８１ （ｉｎ Ｃｈｉｎｅｓｅ）
［３０］　 Ｚｈｏｕ Ｇ⁃Ｆ （周国法）， Ｘｕ Ｒ⁃Ｍ （徐汝梅）． Ｒｅｓｅａｒｃｈ ｏｎ
ｔｈｅ Ｂｉｏｇｅｏｓｔａｔｉｓｔｉｃｓ． Ｂｅｉｊｉｎｇ： Ｓｃｉｅｎｃｅ Ｐｒｅｓｓ， １９９８： ３－
１０ （ｉｎ Ｃｈｉｎｅｓｅ）
作者简介　 夏鹏亮，男，１９８６年生，博士，农艺师．主要从事烟
草植保研究． Ｅ⁃ｍａｉｌ： ２７２８６８４４１＠ ｑｑ．ｃｏｍ
责任编辑　 肖　 红
４５５ 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 应　 用　 生　 态　 学　 报　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 ２６卷