1. 课程导入与目标设定
在地质统计学中，变异函数（variogram）是描述区域化变量空间相关性的核心工具。无论是油气储层建模中的孔隙度、渗透率分布，还是矿产勘查中的品位估算，我们都需要回答同一个问题：相距多远的两个样品点，其属性值才不再具有明显的相关性？变异函数拟合正是把离散的实验变异函数点转化为可参数化的连续模型，而各向异性判断则决定了这种相关性在不同方向上是否存在差异。二者共同构成了空间插值（克里金）与随机模拟的输入基础，直接决定了储量估算、品位估值与不确定性评价的可信度。
本课程围绕一条完整的技术主线展开：从原始样品数据出发，计算实验变异函数，识别并量化各向异性，完成模型拟合与交叉验证，最终把拟合结果交付给储层建模或品位估算流程。本节作为导入，先建立全局工作流的认知框架，并布置一个贯穿全课程的真实数据集任务。
为什么变异函数拟合与各向异性判断是“地基”
克里金插值权重完全由协方差（即变异函数的补函数）决定。若变异函数模型选择不当，或各向异性方向、比率设置错误，插值结果会出现系统性偏差：在储层建模中可能把河道砂体的连续性人为拉长或截断，导致连通性评价失真；在品位估算中可能把矿化带的延伸方向判断错误，造成块段品位高估或低估。更隐蔽的问题是，拟合误差往往不会在插值图上直接暴露，而是通过交叉验证的残差结构、条件模拟的实现变异性间接体现。因此，变异函数拟合不是“画一条好看的曲线”，而是一次有明确物理地质含义的参数估计。
各向异性判断尤其关键。沉积储层常因水流方向形成条带状结构，矿体常因构造控制形成定向延伸。若忽略各向异性，用单一各向同性模型拟合，变程会被“平均化”，既不能反映主方向的连续性，也不能反映次方向的快速变化。判断各向异性的本质，是在多个方位上分别计算实验变异函数，比较其变程、基台与块金效应的差异，再决定是否采用几何各向异性（仅变程随方向变化）或带状各向异性（变程与基台同时变化）模型。
工作流中的关键决策点
上图把整个流程压缩为五个阶段，其中真正决定成败的是三个决策点。第一是方位角分组方案：常用做法是以 22.5° 或 45° 为间隔划分方位带，每个方位带设置角度容差（如 ±22.5°）和滞后距容差（如滞后距的一半）。分组过粗会掩盖方向差异，过细则每个方位带的点对数量不足、实验变异函数噪声过大。第二是各向异性类型选择：若各方向基台接近而变程差异明显，用几何各向异性，通过坐标压缩把椭圆变为圆；若基台也随方向变化，则需带状各向异性模型。第三是拟合权重：滞后距越小、点对越多，实验值越可靠，加权最小二乘通常按点对数量或 1/h² 赋权。
需要特别强调参照约定：本课程统一采用方位角从正北顺时针度量，主变程方向记为方位角 α，次变程方向为 α+90°；各向异性比定义为主变程与次变程之比，比值大于 1 表示主方向连续性更强。后续所有图示与参数表均遵循此约定，切换参照系时会明确说明。
贯穿课程的数据集任务
为使各节知识可落地，本课程布置一个贯穿始终的练习任务。数据集为一份教学用合成数据（非实测），包含 180 个样品点，每个点有平面坐标 X、Y（单位：米）与属性值 Z（模拟品位，单位 g/t）。数据由已知的各向异性球状模型生成：主方向方位角 45°，主变程 320 m，次变程 120 m，块金 0.15，基台 1.0，并叠加了少量随机噪声。学员的任务是仅凭样品点，反演出上述参数，并与“真值”对比评估拟合质量。
建议按以下顺序推进，每完成一步记录一次中间结果：
- 绘制样品点位置图与 Z 值气泡图，检查采样密度与空间覆盖是否均匀；
- 分别沿 0°、45°、90°、135° 四个方位计算实验变异函数，滞后距取 40 m，容差 20 m，最大滞后距取 400 m；
- 比较四个方向的变程，判断主方向是否落在 45° 附近，计算各向异性比；
- 用球状模型拟合主、次方向实验曲线，记录块金、基台、变程三参数；
- 用留一法交叉验证，检查标准化残差均值是否接近 0、方差是否接近 1；
- 把拟合参数与生成真值列表对比，分析偏差来源（分组方案、滞后距、模型类型或噪声）。
本节应达成的目标
学完本节，你应当能够用自己的话说明变异函数拟合与各向异性判断在储层建模、品位估算中分别解决什么问题；能够画出从数据准备到交叉验证的完整工作流，并指出三个关键决策点；能够理解本课程数据集的构成与任务要求，明确每一步要记录什么中间结果。后续各节将沿着这条工作流逐段展开：先讲实验变异函数的计算细节与参数敏感性，再讲各向异性的识别与量化，最后讲模型拟合、验证与向克里金/条件模拟的交付。请先把数据集下载并完成样品点位置图，这是后续所有练习的起点。