地质统计场地调查克里金方法选择标准指南（D5923-18）

📋 概述与适用范围：本标准（D5923-18）由美国材料与试验协会（ASTM）D18委员会下属D18.01分委会制定，首次发布于1996年，2018年完成最新修订。标准旨在为环境与岩土工程场地调查中的数据空间估计提供克里金方法的选择指导。地质统计学在处理土壤、岩石及所含流体的空间变异特征方面具有独特优势，而克里金作为一类最小二乘估计方法，能够有效利用样本的空间相关性进行无偏最优估计。标准明确指出其不以规定统一操作流程为目标，而是基于广泛共识的方法评述，供专业人员在特定项目条件下结合经验判断灵活选用。标准范围涵盖普通克里金、简单克里金、对数正态克里金、通用克里金和指示克里金五种常见形式，不涉及多元、时空及其他低频使用的方法，但此限定不代表对其有效性的否定。标准同时给出了克里金不适用、转而采用地质统计模拟的判别条件，并排除了样条、反距离加权等非地质统计替代方法。标准还引用了ASTM D653（与土壤、岩石及所含流体相关的术语标准）作为基础术语依据，在使用前应确认术语定义的一致性。该指南特别强调，其内容不能替代专业教育与经验积累，所有应用均需结合具体工程背景与判断，标准标题中的“标准”仅指通过了ASTM协商一致审批程序，而不代表对专业服务充分性的唯一衡量尺度。 ⚙️ 方法原理与选择准则：标准要求在选择克里金方法前必须完成三个核心步骤：明确研究目标、开展探索性数据分析以及分析空间变异结构。探索性数据用于识别数据分布特征、趋势成分、异常值及空间离群点，这直接关系到后续对平稳性假设的满足程度与克里金类型的选择。空间变异性的定量描述通常通过变异函数实现，标准虽未详细讲授变异函数理论，但强调其参数（块金、基台、变程等）对不同克里金方法适用性的关键影响。在方法选择上，标准提出了以下基本准则：若数据在空间上呈现均值平稳且无显著趋势，则普通克里金是最为稳健的默认选择；若均值已知且数据满足平稳性，简单克里金可提供更小的估计方差；若数据呈现明显的对数正态分布特征（例如污染物浓度），对数正态克里金可避免有偏估计；若存在明确的空间趋势（如地下水位坡度），通用克里金能够同时估计趋势与残差；若需估计超过阈值概率或识别空间分类边界，指示克里金则更为适合。标准特别指出，当数据高度偏离正态分布或存在强不确定性时，应优先考虑基于模拟的方法，而非简单的克里金插值。整个选择流程强调迭代验证——通过交叉验证和模型诊断来反复确认所选方法的合理性，避免机械套用。标准还提示，真实场地通常具有各向异性，应通过多方向变异函数分析判断是否采用各向同性假设或进行坐标变换。 📊 技术参数与指标：下表归纳了标准涉及的五种克里金方法的核心适用条件与参数特征。表格信息均来源于标准正文中对各方法使用前提的明确说明。

表1 常用克里金方法对比

🟦 方法名称	📏 均值假设	📐 数据分布要求	🎯 趋势处理	⚡ 典型应用
普通克里金	局部平稳，均值未知	无严格分布要求	不处理趋势	土壤属性空间估计
简单克里金	全局平稳，均值已知	无严格分布要求	不处理趋势	已知背景场修正
对数正态克里金	对数服从平稳	正偏态分布，对数正态	不处理趋势	污染物浓度插值
通用克里金	非平稳，含漂移	残差部分要求弱平稳	需指定趋势函数	地下水位、地形面
指示克里金	无均值要求	基于二元指示变量	可处理趋势	污染范围划定

表2 克里金方法选择关键判断指标

🟦 判断条件	📐 推荐方法	📏 参数注意点
数据平稳且无趋势	普通克里金	变异函数块金值不宜过高
数据平稳且均值已知	简单克里金	均值估计准确性至关重要
数据对数正态分布	对数正态克里金	需对原始数据进行对数变换
数据含明显趋势	通用克里金	趋势阶次需合理选择
估计超阈值概率	指示克里金	阈值设置需与工程标准一致

🔬 工程应用与注意事项：本标准在环境与岩土工程中具有广泛适用性，典型场景包括：场地污染调查中污染物空间分布图的生成、采矿区品位估算、地基土层参数空间插值以及地下水水位面的绘制。应用时需特别注意以下质量要点：其一，样本数量与空间布局直接影响变异函数计算的稳定性，标准建议至少应有50个以上采样点才能获得可靠变异函数；其二，数据变换（如对数变换）不能盲目使用，应通过假设检验（如夏皮罗-威尔克检验）确认分布特征；其三，通用克里金中的趋势函数应基于物理背景设定，避免纯粹数学拟合导致的过度参数化。特别警示：当数据存在剧烈空间变异或比例效应时，克里金估计结果可能出现负值或物理意义不合理的情况，此时应考虑转换类型或采用协同克里金等高级技术。标准还强调，克里金方法假设空间相关性是二阶平稳或本征平稳的，当研究区域存在明显不同尺度的变异或断裂时，应首先进行数据分区（如按地层或污染源区域分组），再分别建立局部变异函数模型。在工程应用层面，将克里金结果用于决策时必须伴随不确定性评估，例如通过克里金方差图或模型模拟区间给出估计精度，而非单纯依赖插值曲面。此外，标准明确指出不应以克里金结果作为唯一依据进行风险评估或工程数值计算，需综合地质判断与监测数据验证。对于时间序列或动态数据的空间估计，建议参考第1.2条的说明，考虑时空克里金方法或其他非本指南覆盖的专门技术。 ❓ 常见问题解答：