基于地统计学插值方法的局部奇异性指数计算比较的研究.pdfVIP

第32卷第2期 地 理 科 学 Vol. 32 No. 2 2 0 12 年02 月 SCIENTIA GEOGRAPHICA SINICA Feb . , 2 0 1 2 基于地统计学插值方法的局部奇异性指数 计算比较研究 李晓晖，袁 峰，贾 蔡，张明明，周涛发 （合肥工业大学资源与环境工程学院，安徽合肥 230009 ） 摘要摘要：以铜陵矿集区土壤Pb 元素为例，研究稀疏采样条件下地统计学克里格方法，序贯高斯模拟方法对奇异性 指数计算的影响。研究结果表明，序贯高斯模拟方法强调了短距离范围内的空间不确定性，弥补了克里格方法 平滑效应的不足，对于精细重建土壤元素的空间分布特征具有更好的效果。对于稀疏采样的数据集，较之原始 数据和克里格方法，基于序贯高斯模拟方法获取的奇异性指数能够更精细的刻画局部空间结构，更好的应用于 土壤地球化学异常的提取和识别。 关关 键键 词词：序贯高斯模拟；克里格；奇异性指数；地统计学 中图分类号中图分类号：P628 文章标识码文章标识码：A 文章编号文章编号：1000-0690(2012)02-0136-07 土壤地球化学异常的提取和识别是勘查地球 理。插值方法中，以地统计学为理论基础的克里 [1, 2] [11] 化学和环境地球化学工作重要目标 。异常信息 格插值法由于其无偏最优特性 ，被广泛应用于地 [12] 往往与背景信息相叠加，或是弱异常掩盖在强变 质、环境、生态等行业 。但克里格方法基于滑动 异性的背景中，使得常规统计方法识别异常信息的 加权平均原理，以变异函数为基础对局部范围进 效果并不显著，无法很好的应用于区域地球化学的 行加权平均，因而不可避免的对数据造成平滑效 [3,4] [13] 异常识别 。多维分形理论的提出为地球化学数据 应 。而局部奇异性指数恰恰需要突出数据的差 [1,3~6] 的异常识别和处理提供了强有力的工具 。其中， 异和异常，克里格方法的平滑效应将显著降低数 基于多维分形理论的局部奇异性指数从多维分形 据集的高值和方差从而无法重建原始信息中的高 [14] 的观点解释了尺度独立性，利用局部自相似性来 值异常信息 。序贯高斯模拟方法同样属于地统 度量奇异性，表明了空间变量随尺度变化的统计 计学范畴，与克里格方法相比，其理论上更适于定 [4] [15] 特征 。通过对这种统计特征进行研究，可以有效 量刻画某一属性的非均质性和不确定性 ，得到的 的应用于低缓、微弱、难识别异常信息的获取，并 模拟值不但能保持与样本数据的数学期望、方差 [7] 可以从背景中区分复杂的叠加异常 。近年来，针 和分布函数相同，还能保持与协方差函数或变异