化探数据处理

化探数据处理地调局方法EDA方法将地球化学背景看成平面：①基于元素含量的统计模型，如，（对数）正态分布；②将数据转换为对数，去掉特高和特低含量后进行统计；③按数据服从多重分形分布出发，探索异常下限划分；④基于EDA技术，求上异点作为异常下限的方法。将背景面看成连续曲面：①以构造单元划分子区，分别确定背景和异常下限；②用各种多元统计方法进行“背景校正”，③采用各种数据滤波技术去提取地球化学背景与异常，如，低通滤波方法等；“子区中位数衬值滤波法”将背景面看成不连续曲面：目前常用的化探数据处理方法铂族元素含量的高值区基本上与峨眉山玄武岩分布区相吻合，如图所示，PGE地球化学背景并不是一个连续的曲面。新街大岩子铂族元素地球化学背景面为不完全连续的曲面FC=Cm/A0衬值滤波法示意图FC-衬值Cm-中位数A0-EDA方法上异点值箱形分析图EDA技术经统计，EDA（ExploratoryDataAnalysis）技术获得的Pt、Pd背景值和异常下限（上异点），其结果与基于对数正态分布的数据处理结果类似。基于对数正态分布的背景值和异常下限基于EDA技术水系样品的基岩环境元素样品数背景值C0＋1.5SC0＋1.65SC0＋2S中位数上异点Pd1.665.165.787.531.635.39三叠系Pt12961.694.955.517.071.645.41Pd2.9210.211.5615.482.8713.72石炭系Pt682.457.898.8711.652.398.29Pd2.147.588.6111.562.288.03泥盆系Pt951.494.525.056.541.434.91Pd1.415.656.498.981.225.22基性超基性岩Pt1251.414.374.896.371.34.59基于正态分布的经典数据处理方法与EDA方法处理结果对比表背景和异常的关系正异常负异常背景值：异常下限：异常上限：逐步剔除法选取背景数据范围：xC0SxCSxCAA65.165.121SxxSx22数据取对数计算程序C0-CA.BASCuPbZn0.0270.2510.030.3220.030.5440.030.6570.03190.031.4110.031.9230.032.1450.033.2670.034.3890.035.4940.036.21030.037.1245原始数据输入格式：元素符号、数据从小到大顺序排列输出结果说明:元素符号C0=背景值Ca1=异常上限Ca2＝异常下限背景区:（Max=最大值，Min=最小值）取对数:（X1_ln=总平均值S1_ln=样本方差）X1=总平均值S1=样本方差Cv=变异系数（未取对数）EDA方法化探元素含量等值线图的绘制用Surfer软件绘图地球化学图：等值线间隔约0.1Lg地球化学异常图：2n-1CA121-nCA2n=1,2,3,4,5MAPGIS软件环境下等值线图的制作生成MAPGIS软件能够接受的grd文件在空间分析子系统中产生等值线，并保存为点、线、面文件格式通过坐标转换与相应的地形地质图相套合原始数据转换为带点位坐标的属性数据，通过属性表注释等功能生成实际材料图