变换后数据的因子分区标准化在东昆仑东段地球化学异常圈定

合集下载

拉脊山东段地区Au、Cu地球化学组合异常识别与提取

拉脊山东段地区Au、Cu地球化学组合异常识别与提取姜晓佳;陈鑫;郑有业;高顺宝;欧阳嵩;张永超;郑磊;黄建【摘要】In order to weaken the closure effect of the geochemical data which serve as typical component data,the authors analyzed the spatial distribution of geochemical elements in the Laji Mountains.The authors used the isometric log-ratio(ilr) transformation to "open" the geochemical data of the water sediments,built the combined model using robust principal component analysis (RPCA) to identify the combined geochemical anomalies and employed fractal filtering technique to strengthen weak anomalies and separate anomaly from background.In this way,the authors delineated geochemical anomaly so as to guide mineral exploration better.The result obtained by using RPCA displays two different compositional assemblages: (I) Au-As,probably representing tectonic-altered rock type gold deposit like the east ditch of Jinyuan,and (II) Cu-Ni,likely representing copper-nickel sulfide mineralization such as Lashuixia.The results from S-A filter technique can decompose anomalies further,strengthen weak anomaly and reduce the abnormal area in the strong background and,in combination with the geological survey of the study area and the anomaly map,can better search for tectonic-altered rock type Au-As deposits and Cu-Ni sulfide mineralization.%为削弱拉脊山东段地区化探数据(典型的成分数据)存在的闭合效应,进一步分析该区地球化学元素空间组合分布规律,使用等距对数比变换(ilr)"打开"水系沉积物地球化学数据;采用稳健主成分分析(RPCA)构建组合模型,用于识别组合地球化学异常;通过S-A分形滤波技术强化弱异常并分离异常与背景,圈定地球化学致矿异常,进一步指导矿产勘查.研究结果表明:得到两组与矿化相关的元素组合,第一组为Au-As,与研究区的构造蚀变岩型金(砷)矿(金源东沟等)相关;第二组为Cu-Ni,与区域内铜镍硫化物矿床(拉水峡等)相关.S-A滤波技术可以进一步分解组合异常,强化弱小异常,同时能够缩小强背景下的异常面积,结合研究区地质概况和分解后的异常图可以进一步指导该地区构造蚀变岩型Au(As)矿和铜镍硫化物矿床的找矿勘查工作.【期刊名称】《物探与化探》【年(卷),期】2017(041)003【总页数】9页(P459-467)【关键词】等距主成分分析;稳健主成分分析;多重分形滤波;组合异常提取;拉脊山东段【作者】姜晓佳;陈鑫;郑有业;高顺宝;欧阳嵩;张永超;郑磊;黄建【作者单位】中国地质大学(武汉)资源学院,湖北武汉 430074;中国地质大学(武汉)资源学院,湖北武汉 430074;中国地质大学(武汉)资源学院,湖北武汉 430074;中国地质大学(武汉)地质调查研究院,湖北武汉 430074;中国地质大学(北京)地质过程与矿产资源国家重点实验室,北京100083;中国地质大学(武汉)地质调查研究院,湖北武汉 430074;中国地质大学(武汉)资源学院,湖北武汉 430074;中国地质大学(武汉)地质调查研究院,湖北武汉 430074;中国地质大学(武汉)资源学院,湖北武汉430074;中国地质大学(武汉)地质调查研究院,湖北武汉 430074【正文语种】中文【中图分类】P632由于地球化学元素携带了某种矿化信息, 研究地球化学元素的空间分布特征可直接或间接地揭示矿化溯源问题。

东昆仑布青山地区水系沉积物测量地球化学特征及找矿方向

东昆仑布青山地区水系沉积物测量地球化学特征及找矿方向崔晓亮;刘婷婷;王文恒;景明;白云【摘要】在1∶5万水系沉积物测量的基础上,对青海东昆仑布青山地区的地球化学特征进行了初步分析.对元素分布特征、单元素异常特征、元素相关性分析、异常元素组合特征进行了研究,发现Au、Cu、Co、Ni、Sb为区内具有一定潜力的找矿指标.推断区内的成矿类型为与花岗岩、花岗斑岩有关的斑岩型铜金多金属矿和与基性一超基性岩有关的钴镍多金属矿,划分出马尼特金铜找矿远景区、得力斯坦铜钴镍找矿远景区和哥日卓多金铜找矿远景区.%The preliminary study of the geochemical features in Buqingshan area of East Kunlun in Qinghai Province was based on the 1: 50000 stream sediment survey. The study included element distribution characteristics, element anomalies, element correlation and element anomaly assemblage. The result shows that the study area might be a prospective area for Au, Cu, Co, Ni and Sb. It is inferred that ore deposits in Buqingshan area are probably porphyry Cu-Au polymetallic deposits related to granite and granite porphyry, whereas Co-Ni polymetallic deposits are related to the basic and ultrabasic rocks. The Mnite area is expected to be a prospective area for porphyry Cu-Au polymetallic deposits, the Delisitan area is likely to have Co-Ni polymetallic deposits associated with basic and ultrabasic rocks, and the Gerizhuoduo area is expected to be a prospective area for porphyry Cu-Au polymetallic deposits.【期刊名称】《物探与化探》【年(卷),期】2011(035)005【总页数】6页(P573-578)【关键词】水系沉积物测量;地球化学特征;找矿远景区;东昆仑布青山地区【作者】崔晓亮;刘婷婷;王文恒;景明;白云【作者单位】成都理工大学地球科学学院,四川成都 610059;成都理工大学地球科学学院,四川成都 610059;四川煤田地质局,四川成都 610072;四川华源矿业勘查开发有限公司,四川成都 610072;四川省国土勘测规划研究院,四川成都610031;成都理工大学地球科学学院,四川成都 610059【正文语种】中文【中图分类】P632东昆仑地区位于青海省中西部，以昆中断裂为界，分别归属华北陆块和扬子陆块的古生代—中生代陆缘活动带，南侧与特提斯构造带紧邻，是欧亚板块与印度板块的碰撞拼接部位，与周边地区共同构成青藏高原北部的重要地质单元，特殊的大地构造位置决定了其构造演化的独特性和复杂性［1－6］。

汶川地震前后地球化学和地球物理异常

汶川地震前后地球化学和地球物理异常飞姓名：田仁仁飞摘要摘要：5.12汶川地震造成了惨重的人员伤亡和巨大的财产损失，是地球深部地质作用过程在浅部（地表）响应的直接后果。

本文总结了汶川地震后，在地表发现或观测到的地球化学特征、地球物理特征，以期了解这些特征的变化规律，为今后的中长期地震预报提供有参考价值的基础资料。

1.引言2008年5月12日14时28分，在青藏高原东缘龙门山地区发生了汶川8 级强烈地震，成为继2001年11月14日昆仑山口西8.1 级地震之后在中国大陆发生的又一次强烈地震，造成了惨重的人员伤亡和巨大的财产损失。

据中国地震局汶川地震科学考察队利用川西流动宽频带地震台阵记录到的近场地震波确定了2008年汶川8.0级地震震源深度：18.66±0.49 km。

这是地球深部地质作用过程在浅部响应，导致大量的地表破裂、房屋倒塌、交通和通讯中断，从而造成了大量的人员伤亡和巨大的经济损失。

地震是地应力达到一定极限，从而地壳快速释放能量，在释放能量的过程中会造成振动。

释放能量越大，振动越强烈，在地表的响应越明显。

此次汶川地震发生后，大量的科研人通过实地考察、各地震台站的资料收集，以及模拟实验等研究了汶川地震的产生机制，为今后类似地震的中长期预报提供了大量可供借鉴资料。

汶川地震前前后浅部地球化学特征2.汶川地震长期的观测结果表明，氡等流体地球化学组分存在显著的震前异常变化。

虽然汶川大地震发生前，未见有报道汶川地区的氡等地球化学异常，但是，大地震发生后，很多科研工作者收集了四川、陕西、青海、宁夏、福建、江苏等地震台站测量氡气变化规律，可以较为明显的发现在地震发生时刻及前后具有较大的异常。

刘耀炜等（2009）通过观测汶川地震后氡气观测值的变化，初步分析了汶川地震氡观测值震后效应的特征：认为氡射气具有对地震波响应的特点，也提出了研究地球化学组分变化与含水层动力过程的关系等问题。

谷懿等（2009）运用活动断裂 (带 )剖面上的土壤氡浓度测量，为地震后大成都地区活动断裂的相对活动性提供了科学评价依据。

西昆仑地区水系沉积物地球化学异常识别

西昆仑地区水系沉积物地球化学异常识别罗伟;李佑国;彭静;张小静【摘要】采用子区中位数衬值滤波法和异常下限衬值滤波法对西昆仑地区1∶20万水系沉积物测量数据进行处理，分别圈定了各元素的衬值异常，并在此基础上将两种方法所获衬值异常的叠合部分作为新的异常边界，重新圈定单元素异常和多元素组合异常。

研究表明，该方法能更好地识别异常，根据异常区元素浓集特征并结合区域地质背景，认为区内Au、Ag、Cu、Pb、Zn、Sb、Sn、W具有较明显的成矿潜力与找矿前景，在找矿方面应重点寻找热液型铜多金属矿、热液型铅锌矿、热液型锑矿、花岗岩型钨锡矿、斑岩型铜矿以及岩金矿床。

%Based on the data of the 1 ∶ 200 000 geochemical stream sediment survey, the authors delineated the single element contrast value anomaly by using the method of "subinterval area median contrast filtering" and "outlier contrast filtering" respectively. The o⁃verlap contrast value anomalies by the two methods were taken as a new anomaly border to re⁃delineate single element anomaly and multi⁃element association anomaly. This study indicates that, by the spatial overlay analysis of contrast value anomalies obtained by the two methods, the overlapped portion can be more effectively identified for the anomaly. According to anomaly characteristics of the con⁃centration elements and regional geological background, It is held that this area has obvious potential mineralization and ore prospects of Au, Ag, Cu, Pb, Zn, Sb, Sn, W, and the prospecting work should focus on hydrothermal Cu polymetallic ore deposit, hydrothermal Pb⁃Zn ore deposit,hydrothermal Sb ore deposit, granite⁃type W⁃Sn ore deposit, porphyry Cu ore deposit and rock gold ore deposit.【期刊名称】《物探与化探》【年(卷),期】2016(040)004【总页数】6页(P722-727)【关键词】子区中位数衬值滤波法;异常下限衬值滤波法;水系沉积物异常;异常识别;找矿潜力【作者】罗伟;李佑国;彭静;张小静【作者单位】成都理工大学地球科学学院，四川成都 610059; 四川省地质与矿产勘查开发局化探队，四川德阳 618000;成都理工大学地球科学学院，四川成都610059;四川省地质与矿产勘查开发局化探队，四川德阳 618000;中国人民警察部队黄金第五支队，陕西西安 710000【正文语种】中文【中图分类】P632西昆仑地区是寻找岩浆型与热液型铜多金属矿产的有利地区[1-4]。

分形奇异(特征)值分解方法与地球物理和地球化学异常重建

第29卷第1期地球科学———中国地质大学学报Vol.29　No.1 2004年1月Earth Science — Journal of China University of G eosciences Jan.　2004分形奇异(特征)值分解方法与地球物理和地球化学异常重建李庆谋1,成秋明1,21.加拿大约克大学地球和大气科学系,多伦多,加拿大M3J1P32.中国地质大学地球系统和矿产资源实验室,湖北武汉430074摘要:地球物理和地球化学异常是找矿的重要依据.地球物理和地球化学异常取决于地层、构造在成矿时间上的多样性与空间上排列、叠置的复杂性.地层、构造因素是构造、岩浆、沉积与成矿地球化学等多种动力学过程的综合反映.这些岩石和构造的因素以及动力过程相互渗透和影响决定了最终地质、地球物理与地球化学场.本文提出的在GIS环境下实现的分形奇异值分解(MSVD)异常重建方法,不仅可以提取地球物理和地球化学等异常,而且能够进一步刻画其中的线性和环状构造、细微的局部纹理结构特征.该方法首先对地球物理和地球化学等网格数据进行二维矩阵的奇异值分解,之后用左特征向量矩阵与右特征向量矩阵的直积构造一个正交完备基.地球物理和地球化学二维数据可以投影到该正交基上,其投影系数是矩阵的奇异值.在该正交完备空间的某些子空间上对地球物理和地球化学等数据进行滤波.为了选择子空间,本文定义了上述正交完备基中的能谱密度、能谱半径(或尺度)与能量测度.在此基础上与空间域及频率域类比,探讨了能量测度与能谱密度呈现分形(fractal和bifractal)规律.利用分形关系的间断点,设计分形奇异值重建算子,实现对地球物理和地球化学异常的分解.以加拿大Nova Scotia南部布格重力异常与As地球化学异常为例,采用MSVD方法分解Au、Wu-Sn-U 等已知矿有关的地球化学异常.发现重建异常能很好地用于解释已知矿点的分布规律.重建的地球化学异常显现了地球化学中的线状和环状异常;重建的布格重力异常有效勾勒出原图中不易发现的纹理结构,这些纹理结构可以合理地解释已知矿点在侵入岩体内及其周围的分布规律.应用实例表明,该方法不仅可以从起因复杂的异常中区分出背景、异常场,还可以识别代表了成矿源岩、流体、运移通道、赋存空间等异常因素引起的纹理、结构与构造特征.同时实现了GIS环境下交互可视化的MSVD处理与解释系统,增强了地质异常定量分析的实用性与可操作性.关键词:异常重建;分形奇异值分解;分形奇异值分解图解;地理信息系统.中图分类号:P628 文章编号:1000-2383(2004)01-0109-10 收稿日期:2003-07-30Fractal Singular2V alue(Egin2V alue)Decomposition Method forG eophysical and G eochemical Anomaly R econstructionL I Qing2mou1,CHEN G Qiu2ming1,21.Depart ment of Earth&A t mospheric Science,York U niversity,Toronto,M3J1P3,Canada2.L aboratory of Earth S ystems and Mi neral Resources,Chi na U niversity of Geosciences,W uhan430074,Chi naAbstract:G eochemical and geophysical anomalies are originated from geological processes.These processes involve a great deal of complexity temporally and spatially.It is critical to improve the current anomaly extraction methods from the stand point of the as2 sociation of geophysical and geochemical anomalies for mineral ex ploration.The fractal singular2value2decomposition(MSVD)in GIS environment developed in this study is demonstrated superb in extracting linear and circular geophysical and geochemical anomalies as well as the detailed structural and textural information from2D geochemical and geophysical maps.The MSVD method constructs a self2contained orthogonal basis using the outer product of left and right eigenvector matrixes decomposed from 2D geochemical or geophysical maps.A power2law relationship based on fractal theory has been suggested to associate the s pectrum基金项目:国家“863”计划课题(No.2002AA135090);OMET基金项目.作者简介:李庆谋(1965-),男,教授,1994年毕业于中国地质大学(北京)能源地质系,获博士学位,现在加拿大约克大学作博士后研究.E2mail:qingmou@yorku.ca,qiuming@yorku.cadensity and spectrum radius(or spectrum scale)defined in the paper.Multiple power2law relationships observed between the s pec2 trum density and spectrum radius can help to group singular values and their corres ponding eigenvectors.Each of these groups can be used to reconstruct the geophysical and geochemical maps to reflect decomposed components.The component reconstructed with relative large singular values may correspond to background and those obtained with relatively small singular values may represent anomalies.This method has been demonstrated using datasets from Nova Scotia,Canada.The results obtained for As and other el2 ements from lake sediment sam ples,gravity anomalies and airborne magnetic anomalies have shown that the power2law relationship might exist between s pectrum density and spectrum radius.Several different ex ponents are observed from the datasets which can be based to separate the anomalies from background.K ey w ords:anomaly reconstruction;fractal SVD(MSVD);MSVD plot;GIS.0　引言在找矿实践中,定量地提取、对比、分析、归纳和总结地球物理、地球化学异常有其理论与实践意义(赵鹏大,2002),是将地质异常理论用于找矿实践并使其量化和具可操作性的关键.因为地球化学元素对矿床的直接指示作用,地球化学异常的提取一直且将继续在地质勘探、异常评价中发挥重要作用.地球化学异常的识别与提取也一直是地球化学勘探的主要目标(Grunsky, 1997;Grunsky and Smee,1999;Harris et al., 1999;Cheng et al.,1994).已经使用多年的直方图、Q-Q图、概率图和方形图(box2plot)等可以用于确定异常的阀值并提取地球化学异常(Stanley and Sinclair,1987;Kurzl,1988;Sinclair,1991; Cheng et al.,1996;G arret,1989).这些异常提取的方法没有利用地球化学中的空间相关关系.以利用半变异函数等为代表的空间分析也可以用以区分地球化学异常(G ovtt et al.,1975;Cheng et al., 1994,1999;Panahi et al.,2003).更进一步,可以根据地球化学的尺度独立性(scale invariance)将地球化学场分解并确定异常(成秋明,2001);也可以将空间变化与几何变化相结合,判断地球化学区域变化与局部变化的阀值,从而进行地球化学异常分离(Cheng,1999;G oncalves et al.,2001).地球物理位场处理中也通过局部曲面拟合、最小曲率样条曲面(Mickus et al.,1991;Wang et al.,1997)、空间域与频率域滤波(Pawlowski and Hamsen,1990;Maus and Dimiri,1996;Maurizio et al,1997)、正反演(Sharma et al.,1999;Schwartz and Broome,1994)或结合各种地质特征更为综合性地提取布格异常(Weiland,1989)等方法提取分析地球物理异常.虽然频率域滤波方法是一种有效的异常分离方法,但如何确定滤波器的(振幅谱)形状、截止频率,甚至其相位谱的取舍一直是频率域滤波中的难点. Maurizio et al.(1997)曾经利用地球物理磁场强度随磁性体埋深变化呈现的指数关系确定磁性体的埋深.20世纪70年代随着计算机广泛普及,广泛开展的位场频率域延拓、水平、垂直与任意方向导数、布格重力异常与磁场强度的相互转换(Mesco,1966; Bhatacharyya,1966;Zurflueh,1967;Fedi et al, 1997;Li and Oldenberg,1998)等都可以作为广义的频率域异常提取的经典实例.甚至,同时期开展的以加强、归位微弱反射为目的的频率域地震偏移成像(Gray et al.,2001)也可以认为是广义的异常重建、增强方法.功率谱的尺度不变性(scale invariance)和功率谱的异向尺度不变性(anisotropic scale invariance)分形规律的发现为解决上述难题提供了有效途径(Cheng,1999;成秋明,2000;Cheng and Li,2002).这一特性导致了频率域分形滤波器的设计(Cheng et al.,1999,2001)并且已经成功地用于地球化学异常的分离(Cheng,1995;Cheng et al.,1999, 2001;Xu and Cheng,2001).由于以下的诸多原因,在频率域处理中,截断误差与边界效应是普遍存在的(Li,1996;李庆谋和刘少华,2000;Li and Cheng,2001):(1)任何小范围内测量的地球物理和地球化学场,或内部存在缺失数据的情况均将产生截断误差;(2)功率谱估计要求地球物理和地球化学场具有二阶平稳性(second order stationary)或起码二阶弱平稳(weak stationary).不完全满足条件的物理化学场也将产生功率谱估计误差(Thomson,1982;Li,1996;李庆谋和刘少华, 2000;K laus and Hentschel,2002);(3)频率域滤波器设计中,虽然可以对功率谱截断处进行光滑处理,但是,对于非规则形状的滤波器(Cheng et al.,1999;011地球科学———中国地质大学学报第29卷Xu and Cheng,2001)进行抛光处理将比较困难.付氏变换的基空间是由相对光滑的正弦和余弦函数系组成的.这些函数系的任意次导数仍然是正弦或余弦函数(或仅有一个相位的变化),或者说仅在正弦与余弦函数系之间相互变换.正弦、余弦函数非常光滑平坦,没有明显的变化边界.这决定了基于Fourier变换的频率域处理方法对地球科学中广为存在的离散的点状、线状、环状构造特征不敏感,而这些特征常常对找矿具有特别的意义(Bonham2 Carter et al.,1985).识别这些特点的异常往往要求高密度的空间数据.本文提出的建立在分形与奇异值分解基础上的地质异常分解方法,从一定程度上可克服基于付氏变换方法的不足,能够较有效地提取地球物理和地球化学异常的线性、局部纹理结构和构造特征,而这些线性、局部纹理、构造与金矿等矿产的分布密切相关.奇异值分解方法已经用于医学图像处理.如核磁共振成像(MRI),X射线层析成像(CA T2scan),正电子辐射影像(PET),光子(正演模拟)合成成像(SPECT)和用于重建组织、血流与骨骼图像的超声成像(William et al.,2002;Maintz and Viergever, 1998)等等.奇异值分解也广泛用于基因排列图像中基因模式的研究(Alter et al.,2000;Holter et al., 2001).奇异值分解在地球物理研究中的应用更为广泛.奇异值分解可用于地震剖面的压缩、增强(Vas2 co,1991;Cagnoli and Ulrych,2001)、垂直地震剖面中弱信号的提取(Freire and Ulrych,1988)、探地雷达弱信号的增强(Michelena,1993)以及重磁及电磁勘探的正反演问题(Groom and Bailey,1991).奇异值分解还用于井间层析成像(Michelena,1993)、地震反演问题中的误差分析(Chang et al.,1999)、分辨率控制(Lebrun et al.,2001)和多分量各向异性地震极化滤波(De and Musacchio,2001)等等.在这些与奇异值分解有关的研究与应用中,一个关键而没有完全解决的问题是如何定量估计奇异值在所解问题中的权重.奇异值分解后,左、右特征向量矩阵的直积构成了正交完备基,原始矩阵的能量在这个正交基中的分形分布为异常的重建提供了关键线索.为了定量研究能量分布,首先需要定义在所构造的正交完备空间中的能谱密度半径(或能谱密度尺度)和能量测度.这样定义的能量测度与能谱密度半径有可能会呈现分形规律.在此基础之上,利用最小平方意义下的最优分割算法,在双对数坐标系下,交互地对能量测度随能谱密度尺度变化的序列进行分割图解(称为分形奇异值图解,MSVD plot).在MSVD plot的帮助下,可以迅速、方便地选择奇异值及其对应的子空间,进行分形奇异值分解后场的重建.在GIS环境下,可以迅速比较分形奇异值分解重建结果与其他相关地质现象的关系,进行地质解释.分形奇异值分解方法可以作为基于GIS的空间处理与分析方法.由于奇异值分解过程需要求解二维矩阵的全部特征值与特征向量,只是在近年随着计算机计算能力的巨大进步,这种巨大的计算量才得以在微型计算机上实现.Cagnoli and Ulrych (2001)在M A TL AB环境下,用M宏语言实现了对地震反射剖面和探地雷达(GPR)剖面中弱信号的提取和加强.1　奇异值分解以及能量分形规律1.1　奇异值分解与正交完备基地球物理和地球化学等二维地学数据均可看作二维矩阵A(m,n).根据代数理论,A(m,n)可以分解为左特征向量矩阵、中间奇异值(特征值)对角矩阵和右特征向量矩阵(Ramachadra and Bhimasan2 karam,2000).设A(m,n)的秩为r,则r≤min (m,n)(如果r=min(m,n),则矩阵A是满秩的),矩阵A(m,n)的奇异值分解形式为A m×n=U m×rΣT r×r V T n×r.(1)上式中,U m×r是左特征向量矩阵,V T n×r是右特征向量矩阵的转置矩阵,Σ=diag(λ1,λ2,…,λr)是奇异值对角矩阵.奇异值的平方是特征值,如果用S V代表奇异值,λ代表特征值,则有如下关系S V i=λi.(2)矩阵的迹等于这些奇异值的平方和t race(A(m,n))=6r i=1λi.(3)如果将U和V写成U=(U1…U r)与V=(V1…V r).U和V的直积(outerproduct)张成子空间ΨkΨ(λk)=U k V T k.(4)所有Ψk(1≤k≤r)构成A(m,n)的一个正交完备111　第1期　李庆谋等:分形奇异(特征)值分解方法与地球物理和地球化学异常重建基Ψ=∪k =1,rΨk .方程(1)又可以表示为A (i ,j )=6rk =1S V k ・Ψ(λk ).(5)(5)式意味着,在Ψ空间中,矩阵A (m ,n )可以被完全重建.如果选择部分奇异值,可以重建含有部分信息(能量)的场.这正是奇异值分解重建异常的作法.但是,怎样选取奇异值,选取什么样的奇异值及其对应的子空间往往是依据经验而定(Cagnoli and Ulrych ,2001).本文试图通过研究Ψ空间中能量的分布,利用能量的分形特性给出定量选取奇异值及其对应的Ψk 子空间的标准,用于重建原始场,以便提取地质异常.1.2　能量分形分布本文通过奇异值分解与Fourier 变换的对比,引入在利用奇异值分解所构造的子空间中能量的概念.设f (x )函数的Fourier 变换为F (f )F (f )=6nk =-n[a k cos (2πkf )+i b k sin (2πkf )].(6)方程(6)中频率k 处的能量密度为P f (k )=a 2k +b 2k .(7)对比方程(5)、(6)和方程(7),可以认为原始矩阵投影到第k 个奇异值对应的子空间上的能量密度为奇异值的平方,等于矩阵的第k 个特征值.ED k =S V 2k =λk .(8)方程(8)中ED k 代表矩阵A 在第k 个奇异值处的能量密度,其意义类似于Fourier 变换中的周期(频率倒数).故可将这一能量密度称为能谱半径.那么,由最小奇异值到第k 个奇异值所对应的子空间的总能量可以表示为(求和区间见图1b )E (S V k |S V≤λ12k )=6ri =kλi,1≤k ≤r .(9)这样,一个能谱半径内对应一个能量.也可以表述为在一定的能谱尺度下,能得到一个能量的测度.对秩为r 的矩阵,可以得到长度为r 的能量测度随能谱尺度变化的序列对.与Fourier 变换中的能谱与频率的关系相类比,这样定义的尺度与能量测度之间可能具有分形规律(图1a ).如果具有分形规律,能量测度随度量尺度变化将呈指数规律(Manderbrot ,1983;Donald ,1997;Pecknoald et al .,图1　奇异值分解能量测度定义(a )以及As 的分形奇异值分解图解(b )Fig.1Definition of accumulated energy (b )and MSVD plot(a )1997;Pilkington ,1995;K ovin ,1992;Lewis et al .,1999):E (S V k )=αλβk .(10)方程(10)在双对数坐标系中对应一条以β为斜率的直线.通常在不同的尺度范围内方程(10)会具有不同的指数值或在双对数坐标系中对应多条斜率不同的直线段(图1a ).1.3　分形奇异值图解(MSV D plot)如图1a 和图2所示,以上述定义的能谱密度为X 坐标,以能量测度为Y 坐标在双对数坐标系中绘制能量测度与尺度序列.在这个图形中,用方差最小(L S )最优分割方法,在人机交互帮助下进行能量测度最优分割图解,称为分形奇异值图解(MSVD plot ):在分形奇异值图解中,如果能量测度与能谱密度呈单一分形关系,则能量测度与能谱密度呈简单的直线关系;如果能量测度与能量尺度在不同的尺度范围内呈不同的分形关系,则能量测度与能谱密度呈多个直线关系.用最小均方误差(L S )拟合多段直线,使总误差平方和最小.图1a ,2是地球化学、地球物理场的能量测度与能谱尺度在双对数坐标中的最优线性分割,同时,可以依据经验进行人机交互调整.各个直线段的端点称为能量测度分布拐点(break ).用这些拐点可以选取相应的奇异值以及对应的子空间,并进行部分重建地球物理和地球化学异常.211地球科学———中国地质大学学报第29卷图2　布格重力异常分形奇异值分解图解Fig.2MSVD plot of Bouguer anomaly 通过对Nova Scotia 地区的布格重力异常、航空磁场强度、航空放射性U 、Th 、K 和Au 、Sb 、W 、Sn 等湖相沉积物地球化学场的试验表明,这样定义的能量测度与尺度关系式(10)具有普遍的多段分形关系.2　研究实例2.1　研究区简介试验研究区位于加拿大Nova Scotia 省西南部(图3a ),面积大约25000km 2(图3b ).区内出露古生代浅变质碎屑沉积岩和晚期花岗质杂岩(图3b ).古生代浅变质沉积地层受到了强烈的北东向褶皱作用,形成了一系列北东走向分布且相互平行的褶皱构造.北西向构造以断裂为主要特征,区内发现多处热液型金矿床(Au )和钨-锡-铀矿床.金矿床和矿点主要分布在花岗杂岩体外一定的范围内,而钨-锡-铀矿化主要发生在花岗杂岩体内部(Bonham 2Carter et al .,1988,1991;Bonham 2Carter ,1994;Xu and Cheng ,2001).这些矿床的分布均受侵入岩体、褶皱构造、断裂构造的明显控制;金矿化和钨-锡-铀矿化有关的蚀变还分别受沉积岩岩相变化和花岗杂岩体内部岩相变化的影响.局部的控矿岩相和构造作用往往会造成与矿化有关元素在岩石以及地表介质中呈现与岩相、构造相关的纹理结构变化.2.2　地球物理、地球化学数据处理应用Kriging 方法在ArcInfo 环境下对1948个地表湖泊沉积物样品中砷(As )等元素进行插值,形成了220m ×220m 的网格图.布格重力异常及其U -Th -K 放射性测量图分辨率为1km ×1km.对网格重力异常进行奇异值分解,共得到205个非零奇异值.在这样分解的基础之上,生成了205个特征子空间(构成一个正交完备基).其中的第1(S V 1)、第4(S V 4)、第10(S V 10)和第30(S V 30)个奇异值以图3　研究区位置、岩性单元和矿点分布Fig.3Location of study area ,simplified lithology units and locations of the known mineral de positsore and mineral occurrences(Xu ,2001)311　第1期　李庆谋等:分形奇异(特征)值分解方法与地球物理和地球化学异常重建图4　布格重力异常的部分奇异值、特征向量与正交子空间(原始布格异常见图7a )Fig.4S ome singular values ,eigenvectors and orthogonal subspaces of Bouguer anomaly (Fig.7a )for the original Bouguer anomaly及所对应的特征向量和特征子空间示于图4中.图4中第1行是左特征向量,第2行是右特征向量,第3行是所标奇异值对应的子空间.随奇异值的减小,子空间的频率迅速增加.大的奇异值对应低频成分,小的奇异值对应高频成分,而信号的能量大部分集中在大的奇异值所对应的低频子空间中(图5a ).应注意的是,该图的纵坐标是对数坐标.前几个奇异值已经代表了信号中绝大部分能量.图5b 表示积分能量,前80个奇异值已经拥有了总能量的93%.大奇异值所对应的信号占总信号的大部分能量.较大奇异值描述场中的长周期变化或背景变化.相对较小的奇异值只占少部分能量,描述信号的高频、局部变化或异常场.湖泊沉积物砷元素分布不仅与金矿化关系密切,并且与该区的地方病有极大的关系(Xu and Cheng ,2001).对As 元素含量地球化学测量的分形奇异值分解图解可分为3段(图1a )(其回归方程和相关系数R 标在图中).结果表明能量测度与能谱密度呈多个直线关系.另外,分别对航空磁场强度、航空放射性(U 、Th 、K )、湖泊沉积物中Au 、Sb 、Sn 等元素分别作了分形奇异值分解图解.这些结果都表现了能量测度图5　地球化学As 元素奇异值分解积分能量贡献(a )与能量密度(b )随奇异值的变化Fig.5Accumulated energy (a )and energy spectrum densi 2ty (b )of As versus the number of singular 2value number411地球科学———中国地质大学学报第29卷图6　As地球化学图(对数变换)(a)、As异常图(b)以及解释(c)Fig.6log2transformed As concentration map(a),MSVD anomaly(b)and interpretation(c)与能谱密度之间的多尺度分形特征.2.3　异常分析图6a是As的地球化学图.图6b是用图1a中第Ⅲ段中的奇异值及其对应的子空间重建的As异常图.图6c是解释图.图6c中标示1区异常较图6a中的原生晕异常明显缩小了范围,并能很好地解释已知的几个金矿点.异常带2也有类似的反映.线状异常10,11是所围面状异常的边缘,能很好地解释这些金矿点沿As 异常边缘带分布的规律.图6c中线状的1～7条带与原生晕(图6a)对比有了极大的细化,反应构造变化更明显、清楚.1～3显然是全区北东-南西向构造褶皱的结果.4～7是局部北西-南东构造活动的图7　布格重力异常图(a)、MSVD重建结果(b)以及解释(c)Fig.7Bouguer anomaly(a),MSVD reconstructed ma p(b) and interpretation(c)反映,它们控制了该局部地区的成矿作用.线6的拐弯趋势非常特别,反映全区北东-南西褶皱构造与局部北西-南东构造对As分布的综合影响结果.图7a是研究区的布格重力异常与已知矿点的叠合图.图7b是用图2中第Ⅱ段奇异值与对应的子空间重建的异常场,将岩性地层边界叠加在其上.图7c是解释结果.图7a中布格重力异常平坦宽缓,与已知矿点无法有效对应.图7b的分形奇异值分解重建结果给出了丰富的异常信息.图7c标示1～5处于沉积变质岩区,强烈的北东-南西向构造作用控制了成矿过程,这一结构在图7a中没有反映.标示6～17是岩511　第1期　李庆谋等:分形奇异(特征)值分解方法与地球物理和地球化学异常重建浆侵入岩体与围岩(沉积变质岩)接触带,主要发育Au矿床.分形奇异值重建的异常给出了侵入岩周围详细的结构变化线条.从这个实例来看,在侵入岩体拐弯、分叉与变化复杂的地方,侵入岩与围岩有更大的接触面积,提供了更多的热液活动成矿机会.这些线状、小面状异常纹理并不完全与侵入体的露头吻合,可能反映隐伏侵入体或岩体分支.图7c反映了岩体内部由于岩相相变、热液活动、构造变化等造成的纹理结构(标示18～20),同时较好地刻画了分布在岩体内部的Sn-W-U矿床分布规律(Xu and Cheng,2001).3　讨论起因复杂的成矿作用在时间上具有多期次、空间上呈渐变、穿插、叠合复杂性的特征.所产生的地质、地球物理和地球化学异常极其复杂.从什么角度(如分布密度、几何形状、能量密度结构等),在什么空间(可以是重构的正交空间,如本文)观察这些地质现象,决定了能否抓住异常形成过程中固有的、起主要决定作用的因素,并进而与找矿联系起来.本文就是在这样的总体思路下,尝试了用奇异值(特征值)分解及其能量分形特征选择子空间分解异常.对研究区内的地球物理和地球化学异常的分解结果表明,该方法能较好地刻画二维场内部的异常纹理结构特征.致谢:感谢徐亚光博士对使用Nova Scotia地区的地质、地球物理与地球化学数据所提供的帮助.本文的研究得到了加拿大NSERC(Natural Sciences and Engineering Research Council)基金和OM ET (A T:http://laurentian.ca/geology/OM ET/)基金的部分资助.R eferencesAlter,O.,Brown,P.O.,Boststein,D.,2000.Singular value de2 composition for genome2side expression processing and mod2 eling.Proc.N A,97:10101-10106. Bhatacharyya,B.K.,1966.Continuous s pectrum of the total magnetic field anomaly due to a rectangle prismatic body.Geophysics,31(1):97-121.Bonham2Carter,G.F.,Rencz,A.N.,Harrris,J.R.,et al., 1985.Spatial relationship of gold occurrence with lineament derived from LANDSA T and SEASA T imagery,MegumaGroup,Nova Scotia.In:Bonham2Carter,G.F.,ed.,Pro2 ceedings of the fourth thematic conference of remote sensing for exploration geology,San Francisco.A GU,San Francis2 co,USA,755-768.Bonham2Carter,G.F.,1994.G eographic information systems for geoscientist:Modeling with GIS.Pergamon,Oxford,398. Bonham2Carter,G.F.,Agterberg,F.P.,Wright,D.F.,1988.Integration of geological datasets for golden exploration in Nova Scotia.Phonogram Remote Sensing,54:1585-1592. Bonham2Carter,G.F.,Agterberg,F.P.,Wright,D.F.,1991.Mineral potential mapping with GIS.Unpublished notes fora short course,November11-15,1991,At De pt.of G eolo2gy,University of Ottawa,109.Cagnoli,B.,Ulrych,T.J.,2001.Singular value decomposition and wavy reflections in ground2penetrating radar images of base surge deposits.Journal of A pplied Geophysics,48:175 -182.Chang,X.,Lu,M.,Liu,Y.,1999.Error analysis and appraisal of three general solution in seismic tomography.Acta Geo2 physica S inica,42(5):695-701.Cheng,Q.M.,2001.Multifractal and geostatistic methods for characterizing local structure and singularity properties of exploration geochemical anomalies.Earth Science—Journal of China U niversity of Geosciences,26(2):161-166(in Chinese with English abstract).Cheng,Q.M.,2000,Multifractal theory and geochemical ele2 ment distribution pattern.Earth Science—Journal of Chi2 na U niversity of Geosciences,25(3):311-318(in Chinese with English abstract).Cheng,Q.M.,1999.S patial and scaling modeling for geochemi2 cal anomaly separation.Journal of Geochemical Ex plo2 ration,63(3):175-194.Cheng,Q.M.,Agterberg,F.P.,Ballantyne,B.S.,1994.The separation of geochemical anomalies from background by fractal method.Journal of Geochemical Ex ploration,43(2):91-109.Cheng,Q.M.,Agterberg,F.P.,Bonham2Carter,G.F.,1996.A spatial analysis method for geochemical anomaly separa2tion.Journal of Geochemical Ex ploration,56:183-195. Cheng,Q.M.,1995.The perimeter2area fractal model and its application to geology.M athem atical Geology,27(1):69-82.Cheng,Q.M.,Li,Q.M.,2002.A fractal concentration2area method for assigning color palette for image representation.Com puter&Geosciences,28:567-575.Cheng,Q.M.,Xu,Y.,Grunsky,E.,1999.Integrated spatial and spectrum analysis for geochemical anomaly separation.611地球科学———中国地质大学学报第29卷In:Lippard,J.L.,Naess,A.,Sinding2Lonsen,R.,eds., Proc.Int.Assoc.Mathematical G eology,Trondheim,Nor2 wayⅠ,87-92.Cheng,Q.M.,Xu,Y.,Grunsky,E.,2001.Multifractal power spectrum2area method for geochemical anomaly separation.N atural Resources Research,9(1):43-51.De,F.R.,Musacchio,G.,2001.Polarization filter with sin gular value decomposition.Geophysics,66(3):932-938. Donald,L.T.,1997.Fractal and chaos in geology and geo2 physics.Second edition.Cambridge University Press,Lon2 don,398.Fedi,M.,Quarta,T.,Santis,A.D.,1997.Inherent power2law behavior of magnetic field power spectra from the s pector and grant ensemble.Geophysics,64(4):1143-1150. Freire,S.L.M.,Ulrych,T.J.,1988.Application of singular value decomposition to vertical seismic profiling.Geo2 physics,53(6):778-785.G arret,R.G.,1989.The Chi2square plot:A tool for multivariateoutlier recognition.Journal of Geochemical Ex ploration,32 (1-3):319-341.G oncalves,M.A.,Mateus,A.,Oliveira,V.,2001.G eochemicalanomaly separation by multi2fractal modeling.Journal of Geochemical Ex ploration,72(2):91-114.G ovtt,G.J.S.,G oodfellow,W.D.,Cha pman,R.P.,et al.,1975.Exploration geochemistry distribution of elements and recognition of anomalies.M athem atical Geology,7:415-446.Gray,S.H.,Etgen,J.,Dellinger,J.,et al.,2001.Seismic mi2 gration problems and solutions.Geophysics,66(5):1622-1640.Groom,R.W.,Bailey,R.C.,1991.Analytic investigations of the effects of near surface three2dimensional galvanic scat2 ters on M T tensor decompositions.Geophysics,56(4):496 -518.Grunsky,E.,1997.Numerical techniques and strategies for in2 terpretation of geochemical data current topics in GIS and integration of exploration datasets short course notes.Ex plo2 ration97Workshop,Ottawa,G AC,Canada,1-147. Grunsky,E.,Smee,B.,1999.Differentiation of soil types and mineralization from multi2element geochemistry using mul2 tivariate methods and digital topography.Journal of Geo2 chemist ry Ex ploration,67(1-3):289-301.Harris,J.R.,Wlkinson,L.,Grunsky,E.,et al.,1999.Tech2 niques for analysis and visualization of litho geochemical data with application to Swayze G reenstone Belt,Ontario.Journal of Geochemical Exploration,67(1-3):301-344.Holter,N.S.,Maritan,A.,Cie plak,M.,et al.,2001.Dynamicmodeling of gene expression data.Proc.N atl.Acd.Sci.USA,98:1693-1698.Kurzl,H.,1988.Exploratory data analysis:Recent advances for the interpretation of geochemical data.Journal of Geochem2 ical Ex ploration,30(3):143-163.K laus,H.,Hentschel,K.,2002.Ma pping the spectrum:Tech2 niques of visual representation in research and teaching.Ox2 ford,New Y ork,362.K ovin,G.,1992.Fractal models in the earth science.Elsevier, Amsterdam,408.Lebrun,D.,Richard,V.,Mace,D.,et al.,2001.SVD for mul2 tioffset linarized inversion,resolution analysis multicompo2 nent acquisition.Geophysics,66(3):871-882.Lewis,G.M.,Lovejoy,S.,Schertzer,D.,et al.,1999.The scale invariant generator technique for quantifying anistropic scale puters&Geosciences,25(9):963-978.Li,Q.M.,ankovitch cycles analysis and application by logging data.Chinese Journal of Geophysics,36(5):699-704.Li,Y.G.,Oldenberg,W.,1998.Separation of regional and residual magnetic field data.Geophysics,62(2):431-439. Li,Q.M.,Cheng,Q.M.,2001.Fractal correction of well log2 ging curves.Journal of China U niversity of Geosciences,12(3):272-275.Li,Q.M.,Liu,S.H.,2000.Wavelet scalogram of well logging.Progress of Geophysics,14(1):25-39(in Chinese with English abstract).Maintz,J.,Viergever,M.,1998.A survey of medical image registration.Medical Im age A nalysis,2:1-36. Manderbrot,B.B.,1983.The fractal geometry of nature.Free2 man,New Y ork,468.Maurizio,F.,Tatiana,Q.,Angelo,D.S.,1997.Inherent pow2 er2law behavior of magnetic field power spectra from a spec2 tor and Grant ensemble.Geophysics,62(4):1143-1149. Maus,S.,Dimiri,V.,1996.De pth estimation from scaling spec2 trum of potential field.Geophysical Journal of Internation2 al,124(1):113-120.Mesco,C.A.,1966.Two dimensional filtering and the second derivative method.Geophysics,31(3):606-617. Mickus,K.L.,Aiken.C.L.V.,K ennedy.W.D.,1991.Re2 gional2residual gravity anomaly separation using the mini2 mum2curvature technique.Geophysics,56(2):279-283. Michelena,R.J.,1993.Singular value decomposition for cross2 well tomography.Geophysics,58(11):1655-1661. Panahi,A.,Cheng,Q.M.,Bonham2Carter,G. F.,2003.Modeling lake sediment geochemical distribution using prin2711　第1期　李庆谋等:分形奇异(特征)值分解方法与地球物理和地球化学异常重建。

变换后数据的因子分区标准化在东昆仑东段地球化学异常圈定

变换后数据的因子分区标准化在东昆仑东段地球化学异常圈定中的应用
耿国帅, 杨帆, 郭建娜
1. 中国地质大学(北京) 地球科学与资源学院,北京 100083 2. 中国人民武装警察部队黄金地质研究所,河北廊坊 065000 3. 北京矿产地质研究院,北京 100012 4. 中国地质调查局土地质量地球化学调查评价研究中心,河北廊坊 065000 5. 河北省廊坊
物探与化探,2020,44(1),112-121.DOI:10.11720/wtyht.2020.1014
ห้องสมุดไป่ตู้
市自然资源和规划局,河北廊坊 065000
图7 传统方法与因子分区标准化方法圈定的金异常对比 a—因子分区标准化圈定的金异常点及金矿床点;b—传统方法圈定的金异常点及金矿床点;c—两种方法都圈定出的金异常点及金矿床点;d—只
在其中一种方法中圈定的金异常点

东昆仑东段东昆中构造混杂岩带长石山蛇绿岩年代学、地球化学特征

东昆仑东段东昆中构造混杂岩带长石山蛇绿岩年代学、地球化
学特征
东昆仑东段是中国的一个著名的古生代构造岩石带之一，其地质构造非常复杂，这里的岩石种类很多，其中最典型的岩石类型是长石山蛇绿岩。

长石山蛇绿岩是一种含有大量绿色矿物质的岩石，主要由蛇纹石、辉绿岩及长英质岩石构成。

长石山蛇绿岩的形成年代较为古老，其形成时代主要为古生代晚期，距今已经有近五亿年的历史。

通过对长石山蛇绿岩进行详细的地球化学分析还可以发现，它是一种典型的富硅碱性岩石，它的成分主要由SiO2、Al2O3、Na2O、K2O等组成。

此外，长石山蛇绿岩中还富含一些微量元素，如V、Cr、Ni等，这些元素也是一些重要的地球化学特征之一。

长石山蛇绿岩在东昆仑东段的分布非常广泛，这一区域的地质条件和构造背景是长石山蛇绿岩形成的主要原因。

长石山蛇绿岩的形成和地质演化与构造运动密切相关，这些构造运动包括横向挤压、纵向拉伸、剪切等多种形式。

这些构造运动通过改变地层的位置和形态，影响到了长石山蛇绿岩的地质过程和成因。

总的来说，东昆仑东段是一个非常值得研究的地质地貌区，长石山蛇绿岩作为其一个典型代表岩石，具有重要的地质学和地球化学研究意义。

通过对这一区域的研究，可以进一步了解长石山蛇绿岩的形成演化过程，深入探究其地球化学特征和地质意义。

3-06区域地球化学勘查规范

区域地球化学勘查规范比例尺：1：200000DZ／T 0167－19951 主题内容与适用范围1.1 本规范对区域地球化学勘查（简称区域化探）工作的基本任务，工作性质、设计书的编写，中低山、丘陵区、特殊景观区的野外工作方法，样品加工和样品库，区域化探中的岩石样品采集、野外工作质量检查，样品多元素分析及质量监控。

质量评估，地球化学图和原始图件的编制，异常查证，地球化学图说明书的编写等作了规定，确立了统一标准。

1.2 本规范主要适用于地质矿产行业的区域地质矿产调查中进行的1：200000区域地球化学勘查工作，亦可供比例尺1：500000的同类工作或其他行业进行类似工作时参考使用。

2 引用标准GB9649 地质矿产术语分类代码DZ／T0075DZ／T00113 总则3.1 区域化探属基础地质矿产调查工作之一，应覆盖全部可工作的我国国土面积。

它的主要工作目的是发现由金属或非金属成矿区（带）、矿田和大、中型矿床以及某些地层、构造和火成岩的区域地球学化特征所引起的省的、区域的和局部地球化学异常，并为基础地质研究等领域提供某些基础地球化学资料。

常按国际分幅范围部署测区。

我国常用工作比例尺为1：200000，交通很不发达的边缘地区同时也采用1：500000，它们的采样密度（以水系沉积物测量为例）分别为：0.25～1点／km2和0.04～0.10点／km2。

3.2 按地理特点和区域化探工作条件，可把我国大致分为以下两类地区：a. 中低山、丘陵地区（见图１）；b. 具有各料特殊地理景观条件的边缘地区。

所谓具有特殊地理景观条件的边缘地区主要是指分布在黑龙江、内蒙古东部等地的森林沼泽、内蒙古（大青山南坡除外）、甘肃兰州以西、宁夏、新疆、青海柴达木盆地及其周缘的干旱、半干旱荒漠、青海、西藏、新疆、川西、甘肃祁连山等地的的高寒山区，西藏西部的刘寒湖沼荒漠、广西、贵州、云南等地的岩溶，云南西南及海南岛等地的热带雨林，我国西北地区的黄土高原，以及内地沿海冲洪积平原区等。

青海省东昆仑地区地球化学异常特征及金矿靶区筛选与评价

青海省东昆仑地区地球化学异常特征及金矿靶区筛选与评价安国英【摘要】利用东昆仑地区1∶50万化探扫面成果,研究该区成矿元素的地球化学分布规律,并重点选择与金矿有关的成矿元素Au、Cu、Co、Mo、Pb、Zn、Ag、As、Sb、Hg等进行讨论.从元素地球化学特征入手,分析它们在空间、时间上的分散、富集规律,为区域矿床地球化学预测提供理论依据.以本区两类典型金矿床模型为依据,利用GIS技术进行金矿靶区的筛选与评价,提出了Ⅰ级金矿找矿远景区3个,Ⅱ级远景区6个,Ⅲ级远景区7个,作为下一步矿产勘查工作的重点.【期刊名称】《物探与化探》【年(卷),期】2013(037)002【总页数】7页(P218-224)【关键词】地球化学;金矿;成矿远景;青海东昆仑地区【作者】安国英【作者单位】中国国土资源航空物探遥感中心,北京 100083【正文语种】中文【中图分类】P632东昆仑地区位于青海省中西部，属昆仑山东段，地理坐标为35°20'～36°40'N，93°00'～99°00'E，包括12幅1∶20万区域化探扫面图幅，实际采样面积61 270 km2。

青海省地球化学勘查技术院分别于1992年和1997年两次对该区涉及的39种元素和氧化物进行过1∶50万和1∶100万地球化学图的编制。

笔者在前人工作的基础上，从元素地球化学特征分析入手，对东昆仑地区成矿地质环境进行探讨，并根据地球化学综合异常特征，利用多元地学信息对本区两类典型金矿远景区进行筛选与评价，为进一步的找矿工作提供指导。

1 地球化学特征1.1 元素地球化学丰度特征研究区地处干旱荒漠与高寒山区景观区，化探工作方法以水系沉积物测量为主，统计区内39种指标的含量(表1)。

与全国水系沉积物中各指标含量对比发现，研究区水系沉积物中大部分指标含量低于全国水平。

以全国水系沉积物中指标含量的平均值为标准，计算东昆仑地区各指标的浓集系数，其中浓集系数大于1 的指标为 Sr、Sb、Au、As、Bi、Cd、F、CaO、Na2O、MgO、SiO2;小于0.7 的指标为 Ag、Zr、Hg;浓集系数在0.99 ～0.7 的指标为 Ba、W、Mo、Li、Cu、Ni、Be、P、Y、La、U、Zn、Cr、Nb、Mn、Pb、V、Co、Sn、B、Ti、Th、K2O、Fe2O3、Al2O3等。

青海东昆仑东段区域成矿地质特征

浅谈青海东昆仑东段区域成矿地质特征摘要：本文主要阐述了青海东昆仑东段区域地质矿产概况、区域地球物理化学特征区域成矿地质特征，对该地区开展地质矿产勘查工作具有一定的指导意义。

关键词：青海；东昆仑；东段区域；成矿地质特征中图分类号：f407.1 文献标识码：a 文章编号：1 区域地质矿产概况1.1区域构造体系东昆仑东段地处中国大陆中央造山带西段的东昆仑加里东造山带东段，北与柴达木盆地毗邻，南与可可西里)巴颜喀拉山毗邻，东以鄂拉山断裂为界，西至青藏沿线。

东昆仑东段横跨古亚洲和特提斯两大构造域，分别以昆北、昆中、昆南、北巴断裂为界划分为昆北、昆中、昆南和阿尼玛卿4个三级构造带。

东昆仑东段是晚加里东造山带、晚华力西—印支构造—岩浆活化带(花岗岩及陆相火山岩)，昆北构造带由元古界裂解残块(金水口群)、早古生代裂陷槽沉积(滩间山群)和大量的中—酸性岩浆岩组成；昆中微陆块由前寒武纪基底变质岩系和各时代侵入杂岩组成，该带中nw—nww 向大型剪切带发育，且与金矿关系密切；昆南构造带由中元古代—早古生代洋壳沉积(万保沟群)，晚古生代—早中生代活动陆缘沉积组成，与钴、铜、金矿关系密切；阿尼玛卿造山带是晚古生代—早中生代构造混杂岩带，代表洋壳的蛇绿岩与弧后火山岩及岛弧火山岩在这一个构造带中共存，与铜、金矿关系密切。

1.2区域地质构造演化区域地质构造演化可划分成前寒武纪、早古生代、晚古生代—早中生代、晚中生代—新生代四个发展演化阶段。

前寒武纪新太古代—古元古代形成统一的稳定陆台后，被克拉通化并成为基底；中—新元古代沿昆南构造带发生初始裂解形成昆仑洋；早古生代寒武—奥陶纪进一步拉伸裂解呈多岛洋/裂陷槽，奥陶纪末大洋俯冲、闭合，形成沟—弧—盆体系，志留纪进一步隆升造山；晚古生代阿尼玛卿裂谷开合、华力西期造山活动引发强烈的岩浆侵入活动；早中生代受特提斯洋俯冲、消减、闭合及印支期造山活动影响，昆北、昆中发生大规模岩浆侵入活动，昆南及阿尼玛卿构造带则形成陆缘活动带；晚中生代—新生代全区隆升，处于陆内叠复造山阶段。

奇异性指数理论在东天山覆盖区地球化学数据分析中的应用

奇异性指数理论在东天山覆盖区地球化学数据分析中的应用李逸川;陈志军;徐红燕;李诗珺;丛丽娟【摘要】The singularity indices can be used to identify regional geochemical anomalies and the degree of sin -gularity.Because of its property of scale invariance and self-similarity,the singularity indices play an important role in geo-anomalies extraction in covered area.The basic principles and method of singularity indices are introduced by a case study of Cu geochemical-anomaly identification in Chihu-Sanchakou area in Eastern Tianshan.A compar-ative analysis between traditional anomalies mapping method and singularity indices mapping method is also carried out.The results show that the singularity indices analysis is effective for weak geochemical anomalies identification in covered area.%奇异性指数分析方法可以用于地球化学局部异常的圈定和奇异性程度的度量,其具有尺度不变性和自相似性的特征,使之能够对覆盖区地球化学弱异常信息进行有效的提取.在总结了奇异性指数理论基本原理的基础上,以赤湖-三岔口地区地球化学数据为例,对该区铜元素进行奇异性指数填图,并对奇异性指数填图法与传统异常填图法进行对比分析,结果显示奇异性指数能够有效地在覆盖区提取地球化学元素弱异常信息.【期刊名称】《科学技术与工程》【年(卷),期】2017(017)030【总页数】5页(P9-13)【关键词】东天山;铜镍硫化物矿床;地球化学;奇异性指数【作者】李逸川;陈志军;徐红燕;李诗珺;丛丽娟【作者单位】中国国土资源航空物探遥感中心,北京100083;中国地质大学地质过程与矿产资源国家重点实验室,武汉430074;中国地质调查局地学文献中心,北京100083;中国国土资源航空物探遥感中心,北京100083;中国国土资源航空物探遥感中心,北京100083【正文语种】中文【中图分类】P596东天山地区位于哈萨克斯坦与准噶尔板块的结合部位，其矿产资源十分丰富，是重要的石油、煤、铀、铁、铜、镍、金等成矿远景区和矿业基地。

青海东昆仑地区1：2.5万地球化学测量方法技术及应用成果

沉积物测量不能完全满足地质找矿的需要。青海省１：２．５万地球化学测量以水系沉积物测多点组合的采样方法，其对１：５万水系沉积物测量综合异常具有较好的分
解和重现性，并且新圈出地球化学综合异常。１：２．５万地球化学测量综合异常浓集中心更加明
显、强度更高，对异常源的定位更准。通过近年来在东昆仑地区开展１：２．５万地球化学测量工作，
发现了果洛龙洼大型金矿，按纳格、瓦勒尕中型金矿，浪木日铜镍矿点等一系列矿床点。因此，１：
２．５万地球化学测量工作方法在阿尔金、柴北缘等类似地球化学景观区得到推广应用。
ＳｕｒｖｅｙｉｎＥａｓｔＫｕｎｌｕｎ，ＱｉｎｇｈａｉＰｒｏｖｉｎｃｅ
ＺＨＡＯＪｕａｎ，ＸＵＧｕａｎｇ。，ＹＡＮＧＢａｏｒｏｎｇ。，ＧＵＯＨａｉｍｉｎｇ。，ＬＩＤｅｂｉａｏ，ＭＡＺｈｅｎｇｔｉｎｇ，ＷＥＩＬｉｑｉｏｎｇ，ＸＩＯＮＧＳｈｅｎｇｙｕｎ。
（１．ＱｉｎｇｈａｉＧｅｏｌｏｇｉｃａｌＳｕｒｖｅｙ，Ｘｉｎｉｎｇ８１０００８，Ｑｉｎｇｈａｉ，Ｃｈｉｎａ；２．ＱｉｎｇｈａｉＮｏ．５ＧｅｏｌｏｇｉｃａｌＭｉｎｅｒａｌＥｘｐｌｏｒａｔｉｏｎＩｎｓｔｉｔｕｔｅ，Ｘｉｎｉｎｇ８１０００３，Ｑｉｎｇｈａｉ，Ｃｈｉｎａ；３．Ｎｏ．８ＴｅａｍｏｆＧｅｏｌｏｇｙａｎｄＭｉｎｅｒａｌＲｅｓｏｕｒｃｅｓＢｕｒｅａｕｆｏｒＮｏｎ—ｆｅｒｒｏｕｓＭｅｔａｌｓｏｆＱｉｎｇｈａｉＰｒｏｖｉｎｃｅ，Ｘｉｎｉｎｇ８１００１２，Ｑｉｎｇｈａｉ，Ｃｈｉｎａ）

东昆仑东段杨树沟金矿地质特征及找矿方向

东昆仑东段杨树沟金矿地质特征及找矿方向田永革;陈秉芳;韩晓龙;喇品贤【摘要】The Yangshugou gold deposit is situated in the east section of the Indosinian Kalike-Xiangride Zu-Pb-Zn-Fe-graphite (-Cu-RE-REE) mineralization belt,Eastern Kunlun metallogenic province.The orebodies are strictly controlled by the E-W-and NE-trending fault structures.The gold deposit is of structural altered rock type,with metallogenic epoch of Early Indosinian.The deep part of existed ore belts and the southeast part of the orefield where faults are developed are favorable areas for prospecting.Further exploration should focus on the ore-bearing fracture zones,with engineering and geochemical survey in the southeast and deep of the orefield,to find more gold orebodies.%杨树沟金矿位于东昆仑成矿省伯喀里克-香日德印支期金、铅、锌、铁、石墨(铜、稀有、稀土)成矿带的东段.矿体严格受近东西向、北东向断裂构造控制.成矿类型为构造蚀变岩型金矿床,成矿时代为印支早期.矿区中已发现的金矿带深部及矿区东南部断裂发育地区具有较好的找矿空间,下一步工作应针对含矿构造破碎带,采用探矿工程与化探工作相结合等手段开展矿区东南部及深部找矿工作,寻找金矿体,争取该矿区找矿的更大突破.【期刊名称】《地质与资源》【年(卷),期】2017(026)003【总页数】5页(P248-252)【关键词】东昆仑;杨树沟金矿;构造蚀变岩型;找矿方向【作者】田永革;陈秉芳;韩晓龙;喇品贤【作者单位】青海省青藏高原北部地质过程与矿产资源重点实验室,青海西宁810012;青海省地质调查院,青海西宁810012;青海省青藏高原北部地质过程与矿产资源重点实验室,青海西宁810012;青海省地质调查院,青海西宁810012;青海省青藏高原北部地质过程与矿产资源重点实验室,青海西宁810012;青海省地质调查院,青海西宁810012;青海省青藏高原北部地质过程与矿产资源重点实验室,青海西宁810012;青海省地质调查院,青海西宁810012【正文语种】中文【中图分类】P618.51杨树沟金矿位于青海省都兰县西南约70 km，属于东昆仑成矿省伯喀里克－香日德印支期金、铅、锌、铁、石墨（铜、稀有、稀土）成矿带的东段①易平乾，等.全国矿产资源潜力评价项目“青海省重要矿种区域成矿规律研究”成果报告.青海省地质矿产勘查开发局，2013.②田永革，等.青海省都兰县杨树沟地区金矿预查报告.青海省地质调查院，2015..该成矿带夹持于昆中断裂和昆北断裂带之间［1］，呈东西向展布，长约800 km，宽20～50 km.1990年代后，在成矿带的中部发现评价了以五龙沟金矿［2-3］为代表的一批金、多金属矿床点，2000年后，又在东部的沟里地区［4］发现了一批大中型构造蚀变岩型金矿床，成矿潜力巨大.而杨树沟矿区正处于五龙沟－沟里地区的中间位置，成矿背景也较为相似.本文从区域地质、矿区地质、矿体特征等方面对矿床成因进行探讨，为下一步勘查工作提供参考.杨树沟金矿大地构造位置处于西域板块东昆中陆块岩浆弧带中东段，位于昆中和昆北断裂带之间的华力西岩浆岩带上［1］.区域岩浆活动强烈，岩体出露面积较大，主要以华力西早期及印支期的中－酸性侵入岩为主.构造以北西向、北西西向断裂为主，北东向断裂构造次之.区域上地层出露较全，但面积不大，从古元古代金水口岩群到第四系地层均有不同程度的出露.矿产以金为主，其次为银、铅锌、磁铁矿等，成因类型主要为构造蚀变岩型［2-7］.已发现的矿床大多与区域断裂构造及次级裂隙构造有关，如五龙沟金矿［2-3］、果洛龙洼金矿［8-10］、阿斯哈金矿［11-12］、瓦勒尕金矿、按纳格金矿、洪水河金矿、巴隆金矿等，矿床规模自小型至大型.2.1 地层矿区出露地层较为单一，为古元古代金水口（岩）群片麻岩组及新生代第四系地层（图1）.古元古代金水口（岩）群片麻岩组（Pt1Ja），仅在矿区北西部有小面积出露.该地层呈残留体产于华力西期花岗岩中，呈近东西向展布，长350～400 m，宽40～60 m.岩性为灰黑色角闪斜长片麻岩，与围岩呈侵入接触.第四系地层分为上更新统-全新统（Qp-hal＋pl）与全新统（Qhal），主要分布于矿区中西部山前平地、沟谷及现代河床等地，由冲洪积砂、砾石以及风成黄土等组成.2.2 构造矿区位于东昆仑山北坡断裂隆起带，昆北断裂和昆中断裂带之间，位于哈图弧形构造的西翼转折端［13］.受此影响，本区以褶皱构造不发育，而断裂构造极为发育为特征.总体构造线方向呈近东西向，与区域构造线方向基本一致.野外工作中在矿区共圈定断裂构造破碎带25条，长度一般在100～1000 m，宽2～10 m，最长为2.5 km.按展布方向分为北西向、近东西向、北东向3组.其中以北西向、近东西向为主，是早期破碎带，北东向为后期构造（见表1）.由表1中可见，一般近东西向、北西向断裂构造是区域性断裂，规模较大；北东向断裂构造规模较小；而南北向及其他小型的裂隙构造是大型构造的次级构造，规模更小.从含矿性分析，各类型的断裂构造均有金矿化的显示，这些断裂构造不仅是矿液的重要运移通道，也是矿区主要的金矿存储空间.2.3 岩浆岩矿区岩浆活动强烈，侵入岩出露面积约占矿区面积的3/4以上，主要为华力西期中酸性侵入岩.岩性主要为花岗闪长岩（γδ3c），其次为花岗岩（γ3b）等，岩体主44要呈岩基状产出，形态不规则.局部受动力变质影响，断裂附近和扭动构造发育地区，岩石多具压性破碎构造以及多种围岩蚀变.内生热液金属矿产主要与华力西期中酸性岩及其派生之脉岩关系极为密切.矿区中脉岩也较发育，主要为闪长玢岩脉（δμ），次为二长花岗岩脉（ξγ）、钾长花岗岩脉（γ）、辉绿玢岩脉（βμ）、石英脉，还有少量的花岗细晶岩脉等.局部闪长玢岩脉呈脉状或扁豆状分布在断裂破碎带两边或附近，与金矿化体有一定关系.3.1 地球化学特征3.1.1 1∶5万水系沉积物异常矿区内包含有异常3处，其中两处异常找矿成果较好.异常以As为主元素，Au为组合元素.Au元素峰值7×10－9，均值6.57×10－9；As元素峰值220×10－6，均值42.88×10－6.异常呈椭圆状，形态规整，异常浓度分带不明显，浓集中心较清晰，Au、As元素异常套合紧密.异常为As单元素异常，As元素峰值93.2×10－6，均值40× 10－6.异常呈椭圆状，形态规整，异常浓度分带及浓集中心不明显.3.1.2 1∶1万土壤异常以往工作针对异常开展了1∶1万土壤测量，在异常内圈定出以W、Pb、Au为主元素的综合异常4处（图2）.各元素异常峰值为Cu 482×10－6，W 749×10－6，Pb 307×10－6，Au 155×10－9.各主元素异常均具有十分明显的三级浓度分带，浓集中心明显.其中Au异常与矿区金矿带等套合较好，局部的Pb异常与金矿带套合较好，而Cu、W元素异常与金矿带及Au异常没有明显的相关性.3.1.3 岩石剖面异常针对矿区构造破碎带进行的1∶5000岩石地化剖面工作中，除已有金矿化带中显示有Au异常外，在矿区东部又发现7处Au高值异常段，Au含量为62×10－9～1235×10－9，高值段均与F7、F22、F23等断层破碎带相对应.岩性为碎裂花岗闪长岩、碎裂云英岩等.3.2 矿（化）体特征3.2.1 矿（化）体分布在矿区以Au 0.1×10－6品位圈定出金矿（化）带8条（K1～K8）.矿化带一般长80～800 m，宽0.5～10 m，最大宽度可达18.2 m.其中异常内5条As（Au）异常中3条.金矿（化）带呈近东西向、北东向及北西向展布，均分布于构造破碎蚀变带内.以Au品位1×10－6分别在K1、K2、K4、K5、K7矿化带中圈定出金矿体5条（M1～M5）.矿体长度在80～800 m不等，平均厚度在0.37～2.52 m之间.矿体平均品位Au 1.13×10－6～2.77×10－6，单样最高5.02×10－6.其中M3为矿区主矿体（表2）.矿化蚀变类型有：褐铁矿化、黄钾铁矾化、碳酸盐化、高岭土化、硅化、绿泥石化、碎粒岩化、糜棱岩化，局部见孔雀石化、黄铁矿化.M3金矿体（图1、3）呈条带状产于AS1丁11As异常F1断裂破碎带之K4金矿（化）带中，呈北西向展布，矿体长800 m，平均厚2.15 m，Au品位0.42×10－6～3.97× 10－6，平均品位1.89×10－6，产状215°∠（85～90°）.赋矿岩石为碎裂云英岩及碎裂岩化花岗闪长岩，岩石遭受强烈挤压，呈碎粒状、粉末状.矿化为黄钾铁钒化、褐铁矿化等；蚀变为高岭土化、碳酸盐化、硅化、绿泥石化等.地表以橘红色、褐红色、局部浅灰色为特征.3.2.2 矿石特征矿石类型按容矿岩石划分为碎裂云英岩型、碎裂闪长玢岩型、碎裂花岗闪长岩型金矿石.金矿石均产于构造破碎带中，破碎强烈.矿石金属矿物主要有黄铁矿、褐铁矿、黑钨矿、金红石、锐钛矿、毒砂矿.金矿化与黄铁矿、黑钨矿、毒砂矿以及褐铁矿的含量呈正相关，与金红石和锐钛矿的含量关系不明.脉石矿物为石英、绢云母、碳酸盐岩、绿泥石、绿帘石等.矿石结构主要有自形晶粒结构、半自形粒状结构及土状结构，其中以土状结构和半自形粒状结构最具代表性.矿石构造类型主要有土状构造、稀疏星点状构造以及稀疏细脉浸染状构造.3.2.3 矿床成因分析从矿体展布形态分析，矿（化）体呈脉状、透镜状产于破碎蚀变带中，并严格受近东西向、北东向断裂构造控制，显示断裂控矿特征；矿石及围岩具强烈的黄钾铁钒化、褐铁矿化、黄铁矿化、高岭土化、碳酸盐化、硅化、绿泥石化等矿化蚀变类型，反映出热液蚀变特点.因此，本矿床成因应属构造蚀变岩型金矿.参考区域上产于同一成矿带的五龙沟金矿蚀变绢云母的Ar－Ar坪年龄（236.5±0.5）Ma，Ar－Ar等时线年龄（238.6 ±3.5）Ma［5，9］，应与本区成矿时间接近，为印支早期. 3.3 围岩矿化蚀变金矿体围岩主要为褐红色碎裂云英岩、碎裂岩化花岗闪长岩以及灰黑色碎裂闪长玢岩.金矿体与围岩的界线不甚清晰，局部呈渐变过度接触，但围岩中的金元素含量明显低于蚀变带内.金矿体附近围岩矿化蚀变较强烈，远离矿体则矿化蚀变减弱.矿区围岩矿化类型主要为黄铁矿化、毒砂矿化、黄钾铁矾化、褐铁矿化.从探矿工程中发现，矿化类型、强度等由浅至深具有分带性：褐铁矿化→黄钾铁矾化→硫酸盐矿化.而黄铁矿化主要发育在较深部的原生环境，一般深度大于20 m，与金矿体关系较密切.矿区围岩蚀变具有多阶段性，主要有硅化、云英岩化、碳酸盐化，次为绿帘石化、高岭土化等.从已发现的金矿特征看，含金构造破碎带中蚀变越强的地带越容易成矿，蚀变较弱的地带金矿化也变弱.杨树沟金矿（化）体严格受断裂构造控制，因此，在以后工作中应继续针对矿区断裂构造开展工作.根据目前矿区工作成果，认为该区下一步找矿方向应注重以下几个方面.（1）以往在矿区东部发现有7处Au岩石样高值段，证明与其相对应的F7、F22、F23等断层具有Au矿化显示.从含矿破碎带规模及带中岩石样分析结果分析，这些金矿带应具有较大的规模，下一步应重点对其进行勘查.（2）地表已发现有8条金矿化带、5条金矿体，应对具有找矿前景的矿化体深部进行勘探.如M1矿体的北东部、M3矿体的北西部等.（3）矿区中的M3矿体及F7、F22、F23等含金破碎带，均位于1∶5万水系Au 异常以外.因此下一步应重视对没有Au异常的异常区及其以南区域的勘查，必要时可进行1∶1万岩屑测量及激电测量等，以寻找新的金矿体.通过对杨树沟矿区地质特征等研究，指出了下一步找矿方向，同时总结出该矿床具有以下特点.（1）矿床中含金破碎带矿化很不均匀，无论从金矿体的倾向还是走向统计，金含量变化较大.（2）大多金矿体地表矿化较弱，而深部金矿化较强，说明金矿体地表氧化淋滤较为强烈，下一步应加强金矿体深部的勘查.（3）含金破碎带产状不稳定.如M3金矿体的中部及东部为南倾，而在矿体的北西部则变为北倾斜.综上所述，杨树沟目前已发现8条金矿化带、5条金矿体，加之邻区已发现有大中型金矿床多处，且成矿背景及成矿类型极为相似，因此该矿区具有广阔的找矿前景.相信通过合理的工作方法，必将会取得找矿的重大突破.【相关文献】［1］张雪亭，杨生德，杨站君，等.青海省板块构造研究：1∶100万青海省大地构造图说明书［M］.北京：地质出版社，2007：15－130.［2］李厚民，沈远超，胡正国，等.青海东昆仑五龙沟金矿床成矿条件及成矿机理［J］.地质与勘探，2001，37（1）：65－69.［3］李厚民，钱壮志，刘继庆，等.东昆仑中带找寻五龙沟式金矿的思路［J］.青海地质，1999（2）：47－52.［4］王凤林，肖小强，陈世顺，等.青海沟里地区金矿地质特征及找矿前景［J］.黄金科学技术，2011，19（4）：45－48.［5］张德全，党兴彦，佘宏全，等.柴北缘－东昆仑地区造山型金矿床的Ar－Ar测年及其地质意义［J］.矿床地质，2005，24（2）：87－98.［6］吴庭祥，张绍宁，安汝龙，等.青海东昆仑东段金矿区地层含矿性分析［J］.矿产与地质，2009，23（5）：431－441.［7］李发明，董毅，何财福，等.青海省都兰县沟里地区沟系次生晕找矿实例［J］.矿产与地质，2008，22（5）：456－460.［8］杨宝荣，杨小斌.青海都兰果洛龙洼金矿床地质特征及控矿因素浅析［J］.黄金科学技术，2007，15（1）：26－30.［9］邹定喜，杨小斌，芦文泉.青海都兰果洛龙洼金矿床同位素特征及成因［J］.黄金科学技术，2011，19（2）：26－30.［10］丁清峰，金圣凯，王冠，等.青海省都兰县果洛龙洼金矿成矿流体［J］.吉林大学学报：地球科学版，2013，43（2）：415－426.［11］李碧乐，沈鑫，陈广俊，等.青海东昆仑阿斯哈金矿Ⅰ号脉成矿流体地球化学特征和矿床成因［J］.吉林大学学报：地球科学版，2012，42 （6）：1676－1687.［12］祁月清，刘瑞雪.都兰县阿斯哈金矿地质特征及找矿远景浅析［J］.西部探矿工程，2012，24（6）：140－149.［13］王仁祥，刘冰，王红梅，等.青海省东昆仑昆中断裂带哈图段地质特征及成矿预测［J］.山东国土资源，2012，28（2）：6－10.。

地球化学异常结构模式及应用

Mo-Sn组合异常
Mo、Sn是碱金属元素，它们的多元素组合异常分布较为简单，而且和 Ni-Ti-Co组合不同，它有三个较强的浓集中心，其中强度较大的两个位于区域的东部， Mo和Sn都是高温元素，应该代表一种高温环境，这可能说明了本区的岩浆侵位是由东向西的。
四、因子分析与多元素组合分布特征
二、多元素异常组合的划分
• 化验了Cu、Pb、Zn、Au、Ag、Sb、 As、Bi、Hg、W、Fe203、Mo、Sn、 Ba、Co、Cr、Ni、Sr、Mn、Ti等20 个元素。
• 根据各元素与成矿元素Cu、Au的相关关系、各元素之间的地球化学性质以及它们与地质背景的关系，划分了多元素异常组合。
应用实例
青海德合龙洼铜矿床是位于青海黄南州的一个岩浆热液型矿床，大地构造位置处于位于西秦岭与南祁连两大构造单元的结合部，为秦岭、祁连山、昆仑山成矿带的交汇部位,属于西秦岭成矿区。
图1 德合龙洼矿区区域地质图
一、样品采集
本次研究在德合龙洼现有矿区及其附近，采用1：5000详查，点线距25～ 50m，在测定的采样点周围点线距的 2.5～5m范围内采样，主要采集距地表 20～50cm深处土壤的B层（淋滤层）中的细粒级物质，取样重量不少于500g，在特殊地貌区应根据不同自然地理条件采用不同的采样方法，本区采集植物根以下的残、坡积土，尽量不要带入腐殖质和碎石。共采集样品1228件。
成矿元素异常是指目前已知的达工业品位的成矿元素异常；
指示元素异常包括直接指示元素异常和间接指示元素异常，直接指示元素异常是指主要成矿元素和次要成矿元素及少量伴生元素异常，间接指示元素异常是指与成矿作用有关的重要指示元素异常；
成矿环境元素异常是指那些反映成矿环境和岩性的典型元素异常。

东昆仑造山带东段大格勒花岗岩锆石U-Pb年代学、地球化学特征及其构造意义

第４２卷　第４期２０２０年７月地球科学与环境学报ＪｏｕｒｎａｌｏｆＥａｒｔｈＳｃｉｅｎｃｅｓａｎｄＥｎｖｉｒｏｎｍｅｎｔＶｏｌ．４２　Ｎｏ．４Ｊｕｌｙ２０２０封　铿，李瑞保，裴先治，等．东昆仑造山带东段大格勒花岗岩锆石ＵＰｂ年代学、地球化学特征及其构造意义［Ｊ］．地球科学与环境学报，２０２０，４２（４）：４４２４６３．ＦＥＮＧＫｅｎｇ，ＬＩＲｕｉｂａｏ，ＰＥＩＸｉａｎｚｈｉ，ｅｔａｌ．ＺｉｒｃｏｎＵＰｂＤａｔｉｎｇａｎｄＧｅｏｃｈｅｍｉｃａｌＣｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆＤａｇｅｌｅＧｒａｎｉｔｅｉｎｔｈｅＥａｓｔｅｒｎＭａｒｇｉｎｏｆＥａｓｔＫｕｎｌｕｎＯｒｏｇｅｎｉｃＢｅｌｔ，ＣｈｉｎａａｎｄＴｈｅｉｒＴｅｃｔｏｎｉｃＩｍｐｌｉｃａｔｉｏｎｓ［Ｊ］．ＪｏｕｒｎａｌｏｆＥａｒｔｈＳｃｉｅｎｃｅｓａｎｄＥｎｖｉｒｏｎｍｅｎｔ，２０２０，４２（４）：４４２４６３．收稿日期：２０２００５１０；修回日期：２０２００６２２；网络首发日期：２０２００７２４投稿网址：ｈｔｔｐ：∥ｊｅｓｅ．ｃｈｄ．ｅｄｕ．ｃｎ／基金项目：国家自然科学基金项目（４１５０２１９１，４１４７２１９１，４１８７２２３３）；青海省国土资源厅中国铝业公司公益性区域地质矿产调查项目（中铝基金２００８０１）；中央高校高新技术研究培育项目（３００１０２２７９２０４）；中央高校基本科研业务费专项资金项目（３００１０２２７９１０８）；陕西省自然科学基金项目（２０２０ＪＭ２２９，２０２０ＪＱ３５３，２０１９ＪＭ３１２）作者简介：封　铿（１９９５），男，陕西延安人，理学硕士研究生，Ｅｍａｉｌ：ｆｅｎｇｋｅｎｇ５２１＠１６３．ｃｏｍ。

通讯作者：李瑞保（１９８２），男，山西朔州人，副教授，理学博士，Ｅｍａｉｌ：ｌｉｒｕｉｂａｏ０９７１＠１６３．ｃｏｍ。

犇犗犐：１０．１９８１４／ｊ．ｊｅｓｅ．２０２０．０５００９东昆仑造山带东段大格勒花岗岩锆石犝犘犫年代学、地球化学特征及其构造意义封　铿１，李瑞保１，２，裴先治１，２，李佐臣１，２，裴　磊１，２，刘成军１，２，陈国超１，李小兵１，王　盟１，２，赵少伟１，２，陈有?１，２，周　海１，２（１．长安大学地球科学与资源学院，陕西西安　７１００５４；２．长安大学西部矿产资源与地质工程教育部重点实验室，陕西西安　７１００５４）摘　要：大格勒花岗岩是东昆仑弧岩浆岩的重要组成部分，对解译古特提斯演化具有重要意义。

东昆仑塔格地区变质基性火山岩地球化学及演化特征

作者简介:(1994-),,。

、、。

1地质特征塔格地区火山岩主要为中元古代的一套变质基性火山岩，赋存于狼牙山组地层中[1]。

本次研究区查明的火山岩形态主要呈夹层状或透镜状，主要分布在狼牙山组上部的碳酸岩中。

而狼牙山组的地层由于遭受了多期次的后期构造叠加改造作用，其岩石发生了强烈的变形，从而导致局部可见有发育大量的韧性剪切带[2]。

从而致使其中的火山岩也多发生有变质变形，其部分变为糜棱岩类。

由于岩石成岩时代久远，多期次的构造叠加变形明显，使查明其古火山机构相当困难[3]。

2岩石学特征含绿泥石绿帘石片岩：形态呈细粒鳞片状变晶结构，呈片状构造；其中绿帘石含量在65%左右、石英在15%左右、钠长石及绿泥石均为10%左右；粒度以细粒为主，大小在0.08mm~0.48mm，其中绿帘石呈粒状变晶结构，较均匀的分布在岩石中，具有明显的定向排列特征，石英具有粒状变晶结构，大多数多为粒状集合体。

3岩石地球化学特征（1）岩石化学特征李润根1，宋小坤2，江小强1(1.成都理工大学地球科学学院,四川成都610059;2.四川省众成矿业有限公司,四川成都610045)摘要：通过进行1:5万区调等相关工作，在青海格尔木市塔格地区发现一套主要由变质基性火山岩组成的绿片岩类，本文对其进行了岩石地球化学特征研究。

其特征表明该火山岩为拉斑玄武岩系列；稀土元素特征表明近似于洋岛碱性玄武岩；变质演化总体有变形强，但变质程度低的特征。

关键词：塔格地区；变质基性火山岩；绿片岩类中图分类号：P588.14文献标志码：AGeochemistry and Evolution Characteristics of Metamorphic Basic Volcanic Rocks in the Tage Area,East KunlunLI Rungen 1,SONG Xiaokun 2,JIANG Xiaoqiang 1(1.Chengdu University of technology,Chengdu 610059,China;2.Sichuan Zhongcheng Mining Co.,Ltd.,Chengdu,610045,China)Abstract :A set of greenschist mainly composed of metamorphic basic volcanic rocks was discovered in the tagg area of geermu city,qinghai province through 1:50000district adjustment and other related work.The characteristics show that the volcanic rocks are tholeiite series.The characteristics of rare earth elements are similar to ocean island alkaline basalt.The metamorphism evolution is characterized by strong deformation but low metamorphism.Key words:tag area;metamorphic basic volcanic rock;greenschist2东昆仑塔格地区变质基性火山岩地球化学及演化特征狼牙山组3件绿片岩样品岩石化学特征见表1，由表可知岩石的SiO 2含量为39.84%~48.64%；其Al 2O 3含量为11.56%~13.47%；而TiO 2含量大于1%，其中大部分K 2O的含量大于Na 2O的含量，在尼格里参数中：氧化铝数t=-2~19.7，属于铝过饱和类型，石英数qz=-19.85~13.58，为典型的SiO 2不饱和型，表明岩石中有石英较少，与岩石镜下鉴定结果一致。

元素迁移的分维结构,级序路径及共轭地球化学异常存在的理论依据

均跳跃距离的平方成正比。
１．２扩散方程霹其数学解
对扩散宏观现象数学处理的主要基础是Ｆｉｃｋ定律：
第一作者简舟：）哥永章男３Ｏ岁博士、研究员数学与计算机导向地球化学、矿床地球化学 · 中国科学院留学择忧支持基金、院厂基金和中国科学院矿床地球化学开放研究宣资助项目
维普资讯
地
球
化
学
元由亚晶体、亚颗粒及线状和面状缺陷组成，线状或面状缺陷是对应级别的陷阱。再上一级别是晶体或矿物颗粒规模这一层敬，对应级别的陷阱是颗粒边界。再往上则包括显微裂隙岩石、手标本规模破碎岩石、大断裂地块等一系列层敬，与它们对应的陷阱是显微裂隙、手标本规模破裂和大规模断裂带等。如果把陷阱镶嵌结构中某一个层敬放大或缩小，则上下不同层次的镶嵌图案是充分相似的。这种充满空间的自相似性是数学分维的最主要特征之一。可以认为，多层次陷阱镶嵌结构是一数学分维结构。
关键词歪墨星蔓迁移陷阱坌丝堕塑埋星￡＆主毛－共轭挫彗业主是童
许多金属矿床的形成涉及到矿源岩（层）及微量元素从矿源岩（层）中的初始位置最后成矿位置的迁移问题。前人研究表明，相当大比例的成矿元素是以杂质组分分布于各类造岩矿物中的，成矿元素从造岩矿物中活化出来，以及此后的迁移对矿床的形成是不可缺少的。本文着重从数学观点探讨矿源岩（层）中微量元素在矿物结构不破坏条件下构迁移路径及其总体图象。

东昆仑东段浪木日地区地球化学异常信息提取方法研究及异常查证

东昆仑东段浪木日地区地球化学异常信息提取方法研究及异常
查证
李斌;高强;魏俊浩;赖联新;李笑龙
【期刊名称】《矿产勘查》
【年(卷),期】2024(15)3
【摘要】青海沟里地区位于东昆仑造山带东段,根据水系沉积物地球化学测量已发现了多个矿床(点),具有较好的找矿潜力。

本文选取沟里地区巴加别里赤尔幅(浪木日地区)1∶5万水系沉积物地球化学测量数据为研究对象,利用统计学方法和趋势面分析法两种方法进行异常信息提取。

结果显示,趋势面分析法相较于传统统计法,圈定的异常面积更大、浓集中心更明确、分带特征更明显,且与已知矿床点的吻合度更高,整体上具有强化低背景异常的优势。

同时利用趋势面分析法提取了3个元素组合,并对套合性和相关性突出,且与区内浪木日铜镍硫化物矿床相关的Cu-Zn-Cr-Co-Ni组合异常进行异常剖析。

结合区内地质特征和浪木日矿区矿体分布情况,圈定了哲扎空龙洼东、龙里、瑙格木东和卡鲁南等找矿远景区,为区内下一步铜镍找矿工作提供了方向。

【总页数】17页(P403-419)
【作者】李斌;高强;魏俊浩;赖联新;李笑龙
【作者单位】中国地质大学(武汉)资源学院
【正文语种】中文
【中图分类】P596
【相关文献】
1.东昆仑东大滩地质遥感异常查证与近红外矿谱仪应用研究
2.西昆仑叶尔羌河上游地区区域化探异常查证方法及找矿效果
3.ILR变换后数据的因子分区标准化在东昆仑东段地球化学异常圈定中的应用
4.新疆东天山地区土屋、延东铜矿地球化学特征与异常查证方法
5.东昆仑东段浪木日锰镍多金属矿产勘查中的物探技术应用及研究
因版权原因，仅展示原文概要，查看原文内容请购买。

KNN改进算法在圈定地球化学元素异常区域上的应用

KNN改进算法在圈定地球化学元素异常区域上的应用
曹赛男;张乾荣;刘斌;周仲礼
【期刊名称】《中国科技论文》
【年(卷),期】2016(011)015
【摘要】以青海省大柴旦镇塔塔棱河中下游一带为研究区域,通过开展对1∶10 000土壤地球化学测量,采用KNN分类算法求取Au、Cu、Zn、Pb、As、Sb 6种元素的异常下限值,并作出异常图。

通过后期工程验证,可以看出圈定的异常图对后期找矿具有明显的指示意义,说明运用KNN分类算法求取异常下限和确定异常区域是可行的。

【总页数】4页(P1782-1785)
【作者】曹赛男;张乾荣;刘斌;周仲礼
【作者单位】[1]成都理工大学管理科学学院,成都610059;[2]成都理工大学数学地质四川省重点实验室,成都610059
【正文语种】中文
【中图分类】P595.3
【相关文献】
1.“三江”地区北段区域地球化学元素组合异常提取及其找矿意义 [J], 丛源;陈建平;肖克炎;董庆吉
2.核主成分分析在地球化学元素组合异常圈定中的应用 [J], 黄志强;任鸿飞;梁芳敏;魏继祖;段长生
3.基于实码加速遗传算法的投影寻踪模型在圈定水系沉积物地球化学异常中的应用
——以湖南某铅锌矿床为例 [J], 郝立波;田密;赵新运;赵昕;张瑞森;谷雪
4.KNN改进算法在圈定地球化学元素异常区域上的应用 [J], 曹赛男;张乾荣;刘斌;周仲礼
5.分区上异点校正法在干旱地区1∶5万地球化学异常圈定中的应用:以我国西北某地区为例 [J], 石文杰;魏俊浩;王启;张德才;刘国春;于海涛;陈冲
因版权原因，仅展示原文概要，查看原文内容请购买。