广西六林多金属矿化探数据处理分析及成矿预测
|
摘要
化探数据处理是地球化学勘查的一项重要内容，不同的处理方法直接影响化探找矿的效果。化探数据处理是应用数学方法从地球化学原始数据中提取信息，指示变量（主要是化学元素）与各种地质现象的内在联系。如何科学有效地提取化探异常信息，并从大量异常中进行快速准确地筛选评价，以确定进一步找矿的靶区，则是决定化探工作找矿效果和经济效益的关键。
本文以广西六林多金属矿区为研究区，采用传统迭代法、分形法、多元统计分析法、成矿能量法对1:5万水系沉积物地球化学测量数据进行处理分析，并结合实际地质资料，圈定成矿远景区。
地球化学元素含量数据分布特征表现为不服从正态或对数正态分布，而是呈多峰态、高丛集的多重分形分布特征。采用分形法对研究区地球化学特征进行分析，可以很好地揭示各主要成矿元素空间变化及矿化富集规律。采用分形理论进行分析，在拐点的确定上，即分界点的确定，求和法要比面积法直接一些。而且，求和法分组直观快捷，处理过程相对简单，可优先考虑求和法进行元素数据的处理分析。
通过相关分析、聚类分析与因子分析确定元素共生组合关系，将12种元素分为Pb-Sn-Ag-Sb-Au、W-Mo-Bi、Cu-Ag-Zn、As、F五组。根据元素组合关系得到元素累加组合异常图。Pb-Sn-Ag-Sb-Au元素组合即F1因子组合，分布在研究区的中部，是寻找贵金属和铅多金属矿的有利靶区；Cu-Ag-Zn元素组合即F3因子组合异常区分布在研究区的西南角和中东部，是寻找多金属矿产的有利靶区；F1因子组合与F3因子组合在空间上具有一定的套合关系，具有更为重要的找矿意义；W-Mo-Bi元素组合即F2因子组合，为高温成矿元素组合，与燕山期花岗岩体密切相关。利用多元正态分布理论计算得出的元素组合异常分布相对集中，浓集中心明显，形态较规则，并且套合关系比较明显，其找矿意义比单元素异常更为明确，同时组合异常也更加符合成矿作用和热液活动中元素共生组合的地质特征。
采用成矿能量法对研究区进行化探综合异常的提取与筛选评价，结合本区地质条件，确定了研究区的找矿靶区。E1区可作为寻找热液型铜多金属矿的重点找矿远景区，E2区可作为寻找高温矿物组合—钨钼矿的重点找矿远景区，E3区为寻找多金属矿床的重点找矿靶区；E4区可作为寻找铅锌矿的次级找矿远景区。
Geochemical exploration data processing is an essential part, which may affect thegeochemical exploration effects directly. It is to extract information from the original databy mathematics method, to indicate the interrelationship between chemical elements andgeological phenomena. How to extract geochemical anomaly scientifically and efficiently,and then to select and evaluate them quickly and precisely, to identify the prospect areaand improve the probability of ore, is the key for the exploration effects and economicbenefit.
This paper analyzes the1:50000stream sediment data in the Liulin area, Guangxiprovince and identify the prospect area, by use of traditional iterative method, fractaltheory, multivariate statistical analysis method and metallogenic energy method.
The distribution characteristics of geochemical elements data is not followed by thenormal distribution and lognormal distribution, but show the characteristics of highcluster fractal distribution. Adopting fractal theory to analyze the geochemicalcharacteristics can reveal the spatial distribution and enrichment regularity of the mainore-forming elements fully. In the fractal theory, the content-sum method is more directlythan content-area method in the determination of the cutoff point. And content-summethod groups the data more easily, so it is the better choice to process the data.
By using correlation analysis, Cluster analysis, factor analysis to determine themineral assemblages,12elements are divided into five groups, namely Pb-Sn-Ag-Sb-Au,W-Mo-Bi, Cu-Ag-Zn, As and F. According to the groups, make the elements assembledanomaly maps. Anomalies of Pb-Sn-Ag-Sb-Au, namely factor1, distributing in themiddle part of the area, is the prospect area of precious metals and Pb-bearing polymetaldeposits; Anomalies of F3[Cu-Ag-Zn] distributes at the southwest corner and the middleeast part, which is the target area of searching Cu-bearing polymetal ore; F1coincideswith F3in the middle part, and it is more meaningful for prospecting; Anomalies ofF2[W-Mo-Bi] represents the association of high temperature elements, related with thegranite mass of Yanshanian period closely. The elements assembled anomalies calculatedby using multivariate normal distribution theory, distributes more collective relatively, with clear concentration heart and regular shape. These anomalies which may match withsome others sterically are more meaningful than single element anomaly, and meanwhileassembled anomaly is better in accordance with the geological characteristics of elementassociation in mineralization and hydatogenesis.
In this paper, adopting metallogenic energy method extracts composite anomaly andevaluate. Combined with geological conditions, determine the prospect area. E1area canbe the main prospecting target of copper polymetallic ore. E2area can be priorityprospect area to search wolfram-molybdenum ore. E3area can be considered the majorarea of polymetallic deposits. E4can be treated as the secondary exploration target oflead-zinc deposit.
本文以广西六林多金属矿区为研究区，采用传统迭代法、分形法、多元统计分析法、成矿能量法对1:5万水系沉积物地球化学测量数据进行处理分析，并结合实际地质资料，圈定成矿远景区。
地球化学元素含量数据分布特征表现为不服从正态或对数正态分布，而是呈多峰态、高丛集的多重分形分布特征。采用分形法对研究区地球化学特征进行分析，可以很好地揭示各主要成矿元素空间变化及矿化富集规律。采用分形理论进行分析，在拐点的确定上，即分界点的确定，求和法要比面积法直接一些。而且，求和法分组直观快捷，处理过程相对简单，可优先考虑求和法进行元素数据的处理分析。
通过相关分析、聚类分析与因子分析确定元素共生组合关系，将12种元素分为Pb-Sn-Ag-Sb-Au、W-Mo-Bi、Cu-Ag-Zn、As、F五组。根据元素组合关系得到元素累加组合异常图。Pb-Sn-Ag-Sb-Au元素组合即F1因子组合，分布在研究区的中部，是寻找贵金属和铅多金属矿的有利靶区；Cu-Ag-Zn元素组合即F3因子组合异常区分布在研究区的西南角和中东部，是寻找多金属矿产的有利靶区；F1因子组合与F3因子组合在空间上具有一定的套合关系，具有更为重要的找矿意义；W-Mo-Bi元素组合即F2因子组合，为高温成矿元素组合，与燕山期花岗岩体密切相关。利用多元正态分布理论计算得出的元素组合异常分布相对集中，浓集中心明显，形态较规则，并且套合关系比较明显，其找矿意义比单元素异常更为明确，同时组合异常也更加符合成矿作用和热液活动中元素共生组合的地质特征。
采用成矿能量法对研究区进行化探综合异常的提取与筛选评价，结合本区地质条件，确定了研究区的找矿靶区。E1区可作为寻找热液型铜多金属矿的重点找矿远景区，E2区可作为寻找高温矿物组合—钨钼矿的重点找矿远景区，E3区为寻找多金属矿床的重点找矿靶区；E4区可作为寻找铅锌矿的次级找矿远景区。
Geochemical exploration data processing is an essential part, which may affect thegeochemical exploration effects directly. It is to extract information from the original databy mathematics method, to indicate the interrelationship between chemical elements andgeological phenomena. How to extract geochemical anomaly scientifically and efficiently,and then to select and evaluate them quickly and precisely, to identify the prospect areaand improve the probability of ore, is the key for the exploration effects and economicbenefit.
This paper analyzes the1:50000stream sediment data in the Liulin area, Guangxiprovince and identify the prospect area, by use of traditional iterative method, fractaltheory, multivariate statistical analysis method and metallogenic energy method.
The distribution characteristics of geochemical elements data is not followed by thenormal distribution and lognormal distribution, but show the characteristics of highcluster fractal distribution. Adopting fractal theory to analyze the geochemicalcharacteristics can reveal the spatial distribution and enrichment regularity of the mainore-forming elements fully. In the fractal theory, the content-sum method is more directlythan content-area method in the determination of the cutoff point. And content-summethod groups the data more easily, so it is the better choice to process the data.
By using correlation analysis, Cluster analysis, factor analysis to determine themineral assemblages,12elements are divided into five groups, namely Pb-Sn-Ag-Sb-Au,W-Mo-Bi, Cu-Ag-Zn, As and F. According to the groups, make the elements assembledanomaly maps. Anomalies of Pb-Sn-Ag-Sb-Au, namely factor1, distributing in themiddle part of the area, is the prospect area of precious metals and Pb-bearing polymetaldeposits; Anomalies of F3[Cu-Ag-Zn] distributes at the southwest corner and the middleeast part, which is the target area of searching Cu-bearing polymetal ore; F1coincideswith F3in the middle part, and it is more meaningful for prospecting; Anomalies ofF2[W-Mo-Bi] represents the association of high temperature elements, related with thegranite mass of Yanshanian period closely. The elements assembled anomalies calculatedby using multivariate normal distribution theory, distributes more collective relatively, with clear concentration heart and regular shape. These anomalies which may match withsome others sterically are more meaningful than single element anomaly, and meanwhileassembled anomaly is better in accordance with the geological characteristics of elementassociation in mineralization and hydatogenesis.
In this paper, adopting metallogenic energy method extracts composite anomaly andevaluate. Combined with geological conditions, determine the prospect area. E1area canbe the main prospecting target of copper polymetallic ore. E2area can be priorityprospect area to search wolfram-molybdenum ore. E3area can be considered the majorarea of polymetallic deposits. E4can be treated as the secondary exploration target oflead-zinc deposit.
引文
Agterberg F P,Qiuming C,Wright D F. Fractal modelling of mineral deposits[A]. InternationalSympsium on the Application of Computers and Operations Research in the MineralIndustries[C]. Allege C J. Scaling laws and geochemistry distribution [J].EPSL,1995,132(1~4):1~13.
Canada,Montreal:Canadian Institute of Mining Metallurgy and Petroleum,1993.43~53.
Carlson C A. Spatial distribution of ore deposit [J]. Geology,1991,19(2):107~110.
Kulik D A,Chernovsky M I. Fractal properties of multi-order folding as a tool for exploration ofLow-grade banded iron ores in the kirivoy Rog Basin (Ukraine)[J]. Geologiskhe Rundschau,1996,85(10):3~11.
QiumingC, AgterbergF P, Bonham-CarterG F. A spatial-nalysis method for geochemical anomalyseparation [J]. Journal ofGeo-chemicalExploration,1996,56:183-195.
Sanderson D J, Roberts S A. Fractal relationship between vein thickness and gold grade in drill corefrom La Co-dosera,Spain[J].Economic Geology,1994,89:168~173.
Turcotte D L.Fractals and Chaos in Geophysics [M].2nd ed.CombridgeUK: CambridgeUniversityPress,1996:81-99.
戴慧敏,鲍庆中,宫传东等.因子分析法对内蒙古查巴奇地区水系沉积物地球化学分区的应用研究.现代地质,2010,24(2):245-251.
龚庆杰.一种确定地球化学异常下限的简便方法[J].地质地球化学,2001,29(3):215-219.
郭春影,高帮飞,邢学文.两种容量维方法提取化探数据异常下限效果对比.黄金地质,2008,3(29):13-17.
何明华.MAPGIS数字高程模拟系统在化探数据处理中的应用.矿山测量,2005,2:13-15.
侯春秋.西藏拉萨-泽当地区1:20万铅锌地球化学异常信息处理研究:[硕士学位论文].成都:成都理工大学,2011.
黄瑞.化探数据处理方法研究—以攀西地区铂族元素为例：[硕士学位论文].成都:成都理工大学,2005.
井国正,闫杰等.分形理论在江日呷玛地区地球化学异常圈定中的应用.矿产勘查,2010,4(1):375-379.
李宾,李随民等.分形方法圈定河北省龙关地区化探元素异常.地质调查与研究,2011,35(2)154-161.
李长江,蒋叙良,徐有浪等.浙江中生代热液矿床的分形研究[J].地质科学,1996,31(3):264-273.
李长江.快速逼近潜在超大型矿床的新理论和新方法[J].矿床地质,1998,17(增刊):853-86.
李建东,龚庆杰.分形理论在勘查地球化学中的应用.中国矿业,2006,15(11):91-94.
李世超.基于分形理论的化探数据处理系统GC-GIS的开发与研制:[硕士学位论文].吉林:吉林大学,2007
李随民,姚书振.基于MAPGIS的分形方法确定化探异常.地球学报.2005,26(2):187-190.
李汶奎,杨贵兵,唐文春等.云南省耿马勐简铅锌矿区地球化学特征及找矿前景预测.物探化探计算技术,2010,32(2):200-205.
李佑国.基于”3S”技术的攀西地区铜镍铂族元素矿床找矿靶区筛选:[硕士学位论文].成都:成都理工大学,2007.
连盈盈.盲提取在地球化学异常识别中的应用:[硕士学位论文].成都:成都理工大学,2009.
孟完国,赵鹏大.论地质现象中的分形统计学[J].地球科学,1996,22(1):601~603.
彭涛.基于GIS的化探数据处理方法研究及分析系统开发:[硕士学位论文].湖南:中南大学,2004.
秦长兴,翟裕生.矿床学中若干自相似性现象及其意义[J].矿床地质,1992,11(3):259~265.
申维.分形混沌与矿产预测.北京:地质出版社,2002.
沈步明,沈远超.新疆某金矿的分维特征及其地质意义[J].中国科学(B辑),1993,23(3):297~302.
时艳香,纪宏金,陆继龙等.水系沉积物地球化学分区的因子分析方法与应用[J].地质与勘探,2004,40(5):73-76.
宋运红,贾大成,李健等.河北承德三沟地区化探数据处理与异常特征.吉林大学通报(地球科学版),2008,38增刊:167-169.
王晓.云南省富宁县大湾金矿土壤地球化学异常评价及找矿方向:[硕士学位论文].成都:成都理工大学,2011.
巫晓兵,范良明,毛玉元等.地球化学找矿中的元素成矿能量方法试验研究.新疆地质,1995,13(1):93.
吴锡生,纪宏金,陈明.化探数据处理的发展、现状与趋势.物化探计算技术,1994,16(1)84~88.
吴锡生,孙春岩.金属矿勘查地球化学.2009(未公开发表).
徐锡华.成矿能-克立格双模型在化探数据处理中的特殊功效.地质找矿论丛,2002,17(1).
杨利民,杨自安,罗铁良等.成矿能量在化探综合异常提取中的应用研究.矿产与地质,2008,22(1):74~77
杨自安.西部高寒山区遥感与化探信息综合找矿定位预测研究:[博士学位论文].北京:中国地质大学(北京),2005.
张彪.内蒙古鲁北地区地球化学异常识别:[硕士学位论文].北京:中国地质大学(北京),2008.
张辉善,张杰,史俊波等.青海格尔木红石山地区水系沉积物测量异常特征及优选找矿靶区.物探与化探,2011,35(6):768-772.
张哲儒,毛华海.分形理论与成矿作用.地学前缘,2000,7(1):195-204.
邹林,彭省临,杨自安等.青海阿尔茨托山地区地球化学（异常）场的多重分形研究.中国地质,2004,31(4):436-441.
邹林.青海阿尔茨托山地区化探找矿信息的快速识别与综合评价:[硕士学位论文].湖南:中南大学,2004.
Canada,Montreal:Canadian Institute of Mining Metallurgy and Petroleum,1993.43~53.
Carlson C A. Spatial distribution of ore deposit [J]. Geology,1991,19(2):107~110.
Kulik D A,Chernovsky M I. Fractal properties of multi-order folding as a tool for exploration ofLow-grade banded iron ores in the kirivoy Rog Basin (Ukraine)[J]. Geologiskhe Rundschau,1996,85(10):3~11.
QiumingC, AgterbergF P, Bonham-CarterG F. A spatial-nalysis method for geochemical anomalyseparation [J]. Journal ofGeo-chemicalExploration,1996,56:183-195.
Sanderson D J, Roberts S A. Fractal relationship between vein thickness and gold grade in drill corefrom La Co-dosera,Spain[J].Economic Geology,1994,89:168~173.
Turcotte D L.Fractals and Chaos in Geophysics [M].2nd ed.CombridgeUK: CambridgeUniversityPress,1996:81-99.
戴慧敏,鲍庆中,宫传东等.因子分析法对内蒙古查巴奇地区水系沉积物地球化学分区的应用研究.现代地质,2010,24(2):245-251.
龚庆杰.一种确定地球化学异常下限的简便方法[J].地质地球化学,2001,29(3):215-219.
郭春影,高帮飞,邢学文.两种容量维方法提取化探数据异常下限效果对比.黄金地质,2008,3(29):13-17.
何明华.MAPGIS数字高程模拟系统在化探数据处理中的应用.矿山测量,2005,2:13-15.
侯春秋.西藏拉萨-泽当地区1:20万铅锌地球化学异常信息处理研究:[硕士学位论文].成都:成都理工大学,2011.
黄瑞.化探数据处理方法研究—以攀西地区铂族元素为例：[硕士学位论文].成都:成都理工大学,2005.
井国正,闫杰等.分形理论在江日呷玛地区地球化学异常圈定中的应用.矿产勘查,2010,4(1):375-379.
李宾,李随民等.分形方法圈定河北省龙关地区化探元素异常.地质调查与研究,2011,35(2)154-161.
李长江,蒋叙良,徐有浪等.浙江中生代热液矿床的分形研究[J].地质科学,1996,31(3):264-273.
李长江.快速逼近潜在超大型矿床的新理论和新方法[J].矿床地质,1998,17(增刊):853-86.
李建东,龚庆杰.分形理论在勘查地球化学中的应用.中国矿业,2006,15(11):91-94.
李世超.基于分形理论的化探数据处理系统GC-GIS的开发与研制:[硕士学位论文].吉林:吉林大学,2007
李随民,姚书振.基于MAPGIS的分形方法确定化探异常.地球学报.2005,26(2):187-190.
李汶奎,杨贵兵,唐文春等.云南省耿马勐简铅锌矿区地球化学特征及找矿前景预测.物探化探计算技术,2010,32(2):200-205.
李佑国.基于”3S”技术的攀西地区铜镍铂族元素矿床找矿靶区筛选:[硕士学位论文].成都:成都理工大学,2007.
连盈盈.盲提取在地球化学异常识别中的应用:[硕士学位论文].成都:成都理工大学,2009.
孟完国,赵鹏大.论地质现象中的分形统计学[J].地球科学,1996,22(1):601~603.
彭涛.基于GIS的化探数据处理方法研究及分析系统开发:[硕士学位论文].湖南:中南大学,2004.
秦长兴,翟裕生.矿床学中若干自相似性现象及其意义[J].矿床地质,1992,11(3):259~265.
申维.分形混沌与矿产预测.北京:地质出版社,2002.
沈步明,沈远超.新疆某金矿的分维特征及其地质意义[J].中国科学(B辑),1993,23(3):297~302.
时艳香,纪宏金,陆继龙等.水系沉积物地球化学分区的因子分析方法与应用[J].地质与勘探,2004,40(5):73-76.
宋运红,贾大成,李健等.河北承德三沟地区化探数据处理与异常特征.吉林大学通报(地球科学版),2008,38增刊:167-169.
王晓.云南省富宁县大湾金矿土壤地球化学异常评价及找矿方向:[硕士学位论文].成都:成都理工大学,2011.
巫晓兵,范良明,毛玉元等.地球化学找矿中的元素成矿能量方法试验研究.新疆地质,1995,13(1):93.
吴锡生,纪宏金,陈明.化探数据处理的发展、现状与趋势.物化探计算技术,1994,16(1)84~88.
吴锡生,孙春岩.金属矿勘查地球化学.2009(未公开发表).
徐锡华.成矿能-克立格双模型在化探数据处理中的特殊功效.地质找矿论丛,2002,17(1).
杨利民,杨自安,罗铁良等.成矿能量在化探综合异常提取中的应用研究.矿产与地质,2008,22(1):74~77
杨自安.西部高寒山区遥感与化探信息综合找矿定位预测研究:[博士学位论文].北京:中国地质大学(北京),2005.
张彪.内蒙古鲁北地区地球化学异常识别:[硕士学位论文].北京:中国地质大学(北京),2008.
张辉善,张杰,史俊波等.青海格尔木红石山地区水系沉积物测量异常特征及优选找矿靶区.物探与化探,2011,35(6):768-772.
张哲儒,毛华海.分形理论与成矿作用.地学前缘,2000,7(1):195-204.
邹林,彭省临,杨自安等.青海阿尔茨托山地区地球化学（异常）场的多重分形研究.中国地质,2004,31(4):436-441.
邹林.青海阿尔茨托山地区化探找矿信息的快速识别与综合评价:[硕士学位论文].湖南:中南大学,2004.