化探数据处理化探数据处理地调局方法地调局方法EDA方法方法2021/6/71 将地球化学背景看成平面：将地球化学背景看成平面：基于元素含量的统计模型，如，（对数）正态分布；基于元素含量的统计模型，如，（对数）正态分布；将数据转换为对数，去掉特高和特低含量后进行统计；将数据转换为对数，去掉特高和特低含量后进行统计；按数据服从多重分形分布出发，探索异常下限划分；按数据服从多重分形分布出发，探索异常下限划分；基于基于EDAEDA技术，求上异点作为异常下限的方法技术，求上异点作为异常下限的方法。 将背景面看成连续曲面：将背景面看成连续曲面：以构造单元划分子区，分别确定背景和异常下限；以构造单元划分子区，分别确定背景和异常下限；用各种多元统计方法进行用各种多元统计方法进行“背景校正背景校正”，采用各种数据滤波技术去提取地球化学背景与异常，采用各种数据滤波技术去提取地球化学背景与异常，如，低通滤波方法等；如，低通滤波方法等； “ “子区中位数衬值滤波法子区中位数衬值滤波法” 将背景面看成不连续曲面：将背景面看成不连续曲面：目前常用的化探数据处理方法目前常用的化探数据处理方法2021/6/72铂族元素含量的高值区基本上与峨眉山玄武岩分布区相吻合，铂族元素含量的高值区基本上与峨眉山玄武岩分布区相吻合，如图所示，如图所示， PGEPGE地球化学背景并不是一个连续的曲面。地球化学背景并不是一个连续的曲面。新街新街大岩子大岩子铂族元素地球化铂族元素地球化学背景面为不完学背景面为不完全连续的曲面全连续的曲面2021/6/73FC=CFC=Cm m/A/A0 0衬衬值值滤滤波波法法示示意意图图FC-FC-衬值衬值C Cm m- -中位数中位数A A0 0-EDA-EDA方法方法上异点值上异点值2021/6/74箱形分析图箱形分析图EDA技术技术2021/6/75 经统计，经统计，EDAEDA（Exploratory Data Analysis）技术获技术获得的得的PtPt、PdPd背景值和异常下限（上异点），其结果与基背景值和异常下限（上异点），其结果与基于对数正态分布的数据处理结果类似。于对数正态分布的数据处理结果类似。基于正态分布的经典数据处理方法与基于正态分布的经典数据处理方法与EDA方法处理结果对比表方法处理结果对比表 2021/6/76背景和异常的关系背景和异常的关系正异常正异常负异常负异常背景值：背景值：异常下限：异常下限：异常上限：异常上限：逐步剔除法选取背景数据范围：逐步剔除法选取背景数据范围：数据取对数数据取对数2021/6/77计算程序C0-CA.BAS原始数据输入格式：元素符号原始数据输入格式：元素符号、数据从小到大顺序排列、数据从小到大顺序排列 输出结果说明输出结果说明: :元素符号元素符号C0= C0= 背景值背景值 Ca1= Ca1= 异常上限异常上限 Ca2Ca2异常下限异常下限 背景区背景区: :（Max= Max= 最大值，最大值，Min= Min= 最小值）最小值） 取对数取对数: :（X1_ln= X1_ln= 总平均值总平均值 S1_ln= S1_ln= 样本方差）样本方差）X1=X1=总平均值总平均值 S1=S1=样本方差样本方差Cv= Cv= 变异系数（未取对数）变异系数（未取对数） 2021/6/78EDA方法方法2021/6/79化探元素含量等值线图的绘制化探元素含量等值线图的绘制用用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,52021/6/710MAPGISMAPGIS软件环境下等值线图的制作软件环境下等值线图的制作生成生成MAPGISMAPGIS软件能够接受的软件能够接受的grdgrd文件文件在空间分析子系统中产生等值线，并保存在空间分析子系统中产生等值线，并保存为点、线、面文件格式为点、线、面文件格式通过坐标转换与相应的地形地质图相套合通过坐标转换与相应的地形地质图相套合原始数据转换为带点位坐标的属性数据，原始数据转换为带点位坐标的属性数据，通过属性表注释等功能生成实际材料图通过属性表注释等功能生成实际材料图2021/6/711部分资料从网络收集整理而来，供大家参考，感谢您的关注！