2020地学计算第四章.pptVIP

三种常用模型比较 0.95 2020/1/29 华中农业大学 资源与环境学院 21 有基台值模型 — 线性有基台值模型 C0 ：块金常数 C0+C ：基台值 C ：拱高 A ：常数，表示直线斜率 当 C0=0,C=1 时，称为标 准指数函数模型 2020/1/29 华中农业大学 资源与环境学院 22 有基台值模型 — 纯块金效应模型 2020/1/29 华中农业大学 资源与环境学院 23 无基台值模型 —— 线性无基台值模型 2020/1/29 华中农业大学 资源与环境学院 24 无基台值模型 —— 幂函数值模型 2020/1/29 华中农业大学 资源与环境学院 25 无基台值模型 —— 对数值模型 2020/1/29 华中农业大学 资源与环境学院 26 套合模型 ? 在实际中，有时区域化随机变量 Z(x) 的变化 相当复杂，往往包含各种 尺度 及各种层次 的变化，反映在变异函数 r(h) 上，就是单一 的模型结构不能将其合理表达，而是多层 次的结构相互叠加在一起，地统计学上称 为套合。所谓 套合结构 ，就是把分别出现 在不同距离 h 上或不同方向上同时起作用的 变异性组合起来，对全部有效的结构信息， 作定量化的概括，以表示区域化变量的主 要特征。 华中农业大学 资源与环境学院 27 2020/1/29 套合模型 ? 土壤是一个不均与、具有高度空间异质性的复合 体，它与土壤母质、气候、水文、地形和生物等 因素有关，分析土壤空间变异的因素，可将其变 异分为系统变异（土壤形成因素相互作用造成） 和随机变异（可以观测到的，但与土壤形成印务 无关且不能直接分析的）两大类。如由 h 分开的两 个点 x 和 x+h 的土壤某一性质 Z(x) 和 Z(x+h) 。当 h 趋 近于 0 时，可以认为两点间的差异完全是由取样和 测定误差造成，当 h 逐步增大，如 h1m ，差异可 能还要加上诸如水分等因素，当 h100m 时，在新 的变异要考虑地形的作用。 华中农业大学 资源与环境学院 28 2020/1/29 套合模型 ? ? 当 h 一定时，变异函数 r(h) 应包含小于 h 的所有影 响因素，因此，绝大多数变异函数都由下面两个 变异函数组成： r(h)=r 0 (h)+r 1 (h) ，即一个代表纯块金方差，一个 代表空间相关的方差。一般情况下，套合模型可 以用放映各种不同尺度变化的多个变异函数之和 表示，即： n r ( h ) ? r 0 ( h ) ? r 1 ( h ) ? ? ? r n ( h ) ? ? r i ( h ) i ? 0 r i (h) 可以是相同的或不同的理论模型 2020/1/29 华中农业大学 资源与环境学院 29 套合模型 ? 如，区域化变量 Z(x) 的变异性由 r 0 (h),r 1 (h) 和 r 2 (h) 组成，其中 2020/1/29 华中农业大学 资源与环境学院 30 套合模型 ? 三者组成的套合模型为： 2020/1/29 华中农业大学 资源与环境学院 31 套合模型 2020/1/29 华中农业大学 资源与环境学院 32 套合模型实例 ? 某区域镍蕴藏量的变异函数计算值 2020/1/29 华中农业大学 资源与环境学院 33 套合模型实例 ? 镍含量的变异函数曲线图 可以看出基台值在 2.55 左右 对头两点作切线与纵坐标 相交的交点就是块金值， 大约在 0.4 左右 0.4/2.55=0.156 则大约 16% 的误差是随机现 象造成的，而空间自相关 现象造成的误差为 84% 2020/1/29 华中农业大学 资源与环境学院 34 套合模型实例 C 0 =0.4 C 1 =1.55 a 1 =14 C 2 =0.6 a 2 =50 从图中可看出，理 论值与实际值差异 较大，尤其是在 15 到 40m 之间，因此， 需进行反复修改 2020/1/29 华中农业大学 资源与环境学院 35 套合模型实例 C 0 =0.4 C 1 =1.15 a 1 =12 C 2 =1 a 2 =60 从图中可看出，理 论值与实际值差异 拟合较好 2020/1/29 华中农业大学 资源与环境学院 36 变异函数理论模型的最优拟合 ? 要使变异函数的理论模型真实地描述变量 的变化规律，在建立理论模型过程中