基于MAPGIS的含量—面积法确定异常下限——以个旧西区为例-论文.pdf

《基于MAPGIS的含量—面积法确定异常下限——以个旧西区为例-论文.pdf》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、第35卷第4期物探化探计算技术2013年7月文章编号：1001—1749(2013)04一0477—04基于MAPGIS的含量一面积法确定异常下限——以个旧西区为例杨震，梅红波(中国地质大学资源学院，武汉430074)摘要：含量一面积法是基于分形理论，一种研究含量与曲面面积之间关系的方法。这里以云南省个旧西区某地Sn元素土壤化探数据为例，利用MAPGIS平台进行化探数据处理，研究其含量与面积的关系，并确定异常下限，圈定异常范围。同时与传统方法进行对比分析，结果表明，基于分形含量～面积法圈定的异常区域更加有效、合理。关键词：分形；含量一面积法；MAPGI

2、S平台；异常下限中图分类号：P632文献标识码：ADOI：10．3969／j．issn．1001—1749．2013．04．18区sn元素的异常下限值，同时与传统方法进行对0前言比分析，探讨其可行性与合理性。元素地球化学异常下限值的确定，是地球化学I含量一面积分形模型原理研究中的重要问题之一[1]，通常元素地球化学异常下限值的确定应用于金属矿产的地球化学勘探[1]，分形是指局部与总体具有某种相似性的形状，为地质找矿提供科学依据。长期以来，人们都采用或者是在不同尺度上看起来基本相似的形状，即具经典的统计学方法，以样品数据呈正态分布或者近有自相似性(sel

3、f—similarity)，而自相似性就是局似正态分布为前提，通过计算其统计学参数(如均部与整体在形态、功能和信息等方面具有统计意义值和标准偏差)，对数据进行筛选和评价。我国目上的相似i7。。定量描述这种自相似性的参数称为前采用的地球化学勘探标准中规定，均值与二倍的“分维数”，用D表示。分形分布的特点要求大于标准偏差作为异常下限[2~4]。这种方法忽视了元等于某一尺度的数目与物体大小之间存在幂指数素的空间分布信息，且仅仅适用于元素含量在空间关系，由此可得到含量一面积分形统计模型，即上呈理想情况下的正态分布，但是有很多元素，特A(，-)一Cr—D，r>0

4、(1)别是微量元素并不服从正态分布，而是呈明显的正式中，．表示元素含量值；A(r)表示元素含量大于向偏斜，或表现为一种幂型的拖尾分布[5]。李长江等于r的面积值；C>0称为比例常数；D>0称为等[6]研究揭示的地球化学景观，可能是一个具有低分维数。维(D=2．9)吸引的混沌系数，具有典型的分形特分形方法的关键在于确定分维数，具体思路是征，因此可以作为分形来处理。将公式A(r)=Cr加二端同时取对数，得出线性作者在本文选用分形中含量一面积法，并结合回归方程(2)。MAPGIS平台的属性统计分析功能，确定出研究lgA(r)一一Dlgr+lgC(2)收稿日期

5、：2012一08—27改回日期：2013～03—26478物探化探计算技术35卷将数据[A(r。)，A(rz)，⋯，A(^)]和(r1，r2，文件，并将其命名为Sn．txt。⋯，rn)分别代入公式(2)中，在GRAPHER软件(2)在MAPGIS平台上，首先打开空间分析模中，采用最小二乘法进行直线拟合，从而求出斜率块中的DTM子模块，对刚生成的Sn．txt文件进D的估计量。一般情况下，散点会分别在两段或者行网格化。网格化方法选用泛克里格法，输出网格多段直线上，这时可以采用分段拟合方式。由于人文件Sn．Gr&然后打开“平面等值线图绘制”菜为确定分段拐点所

6、造成的误差较大，为了提高精单，选择网格数据Sn．Grd，设置相应的参数，从而度，可以采用最优法确定分界点。散点分布在两段生成一个区文件Sn．WP。直线上的情况下如何确定分界点，基本思路是找出(3)打开MAPGIS中的空间分析子模块，装入合适的分界点‰，使二个区间拟合的直线与原始区文件Sn．WP；打开“属性分析”菜单下的“双属数据之间的残差平方和Ei(i一1，2)在各区间的总性分类统计”子菜单，其中分类属性字段选择起始和以式(3)表示为最小。值，分类方式选择分段方式，统计属性字段选择面i0积，统计方式选择累计方式；然后将分析结果保存E—E1+E2一∑11

7、9A(n)+D1lgri—lgCl]2+i一1为Zn．wB；最后在MAPGIs库管理模块下的属月性库管理子模块中导人Zn．WB，点击“输出属性”∑[1醇(t)+D2lgn—lgC2]2(3)20+1菜单，导出Excel能够识别的数据文件Zn．DBF。其中‰为分界点；D。和D。分别为相应区间的(4)在Excel软件中，打开Zn．DBF，逐个累加斜率，即分维数。面积值，其原因是采用MAPGIS双属性分类统计为了检验回归方程的显著性，需对每个回归方法统计出来的面积，只能表示二个分段点之间的面程进行方差分析检验和相关系数检验。该方法也积，并不能表示大于等于某一

8、分段点的所有面积之可以推广到散点分布在多段直线上的情况。和。这样就町以得到一系列不同的r值及其