变差函数球状模型的自动拟合与实现-物探与化探.PDF

　第３４卷第２期 物　探　与　化　探 Ｖｏｌ．３４，Ｎｏ．２　 　２０１０年４月 ＧＥＯＰＨＹＳＩＣＡＬ＆ＧＥＯＣＨＥＭＩＣＡＬＥＸＰＬＯＲＡＴＩＯＮ Ａｐｒ．，２０１０　 变差函数球状模型的自动拟合与实现 张小艳，谭勇 （西安科技大学计算机学院，陕西 西安　７１００５４） 摘 要：根据变差函数的随机性和空间结构性，综合利用变差函数计算方法和加权线性规划拟合方法，分别拟合出各 主要方向上的模型参数，再根据各向异性情况进行结构套合，实现了变差函数的计算及球状模型的自动拟合。针 对样本中存在特异值的情况，算法中还提供了相对变差函数方法，有效地抑制了特异值对变差函数的影响，保证了 球状模型拟合的精度。本算法在ＶＣ＋＋６．０中实现，利用拟合出的模型，对样本区域进行插值得到网格文件，调用 Ｓｕｒｆｅｒ８．０绘制了等值线图。通过交叉验证和综合法验证，表明拟合度较高。 关键词：变差函数；相对变差函数；球状模型；自动拟合；线性规划 中图分类号：Ｐ６３１　　　文献标识码：Ａ　　　文章编号：１０００－８９１８（２０１０）０２－０２５３－０５ 　　地质统计学经近半个世纪的发展，变差函数理 金法的球状模型插值精度最好，且优于常规插值方 ［６］ 论模型的拟合问题已得到了深入的研究，产生了多 法 ，球状模型的表达式为 ［１］ 种拟合方法。但这些方法通常采用人工拟合 ，在 ｒ（ｈ）＝　　　　　　　　　　　　　　　　　 拟合过程中，人为因素起主导作用，需要反复调整参 ０， ｈ＝０； 数，进行多次试验才能拟合出一条较理想的变差函 ３ ３ ｈ １ ｈ Ｃ＋Ｃ · － · ，　０＜ｈ ａ；（１） ０ ( ２ ａ ２ ３) ≤ 数曲线。这样导致拟合的自动化程度不高，降低了 {Ｃ＋Ｃ ａ ｈ＞ａ。 基于地质统计学商用软件的通用性和简便性。 ０ 目前已提出的典型拟合方法有：加权回归多项 式中，Ｃ为块金效应值，表示ｈ很小时两点间变量 ０ ［２］ ［１］ 值的变化，Ｃ为基台值，反映变量在研究范围内的变 式拟合法 、线性规划拟合法 、目标规划拟合 法［３］ ［４］ 异程度， 、加权线性规划拟合法 。其中，加权线性规 ａ为变程，ｈ为滞后距离。 划拟合方法综合了前两种方法的优点，既考虑了不 １．２　球状模型的拟合 同滞后距离下的数据对数目对实验变差函数计算结 拟合的主要任务是通过一组试验变差函数值， 果的可靠性影响，又保证了变差函数的成功拟合，而 确定出该模型的参数 Ｃ、Ｃ、ａ，对式（１）来说，主要 ０ 目标规划则缺乏简便性，不便于用计算机进行自动 对 范围内的变差函数表达式 化拟合。为此，笔者利用变差