3. 网格采样类型与适用条件

网格采样并不是“画格子、取土样”这么简单。格子的形状、布设方式以及是否叠加分层信息,直接决定了样本对目标变量的代表性。本节围绕三类常见网格——规则网格、系统随机网格、分层随机网格——讲清它们的构造方式、优缺点和适用条件,并用一块农田与一块污染场地的对比案例,说明如何依据目标变量和空间尺度做出选择。

三类网格的构造逻辑

规则网格是最基础的形式:在采样区内按固定间距布设节点,形成正方形、矩形或三角形(六边形)网格。正方形网格实现简单、坐标易算,适合用GPS或RTK直接放样;三角形网格在同等点密度下对空间连续变量的覆盖更均匀,常用于土壤属性制图。规则网格的最大优点是空间覆盖均匀、可重复性强、便于插值成图;缺点是当采样区存在明显的周期性条带(如耕作垄向、施肥行)时,固定间距可能与条带周期“共振”,导致系统性偏差。

系统随机网格保留规则网格的骨架,但给每个节点加一个随机偏移量,或随机抽取一个起点后再按固定间距布点。它既保持了空间均匀性,又打破了规则网格可能带来的周期性偏差,同时满足经典统计对“随机性”的要求,便于计算均值、方差等统计量。代价是点位坐标不再整齐,放样和复核稍麻烦,且随机偏移过大会破坏网格的均匀性。

分层随机网格先依据辅助信息把采样区划分为若干相对均质的子区(层),再在各层内独立布设随机或系统随机点位。分层依据可以是土壤类型图、土地利用、地形高程、历史监测数据或污染源分布。它的核心价值在于把样本量按层内变异大小分配,变异大的层多布点、变异小的层少布点,从而在总样本量不变时降低估计方差。缺点是分层边界本身可能不准,层划分过细会导致每层样本不足。

三类网格采样布设方式对比 规则网格 等间距 · 覆盖均匀 · 易放样 系统随机网格 规则骨架 + 随机偏移 · 破周期性 分层随机网格 先分层 · 层内随机 · 按变异配点 适用:均质农田、属性制图 空间连续、无强周期条带 适用:需统计推断、存在条带 兼顾均匀性与随机性 适用:异质场地、污染调查 有辅助分层信息、变异差异大

如何根据目标变量与空间尺度选择

选择网格类型,先问两个问题:目标变量是什么性质,采样区在什么尺度上变化?

若目标是土壤养分、有机质、pH等连续渐变型变量,且农田地块平整、耕作管理一致,空间变异以渐变为主,规则网格或系统随机网格即可胜任。此时关键是网格间距要与变异尺度匹配:间距应小于目标变量空间相关距离的一半左右,否则插值图会出现“块状失真”。例如某地块实测(示例假设)养分空间相关距离约60 m,则网格间距宜取25~30 m,而不是50 m以上。

若目标是污染物的点源扩散、历史堆存或事故泄漏,空间变异呈“热点+背景”的强异质性,规则网格容易漏掉小范围高浓度热点。此时应优先采用分层随机网格:以污染源、地下水流向、历史工艺分区作为分层依据,在疑似源区和下游羽流区加密布点,在背景区降低密度。分层随机网格还能配合序贯采样——先粗网格筛查,发现异常后在异常点周围加密,逐步逼近热点边界。

空间尺度还影响样本量分配。大尺度区域调查(如县域尺度)常用规则网格保证覆盖;中尺度场地调查(如几公顷的污染场地)适合分层随机;小尺度精细调查(如源区几十平方米)则可在分层基础上采用三角形网格或加密系统随机点。

案例:同一块地,两种目标两种选法

假设有一块约8 hm²(约120 m×670 m)的农田,同时承担两项任务:一是绘制养分分布图指导变量施肥,二是排查历史上游某电镀厂可能造成的重金属污染。两项任务的目标变量和空间结构完全不同,网格选择也应不同。

任务一(养分制图):目标变量连续、变异尺度较大,选规则网格。按30 m间距布设,约需 (120/30+1)×(670/30+1) ≈ 5×23 = 115 个点位;若人力有限,可改为系统随机网格抽取约60~70点,既保证覆盖又满足统计要求。采样深度按耕作层0~20 cm,每点用土钻取5~8个子样混合成1个复合样,记录RTK坐标。

任务二(重金属排查):目标变量异质性强,选分层随机网格。以“距污染源距离+地下水流向”分为三层:源区下游0~100 m(高疑)、100~300 m(中疑)、300 m以外(背景)。按变异大小分配样本量,示例假设三层分别布设20、15、10个点,层内采用系统随机布设,源区层间距加密至15 m。先做粗筛,若某点超标,则以其为中心在10 m半径内追加4~6个点刻画热点边界。

判断要点:连续渐变变量看“相关距离”,选间距匹配的规则或系统随机网格;热点型变量看“源与路径”,选能按疑似的污染结构分配样本的分层随机网格。两者都要求网格间距与目标变量的空间尺度相匹配,而不是越密越好。

最后需要强调,网格类型不是一成不变的。实际项目中常采用“分层+系统随机”的混合方案:先按辅助信息分层,层内再用系统随机布点,既控制异质性又保持均匀覆盖。无论选哪种,都应在采样前用历史数据或踏勘信息估算变异尺度,采样后用半方差函数或插值交叉验证检查网格是否足够密——若验证误差偏大,说明网格间距仍大于变异尺度,需要加密或调整分层方案。