2. 土壤空间变异基础与采样设计关联

土壤并不是一块均匀的“面团”，即便在同一块田里，pH、有机质、有效磷、重金属含量也会随位置和深度发生变化。这种变化就是土壤空间变异。采样设计的本质，是用有限的采样点去“估计”这种变异。因此，理解变异的三个核心指标——变异系数、空间自相关和半方差函数——直接决定了网格布多密、分层怎么分。

一、变异系数：先判断“波动有多大”

变异系数（CV）是标准差与均值之比，用于比较不同量纲属性的离散程度。经验上常按 CV 划分：CV < 10% 为弱变异，10%～100% 为中等变异，> 100% 为强变异。例如某地块有机质均值 22.5 g/kg、标准差 3.6 g/kg，则 CV ≈ 16%，属中等变异；而有效磷均值 18 mg/kg、标准差 25 mg/kg，CV ≈ 139%，属强变异。

CV 的意义在于给出“采样密度下限”的直觉：弱变异属性少量点即可代表整体，强变异属性则必须加密。但 CV 只反映离散程度，不反映空间结构——两个 CV 相同的属性，可能一个是随机噪声，一个是明显的斑块分布，采样策略完全不同。这就引出了空间自相关。

二、空间自相关与半方差函数：判断“波动有没有规律”

空间自相关指距离越近的样点，属性值越相似。常用 Moran's I 或半方差函数刻画。半方差函数 γ(h) 定义为距离 h 的两点属性差平方的一半的期望值。随 h 增大，γ(h) 通常先上升，达到一个平台后趋于稳定，这个平台叫基台值，对应的最小距离叫变程，而 γ(0) 处的截距叫块金值。

三个参数各有采样含义：变程是空间自相关的有效范围，采样间距应小于变程的一半，否则相邻点已“不相关”，插值会失真；块金值/基台值反映随机变异占比，比值高说明小尺度随机波动强，需要更密的网格；基台值则对应总体变异幅度。下图展示了同一地块上两种属性的半方差曲线差异及其对布点密度的直接含义。

三、变异信息如何影响网格密度与分层策略

把上述指标落到设计上，可归纳为两条主线。第一条是网格密度：以变程为基准，常规做法是网格边长取变程的 1/2～1/4。若某属性变程约 40 m，则 10～20 m 网格较合理；若变程仅 18 m，则需 5～9 m 网格，否则插值图会出现虚假平滑。第二条是分层策略：当变异主要来自垂直方向（如耕作层与犁底层、不同母质层），应按发生层或固定深度分层采样；当变异主要来自水平方向（如施肥斑块、污染羽流），则应在同一深度内加密平面网格，而非盲目增加层数。

判断变异方向的一个实用做法，是分别计算平面半方差和垂直半方差，比较二者基台值。若垂直基台值显著大于平面，说明层间差异主导，分层优先；反之则平面加密优先。这一步能避免“层分了很多、平面点却不够”的常见浪费。

四、案例：用变异信息指导初步采样设计

假设某蔬菜基地约 4 hm²（200 m × 200 m），拟调查 pH、有机质、有效磷和镉。初步踏勘采集 30 个点做半方差预分析（示例数据，用于说明流程）：pH 的 CV 约 6%、变程约 55 m、块金比 25%；有机质 CV 约 16%、变程约 40 m、块金比 35%；有效磷 CV 约 139%、变程约 18 m、块金比 60%；镉 CV 约 45%、变程约 30 m、块金比 40%。

据此可做差异化设计：pH 变程大、块金低，20 m 网格（约 100 点）已足够；有机质取 15 m 网格；有效磷变程小且块金比高，若全场按 9 m 网格需约 500 点，成本过高，可采用“两级布点”——先按 20 m 网格做初步调查，对有效磷高值区（如棚室周边）再局部加密到 8～10 m；镉作为污染指标，按 15 m 网格并在灌溉入口、道路两侧增设判断性样点。分层上，耕作层 0～20 cm 全场必采，犁底层 20～40 cm 仅在怀疑累积或背景调查时按 1/4 网格抽采。

概括来说，变异系数回答“要不要密”，空间自相关和半方差回答“密在哪里、往哪个方向密”，分层策略则回答“密在平面上还是剖面上”。把这三者与成本约束结合，才能得到既科学又可执行的初步采样方案。