用户名: 密码: 验证码:
地球化学矿致异常非线性分析方法研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
本论文选题来源于中国地质调查局地质调查工作项目“火山岩盖层下找矿技术攻关研究”子课题:地球化学矿致异常定量评价方法研究。由于火山岩盖层对矿体的封闭作用,对成矿元素的向上扩散和微渗透作用具有很强的阻滞能力,使传统地球化学异常的评价解释变得非常困难。要突破传统,就需要应用新理论、新方法,除了要加强深穿透地球化学等新的化探技术方法研究以外,从数学地质和地质统计学角度,应用现代非线性科学的最新成果,对化探异常的再评价,再解释,进一步挖掘化探数据中蕴含的成矿信息无疑也是一条重要的途径。
     本论文研究区是西藏洞嘎普铜矿勘查区,论文从地球化学矿致异常定量评价方法研究出发,利用多种现代数学方法和非线性分析方法,挖掘化探数据蕴含的成矿弱异常信息。首先利用人工免疫算法构建了地球化学数据降噪模型,对地球化学元素数据进行降噪处理,以此消除特高值带来的屏蔽及剔除假异常;其次,根据盲源分离算法建立反演化探数据元素组合模型,以此确定地球化学成矿元素组合;然后,利用分形含量面积法确定地球化学单元素及元素组合的异常下限,再利用分形含量梯度法圈定异常浓集中心,进而确定异常分带性;最后,将元素分带特性研究与研究区地质特征相结合,对比单元素异常图及组合异常图,对研究区的地球化学元素作出异常分类和异常评价解释。本论文的研究将为火山岩盖层下的化探找矿技术方面提供新的技术方法。
     本文的创新点在于:
     (1)论文在地球化学异常识别方面采用了多种非线性分析方法,克服了传统方法中将地球化学数据看成是符合正态分布或对数正态分布而带来的不合理性,以及传统上将地球化学元素分布曲面看成是连续光滑曲面的不科学性。
     (2)采用人工免疫理论对地球化学原始数据进行降噪处理,剔除了特高值带来的屏蔽和假异常,使数据更具有效性。
     (3)由于地球化学元素分布曲面是分形曲面,因而本论文通过分析多种分形方法,使用了将陈秋明的含量-面积模式及谢和平的投影覆盖法结合并改进的分形含量面积法来确定异常下限,并用分形含量梯度法确定异常浓集中心,以此圈定异常分带。
     (4)论文首次尝试将盲源分离技术用于地球化学领域中,即将盲源分离理论中的FastICA算法应用于地球化学数据寻求区域地球化学矿化元素组合之中,形成一套表生地球化学自动求取元素组合的方法模型,从而反演出地球化学各元素组合特征。
     (5)论文开展了非线性科学与高新信息处理技术相结合的地球化学数据处理新方法的研究,为多尺度地球化学异常圈定及深部成矿预测理论提供了新的思路,研究选题不仅具有新颖性和实用性,而且所提出的具有扎实工作基础的研究方案将为现有的化探工作提供有益补充。
     通过研究,在MAPGIS下建立了地球化学异常圈定系统,对西藏自治区洞噶普铜矿的1:1万土壤地球化学测量资料进行处理后,绘制了测区地球化学元素的异常图和组合异常图,并进行了异常解释评价。经过验证于测区内共划分出4个异常,其中甲(1)类异常2个,编号为Ht-甲(1)-1,Ht-甲(1)-2;甲(2)类异常1个,编号为Ht-甲(2)-1;乙类异常1个,编号为Ht-乙-1。所划的异常整体分布反映的元素套叠状况、异常强度均反映了测区的实际状况。异常与目前发现的矿体、矿化体吻合。异常完整、闭合好,表明异常属于矿致异常。其中Ht-甲(1)-1异常与Ⅲ号强矿化带吻合,Ht-甲(1)-2异常与Ⅳ号矿体吻合,Ht-甲(2)-1异常位于③号韧性剪切带上盘,与2004年发现的Ⅴ号矿体关系密切。
This article is from the Geological Survey projects of China Geological Survey, sub-topics of‘the exploration technology research Under the volcanic cover’: quantitative evaluation method of Geochemical ore-forming anomaly. Since the closure of volcanic cover is that the body effect on the forming elements of the upward diffusion and micro-infiltration block has a strong ability to evaluate the traditional interpretation of geochemical anomalies become very difficult. To break through the traditional, you need to apply new theories, new methods. n addition to strengthening the deep-penetrating geochemical exploration of new technology and methods of research other than geology and geostatistics from the mathematical point of view, the application of the latest achievements of modern nonlinear science, the re-evaluation of geochemical anomalies, to explain, and further excavation geochemical metallogenic information contained in the data is undoubtedly an important way.
     This study area is the Dong Gapu copper exploration area in Tibet. From quantitative evaluation method of Geochemical ore-forming anomaly, making use of a variety of modern mathematical methods and nonlinear analysis, mining ore weak anomalies contains in geochemical data. First use of artificial immune algorithm to construct a noise model of geochemical data, geochemical elements data on noise reduction, in order to eliminate high-value brought by special shielding and removing false anomaly. Secondly, based on blind source separation algorithm for establishing anti-evolution geochemical data elements combined model, in order to determine combination of geochemical ore-forming elements. Then, using the fractal concentration-area method to determine single-element geochemistry and elemental composition of the anomaly threshold, re-use content gradient method of fractal delineate abnormal concentration center and to determine the zoning exception. Finally, the elemental characteristics of research and study area with geological characteristics of the combination, compared to single-element anomaly maps and the composite anomaly map of the study area of the geochemical elements to explain the anomaly classification and anomaly evaluation. This study will provide new technical methods for geochemical exploration technical aspects under the volcanic cover.
     Innovation of this article is:
     (1)The article in terms of geochemical anomaly identified by a variety of nonlinear analysis method to overcome the traditional methods of geochemical data will be seen as consistent with normal or lognormal distribution brought about irrationality, and the distribution of surface geochemical elements traditionally seen as a continuous smooth surface of the unscientific nature.
     (2)Artificial immune theory of noise on the geochemistry of the raw data processing, removing the special shielding to bring high-value and false exceptions, to make the data more effective.
     (3) The surface geochemical element distribution is fractal surface, and thus this article, through analysis of various fractal methods, the use of the content-area model of Chen Qiuming and the Projective covering method of Xie Heping and improved method combined with the fractal method to determine anomaly threshold, and the gradient method used to determine the fractal content of abnormal concentration center, as delineated anomalous zone.
     (4) The article first attempt to blind source separation technique for geochemical field, blind source separation theory about the FastICA algorithm is applied to regional geochemical data for mineral elements in combination, forming a supergene geochemical method of automatic model to strike a combination of elements, which combined the elements of anti-geochemical characteristics.
     (5) The article carried a nonlinear information processing science and high technology combined geochemical data processing method of research for multi-scale delineation of geochemical anomalies and deep metallogenic prediction theory provides a new way of thinking, research topics not only novel and practical, but also has a solid work proposed research program based on the existing geochemical exploration will provide a useful supplement.
     Through research, established geochemical anomaly delineated system under the MAPGIS . after the data processing of 1:1 million soil geochemical survey of Dong Gapu copper in Tibet Autonomous Region, drawing elements of the surveyed area geochemical anomaly maps and the combination of anomalies , and had anomaly interpretation and evaluation. Validated in the survey area was divided into four abnormal. There have two A (1) class anomaly , numbers are Ht-A (1) -1 and Ht-A (1) -2; one is A (2) class anomaly , number is Ht-A (2) -1; one is anomaly B , number is Ht-B -1. The overall distribution of anomalies by overlaying the elements reflected in the situation, unusual intensity reflects the actual situation of the surveyed area. Anomalies and ore bodies found so far, consistent with mineralization. anomaly complete, close is good, that anomaly is an anomaly of mineralization. Where Ht-A (1) -1 andⅢof ore-forming anomaly with strong fit, Ht-A (1) -2 anomaly consistent withⅣorebody, Ht-A (2) No -1 anomaly in the ductile shear belt drive③, found in 2004 closeⅤorebody.
引文
[1]赵鹏大,陈永清,刘吉平等.地质异常成矿预测理论与实践[M].北京:中国地质大学出版社,1999
    [2]孙霞,吴自勤,黄畇.分形原理及其应用[M].合肥:中国科技大学出版社,2003
    [3]王兴元.复杂非线性系统中的混沌[M] .北京:电子工业出版社,2003
    [4]申维.分形混沌与矿产预测[M] .北京:地质出版社,2002
    [5]赵鹏大,池顺都.初论地质异常[J] .地球科学-中国地质大学学报,1991,16(3):241-248
    [6]尹冰川,谢学锦.地球化学省的概念、特征及其与成矿省的关系[J] .长春地质学院学报,1994,24(1):37-43
    [7]谢学锦,刘大文,向运川,严光生.地球化学块体——概念和方法学的发展[J.中国地质,2002,29(3):225-233
    [8]谢学锦.勘查地球化学的现状与未来展望[J] .地质论评,1996,42(4):346-356
    [9]谢学锦.区域化探[M] .北京:地质出版社,1979
    [10]李长江,麻士华.矿产勘探中的分形、混沌与ANN[M] .北京:地质出版社,1999
    [11]Allegre C J.Lewin E.Scaling laws and geochemistry distributions[J].Earth and Planetary Letters,1995,132:1-13
    [12]赵鹏大,陈永清.地质异常矿体定位的基本途径[J] .地球科学-中国地质大学学报,1998,23(2):111-114
    [13]王学求,谢学锦.金的勘查地球化学—理论与方法-战略与战术[M] .济南:山东科技大学出版社,2000
    [14]孟宪国.试论非线性科学在数学地质中的地位和作用[J].地球科学进展,1993,8(1):66-71
    [15]程志中,谢学锦.岩石中元素背景值变化对地球化学成矿预测的影响[J].中国地质,2006,33(2):411-417.
    [16]郝立波,李巍,陆继龙.确定岩性复杂区地球化学背景与异常的方法[J] .地质通报, 2007, 26(12):1531-1535.
    [17]郭科.复杂地质地貌区多尺度地球化学异常识别的非线性研究[D].成都理工大学,2005.
    [18]陈聆,郭科,唐菊兴等.人工免疫系统在地球化学异常中的应用[J].地球物理学进展,2009, 24(3):1136-1141.
    [19]Jerne N K .Towards a Network Theory of Immune System.Annual Immunology,125C:1974:373-389
    [20]Farmer J D.,Packard,N H.&Perelson,A.S. The Immune System,Adaptation, And Machine Learning , physica 22D,1986.187~204
    [21]Perelson A S. Immune Network Theory. Immunological review .1986(10):5-36
    [22]Bersini H, Varela F J. Hints for Adaptive Problem Solving Gleaned from Immune System .In Proc The First Workshop Parallel Problem Solving From Nature, 1990.343-354
    [23]Varela F J, Stewart J. Dynamics of a Class of Immune Network.. Global Stability of ldiotype Interactions. J . Theoretical Biology,1990(144):93-101
    [24]肖人彬,王磊.人工免疫系统:原理、模型、分析及展望[J].计算机学报,2002
    [25]吴作顺等.负选择模型中初始检测器的一个生成算法[J].电子学报,2003
    [26]邵峰晶,于忠清.数据挖掘—原理与算法[M] .中国水利水电出版社,2003
    [27]Mehmed Kantardzic著,闪四清,陈茵,程雁等译.数据挖掘—概念、模型、方法和算法.清华大学出版社,2003
    [28]刘越江,黄今慧.数据挖掘中的数据预处理技术[J].科技情报开发与经济,2003
    [29]黄守坤.异常数据挖掘及在经济欺诈发现中的应用[J].统计与决策,2003
    [30]蒋良孝,蔡之华.时序数据库中的数据挖掘研究[J] .微机发展,2003
    [31]易顺明.黄土地区重点城镇滑坡灾害预测预报研究.四川大学博士后学位论文,2000
    [32]马建仓,牛奕龙,陈海洋.盲信号处理[M] .北京:国防工业出版社,2006.
    [33]边肇祺,张学工.模式识别[M] .北京:清华大学出版社,1999.
    [34]张贤达,保铮.盲信号分离[J] .电子学报,2001,29(12A):1766-1771.
    [35]Jutten C, Herault J. Blind separation of sources, PartⅠ:An adaptive algorithm based on neuromimatic architecture, Signal processing, 1991,24(1):1-10
    [36]TONG L,LIU R,Soon V C,et al.Indeterminacy and identifiabililty of blind identification[J].IEEE Trans on Circuits and Systems,1991,38(5):499—506.
    [37]Tong L, Li R W, Soon V C, Indeterminacy and identifiability of blind identification. IEEE Trans. On Circuits and Systems, 1991,38(5):499-509
    [38]Burel G. Blind separation of sources: A nonlinear neural algorithm. Neural Network, 1992,5:937-947
    [39]CARDOSO J F. Blind beamforming for non-gaussian signals [J]. IEE Proceedings-F, 1993,140(6):224-230.
    [40]Comon P. Independent component analysis, a new concept. Signal Processing, 1994,36:287-314
    [41]Bell A J, Sejnowski T J. An information-maximization approach to blind separation and blind deconvolution. Neural Computation, 1995,7(6):1004-1034
    [42]HYVARINEN A,OJA E.A fast fixed point algorithm for independent component analysis[J].Neural Computation.1997,9(7):1483-1492.
    [43]AMARI S.Natural gradient works efficiently in learning[J].Neural Computation,1998, 10(2):241-276.
    [44]Hyvarinen S, Karhunen J, Oja E. Independent Component Analysis, New York: John Wiley, 2001
    [45]Tan Y, Wang J, Zurada J. Nonlinear blind source separation using a radial basis function network. IEEE Trans. On Neural Networks, 2001,12(1):124-134
    [46]高颖.盲源分离理论及其在地球物理勘探中的应用[D] .吉林大学博士学位论文.2008.06
    [47]Carlson C A.Spatial distribution of ore deposit[J].Geology,1991,19(2):107-110
    [48]Blenkinsop T.The fractal distribution of gold deposit:two examples from the Zimbabwe Archaen Craton[A].Barton C C.Fractal in Petroleum Geology and Earth Processes[M].New York and London:Plonum Press,1994,247-258
    [49]Agterberg F P,Qiuming C,Wright D F.Fractal modelling of mineral deposits[A].International Sympsium on the Application of Computers and Operations Research in the Mineral Industries[C].Canada,Montreal:Canadian Institute of Mining Metallurgy and Petroleum, 1993,43-55
    [50]Sanderson D J,Roberts S A.Fractal relationship between vein thickness and gold grade in drill core from La Codosera,Spain[J].Economic Geology,1994,89:168-173
    [51]Allege C J.Scaling laws and geochemistry distribution[J].EPSL,1995,132(1-4):1-13
    [52]Mandelbrot B. B. The fractal geometry of nature[M]. San Francisco:Freeman, 1982
    [53]Qiuming C,Agterberg F P,Ballantyne S B.The separation of geochemical anomalies from background by fractal methods[J].Journal of Geochemical Exploration, 1994,51:109-130
    [54]Qiuming C,Agterberg F P,Bonham-Carter G F.A spatial analysis method for geochemical anomaly separation[J].Journal of Geochemical Exploration,1996,56:183-195
    [55]秦长兴.矿床学中若干问题的分形几何学研究[D].武汉:中国地质大学档案馆,1989
    [56]秦长兴,翟裕生.矿床学中若干自相似性现象及其意义[J].矿床地质, 1992,11(3):259-265
    [57]程小久,卢建杭,宋亮明.铅锌品位分维D值的意义和计算程序[J].地质与勘探,1994,30(5):32-35
    [58]沈步明,沈远超.新疆某金矿的分维特征及其地质意义[J].中国科学(B辑),1993,23(3):297-302
    [59]申维.初论非线性地质学与地质异常[J]. 1997
    [60]申维.成矿预测中分形模型分维数估计的新方法[J].长春地质学院学报,1997,1:58-70
    [61]孟宪国,赵鹏大.论地质现象中的分形统计学[J].地球科学,1996,22(1):601~603
    [62]赵鹏大,周有武.金矿化变异的数学地质研究[M].武汉:中国地质大学出版社,1992,1-4
    [63]赵鹏大,池顺都,陈永清.查明地质异常:成矿预测的基础[J].高校地质学报,1996, (2): 360-373
    [64]连长云.大兴安岭西坡地区与超大型矿床有关的非线性研究[D].长春:长春科技大学,1997,4-8
    [65]毛华海.微细浸染型金矿床地球化学,热力学,动力学及分形研究[D].贵阳:中国科学院地球化学研究所教育处,1997
    [66]王祖伟.蚀变构造型金矿床地质地球化学及资源潜力[D].贵阳:中国科学院广州地球化学研究所教育处,1998
    [67]李长江,蒋叙良,徐有浪等.浙江中生代热液矿床的分形研究[J].地质科学,1996,31(3):264-273
    [68]李长江.快速逼近潜在超大型矿床的新理论和新方法[J].矿床地质,1998,17(增刊):853-860
    [69]周永章,Edward H G, Guha J.地质热场中微量元素迁移的方向性和分维结构图象[J].中国科学(B辑),1994,24(12):1308-1313
    [70]武生智等.沙粒粗糙度和粒径分布的分形特性[J].兰州大学学报, 1999,35(1),
    [71]刘英俊,邱德同.勘查地球化学[M].北京:科学出版社,1987.
    [72]吕国安.成矿区带地球化学异常评价方法[M].北京:冶金工业出版社,2002.
    [73]郭科,陈聆,唐菊兴,吴艳蕾.分形含量梯度法确定地球化学浓集中心的新探索[J].地学前缘, 2007,14(5):285-289
    [74]张发启.盲信号处理及应用[M].西安:西安电子科技大学出版社2006.
    [75](英)肯德尔(M. Kendall).多元统计[M] .北京科学出版社,1983.
    [76]胡增辉,炬波,梁甸农.基于盲源分离的近场源参数估计[J]. 2010,26(6).
    [77]唐江波,曾庆宁,郭威.基于ICA的快速定点算法[J].桂林电子科技大学. 2007,(08)
    [78]郭武,朱长仁,王润生.一种改进的FastICA算法及其应用[J].计算机应用. 2008,(04)
    [79]郭科,刘诚.基于独立分量分析的地震信号多次波盲分离方法研究[J].成都理工大学学报2008.05.
    [80]张济忠.分形[M].清华大学出版社,1995
    [81]申维,赵鹏大.分形统计模型的理论研究及其在地质学中的应用[M].地质科学,1998, 33(2):234-243
    [82] Crovelli R A, Schmoker J W, Balay R H. U.S. Department of the interior U.S. geological survey: Fractal lognomal percentage analysis of U.S. geological survery’s 1995 national assessment of conventional oil and gas resources. Nonrewable Resources, 1997, 6(1):43~51
    [83]申维.多重分形模型的理论研究及其在地质学中的应用[M].中国数学地质(8),北京,地质出版社,1997.117~122
    [84] John N L and Kotz S. Continuous Univariate Distributions-I. John Wiley & Sons. Inc.New York,1970. 300p
    [85]Mandelbrot B B. The fractal geometry of nature[M].San Francisco:Freeman, 1982.
    [86]赵鹏大,周有武.金矿化变异的数学地质研究[J].武汉:中国地质大学出版社,1992. 1~4
    [87]申维,基于奇异值分解法的含量-面积法对化探异常的确定[J].《地质通报》,08年5月,27卷5期
    [88]上海交通大学,集美大学编.高等数学——及其教学软件(下册)第二版[M].北京:科学技术出版社,2003.
    [89]Rockafellar P.T.,Converx Analysis[M],princeton University press,1972.
    [90]Aubin J P ,Gradients Generalises Clarke[J],Ann.Sci.Math.Quebec,1978b,2:197-252.
    [91]Aubin J P.,Contingent Derivatives of Set-valued Maps[J],Math.Research Center,Univ.of Wisconsin,Madison,Tech.Red,1980b,2044.
    [92]Ioffe A D. , Nonsmooth Analysis, Differential Calculus of Nondifferent Mappings[M],Trans.Am.Math.Soc.,1981,266:l-56.
    [93]Loffe A D.,Approximate Subdifferentials of Nonconvex Functions[J]. Cahiers Math,de La Decision,Univ,de Paris-Dauphine,no.8120.
    [94]Clarke F H.,Optimization Nonsmooths Analysis [M].New York:JohnWiely,1983,227-231
    [95]王辉,刘天放.利用重采样时差梯度法定量估算小断层落差[J].煤田地质与勘探,2000,(28):59-61.
    [96]姚长利,黄卫宁,管志宁.综合利用位场及其垂直梯度的快速样条曲化平方法[J].石油地球物理勘探,1997,32(2):229-236.
    [97]岳中琦.多层与梯度非均匀材料弹性力学问题解析解的简明数学理论[J].岩石力学与工程学报,2004, 23(17):2545-2554.
    [98]王家华,高余海,周叶.克里格地址绘图技术——计算机的模型和算法[M] .北京:石油工业出版社,1999.
    [99]吕国安.成矿区带地球化学异常评价方法[M] .北京:冶金工业出版社,2002
    [100]程其襄.实变函数与泛函分析基础[M.北京:高等教育出版社,2001.
    [101]董连科,王晓.函数的H-导数和分形动力学[J] .高压物理学报,1977,12:13-20.
    [102]Besicovitch A S.Proceedings and the Combridge Pholosophical Sciety,1945,14:101

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700