化探数据处理

1、化探数据处理化探数据处理地调局方法地调局方法EDA方法方法 将地球化学背景看成平面：将地球化学背景看成平面：基于元素含量的统计模型，如，（对数）正态分布；基于元素含量的统计模型，如，（对数）正态分布；将数据转换为对数，去掉特高和特低含量后进行统计；将数据转换为对数，去掉特高和特低含量后进行统计；按数据服从多重分形分布出发，探索异常下限划分；按数据服从多重分形分布出发，探索异常下限划分；基于基于EDAEDA技术，求上异点作为异常下限的方法技术，求上异点作为异常下限的方法。 将背景面看成连续曲面：将背景面看成连续曲面：以构造单元划分子区，分别确定背景和异常下限；以构造单元划分子区，分别确定背景和异

2、常下限；用各种多元统计方法进行用各种多元统计方法进行“背景校正背景校正”，采用各种数据滤波技术去提取地球化学背景与异常，采用各种数据滤波技术去提取地球化学背景与异常，如，低通滤波方法等；如，低通滤波方法等； “ “子区中位数衬值滤波法子区中位数衬值滤波法” 将背景面看成不连续曲面：将背景面看成不连续曲面：目前常用的化探数据处理方法目前常用的化探数据处理方法铂族元素含量的高值区基本上与峨眉山玄武岩分布区相吻合，铂族元素含量的高值区基本上与峨眉山玄武岩分布区相吻合，如图所示，如图所示， PGEPGE地球化学背景并不是一个连续的曲面。地球化学背景并不是一个连续的曲面。新街新街大岩子大岩子铂族元素地球

3、化铂族元素地球化学背景面为不完学背景面为不完全连续的曲面全连续的曲面FC=CFC=Cm m/A/A0 0衬值滤波法示意图衬值滤波法示意图FC-FC-衬值衬值C Cm m- -中位数中位数A A0 0-EDA-EDA方法方法上异点值上异点值箱形分析图箱形分析图EDA技术技术 经统计，经统计，EDAEDA（Exploratory Data Analysis）技术获得技术获得的的PtPt、PdPd背景值和异常下限（上异点），其结果与基于背景值和异常下限（上异点），其结果与基于对数正态分布的数据处理结果类似。对数正态分布的数据处理结果类似。基于对数正态分布的背景值和异常下限基于对数正态分布的背景值和异

4、常下限 基于基于 EDA 技术技术 水系样品的水系样品的 基岩环境基岩环境 元元素素 样样 品品 数数 背景值背景值 C01.5S C01.65S C02S 中位数中位数 上异点上异点 Pd 1.66 5.16 5.78 7.53 1.63 5.39 三叠系三叠系 Pt 1296 1.69 4.95 5.51 7.07 1.64 5.41 Pd 2.92 10.2 11.56 15.48 2.87 13.72 石炭系石炭系 Pt 68 2.45 7.89 8.87 11.65 2.39 8.29 Pd 2.14 7.58 8.61 11.56 2.28 8.03 泥盆系泥盆系 Pt 95 1.

5、49 4.52 5.05 6.54 1.43 4.91 Pd 1.41 5.65 6.49 8.98 1.22 5.22 基性基性 超基性岩超基性岩 Pt 125 1.41 4.37 4.89 6.37 1.3 4.59 基于正态分布的经典数据处理方法与基于正态分布的经典数据处理方法与EDA方法处理结果对比表方法处理结果对比表 背景和异常的关系背景和异常的关系正异常正异常负异常负异常背景值：背景值：异常下限：异常下限：异常上限：异常上限：逐步剔除法选取背景数据范围：逐步剔除法选取背景数据范围：xC 0SxCSxCAA65. 165. 121SxxSx22数据取对数数据取对数计算程序C0-CA.

6、BASCu Pb Zn 0.027 0.25 1 0.03 0.32 2 0.03 0.54 4 0.03 0.65 7 0.03 1 9 0.03 1.4 11 0.03 1.9 23 0.03 2.1 45 0.03 3.2 67 0.03 4.3 89 0.03 5.4 94 0.03 6.2 103 0.03 7.1 245 原始数据输入格式：元素符原始数据输入格式：元素符号、数据从小到大顺序排列号、数据从小到大顺序排列 输出结果说明输出结果说明: :元素符号元素符号C0= C0= 背景值背景值 Ca1= Ca1= 异常上限异常上限 Ca2Ca2异常下限异常下限 背景区背景区: :（M

7、ax= Max= 最大值，最大值，Min= Min= 最小值）最小值） 取对数取对数: :（X1_X1_lnln= = 总平均值总平均值 S1_S1_lnln= = 样本方差）样本方差）X1=X1=总平均值总平均值 S1=S1=样本方差样本方差CvCv= = 变异系数（未取对数）变异系数（未取对数） EDA方法方法化探元素含量等值线图的绘制化探元素含量等值线图的绘制用用SurferSurfer软件绘图软件绘图地球化学图：等值线间隔约地球化学图：等值线间隔约0.10.1LgLg地球化学异常图：地球化学异常图：2 2n-1n-1C CA A1 1 2 21-1-n n C CA A2 2 n=1,2,3,4,5 n=1,2,3,4,5MAPGISMAPGIS软件环境下等值线图的制作软件环境下等值线图的制作生成生成MAPGISMAPGIS软件能够接受的软件能够接受的grdgrd文件文件在空间分析子系统中产生等值线，并保存在空间分析子系统中产生等值线，并保存为点、线、面文件格