资源描述:
书书书 第 3 4卷第 6期物 探 与 化 探V o l . 3 4 , N o . 6 2 0 1 0年 1 2月G E O P H Y S I C A L&G E O C H E M I C A LE X P L O R A T I O ND e c . , 2 0 1 0 区域化探数据处理中几种异常下限确定方法的对比 以内蒙古查巴奇地区水系沉积物为例 戴慧敏1, 宫传东2, 鲍庆中1, 孙中任1, 尤宏亮1, 金鑫1, 高飞1 ( 1 . 沈阳地质矿产研究所, 辽宁 沈阳 1 1 0 0 3 4 ; 2 . 沈阳市东陵区安全生产监督管理局, 辽宁 沈阳 1 1 0 0 0 0 ) 摘 要 采用传统统计、 稳健估计、 8 5 %累计频率和多重分形确定异常下限的计算方法, 对查巴奇地区水系沉积物数 据进行统计计算, 对取得的异常下限结果进行了对比。同一元素采用不同方法确定的异常下限存在明显的差异, 采用 8 5 %累计频率法除铅异常下限略低于多重分形法得到异常下限外, 其他 1 1个元素异常下限最低; 各个元素采 用其他三种方法确定的异常下限没有固定的规律性。这种结果与元素测量数据的分布特征及使用的统计方法密 切相关, 通过与研究区地质特征对比认为, 多重分形法确定的异常下限更适合于本研究区, 累计频率法确定异常下 限虽然较少遗漏低缓弱异常, 在一定程度上为地质背景的地球化学反映。因此, 在实际工作中要结合实际地质情 况, 在充分研究地质背景的基础上, 合理地确定了使用的计算方法。 关键词 区域化探数据; 异常下限; 传统统计法; 稳健估计法; 多重分形法; 8 5 %累积频率法 中图分类号P 6 3 2 文献标识码A 文章编号 1 0 0 0- 8 9 1 8 ( 2 0 1 0 ) 0 6- 0 7 8 2- 0 5 勘查地球化学方法找矿过程中进行地球化学异 常筛选, 所用的各种方法的一个共同目, 就是确定异 常下限, 以便研究区根据样品的浓度值进行分类。 随着计算机技术的发展及人们对地球化学数据分布 特征的认识, 发展了多种计算异常下限的方法, 如移 动平均法、 趋势面法、 克力格法、 概率格纸法和均值 加标准离差法等。这些方法的依据是假定地球化学 数据分布是服从正态分布的, 当数据不服从正态分 布时, 前人提出通过一些转换方法, 如对数正态变 换、 平方根变换等, 使其尽可能地服从正态分布[ 1 ]。 地球化学取样和对样品进行的各种化学分析结果常 具有不确定性, 而且元素在地壳中的分布本身就具 有不均匀性和区域随机性, 即区域地球化学数据服 从“ 重尾分布” , 这种分布也称分形或多重分形特 征[ 1 , 2 ], 这些数据即使经过变换处理, 仍然不能服从 正态分布。因此, 适应偏态数据分布的一些更合适 及解决不同地质单元间背景差异的统计方法便应用 于地质地球化学数据处理中, 同时, 这些方法在计算 过程和结果上在某些条件下有效性相对弱, 如对数 据有充分的了解且足以确定其分布类型时, 传统统 计方法相对稳健估计方法具有更强的针对性[ 3 ]。 笔者将应用于服从正态分布的传统统计方法与适用 于偏态分布的稳健估计法、 多重分形法、 8 5 %累计频 率法计算的异常下限进行对比, 以探讨不同方法获 得的异常下限之间的区别及意义。 1 研究区成矿地质概况 研究区位于内蒙古自治区与黑龙江省接壤部 分, 属于大兴安岭山地与松嫩平原过渡区, 伊尔施早 华力西地槽褶皱带与东乌珠穆沁旗晚华力西地槽褶 皱带接触部位, 被大兴安岭中生代火山岩带切割, 位 于兴安地槽褶皱系。大兴安岭火山喷发带东西两侧 分别由嫩江八里罕断裂和乌奴耳断裂隔开, 断裂 两侧分别为嫩江断陷盆地和海拉尔盆地。研究区位 于嫩江断裂带的西侧, 发育各个时期的 N E压性 压扭性构造和 N N W 向脆性断裂, 东侧为隆起带与 断陷盆地相联结处, 为嫩江八里罕断裂带, 沿嫩江 河谷延伸。区内地质构造复杂, 发育前中生代基地 构造和晚中生代断陷火山盆地, 晚侏罗世的侵入岩、 光华期的火山岩、 早燕山期和中燕山期的变质岩均 有分布( 图 1 ) , 其中地层由老至新为古生界、 中生界 及新生界。同时, 研究区在内蒙兴安岭成矿带( Ⅱ 级) 梨子山鄂伦春大兴安岭华力西期燕山 期铁铜钼金铅锌钨成矿带( Ⅲ级) 多宝山阿荣 收稿日期 2 0 0 9- 1 0- 1 6 基金项目 中国地质调查局矿调项目 内蒙古查巴奇地区矿产远景调查 ( 矿调[ 2 0 0 4 ] 1 0 2号) 6期戴慧敏等 区域化探数据处理中几种异常下限确定方法的对比 Q 第四系地层; K 1g 甘河组玄武岩; K1g n 2光华组二段安 山岩英安岩等; K 1g n 1光华组一段沉凝灰岩、 凝灰砂砾岩; T 1l 老龙头组变安山岩; P3l n 林西组砂岩; J3k γ 燕山期碱 长花岗岩; J 3η γ 燕山期二长花岗岩; J1δ o 燕山期中粒石英 闪长岩; J 1γ o 燕山期中粒英云闪长岩; 1 脉岩( γ π 花岗斑 岩脉; δ μ 闪长玢岩脉) ; 2 实测断裂; 3 推测断裂 图 1 查巴奇地区地质概况 旗华力西期燕山铜钼金成矿带( Ⅳ级) 阿荣旗 扎兰屯铜钼金成矿远景区( Ⅴ级) 内。 2 异常下限确定方法 2 . 1 传统统计方法 传统的计算异常下限的方法是建立在数据符合 正态或者对数正态分布基础上的, 故如果测量数据 符合正态分布, 求出背景值和标准离差计算异常下 限。对于测量数据不符合正态分布的情况, 计算异 常下限的首要步骤是对所取得的化探数据进行常规 的数据处理,即进行离群点( 最高值、 最低值) 的迭 代处理, 一般采用 > X+ 3 S及 < X- 3 S迭代剔除, 直 至无离群点数值可剔除为止, 即所有数据全部分布 在 > X+ 3 S 与 < X- 3 S 间, 形成背景数据, 以背景值 加 2倍标准离差为异常下限。 2 . 2 稳健估计法 稳健估计法是对偏态分布的数据建立更优的数 学模型以使估计更为稳健, 其性能对模型的微小变 化反映不敏感, 即当样本中混入少量的异常值时, 对 统计量的值影响不大。稳健估计的原则是要充分利 用观测数据( 或样本) 中的有效信息。一个稳健的 统计过程需要具备以下特征 在假定分布形式下, 该 过程应该具有最优性或接近最优性; 当实际分布与 假定分布有少许出入时, 该过程应该具有一定程度 的稳健性, 即只允许对结果有稍微的不利影响; 当实 际分布与假定分布相差很悬殊时, 该过程计算的结 果不至于产生“ 灾难” 性的影响。在稳健估计学中, 位置( l o c a t i o n ) 及尺度( s c a l e ) 的估计是两个基本问 题, 相当于传统统计学中的均值与均方差。位置估 计一般包括 a 切尾均值法、 H u b e r 估计和 A n d r e w s 估 计, 尺度估计包括中位绝对中位差法和 H u b e r 法, 这 里使用 5 %切尾均值法进行位置估计采用中位绝对 中位差法进行尺度估计, 计算原理及方法如下 ( 1 ) 位置估计 X= 1 ( 1-2 a ) n n i=1 a ixi, 0<a<0 . 5, 当 a n + 1< i < n - a n时, a i= 1 ; 当 i =a n+ 1或 i =n - a n 时, a i= p ; 当 i <a n+ 1或 i >n-a n时, ai= 0 。 式中, X为位置估计; x i为元素分析结果; p = 1+a n - a n ; a n 为 a n 的整数部分; n 为样品总数; a = 5 %。 ( 2 ) 尺度估计 S=m e d [ x i-m e d ( xi) / 0 . 6 4 5 ] S 为相对于中位数偏离的绝对值的中位数, 即为尺 度估计, x i为元素分析结果。根据上述方法所得位 置及尺度, 以各元素位置( 背景值) 加 2倍离散尺度 ( 均方差) 作为异常下限。 2 . 3 多重分形法 分形理论以自然界和社会活动中广泛存在的无 序( 无规则) 而具有自相似性的系统为研究对象, 提 取出确定性、 规律性的参量, 揭示复杂事物中新的深 刻而定量的规律。目前利用分形技术进行地球化学 异常下限确定的方法主要有含量周长法、 含量 面积法、 含量距离法、 含量频数法等, 这里采用 含量频数法。 设分形模型 N ( r )=C r - D, r>0, ( 1 ) 式中, r 为元素含量对数, 又称特征尺度; C为比例常 数, 且 C> 0 ; D为元素分维, 且 D> 0 ; N ( r ) 为含量的 387 物 探 与 化 探3 4卷 对数大于等于 r 的频数, 又称尺度。 如果元素服从多重分形分布, 则存在如下关系 N ( r ≤ T ) ∞r - D1; N ( r ≥ a ) ∞r- D2, ( 2 ) 式中, T为区域异常下限。该式表示区域上元素分 布在空间上存在层次性特征, 即存在背景分形分布 和异常分形分布。具体方法是, 将研究区的每个元 素进行含量排序, 按照 0 . 1 l g X间隔统计样品的频 数, 在含量和频数的双对数坐标中绘制散点图( 图 2 ) , 用拟合方法求出无标度区元素的分维, 式( 2 ) 就 是含量频数分形计算公式。 图 2 A u及 Z n含量频数双对数曲线 2 . 4 8 5 %累积频率法 累计频率就是一个数值的频率和比它的频率高 的数值的频率的总和。设 x 1<x2< <xm是不重 复的样本值, m <n 。把样本值小于或等于某个样 本的数据 x i的频率累加, 就得到小于或等于 xi的累 积频率。笔者使用累积频率 8 5 %时的 x i值作为异 常下限, 得到查巴奇地区水系沉积物的异常下限 ( 表 1 ) 。 表 1 查巴奇地区水系沉积物数据采用不同方法确定的异常下限值对比 统计方法A uA gA sB iC uH gM oP bS bS nWZ n 传统方法2 . 81 7 03 60 . 9 53 53 543 82341 1 0 稳健估计法1 . 61 7 0 . 72 0 . 2 60 . 5 83 4 . 6 82 1 . 2 22 . 7 54 8 . 0 21 . 23 . 4 53 . 5 51 1 3 . 7 多重分形法3 . 01 4 02 813 81 833 01 . 32 . 53 . 61 0 0 8 5 %累积频率法1 . 3 11 2 12 00 . 4 12 8 . 31 6 . 22 . 2 93 11 . 0 32 . 3 12 . 3 58 3 注 含量单位为 w ( A u 、 A g ) / 1 0 - 9, w ( 其他元素) / 1 0- 6 3 对四种方法确定异常下限的认识 应用四种方法对查巴奇地区水系沉积物数据计 算确定异常下限。从表 1中看到, 同一元素采用不 同的方法确定的异常下限存在明显的差异 只有铅 在使用多重分形法得到的异常下限略低于 8 5 %累 计频率法外, 其他所有元素使用 8 5 %累计频率法得 到的异常下限最低; 各个元素间在使用传统统计法、 多重分形法和稳健估计法得到的异常下限没有固定 的规律, 这与元素测量数据的分布特征及统计方法 密切相关。地球化学场元素的分布形式可归纳为两 大类 正态和多模式分布、 分形和多重分形分布, 正 态和多模式分布对应于混合作用占主导地位的地球 化学过程, 而分形和多重分形分布对应于分异作用 占主导地位的地球化学过程[ 4 ]。传统统计方法的 前提是数据必须处于正态分布或对数正态分布, 因 此, 对于其他分布形式的数据则采取剔除离群点数 据直至使数据服从正态分布。离群点数据(>X+ 3 S 及 < X- 3 S ) 不参与统计, 其确定的异常下限值并 非由全部数据而只是由背景数据决定, 这种方法对 度量一般值效果较好; 稳健估计法目标旨在地球化 学异常评价中对偏态分布的数据总体建立较优的数 学模型, 以使特异值的影响相对减小, 但其 5 %切尾 均值法也存在类似传统方法的缺陷性, 从计算式中 看到, 当 i < a n+ 1或 i >n-a n时, a i= 0 , 即元素含 量排序后, 这两个区间的数据不参与统计; 同时在确 定异常下限值使用的参数与传统方法也不同, 传统 统计方法使用剔除离群数据的 X 2 S ( 或 2 . 5 S ) 作 为异常下限值, 稳健估计法是剔除“ 重尾” 数据后使 用位置( 背景值) 加 2倍离散尺度( 均方差) 作为异 常下限, 在剔除离群点数据或“ 重尾” 数据过程中, 两种方法剔除数据量也是不同的。因此, 稳健估计 法确定的某些元素异常下限比传统方法高, 而某些 元素比传统方法低。 综上所述, 基于正态分布及考虑数据偏态分布 的异常下限计算方法中, 传统方法和稳健估计法原 理上有雷同点, 分形理论确定异常下限的过程充分 利用了所有测量数据, 不像传统的只是一个数值; 8 5 %累计频率法同样充分利用了景观区的全部数据 信息, 将景观区 1 5 %的数据量作为异常数据来评 487 6期戴慧敏等 区域化探数据处理中几种异常下限确定方法的对比 价, 从表中可以看出, 除 P b 外的1 1种元素累积频率 和分形方法确定的异常下限差别明显。元素含量在 地质体中的分布及统计方法的原理差异双重因素决 定了四种异常下限值的不同。以 A u为例, 8 5 %累 计频率和稳健估计法确定异常下限的结果比较接 近, 多重分析法和传统统计法确定异常下限相对接 近, 利用四种方法确定异常下限分别作异常图( 图 3 ) 。从图可以看出, 稳健估计法和 8 5 %累频法做出 的异常面积明显大于其他两种方法确定的异常, 位 置及形态基本上覆盖林西组浅变质岩, 与矿床分布 位置吻合较好, 更低程度地遗漏低缓的弱异常, 但这 两种方法圈定的异常面积较大, 给异常查证带来了 难度; 传统方法和多重分形法确定异常下限得到的 异常图异常比较零星、 分散, 对比地质图可以看出, 异常主要分布于林西组砂岩中, 研究区北部异常点 较多, 南部则呈面积分布。 图 3 查巴奇地区四种异常下限 A u异常对比 4 结论与讨论 通过查巴奇水系沉积物异常下限计算, 对比传 统统计法和稳健估计法, 可以得出 对数据总体分布 有充分的了解且足以确定其分布类型时, 稳健估计 方法不如传统统计方法具有更强的针对性, 有效性 相对较差, 有些元素运用稳健估计法确定的异常下 限比传统方法大, 有些却比传统方法小; 分形理论确 定异常下限时一般分为两段,不像传统的方法只是 一个数值, 从方法原理看, 多重分形法和累计频率法 利用了所有测量数据, 其计算的异常下限也更具合 理性, 但是从本研究区来看, 累计频率法确定的异常 下限相对于其他方法最低, 一定程度上为地质单元 的地球化学背景的反映, 给异常查证带来困难。综 合各个元素四种异常下限对比及地质构造分析, 研 究区采用多重分析法确定异常下限更合理。不同的 元素使用不同的方法确定的异常下限值不同, 某一 方法对于某种元素而言, 确定的异常下限值是最低 的, 能将一些弱异常反映出来, 但相反也会给异常查 证带来困难; 同时, 不同的地质单元一般具有不同的 地球化学背景, 在分布有多个地质单元的地区开展 化探工作时, 还要考虑地质单元的地球化学背景差 异而不能在全区采用统一的异常下限, 如采用移动 平均法、 子区中位衬值滤波法等确定异常下限。因 此, 在实际工作中计算异常下限首先要考虑地球化 学数据的分布模式, 其次要与实际地质情况相结合, 充分研究数据分布模式和地质背景的基础上, 才能 确定使用哪种计算方法更合理。 参考文献 [ 1 ] 李蒙文, 战明国, 赵财胜. 稳健估计方法在内蒙古新忽热地区 水系沉积物测量异常评价中的应用[ J ] . 矿床地质, 2 0 0 6 , 2 5 ( 1 ) 2 7- 3 5 . [ 2 ] 孙忠军. 矿产勘查中化探异常下限的多重分形计算方法[ J ] . 物探化探计算技术, 2 0 0 7 , 2 9 ( 1 ) 5 4 . [ 3 ] 王星. 非参数统计[ M ] . 北京 中国人民大学出版社, 2 0 0 5 2 8 . [ 4 ] A l l e g r eCJ ,L e w i nE . S c a l i n gl a w s a n dg e o c h e m i c a l d i s t r i b u t i o n s [ J ] . E a r t ha n dP l a n e t a r yL e t t e r s , 1 9 9 5 , 1 3 2 1 . 587 物 探 与 化 探3 4卷 AC O MP A R I S O NO FS E V E R A LT H R E S H O L DD E T E R MI N A T I O NME T H O D SI N G E O C H E MI C A LD A T AP R O C E S S I N GAC A S ES T U D YO FS T R E A M S E D I ME N T SI N C H A B A Q I A R E AO FI N N E RMO N G O L I A D A I H u i m i n 1, G O N GC h u a n d o n g2, B A OQ i n g z h o n g1, S U NZ h o n g r e n1, Y O UH o n g l i a n g1, J I NX i n1, G A OF e i1 ( 1 .S h e n y a n gI n s t i t u t e o f G e o l o g ya n dM i n e r a l R e s o u r c e s ,S h e n y a n g 1 1 0 0 3 4 ,C h i n a ; 2 .P r o d u c t i o nS a f e t y S u p e r v i s i o na n dA d m i n i s t r a t i o no f D o n g l i n g D i s t r i c t ,S h e n y a n g ,S h e n y a n g 1 1 0 0 0 0 ,C h i n a ) A b s t r a c t T h ea u t h o r s m a d es t a t i s t i cc a l c u l a t i o no f s t r e a ms e d i m e n t d a t au s i n gs u c hm e a n s a s t r a d i t i o n a l s t a t i s t i c s ,r o b u s t e s t i m a t i o n , c u m u l a t i v ep e r c e n t a g eo f 8 5 % a n dm u l t i f r a c t a l s t a t i s t i c a l m e t h o di nC h a b a q i a r e aa n dc o m p a r e dt h et h r e s h o l d so b t a i n e db yt h ef o u r m e t h o d s .T h e t h r e s h o l d s o f t h e s a m e e l e m e n t s o b t a i n e db y t h e f o u r m e t h o d s a r e o b v i o u s l y d i f f e r e n t f r o me a c ho t h e r .T h e t h r e s h o l d s c a l c u l a t e db y t h e c u m u l a t i v e p e r c e n t a g e s t a t i s t i c s a r e l o w e s t e x c e p t f o r l e a d ,a n dt h e t h r e s h o l d s c a l c u l a t e db y t h e o t h e r t h r e e m e t h o d s f a i l t os h o wo b v i o u s r e g u l a r i t y . T h er e s u l t s a r ec l o s e l yr e l a t e dt ot h ed i s t r i b u t i o nc h a r a c t e r i s t i c s o f t h eg e o c h e m i c a l d a t aa n dt h es t a t i s t i c a l m e t h o d s u s e d .C o m p a r i n gt h ea n o m a l yc h a r a c t e r i s t i c s w i t ht h eg e o l o g i c a l c h a r a c t e r i s t i c s ,t h ea u t h o r s h a v ef o u n dt h a t ,a l t h o u g hw e a k a n dg e n t l eg e o c h e m i c a l a n o m a l i e s w i l l b em o s t l yd e t e c t e db yc u m u l a t i v ep e r c e n t a g eo f 8 5 %,t h el a r g e r a n o m a l ya r e a s w i l l b r i n gd i f f i c u l t i e s t oa n o m a l y i n s p e c t i o ni nt h e f i e l d .I t i s t h e r e f o r e t h o u g h t t h a t t h e t h r e s h o l d s o b t a i n e db y m u l t i f r a c t a l s t a t i s t i c a l m e t h o da r e s u i t a b l ef o r t h es t u d ya r e a .Ac o m p a r a t i v es t u d yo f t h et h r e s h o l d s h a s l e dt h ea u t h o r s t ob e l i e v et h a t t h er e a s o n a b l em e t h o dm u s t b ec o m b i n e dw i t ht h ea c t u a l g e o l o g i c a l c o n d i t i o n s a n db a s e do nat h o r o u g hs t u d yo f t h eg e o l o g i c a l b a c k g r o u n do f t h es t u d ya r e a . K e yw o r d s C h a b a q i a r e ao f I n n e r M o n g o l i a ;r e g i o n a l g e o c h e m i c a l d a t a ;t h r e s h o l d s ;t r a d i t i o n a l s t a t i s t i c s ;r o b u s t e s t i m a t i o n ;m u l t i f r a c t a l s t a t i s t i c a l m e t h o d s ;c u m u l a t i v ep e r c e n t a g eo f 8 5 %;s t r e a ms e d i m e n t s 作者简介 戴慧敏( 1 9 7 9- ) , 地球化学专业硕士, 工程师。主要从事地球化学研究 檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵 。 上接 7 8 1页 A b s t r a c t A f t e r t h e 8 . 0m a g n i t u d e e a r t h q u a k e t h a t o c c u r r e di nWe n c h u a nC o u n t y ,S i c h u a nP r o v i n c e , t h e a u t h o r s e s t a b l i s h e dm o n i t o r i n g p o i n t s o f a f t e r s h o c ki nN i n g q i a n g C o u n t y ,S h a a n x i P r o v i n c e .T h e r e l a t i o n s h i pb e t w e e nt h e a f t e r s h o c k s a n dt h e R n H g c o n t e n t c h a n g e i n s o i l g a s w a s o b s e r v e do nt h es p o t a n da t r e g u l a r i n t e r v a l s i n1 7d a y s .T h er e s u l t s s h o wt h a t t h es o i l g a s m e r c u r yc o n t e n t i s g e n e r a l l y m o r e t h a n 2 0 0n g/ m 3i nt h e h a n g i n g w a l l o f t h e e a r t h q u a k e f a u l t , t h e c o n t e n t v a l u e i s t w o t i m e s h i g h e r t h a nt h a t o f t h e f o o t w a l l ( I nf r a c t u r e s a l o n gt h es e c t i o n , m e r c u r ya n o m a l yv a l u ec a nr e a c h5 6 0n g / m 3) , b u t r a d o nc o n t e n t v a l u es h o w sn od i s t i n c t c h a n g e .A l o n gt h e s t r i k eo f t h e s e i s m o g e n i c m a i nf a u l t ,t h e c l o s e r i t i s t o t h e e p i c e n t e r ,t h e h i g h e r t h e R n ,H g b a c k g r o u n dv a l u e s a r e .I f m e r c u r y c o n t e n t i s m o r et h a n2 0n g / m 3i nt h eb a c k g r o u n ds e c t i o n ,t h e r ei s R na n o m a l ys i m u l t a n e o u s l y ,w h i c hI n d i c a t e s t h a t t h e r em i g h t o c c u r h i g h e r p r o b a b i l i t y o f a f t e r s h o c kw i t hh i g h f r e q u e n c ya n do v e r 4 . 0m a g n i t u d ew i t h i n3 0 0k ma n d4~ 4 8h o u r s ;i f t h i n g s a r en o t l i k et h a t ,t h e a f t e r s h o c kw i l l n o t o c c u r . K e yw o r d s a f t e r s h o c k s s u r v e i l l a n c e ; s o i l g a s ; m e r c u r ya n dr a d o nc o m b i n e dm e a s u r e m e n t ;N i n g q i a n gi nS h a a n x i 作者简介 杨少平( 1 9 5 3- ) , 男, 陕西华阴人, 1 9 7 8年毕业于北京大学地质地理系地球化学专业, 1 9 9 4年获中国地质大学( 武 汉) 地球化学硕士学位, 研究员, 长期从事勘查地球化学基础理论和方法技术研究工作。 687
展开阅读全文