GS_变异函数摘要.pptVIP

变异函数的计算与拟合 人工拟合方法：观察-调整； 最小二乘方法：确定模型-用最小二乘法求最佳参数； 加权多项式回归方法：确定模型-用对应滞后距的点对数N(hi)为权系数，在最小二乘意义下拟合。 变异函数的计算与拟合 加权多项式回归拟合单个球状模型 模型选择：单个球状模型 模型变换:变成多项式形式 加权处理：每个点对计算结果乘以权系数 最小二乘计算： 变异函数的计算与拟合 加权多项式回归方法 计算示例： b0=1.02;b1=2.28;b2=-0.01 C0=1.02;a=8.312;C=12.634 C0+C=13.654 h r(h) 1.5 4.1 2.3 6.3 3.2 8.0 3.9 9.5 4.7 10.7 5.4 11.7 6.2 12.5 7 13.0 7.7 13.3 8.5 13.5 9.2 13.7 变异函数的计算与拟合 加权多项式回归拟合二级套合球状模型 模型选择：二级球状模型 模型变换:变成多项式形式 加权处理：每个点对计算结果乘以权系数 最小二乘计算：分段用最小二乘拟合计算 变异函数的计算与拟合 加权多项式回归拟合二级套合球状模型 模型选择：二级球状模型 模型变换:变成多项式形式 加权处理：每个点对计算结果乘以权系数 最小二乘计算：分段用最小二乘拟合计算 变异函数的计算与拟合 加权多项式回归方法 计算示例： C0=0.3353;a1=2.198;C1=0.8665 a2=8.1932;C2=1.43 h r(h) Ni 0.5 0.73 19 1 1.16 18 1.5 1.44 18 2 1.65 17 2.5 1.79 17 3 2 16 3.5 2.07 15 4 2.18 15 4.5 2.25 14 5 2.3 14 5.5 2.44 13 6 2.53 13 6.5 2.51 12 7 2.59 12 7.5 2.64 11 变异函数的计算与拟合 加权多项式回归方法 计算示例： C0=0.273;a=7.383;C=5.4733 h r(h) 0.9 1.3 1.9 2.2 3.1 3 4 4.5 4.8 5.8 6.2 4.7 7.9 6.1 9.4 5 变异函数的计算与拟合 人工拟合方法：观察-调整； 最小二乘方法：确定模型-用最小二乘法求最佳参数； 加权多项式回归方法：确定模型-用对应滞后距的点对数N(hi)为权系数，在最小二乘意义下拟合。 区域化变量结构分析 实验变异函数-用理论函数似合-区域化变量分析 结构分析-构造变异函数模型对全部有效结构信息作定量化的概括，以表征区域化变量的主要特征。 主要使用方法套合结构法（nest structure); 套合结构是把分别出现在不同滞后距h上和不同方向a上同时起作用的变异性组合起来。形成一个统一的描述区域化变量特征的变异函数 套合结构表示为： γ(h)=γ0(h)+γ1(h)+γ2(h)+…+ γ(h) +… 区域化变量结构分析一个方向上的套合结构 套合结构表示为： γ(h)=γ0(h)+γ1(h)+γ2(h)+…+ γi(h) +… 套合结构中每个变异函数模型代表一个特定尺度上的变异性，可以是用不同模型表示变异函数； 例如：某区域化变量在某一方向上的变异性可用三个层次的变异性表示： γ(h)=γ0(h)+γ1(h)+γ2(h) 区域化变量结构分析不同方向的套合结构 概念： 各向同性：区域化变量在各个不同方向上性质相同； 各向异性：不同方向上性质不同； 各向异性表现为不同方向的变异函数的差异； 各向异性的分类及特征： 几何异向性和带状异向性 区域化变量结构分析不同方向的套合结构 几何异向性：（geometric anisotropy) 区域化变量在不同方向上表现出变异程度相同（基台相同）而连续性不同（变程不同）的异向性 为什么称为几何异向性？ 这种异向性可以通过简单几何变换变为各向同性； 各向异性比K: 变程比a1/a2：表示不同方向的变异性之差异大小 各向异性表现为不同方向的变异函数的差异； 区域化变量结构分析不同方向的套合结构 二维几何异向性变换： 水平方向球状模型1：参数C，a2; 垂直方向球状模型2：参数C，a1; 各向异性比K=a2/a1； 区域化变量结构分析不同方向的套合结构 带状异向性：（zonal anisotropy) 区域化变量在不同方向上表现出的变异性之差不能简单用变程比得到，即无法用几何变换将它们变为各向同性 为什么称为带状异向性？ 不同的变异性只局限在该方向带中。 例：多层状矿区，由于矿层及夹层组分变化显著，垂直方向的变异比平面上的变异大。垂直方向含有重力作用因素，矿化作用包括成矿元素的自然扩散和重力双重作用。 平面γ1(h) γ(hw)=γ1(h)+ γ2(hw) 区域化变量结构分析几