联亩他·贾第卷第期年月区域化探异常研究方法探讨李武俊唐开金陕西省地质矿产局物化探队摘要本文对区域化探异常图示、异常评序、可能成矿元素、综合解释及成矿预测等方面进行了探索性研究,这对异常深入的刽析研究,确定异常检查范围及成矿预测都十分有益。数据处理系统是在无一绘图仪的微机上实现数据处理、成图及各类处理分析的程序软件。关键词变异度叠加强度离散度异常结构找矿模式地球化学图及说明书是区域化探扫面的重要成果,是地质找矿、成矿预测、基础地质、理论地球化学、农、林、牧、环境保护等各个领域的基础资料。近儿年来,笔者在区域化探资料研究中,除用常规的方法外,在数据处理、异常图示、异常评序、可能成矿元素的探讨、综合异常解释及成矿预测方面作了一些探索性研究,取得了较好效果。区域化探的数据处理根据《区域化探全国扫面工作方法若干规定》,我们在微机上设计了万区域化探数据处理系统但无一绘图仪,该系统的内容如图所示也适用于万区域化探,采用微机处理与手工相结合的方法来完成地球化学编图、各类计算及综合研究。经新狠、陕西共十个万图幅的全部数据处理结果表明,能达到万区域化探编图要求。在微机上利用区域化探数据处理系统处理数据有以下几方面的优点各元素数据图经检查无误后即可入库进行各类处理,实现了数据处理一入库一校公‘处理一成图的作业流程,避免了等待整图幅的所有数据整理一校核完成后去外地计算作图,多数人等待的弊病。原需用手工计算的色区等值线、元素化合物地球化学参数、全域及各地质单元李武俊女,岁,年毕业于成都地质学院工程师参与主编的《航磁么等值线图》获科技成果三等奖并今与编写过多份物化探研究报告店开金男,“岁,年毕业于成都地质学院,高级工程师。曾编写过多份物化探普、详查报告,陕西地质第卷样品编码图原始数据图一二止厂质量评序竺始数据入库数据检索隆陇隆曦陌味博圈图国国圆图国圈图口逐步回方归程回代求测值预累加累乘晕计算异构结计常图算多元统计分析成元矿素计参数统直方图参数统计下鑫移动平均分析、…今一色区计算一’一布打印地化字符图一下否一程手工勾绘等值线图统计剩余异常值计算一二—统计翼薰篡异异异常常常评评评序序序山冈因回圆向圈甲日圆去自周回,,‘‘匕图数据整理及微机处理流程示意图!沮。的参数统计表、各地球化学分区中主要元素地球化学参数对比表等,现全部用微机计算,大大减轻了手工劳动,节省人力和时间。各类综合研究的数据处理方便,各类研究方法处理易于实施、修改,研究的方法获得的信息量大,为深入研究异常创造了条件。利用我队微机及在外单位处理一个完整的万图幅化探资料的经济效益对比如多表第期区域化探异常研究方法探讨表万图幅化探资料经济效益对比表,项目本队外单位需用人数各类参数计算综合研究处理人年机器自动处理随时可以计算元一人年元元一元二元人年手工计算去外单位计算。。元含差旅费无用用费费时约机节从表中可以明显看出,用本队微机处理数据不仅减轻了大量的手工劳动,而且可节约万元的费用。可能成矿元素的选择方法及效果区域化探样品分析元素多,化探工作者最关心的是哪些元素可能成矿解决这个问题的方法很多,现仅将我们实际工作中探索的方法介绍如下。在地球化学找矿中,通常根据元素含量高低、变异程度的强弱、离散程度大小来推断元素成矿的可能性,依据这一原则,我们选择如下参数变异度二。,叠加强度二天。离散度二。,综合指标变异度叠加强度离散度式中,、,、,分别为全图幅某元素的网格化原始数据经迭代剔除后计算的变异系数、平均值、标准离差。。、。、。分别为全幅图某元素网格化原始数据未经迭代剔除的达异常值的数据的变异系数、平均值、标准离差。综合指标值越大,元素成矿的可能性就越大,以此值来选择可能成矿的元素。在陕西的万区域化探编图的资料研究中,经多个图幅试用,取得了较好的效果。如图,是一个图幅个元素以综合指标的大小排序结果,其中为已知的、矿床,为大型金矿床,为大型的铝矿床,这三个元素的综合指标值都远远大于其它元素,其效果是比较明显的。我们认为,在异常带、成矿远景区、特别是工作程度较低的地区,应用此方法确定可能成矿的元素更为优越。异常评序区域化探扫面工作中获得大批区域化探异常,短期内如何从众多异常中获得找矿成效,关键是如何筛选出最有找矿意义的异常。而异常评序是目前筛选异常的重要途径之一,评序陕西地质第,卷院合指标尸卜劫‘一,心‘一人己图陕西相万区化成矿元紊选实例一!∀方法较多,口前尚未见到统一模式,现就我们应用的方法作一简介。单元素异常评序目前在于将研究区内某元素异常,按其找矿意义进行排序。勘查地球化学的实践证明,一定规模的矿床所引起的异常总是具有一定的面积、且浓度分带好为特征。据此,我们采用值为元素异常面积乘以衬度对单元素异常进行评序。综合异常评序以往常采用直观法分析综合异常,并结合地质、构造及矿产资料来评价综合异常,此方法往往有使用者的主观性,为克服主观性、较真实地反映各元素的地球化学特征,采用了综合异常登记卡片中的每个单元素异常的面积、衬度、浓度分带三个参数,按下列公式矿化度一凰”‘’署丹平均规模艺平均衬度艺,第期区域化探异常研究方法探讨—一一一一一一一一一一一一」一一一一,一一一一一一一平均浓度=艺LD、/N式中N、S、T、L。分别为每个综合异常中单元素异常数、面积、衬度及浓度值。以此求得每一个综合异常的上述四个参数,再将每一个参数由大至小排序,然后将矿化度序号、平均规模序号、平均衬度序号、平均浓度序号相加求得每一个综合异常的综合序号,按综合序号由小至大将异常排序(见表2)。经多个1:20万图幅应用,均能起到筛选异常、指导异常检查工作之目的。衰2区城化探异常评序结果Table2.Theresultsofevaluationsequeneeofthereg三onallygeoehemiealProsPeetionanomalies.序号}异常编号霸不赢卿画奎颤全巨鲤翌叁竺元素顺序HS480O00004。002.003。007,005002.004。OU13。0015.00乙丙丙丙丙丙丁乙丁了乙丙乙丙丁HS39HS44HS426.002.003.0018.00000000000000196C匀80nonno‘:60口nJ八0On0OC0.…飞0九O只...工,1HS47HS58HS38HS54一HS45一aHS51HSS连一2HS45一bHS555.004008.001l.OU9.0013.0010.007.0014.0011,007.0027.0029.0013。009.001户.0021.004.0016.005.0013.0012.0034。0036.0037。0038.0D38.0047.0048.00HS4x12.002.0014.0016。0054.002345归了s。比1112比1415.HS4317.0014.0013_0015.0059.00MoCdVNiUZnAgHgCuBaAuAsVMoCdAgBaZnAsSbCuVAgSbBaAsHgCrCdCoZnCuTiSbWBaCdAsCrVNbHgPbZnAgAgVCdSbAuZnMoHgMoUAgNbAsZrVCuTiHgCdAgVMoWAsSbZnMoAgCrCdBaHgCoMoZnUNiCdCuAuNbAsCrAuHgBaNbUCdCrZnNiVAuAgAuWHgMoCrTiCtzUVCdAgZnNiNbZrCoW贪犷手公货愁万‘NbZrCrWCOSbPbNiCrCdTiWA比Zn飞『CuZrWSbMoAgCdAsNiAuZnBaCrCdAuZnWAgMoBaPbAsSbCdCrAsMo人uFbAuPbZ吐SbAgZrNbTiZnCuVCoWAsSbMoCdBaAgUPbZnZrNbWMoCoPb日00U‘U八甘UO八U0000001八U00626570”的8286579500000000000000000014171820192421222300000 UnU00000UU0CUCUO1512[719邪2022咒.240OCnUU八U00;:al‘OqJl‘.上,L01U八UUo000n八0..…68204、1102200000000000000000018191615222023叭24丁丙丁丙丁丁丁丁丁HS40HS54一3HS50HS52HS49HS56HS53HS57I江546比1718192021丝邓2’陕西地质第9卷4地球化学异常结构分带图综合异常图是区域化探资料整理及其研究中的重要图件之一,目前一般多将几个单元素异常放在同一张图上构成综合异常图,这类图对三、四个元素还较清晰,但元素多时,图面负担太重,许多元素的不同线条、花纹交叉、重叠在一起,使人难以清晰辨认。我们在工作实践中探索出一种剖析综合异常的方法一地球化学异常结构分带图。该图属多元素综合异常图,能清晰、客观的反映元素的组合、分布及其水平分带特征,准确地划分组合异常的界域。4.1编图的方法在研究区内,将所有元素数据进行点群分析和因子分析,按区内成矿规律、矿产类型对元素进行分组,同一组元素编制一张地球化学结构分带图。其编制方法步骤简述如下:(1)用迭代剔除法求各元素的异常下限(T二X+25)。(2)求取以异常下限为分母的各元素每一点的衬度值,这样每一个分析点上的衬度值比例尺1:200000图3SnBIWPbBeBaCd异常结构分带图F19.3.Thezonationaldia又ramshowingtheanomaliesofSn,Bi,Mo,W,Pb,Be,BaandCd.第l期区域化探异常研究方法探讨63大于或者等于1者均为异常点。(3)计算异常结构值为了便于勾绘异常结构分带图,将异常范围内每个点上出现的元素组合,同一个由“1”和“O”组成的多位数来表示,这个多位数称为异常结构值。现举例说明它的求法:假设将Sn、Bi、M。、W、Pb、Be这六个元素编一张地球化学异常结构分带图,按下列顺序给定值。元素SnBiMoWPbBe定值210100100010000100000Sn元素占个位,将某点Sn元素的异常衬度值》1者给1,1者给。;Bi元素占十位,’衬度》1者,在十位给1,1者给。;MO元素占百位,“一依次类推。如某点出现Sn、M。、Pb元素的异常组合,则该点的异常结构值为10101,这样就求得区内每一点的异常结构值。以上均由计算机处理。(4)勾绘异常结构分带图1:200000图4SnBIMoWPbBeZnBaCd合综异常困F19。4.ComprehensiveanomalymapofSn,Bi,Mo,W,Pb,Be,Zn,BaandCd.“陕_西地质第9卷由结构值的求法可以看出,该值仅表示某点异常元素的组合,不具强度概念,故该图不能由计算机自动成图,只能按结构值代表的元素组合来手工勾绘成图。4.2应用效果利用xx地区化探扫面的七个元素的分析数据,运用上述方法编制了地球化学异常结构分带图(图3),将该图与相同异常的一般综合异常图(图4)对比,并结合地质、矿产及异常检查结果进行了分析研究,认为结构分带图有以下几个特点:(1)这种图件在微机上进行计算,最终输出结构值,手工勾绘异常结构分带图,方法简便,容易在一般微机上实现。(2)能体现区域异常组分分带特征,指明结构中心,这对认识异常性质、确定矿化类型都有意义,并对区域异常的综合解释推断有着特殊的作用。例如某地区Ba等九个元素的地球化学异常结构分带图(图3)中,可以清楚看到,1号异常表现出似同心圆状的结构分带特征(经地表检查见富锡矿),而2号异常表现出不同元素的结构分带特征,这两个异常虽然相邻,却反映了两种不尽相同的地球化学作用。(3)能较为准确、客观的圈出综合异常的范围,在多元素综合异常图上确定综合异常的范围是一项较困难的工作,其结果往往带有工作者的主观因素在内,该图较为客观的解决了这个问题。5异常解释及找矿预测6.1异常解释在xx地区两个1:20万的综合异常图上,出现A、B、C、D四个异常带,均有15种以上元素异常集中呈带状分布,根据异常带的元素组合及分布特征,很难确定哪些元素的异常具有找矿意义,