自1979年区域化探全国扫面计划开始和自1999年多目标区域地球化学调查开始至今，我国已获得覆盖陆地国土7百余万平方千米的水系沉积物和土壤样品的39项元素分析数据(Xie et al., 1997; 谢学锦和周国华，2002; 杨忠芳等，2007; 谢学锦等，2009; Gong et al., 2013; Xie and Cheng， 2014; Li et al., 2014)。所取得的海量地球化学数据在我国矿产资源勘查、基础地质研究、环境、生态及农业等研究方面发挥了重要的作用(Cheng et al., 2014; Yang et al., 2014a; Yang and Badal， 2013; Deng et al., 2010; Gong et al., 2008)。全国矿产资源潜力评价项目自2006年执行以来，全国化探工作组与大区、省级化探工作者共同补充完善了我国区域化探数据库，制定了化探资料应用技术流程(向运川等，2010)。这一成果为海量区域地球化学数据在基础地质、矿产勘查以及环境、生态等方面的应用提供了可靠的数 据。

区域地球化学数据是地表物质化学组成的直接反映，具有不同化学组成的地表物质其差异必将反映在区域地球化学数据之中，因此利用区域化探数据的空间分布特征可以推断地表地质体的空间分布规律。花岗岩类侵入体不仅是地表出露的主要岩石类型之一，而且与许多大型、超大型矿床关系密切(Hua et al., 2004; Yang et al., 2014b; Deng et al., 2014a; Wang et al., 2014a)。花岗岩类侵入体的主要地球化学特征在氧化物方面表现为富K₂O贫Fe₂O₃(区域化探全国扫面计划所分析的7项氧化物不包括FeO)，在微量元素方面表现为富高场强元素贫相容元素等(Yan and Chi， 2005; 史长义等，2005)。本文研究主要是基于花岗岩类侵入体的这一地球化学特征，依据区域化探数据在模型区构建推断花岗岩类侵入体的地球化学模型，进而将其应用到研究区。 2 模型区与研究区

根据全国矿产资源潜力评价化探资料应用研究成果，全国共划分出5个地球化学域和25个地球化学省(表 1、图 1)。全国5个一级地球化学域按照编号分别为1古亚洲地球化学域、2秦祁昆地球化学域、3特提斯地球化学域、4扬子地球化学域和5滨太平洋地球化学域。全国25个二级地球化学省编号采用2位编码，第1位为地球化学域编码，第2位为地球化学省编码，编码自北向南、自西向东顺序编号(图 1)。

图 1

Fig. 1

图 1 全国地球化学域和地球化学省划分方案 地球化学图采用GeoExpl软件中累频19级方法绘制，色标从深蓝色渐变至深红色代表从低值渐变至高值.此图为全国矿产资源潜力评价化探资料应用研究成果 Fig. 1 Classification of geochemical domain and geochemical province in Chin l and surface

表 1

Table 1

表 1(Table 1)

表 1 全国二级地球化学省划分方案Table 1 Classification of geochemical province in China l and surface 编码 地球化学省名称 11 阿尔泰地球化学省 12 准噶尔地球化学省 13 天山地球化学省 14 内蒙古高原地球化学省 15 阴山地球化学省 16 浑善达克沙地地球化学省 17 大兴安岭地球化学省 21 祁连山-昆仑山地球化学省 22 秦岭地球化学省 23 伏牛山-大别山地球化学省 31 阿尼玛卿山地球化学省 32 巴颜喀拉山-无量山地球化学省 33 唐古拉山地球化学省 34 冈底斯山-念唐-腾冲地球化学省 35 喜马拉雅山地球化学省 36 大雪山-哀牢山地球化学省 41 四川盆地地球化学省 42 大凉山-乌蒙山-大娄山地球化学省 43 苗岭-六诏山地球化学省 51 长白山-小兴安岭地球化学省 52 东北-华北平原地球化学省 53 黄土高原地球化学省 54 山东丘陵地球化学省 55 巫山-雪峰山地球化学省 56 华南丘陵地球化学省

表 1 全国二级地球化学省划分方案 Table 1 Classification of geochemical province in China l and surface

图 1全国地球化学分区图中底图分别采用K₂O和Fe₂O₃的地球化学图来表示，其目的有二。一是全国地球化学分区图的划分必然参考了K₂O和Fe₂O₃的地球化学图，据此可以辨析地球化学分区的合理性；另一是花岗岩类侵入体具有富K₂O贫Fe₂O₃的特征，据此可以推断可能的花岗岩类侵入体分布区。

对比图 1中的K₂O和Fe₂O₃的地球化学图可以看出存在明显富K₂O贫Fe₂O₃的地区有三处：大兴安岭地球化学省(图 1中17)、华南丘陵地球化学省(图 1中56)和云南西部的腾冲地球化学亚省(图 1中34南段)-哀牢山地球化学亚省(图 1中36南段)。

在大兴安岭地球化学省由于受森林景观影响目前对其基础地质研究程度不很高，不考虑将其作为已知区开展建模研究。在华南丘陵地球化学省目前对其基础地质研究程度高，尤其是对花岗岩类侵入体的空间分布也已基本查明，但花岗岩类侵入体在空间上呈零散片状分布，规律性不很明显，不易选为模型区。在云南西部的腾冲地球化学亚省和哀牢山地球化学亚省，花岗岩类侵入体分布相对比较集中，基本呈弧形带状和条带状分布，规律性比较强(图 2a)。因此本文选择云南西部的腾冲地球化学亚省和哀牢山地球化学亚省为模型区来构建花岗岩类侵入体的区域地球化学推断模型，进而选择华南丘陵地球化学省的中段为模型应用研究区以检验模型的适用性。

云南西部的腾冲地球化学亚省和哀牢山地球化学亚省是特提斯地球化学域的南段，即三江特提斯造山带和成矿带(Deng et al., 2013a，2014b)。三江特提斯成矿带以大型、超大型Au、Cu、Pb-Zn等矿床而备受关注，尤其是哀牢山金成矿带(Yang et al., 2010; 杨立强等， 2011a，b; 邓军等， 2012，2013b; Wang et al., 2013b，2014b; 张静等，2010; Zhang et al., 2013)。在澜沧江成矿带发育有云南著名的临沧花岗岩体和勐海花岗岩体(Gong et al., 2011)，在腾冲地区发育有东河-明光花岗岩带、古永花岗岩带和槟榔江花岗岩带(江彪等，2012; Wang et al., 2013b)。 3 花岗岩类侵入体地球化学推断模型 3.1 模型区地球化学特征

基于花岗岩类侵入体富K₂O和高场强元素贫Fe₂O₃和相容元素的地球化学特征，在全国矿产资源潜力评价成果之全国地球化学图中直接裁剪出云南西部地区部分高场强元素和相容元素(即在基性岩中富集和在花岗岩中贫化的元素)的地球化学图分别如图 2和图 3所示。

图 2

Fig. 2

图 2 全国侵入岩及高场强元素地球化学图中云南西部地区视图 图(a)侵入岩分布图引自全国矿产资源潜力评价研究成果，粉红色区域为花岗岩类侵入体分布区.绘图数据来源于全国区域化探数据库Fig. 2 Map of intrusions and geochemical map of high field strength elements in the western of Yunnan Province，China

图 3

Fig. 3

图 3 全国相容元素地球化学图中云南西部地区视图 图区范围与图 2一致.绘图数据来源于全国区域化探数据库 Fig. 3 Geochemical map of compatible elements in the western of Yunnan Province，China

从图 2可以看出，高场强元素Th、U、La、Y、Zr在澜沧江带均呈带状显著富集，在腾冲地区均呈弧形带状显著富集，其富集区与该区花岗岩类侵入体的空间分布区相吻合。

在图 3中，相容元素Co、Ni、V、Cr、Ti、Fe₂O₃在澜沧江带均出现贫化现象，但在其东侧出现带状富集现象；在腾冲地区相容元素均出现显著贫化现象，但在保山地区及其南部相容元素均出现显著富集的现象。

将图 2与图 3对比可以发现，相容元素Co、Ni、V、Cr、Ti、Fe₂O₃的贫化区均位于花岗岩类侵入体的分布区。该区花岗岩类侵入体的带状空间分布决定了高场强元素与相容元素的带状分布特征，这两类元素的富集与贫化区在云南西部地区似乎具有相间排列的规律。 3.2 模型构建

在云南西部模型区内上述高场强元素和相容元素富集与贫化具有相间排列的规律主要取决于该区花岗岩类侵入体的带状空间分布。区域化探全国扫面计划规定分析测试项目为39项，为了更好地揭示元素含量与矿床及其它地质体的关系，许多研究者提出了新的数据处理方法技术(Hao et al., 2014; Ma et al., 2014; Deng et al., 2011; Sun et al., 2010; 韩东昱等，2004)。本文为揭示反映花岗岩类侵入体的分析指标，此处利用因子分析首先查明指示花岗岩类侵入体的分析指标，进而依据因子载荷矩阵建立回归分析模型获得指示花岗岩类侵入体的综合指标。

在全国区域化探数据库中检索出云南西部模型区的39项分析测试数据，对模型区39项指标进行因子分析，采用主成分分析法利用特征根大于1提取公因子并采用方差极大法进行因子旋转获得旋转因子载荷矩阵(表 2)。

表 2

Table 2

表 2(Table 2)

表 2 云南西部地区旋转因子载荷矩阵 Table 2 Rotated component matrix of geochemical survey data in the western of Yunnan Province，China 变量 Ag As Au B Ba Be Bi Cd Co Cr Cu F Hg La Li Mn Mo Nb Ni P F1 0.04 0.15 0.02 -0.09 0.02 -0.01 0.00 0.02 0.92 0.81 0.58 0.07 0.02 0.12 0.05 0.69 0.30 0.57 0.80 0.49 F2 0.00 0.01 0.00 -0.18 0.35 0.53 0.02 -0.01 -0.08 -0.05 -0.11 0.20 -0.02 0.67 -0.09 -0.07 0.03 0.32 -0.05 0.10 F3 0.42 0.19 0.03 0.01 0.04 -0.01 0.03 0.92 0.01 0.03 0.02 -0.01 -0.09 -0.02 0.00 0.07 0.05 -0.01 0.02 0.00 F4 -0.09 0.23 -0.04 0.51 0.02 0.50 -0.02 0.04 -0.01 -0.05 -0.09 0.61 -0.02 -0.03 0.70 0.15 0.21 0.06 0.02 0.13 F5 0.00 0.06 0.03 0.02 -0.05 -0.07 -0.02 0.00 0.15 0.24 0.02 0.11 -0.03 0.04 -0.10 -0.09 -0.11 -0.20 0.23 0.03 F6 0.21 0.10 -0.01 0.21 0.00 0.27 0.77 -0.04 -0.01 0.00 0.25 0.21 -0.01 0.02 -0.08 -0.06 -0.02 0.06 -0.02 -0.02 F7 -0.42 -0.56 0.04 0.05 0.01 0.03 -0.18 0.02 -0.05 0.01 -0.37 -0.05 -0.64 0.02 -0.17 -0.15 -0.13 0.07 -0.02 -0.01 F8 0.04 0.04 -0.03 -0.01 0.77 0.00 0.04 0.02 0.05 -0.06 0.06 0.21 -0.01 0.40 -0.03 0.15 0.08 -0.05 -0.05 0.59 F9 0.05 0.27 -0.03 0.36 0.07 -0.33 -0.02 -0.01 0.03 0.10 -0.06 -0.08 -0.15 0.07 -0.06 0.16 0.09 -0.38 0.16 -0.03 F10 -0.30 -0.26 -0.81 0.05 0.04 0.00 -0.03 0.07 0.01 0.06 -0.18 -0.08 0.25 0.03 0.02 -0.12 -0.37 -0.12 0.02 -0.08 变量 Pb Sb Sn Sr Th Ti U V W Y Zn Zr SiO2 Al2O3 Fe2O3 K2O Na2O CaO MgO F1 0.01 0.05 0.00 0.21 -0.12 0.93 -0.14 0.93 -0.02 0.11 0.08 -0.06 -0.72 0.37 0.95 -0.27 -0.12 0.17 0.30 F2 0.05 0.00 0.04 -0.09 0.87 -0.02 0.77 -0.05 0.14 0.86 0.00 0.77 -0.29 0.58 0.00 0.61 0.18 -0.12 -0.03 F3 0.83 0.23 0.00 0.08 0.06 0.00 -0.02 0.03 0.02 0.04 0.96 -0.02 -0.02 -0.01 0.02 -0.02 0.01 0.00 0.02 F4 -0.04 0.10 0.12 0.01 -0.07 -0.10 0.11 -0.01 0.35 0.04 0.03 -0.30 -0.21 0.34 0.02 0.30 0.06 -0.02 -0.01 F5 0.02 0.02 -0.02 0.23 0.00 -0.04 -0.02 0.08 0.05 -0.01 0.00 -0.06 -0.30 -0.12 0.04 -0.14 0.18 0.84 0.83 F6 0.06 -0.05 0.80 -0.01 0.09 -0.01 0.05 -0.02 0.42 0.01 -0.01 0.04 0.06 -0.11 0.02 -0.01 0.03 -0.02 0.01 F7 -0.22 -0.64 0.07 -0.02 0.00 0.03 0.03 -0.02 0.11 -0.05 -0.04 0.05 0.03 0.01 -0.03 0.03 0.07 -0.04 0.03 F8 0.02 -0.02 0.00 0.68 0.00 0.08 -0.08 0.08 -0.11 -0.02 0.04 0.10 -0.21 0.15 0.10 0.21 0.13 0.03 0.08 F9 0.01 0.10 0.01 -0.35 -0.08 -0.09 -0.28 -0.05 0.02 0.04 -0.02 0.08 0.04 -0.06 -0.01 -0.27 -0.75 -0.11 -0.03 F10 -0.13 -0.06 0.07 -0.03 -0.01 -0.02 -0.09 -0.03 -0.07 -0.01 0.00 0.04 0.04 0.02 -0.02 0.01 0.05 -0.05 0.06

表 2 云南西部地区旋转因子载荷矩阵 Table 2 Rotated component matrix of geochemical survey data in the western of Yunnan Province，China

在表 2中，F1主要代表Fe₂O₃、Ti、V、Co、Cr、Ni、SiO₂、Mn、Cu、Nb计10项指标，主要为相容元素组合(除SiO₂、Nb外)。F2主要代表Th、Y、U、Zr、La、K₂O、Al₂O₃、Be计8项指标，主要为高场强元素组合(除K₂O、Be、Al₂O₃外)。F3主要代表Zn、Cd、Pb元素组合，F4主要代表Li、F、B、Be元素组合，F5主要代表CaO、MgO元素组合，F6主要代表Sn、Bi元素组合，F7主要代表Hg、Sb、As元素组合，F8主要代表Ba、Sr、P元素组合，F9代表Na₂O，F10代表Au。

在F1中因SiO₂具有负相关性，因此在回归计算F1因子得分时只采用排在SiO₂前的Fe₂O₃、Ti、V、Co、Cr、Ni计六项指标，回归系数采用旋转因子系数，即

F1=0.95Fe₂O₃+0.93V+0.93Ti+0.92Co+0.81Cr+0.80Ni。

该综合指标代表相容元素组合，在花岗岩类侵入体中应该具有较低数值。

在F2中8项指标均为正相关，因此在回归计算F2因子得分时全部采用这些指标，回归系数采用旋转因子系数，即

F2=0.87Th+0.86Y+0.77U+0.77Zr+0.67La+0.61K₂O+0.58Al₂O₃+0.53Be。

该综合指标主要代表高场强元素组合(属于不相容元素组合)，在花岗岩类侵入体中应该具有较高数值。

为了突出强化差异，选取在花岗岩类侵入体中具有较高值的F2与具有较低值的F1进行比值，获得F2/F1综合指标。利用GeoExpl软件绘制F2/F1累频19级等值线图(图 4a)，在此基础上选择累频85%、92%、98%值分别作为其异 常外、中、内带起始值绘制F2/F1的异常图(图 4b)。图 4c为F2/F1的异常等值线与花岗岩类侵入体的空间套合关 系。

图 4

Fig. 4

图 4 F2/F1等值区及其异常与花岗岩类侵入体的关系 (a)为F2/F1等值区图，(b)为F2/F1异常等值线图，(c)为侵入体图叠加F2/F1异常等值线图，其中粉红色区域为花岗岩类侵入体分布区，F2/F1异常外带为黑色线，异常中带为浅红色线，异常内带为棕红色线 Fig. 4 Geochemical map and geochemical anomaly map of F2/F1 which overlays on the granitic intrusion map in the western of Yunnan Province，China

从图 4可以看出，在花岗岩类侵入体分布区F2/F1显著富集。F2/F1的异常区与该区花岗岩类侵入体的空间分布区吻合较好。这表明基于因子分析结果利用代表高场强元素组合与代表相容元素组合的公因子比值所构建的模型可以较好的反映花岗岩类侵入体的空间分布范围。 4 地球化学推断模型应用

利用在云南西部地区构建的区域地球化学推断花岗岩类侵入体模型所确定的参数，即

F1=0.95Fe₂O₃+0.93V+0.93Ti+0.92Co+0.81Cr+0.80Ni

F2=0.87Th+0.86Y+0.77U+0.77Zr+0.67La+0.61K₂O+0.58Al₂O₃+0.53Be

将其推广到全国范围进行推演。

基于全国区域化探上述F1与F2所涉及的14项元素含量数据计算F2/F1值。然后绘制全国F2/F1累频19级等值线地球化学图及其外、中、内三级异常图。在全国F2/F1等值区图及其异常图中直接裁剪出华南中部地区F2/F1的等值区图(图 5a)和异常图(图 5b)。将图 5b中异常等值线叠加到华南中部侵入体图中获得图 6。

图 5

Fig. 5

图 5 华南中部F2/F1等值区图及其异常图 (a)为F2/F1等值区图，(b)为F2/F1异常等值线图，图中地名与图 6中相同 Fig. 5 Geochemical map and geochemical anomaly map in the middle part of South China

图 6

Fig. 6

图 6 华南中部F2/F1异常等值线与侵入岩体空间关系 该区从全国图件中直接裁剪获得，图中粉红色区域为花岗岩类侵入体分布区，F2/F1异常外带为黑色线，异常中带为浅红色线，异常内带为棕红色线 Fig. 6 Geochemical anomaly map of F2/F1 which overlays on the granitic intrusion map in the middle part of South China

从图 6可以看出，华南中部地区F2/F1异常区基本位于该区花岗岩类侵入体的空间分布区内，异常区形态和规模与侵入体的形态和规模吻合程度很高，尤其在韶关至桂林、贺州一带，花岗岩类侵入体边界与F2/F1勾绘的异常区边界几乎一致。此外还发现在局部地区F2/F1异常区范围明显大于实测花岗岩类侵入体的分布范围，如清远-惠州一带和茂名-阳江一带，据此推断这两处可能存在尚未证实的花岗岩类侵入体。

上述基于云南西部模型区建立的区域地球化学推断花岗岩类侵入体模型在华南中部地区得到有效检验。这一结果表明，在应用区域地球化学资料推断花岗岩类侵入体(或其它地质体)时只要方法技术合理、参数选择合适推断出的地质体有较强的可信性。这一方法技术在基础地质研究、矿产资源勘查及环境、生态等领域研究中将具有潜在应用价值。 5 结论

全国矿产资源潜力评价化探工作组编制了全国二级地球化学省划分方案，将全国共划分出5个地球化学域和25个地球化学省。在云南西部地区花岗岩类侵入体具有显著富K₂O、Th、U、La、Y、Zr而贫Co、Ni、V、Cr、Ti、Fe₂O₃的特征。

选择云南西部地区利用因子分析数据处理技术构建了推断花岗岩类侵入体的区域地球化学推断模型。基于构建的推断模型在全国范围内进行推演并选择华南中部地区进行验证和应用，结果发现模型指标异常区与实测花岗岩类侵入体在空间形态和规模上相吻合，这不仅证实了所建模型的可行性，而且该模型在基础地质研究方面具有潜在应用价值。

致谢

研究工作得到了全国矿产资源潜力评价工作组，尤其是全国化探组及大区级、省级化探组成员的指导与帮助，在此致以诚挚的感谢！