自热液矿床原生晕典型分带序列被提出以来(Beus et al.,1997)，该方面的研究愈加成熟、完善，并形成了一套完整的原生晕找矿方法，并在“攻深找盲”“探边摸底”等方面取得了较好的找矿效果(邵跃，1997；刘崇民，2006；刘崇民等，2010；李惠等，2010a;李惠等，2010b)。因此，在热液矿床中选用原生晕找矿方法，可为下一步找矿工作提供一定的参考依据。

玛日当沟钨矿床大地构造位置位于西秦岭造山带与东昆仑造山带东部衔接部位之青海南山造山带(III1，图1)内，属于布喀达坂—青海南山成矿省同德—泽库印支期汞、砷、铜、铅、锌、金(锑、钨、铋、锡)成矿带，是青海省重要的金及多金属成矿带之一。该矿床受北北西向断裂控制(图2)，目前已达大型规模，且在该矿床的南西部(上盘)发育加吾金矿床。前人主要对该矿床南西部的加吾金矿床地质特征、成矿条件等进行了研究(李学虎等，2010；李学虎等，2013，袁士松等，2011；畅捷，2015；范增平等，2015；郝骐，2015；王斌等，2015；张嘉升，2015)，虽二者相距较近，且同属一个构造带及成矿带内，但玛日当沟钨矿床的研究程度相对较低，仅开展了地质特征及成矿特征等方面的研究。虽然将原生晕方法应用于钨矿的找矿研究在国内外鲜有报道，但热液成因的钨矿床在形成过程中，其成矿元素及伴生元素在运移时的赋存形式和运移方式亦会存在差异，且会发生先后沉淀现象，从而在空间上也会呈现出异常分带现象。而通过对异常分带的解读，亦能发现一些深部找矿的线索。为此，该文旨在通过对该钨矿床成矿元素的原生晕分带序列进行定量研究，并建立深部矿体预测模型，为研究区攻深找盲提供方法依据。

1 矿床地质特征

1.1 矿区地质特征

研究区内出露地层(图2)较简单，主要为下—中三叠统隆务河组第二岩段(T1-2l2)砂质板岩夹薄层碳质板岩及砂岩，下—中三叠统隆务河组第三岩段(T1-2l3)炭质板岩夹变质砂岩，新近系咸水河组上段(Nx2)紫红色夹青灰色砂质粘土岩、含砾砂岩，以及下更新统共和组(Qp1g)河—湖沉积、洪积、风沙堆积、冰水堆积的砾石、泥土层。矿体多产于北东向断层下盘的下一中三叠统隆务河组第二岩段(T1-2l2)砂岩破碎带中，二者产状近一致。隆务河组发生了浅变质作用，主要变质岩性为板岩及变砂岩，地层中褶皱变形较强，以不同形态的褶皱构造、断裂构造以及板岩发育透入性板理和砂岩发育非透入性间隔劈理为特征，总体具有强变形弱变质的特征，属于较强应力条件下的低绿片岩相区域变质作用的产物。

  

图1 青海省东部构造分区略图(据潘桂堂等，2002)

 

I1—北祁连造山带；I2—中祁连造山带；I3—拉脊山造山带；II1—柴达木北缘过渡带微地块；II2—柴达木北缘地块：II3—柴北缘造山带；II4—柴达木中间地块；II5—祁漫塔格微造山带；II6—东昆仑北坡过渡带；II7—柴达木南缘微地块；III1—青海南山造山带；III2—西倾山中间地块；III3—阿尼玛卿山微造山带；III4—北巴颜喀拉造山带；III5—南巴颜喀拉造山带；III6—通天河造山带

区内仅南部发育少量花岗斑岩脉、闪长玢岩脉。由于本区经历的构造活动持续时间长，不同期次、强度及类型的构造作用发生了复杂的复合现象，导致区内构造较为复杂。其中，断裂构造较发育，褶皱次之；断裂构造以北北西向压扭性断裂为主，近东西向断裂相对较少。矿体及岩脉的空间分布严格受断裂构造控制，其中北北西向断裂与地层、斑岩脉的走向相同，产状基本一致，局部斜切地层；断裂两侧岩石破碎，产状紊乱，发育牵引褶曲，且多被花岗斑岩脉、石英脉、矿(化)脉所充填，是该区的主要控矿、容矿断裂。近东西向断裂中多充填有石英脉、闪长玢岩脉及少量花岗斑岩脉。区内褶皱构造以北北西—南南东向及近东西向为主，局部被北北西向压扭性断裂截切。

1.2 矿体特征

玛日当沟钨矿床目前由22个槽探工程及9个钻探工程控制。地表控制钨矿化带长度847 m，厚1.9～88.95 m，总体走向198°，倾向北东，倾角50°～70°，主要赋存在北北西向断层下盘的隆务河组第二岩段(T1-2l2)砂岩破碎带中。钨矿化主要为白钨矿化，黑钨矿化较少。其中，白钨矿主要有3种赋存状态：(1) 赋存在顺层或沿节理充填的石英细脉中(图3(a))，这种赋存形式的矿化是该区最主要的矿石类型；(2) 砂岩的裂隙面中(图3(b))；(3) 呈细粒星点状分布在砂岩中，较少见。

根据该区矿体中钨矿的总体赋存状态确定该区矿石类型为石英脉型，钻孔中主要表现为网脉及单脉状含钨石英脉充填于砂岩节理、裂隙中，厚度多为0.2～15 cm，钨矿品位和石英网脉的发育程度呈正比。围岩主要为碎裂长石石英砂岩、含炭粉砂质板岩、粉砂岩，长石石英砂岩中的裂隙发育程度及含矿性要好于含炭粉砂质板岩及粉砂岩。围岩蚀变主要为硅化、高岭土化、碳酸盐化、绢云母化及绿泥石化；金属矿化主要为毒砂、白钨矿、黄铁矿，少量黑钨矿、锡石、黄铜矿、方铅矿化，氧化带见有褐铁矿化及少量铜蓝、孔雀石化。白钨矿主要呈片状、块状、星点状产出，粒度多为0.2～5 mm。矿体的WO3品位为0.079%～0.239%，平均0.139%。

1.3 成矿阶段划分

矿石中主要矿石矿物为毒砂、白钨矿、黄铁矿、黑钨矿、锡石、黄铜矿、闪锌矿等，脉石矿物主要有石英、方解石等。矿石构造主要为脉状、网脉状、块状构造，矿石的结构主要有自形晶、半自形晶粒状结构、他形粒状结构、碎裂结构、交代残余结构和环状胶状结构等。金属矿物毒砂、白钨矿、黄铁矿、黑钨矿等多呈团块(斑)状、针柱状、稀疏浸染状或星点浸染状与脉状、网脉状石英伴生，属成矿热液沿构造裂隙充填而成。

矿床(体)的形成会因其多期、多阶段的成矿过程而发生原生晕叠加的现象，各元素含量单一升降的规律就可能会发生变化。而基于元素分带指数值，利用前缘晕组合元素的分带指数累乘值与尾晕组合元素的分带指数累乘值之比，可建立有效的深部矿体预测模型(Beus and Grigorian,1997)。

经野外观察和室内研究，结合矿石结构构造、蚀变特征及矿物共生组合特征，将成矿阶段划分为3个阶段：早期石英—黑钨矿阶段(早期成矿阶段：自形石英—黑钨矿—锡石等)、白钨矿—多金属硫化物阶段(主成矿阶段：白钨矿—黄铜矿—自形毒砂—自形黄铁矿等)、晚期硫化物阶段(成矿晚期阶段：他形黄铁矿—闪锌矿—方铅矿等)。

  

图2 玛日当沟钨矿床地质图

 

1—下更新统共和组河—湖沉积层；2—新近系咸水河组上段；3—下-中三叠统隆务河组第三岩段；4—下-中三叠统隆务河组第二岩段；5—花岗斑岩脉；6—闪长玢岩脉；7—石英脉；8—实测/推测的性质不明断层；9—实测逆断层；10—断层破碎 带；11—地层倾向及倾角；12—倒转地层倾向及倾角；13—金矿(化)体；14—铅矿体；15—钨矿体；16—钻孔及编号

  

图3 玛日当沟钨矿典型矿石照片

 

(a)-毒砂-黑钨矿-白钨石英脉；(b)-砂岩节理中发育的星点状钨矿化；(c)-白钨矿-毒砂矿化(10×10(-))；(d)-白钨矿不完全 交代黑钨矿(10×2.5(+))；Apy-毒砂；Sh-白钨矿；Wol-黑钨矿

2 原生晕分带特征

1）危害症状。苹小食心虫多在果实胴部为害，幼虫蛀果后在果皮下浅层果肉，一般不深入果心，被害部位果皮变褐、干裂，形成直径1 cm左右近圆形干疤，稍凹陷，干疤上有少量虫粪。脱果孔在蛀入孔或果实虫孔疤边缘。

总体上，该钨矿床原生晕轴向分带序列与多金属矿床标准轴向分带序列的元素分布相似，且较能客观地反映出各元素在深部的相对浓集位置，故对深部矿体资源潜力评价具有重要的指导作用。

2.1 成晕元素分带

从表2中可以看出，前4个主因子的累计方差贡献达到70.70%，可认为这4个因子已经包含了原始变量的大部分信息。初始因子载荷中各因子载荷的分异度不高，特征表现不明显。经正交旋转后，因子载荷的分异变得相对清晰：F1为Ag-Pb-Zn因子组合；F2为Cu-W组合；F3为Au-Sb-Bi因子组合；F4为Mo因子。

依据上述参数，绘制了各元素异常的剖面图(图4)。从0号勘探线剖面图上可以看出：(1)Au、Ag、Mo、W元素异常形态相似，且与矿体形态相似，反映了矿体向东倾伏的特点，在矿体的中下部具有较强的异常，内、中、外带发育(W元素含量因测试限值不超过200，其内带不发育)，浓度分带清晰完整，Au、Mo、W异常规模向下变大，富集趋势明显；(2)Sb元素异常在矿体的中部、中上部相对发育，表现为矿体中上部的中内带异常；(3)Bi异常主要分布在矿体下部，表现为矿体下部的中外带异常，且向下具有明显富集增大的趋势；(4)Cu、Pb、Zn元素异常在矿体的中(上)部具有相对较强的异常，表现为矿体中上部的中内带异常。

总体上，Ag、Cu、Pb、Zn、Bi异常向上有尖灭闭合的趋势，Au、Sb、W、Mo异常向上封闭趋势较弱。对比元素异常及矿体的空间分布特征，初步判断矿体的前缘指示元素为Sb，近矿指示元素为Cu、Mo、Au、Pb、Ag、Zn，尾晕元素为W、Bi。

为了开展地球化学原生晕分带研究，分别对I号脉0线上的3个钻孔(ZK002、ZK004、ZK008)进行采样(ZK012刚施工)，共采集原生晕样品510件；采样方法为除第四系外，基岩部分每5m或小于5m(视地质界线及矿化蚀变情况而定)组合采集1件样品，连续采样。样品测试由武警黄金第六支队实验室承担(乙级资质)，定量分析Au、Ag、Cu、Pb、Zn、Sb、Bi、W、Mo等元素，并以元素含量作为进一步定量研究原生晕分带的基础。元素检出限及分析方法分别为：Au，0.3×10-9，原子吸收分光光度计(Z-2000/HLS-039)；Ag，0.03×10-6，平面光栅摄谱仪(WP1/HLS-006)；Cu，1.5×10-6，X荧光光谱仪(AxiosMAX型-PW4400/HLS-071)；Pb，5×10-6，X荧光光谱仪(AxiosMAX型-PW4400/HLS-071)；Zn，10×10-6，X荧光光谱仪(AxiosMAX型-PW4400/ HLS-071)；Sb，0.2×10-6，原子荧光分光光度计法(AFS-9800/ HLS-072)；Bi，0.1×10-6原子荧光分光光度计法(AFS-9800/ HLS-072)；W、Mo，0.5×10-6，电感耦合等离子体质谱法(ICP-MS)。

2.2 矿体原生晕的元素组合特征

选取Au、Ag、Cu、Pb、As、Sb、Mo、W等元素的含量进行处理。在元素含量通过KMO和BartIett’s球度检验后对其进行因子分析，确定元素的共生组合特征(表2)。

线下：由于线下实地调研的原因，原本是单选题的此题被很多参与者当成了多选题进行回答，现就以多选题形式对此题进行分析。有52.63%的人认为对于南京城墙应该以保护为主开发为辅，21.05%的人认为应与周围环境相协调，另外还有21.05%的人认为应该注重宣传。说明绝大多数人对于历史遗迹的态度还是以保护为主。

 

表1 玛日当沟钨矿床原生晕元素分带参数

  

含量分级Au/10-9ωB /10-6AgCuPbZnSbBiWMo外带12.00～24.000.75～1.50120～24015～3091～1820.75～1.507.5～1570～1402.25～4.5中带24.00～48.001.50～3.00240～48030～60182～3641.50～3.0015～30140～2804.5～9.0内带>48.00>3.00>480>60>364>3.00>30>280>9.0

在研究成晕元素分带之前，要合理地确定原生晕的异常下限。由于此次研究的样品均采自矿区内，背景值普遍较高。该文首先利用SPSS软件对原始数据进行概率分布分析，对于近似满足对数正态分布的元素(Au、Ag、Cu、Pb、Sb、Bi、Mo)，将其对数均值的指数作为异常下限；对于取对数后满足双峰分布的元素(W)进行概率拟合计算，求出两条概率分布曲线的参数，将第一条曲线的均值的指数作为异常背景值，两条曲线交点的指数作为异常下限(蒋敬业等，2006)。将原生晕划分为3个浓度带(外带、中带、内带)，以1倍、2倍、4倍的异常下限圈出原生晕的外带(Ca)、中带(2Ca)、内带(4Ca)(表1)。

 

表2 各元素旋转因子载荷表

  

元素旋转因子载荷F1F2F3F4Au-0.030.340.740.18Ag0.810.320.230.07Cu0.240.760.17-0.35Pb0.89-0.090.050.15Zn0.500.34-0.05-0.53Sb0.19-0.260.73-0.30Bi0.210.240.620.48W0.020.790.040.18Mo0.170.030.030.67累计方差贡献/%20.6639.3156.5470.70

以特征值λ>1为选取因子个数的标准。

选取(Cu×Sb×Mo)D/(Ag×W×Bi)D作为深部矿体预测模型的评价指标。该指标自矿体的上部至下部(由浅至深)明显降低(图5)：矿体上部(3540 m标高)4.26→矿体中部(3350 m标高)0.71→矿体下部(3100 m标高)0.39，表明该指标随着矿体深度的增加而发生明显的减小，可作为深部预测的有效指标。

1977年，文艺复兴史学家Joan Kelly-Gadol(1928—1982)提出了一个问题：“妇女有一个文艺复兴吗？”经过对历史资料的研究，她得出结论，在文艺复兴时期的欧洲，妇女仍旧只是男人的附庸，她们的处境并没有得到真正意义上的改善[1]。

2.3 矿床原生晕轴向分带序列

原生晕分带可反映含矿溶液的运移方向，其分带序列对判定矿体剥蚀程度及“攻深找盲”具有重要的意义。该文采用应用较为广泛的分带指数方法计算分带序列。该方法以元素异常的线金属量为基础，对线金属量进行标准化处理后，将其总和起来，可减小原始数据的偶然误差；该方法较为实用，己在诸多矿区中取得了较好的找矿效果(Chen et al.,1998;李强等，2005；晁会霞等，2006；陈伟军等，2007；杨怀辉等，2007；陈永清等，2010)。

如表3所示，基于分带指数法所获得的分带序列为Cu→Sb→Mo→Au→Pb→Zn→Ag→W→Bi。0号勘探线上Cu、Sb、Mo位于分带序列的上部，Au、Pb、Zn、Ag基本位于中(上)部，W、Bi位于序列尾部，除上部序列中有高温元素Cu、Mo外，总体上与Beus and Grigorian(1997)总结的多金属矿床标准轴向分带序列[(Sb,As1,Hg)-Cd-Ag-Pb-Zn- Cu-Bi-Mo-Sn-As2-W]基本一致。不同的是0号勘探线上Cu、Mo高温元素与前缘晕元素Sb共存，存在反分带现象，其他为正常序列。

水利水电工程项目建设需要考量项目建设成本问题，由于水利水电工程项目普遍存在建设周期较长、建筑成本投入较大等特点，建设成本问题也是工程项目建设管理需要关注的重点问题。项目建设过程中，加强施工技术管理，能够保证各项施工技术有效落实到工程项目中，从而实现建设资源的合理配置和使用，提升项目管理的有效性和科学性，避免出现施工技术落实不到位影响工程质量，同时也能够有效避免施工成本的过度投入对项目管理效能造成负面影响。强化施工技术管理，能够实现对成本流向的合理监控，保证工程项目的合理投入，同时施工技术的积极创新，也能够从缩短施工周期、减少施工资源投入等方面实现有效的成本管理。

3 深部矿体预测模型

随着高职教育改革的不断深入，各种新的教学模式不断研究和推广，项目教学、任务教学、案例教学的广泛应用使跨境电商英语写作课程与传统写作课程有了很大区别，提高了课堂气氛的活跃性，帮助学生克服了因商务知识不足而造成的写作障碍，提高了学习兴趣。本文主要探讨跨境电商英语案例分析库的建设思路及在跨境电商英语写作教学中的应用。

(1)对于“家长放任型”的学生，教师需要对他们施以情感教育来抚平他们心灵的创伤，为他们创造温暖的大家庭氛围。这就要求我们教师从爱出发从严要求。

  

图4 I号脉0号勘探线成矿成晕元素浓度分带

从因子分析的元素组合来看，F1为Ag-Pb-Zn因子组合，以中温元素组合为主，累计方差贡献为20.66%，表明该区可能存在Ag-Pb-Zn等矿化(该矿床南西部的加吾金矿床中发育银、铅、锌矿(化)体)，且该因子组合与晚期硫化物阶段中的矿物组合(他形黄铁矿—闪锌矿—方铅矿等)相对应。F2以高温元素组合为主，与白钨矿—多金属硫化物阶段的矿物组合(白钨矿—黄铜矿—自形毒砂等)相对应；F3为Au-Sb-Bi因子组合，低、中、高温元素均有，鉴于三者的地球化学特征差异，认为这是由于不同期次或不同阶段的热液活动叠加的结果；F4为Mo因子，属高温元素。Bi、W、Mo、Cu高温元素分别位于不同的因子中，结合该区的成矿地质特征以及对F3因子的分析，认为该区成矿具有多期多阶段的特征。

基于最小二乘法对标高(y)及分带指数累乘比进行直线拟合，得出拟合函数为：y = 165.89*x + 3321.2，x= ln[(Cu×Sb×Mo)D/(Ag×W×Bi)D]，R2 = 0.8761。由此可以看出：异常中心(标高)与分带指数累成比之间具有较好的线性关系。因此，该指标可用于预测深部矿体。若该评价指标随着深度发生线性降低之后，在一定深度又突然升高，则推测深部可能出现了盲矿体(Chen et al.,1998;陈永清等，2010；Chen et al.,2008)，因为这种突然升高，很可能是深部矿体或盲矿体的前缘晕叠加到上部矿体的尾晕所造成的(陈永清等，2010；Li et al.,1995;Li et al.,199;李惠等，1999；孙华山等，2008)。

 

表3 玛日当沟钨矿床原生晕轴向分带序列

  

元素线金属量(初始)/10-6·mZK002ZK004ZK008标准化系数线金属量(标准化后)/10-6·mZK002ZK004ZK008分带指数ZK002ZK004ZK008Au8.588 324.771 1411.596 3510008 588.324 771.1411 596.350.0590.0140.024Ag150.41421.50587.3910015 041.1542 150.2758 738.630.1030.1270.121Cu32 070.9265 120.7091 490.51132 070.9265 120.7091 490.510.2190.1970.189Pb2010.595492.487356.001020 105.9254 924.7973 560.000.1370.1660.152Zn14 608.2447386.6655631.90114 608.2447 386.6655 631.900.1000.1430.115Sb261.18321.69416.3210026 117.8232 168.5541 632.080.1780.0970.086Bi1 013.704 153.217 311.921010 137.0441 532.0973 119.180.0690.1260.151W11 158.0426 851.0554 243.95111 158.0426 851.0554 243.950.0760.0810.112Mo868.301 590.442 368.99108 682.9815 904.3723 689.860.0590.0480.049ΣM146 510.42330 809.62483 702.45

  

图5 玛日当沟钨矿床深部矿体预测模型

4 结论

通过对玛日当沟钨矿床地质特征及原生晕轴向分带特征的研究，得出以下几点认识：

(1) 通过因子分析，得出4个因子组合：Ag-Pb-Zn、Cu-W、Au-Sb-Bi、Mo，元素组合特征与该区矿化组合特征相对应。

(2) 基于原生晕元素异常特征及原生晕分带序列分析结果，确定了该矿床原生晕轴向分带序列为Cu→Sb→Mo→Au→Pb→Zn→Ag→W→Bi。除上部序列中有高温元素Cu、Mo外，总体上与Beus and Grigorian(1997)总结的多金属矿床标准轴向分带序列基本一致。

(3) 建立了深部矿体预测模型，其预测评价指标(Cu×Sb×Mo)D/(Ag×W×Bi)D的自然对数与深度之间具有较好的线性关系，该指标的比值随深度的增加而发生线性降低，可作为深部矿体预测的有效指标。

致谢 研究工作得到了牟长贤工程师及陈俊瑾工程师的热情帮助，武警黄金第六支队实验室完成了研究样品的测试分析工作，审稿老师对文章进行了细心修改，并提出了建设性意见和建议，在此深表感谢。

全内脏反位(Situs Inversus Viscerum)又称“镜面人”或“镜像人”，是指人的心脏、 肝脏、脾脏、胆囊、肾脏、胃、肠等器官的解剖位置与正常人完全相反。此种病例极其罕见，最近本院脊柱外科收治1例“镜面人”合并腰椎间盘突出症病例，采用经Quadrant通道系统椎间孔腰椎椎间融合术(Mis-TLIF)治疗，现将诊治过程报道如下。

参考文献

Beus A A， Grigorian S V. 1997. Geochemical exploration methods for mineral deposits [J]. Applied Publishing Ltd., Wilmette Illinois, U.S.A., 287.

Chen Yongqing, Huang Jingning, Liang Zhen. 2008. Geochemical characteristics and zonation of primary halos of Pulang porphyry copper deposit, northwestern Yunnan Province, southwestern China [J]. Journal of China University of Geosciences, 19(4):371-377.

Chen Yongqing, Zhao Pengda. 1998.Zonation in primary halos and geochemical prospecting pattern for the Guilaizhuang gold deposit, Eastern China [J]. Natural Resources Research, 7(1): 37-44.

Li Hui, Wang Zhinong, Li Fuguo. 1995.Ideal models of superimposed primary halos in hydrothermal gold deposits [J]. Journal of Geochemical Exploration, 55(1-3): 329-336.

Li Yinggui, Cheng Hangxin, Yu Xuedong, Xu Waisheng. 1995. Geochemical exploration for concealed nickel-copper deposits [J]. Journal of Geochemical Exploration, 55(1-3): 309-320.

畅捷. 2015.青海省加吾金矿构造地球化学特征及深部成矿潜力分析[D]. 长安大学, 1-125.

晁会霞, 杨兴科, 王磊, 梁广林, 孙继东, 马立成. 2006.新疆梧南金矿床原生晕特征与深部预测[J]. 地质与勘探, 42(1): 72-76.

陈伟军, 刘红涛. 2007. 内蒙红花沟金矿床地质特征与原生晕深部预测[J]. 金属矿山, (4): 44-48.

陈永清, 韩学林, 赵红娟, 程志中, 唐 宇, 陈 武. 2010.内蒙花敖包特Pb-Zn-Ag多金属矿床原生晕分带特征与深部矿体预测模型[J]. 地球科学--中国地质大学学报, 36(2): 236-246.

范增平, 王根祥, 赵向红, 段留安. 2015.青海省加吾金钨矿地质特征及钨成矿评价[J]. 矿物学报, 35(S1): 475.

郝骐. 2015.新构造运动对加吾金矿床赋存状态的影响研究[D]. 长安大学, 1-82.

蒋敬业, 程建萍, 祁士华, 向 武. 2006.应用地球化学[M]. 武汉: 中国地质大学出版社, 1-180.

李惠, 禹斌, 李德亮, 马久菊, 许丽云, 魏江, 赵佳祥, 张 强, 司淑云, 杨亚娟, 贺容华. 2010a.化探深部预测新方法综述[J]. 矿产勘查, 1(2): 156-160.

李惠, 张国义, 禹斌, 李德亮. 2010b.构造叠加晕找盲矿法及其在矿山深部找矿效果[J]. 地学前缘, 17(1): 287-293.

李惠, 张文华, 刘宝林, 常凤池. 1999.金矿床轴向地球化学参数叠加结构的理想模式及其应用准则[J]. 地质与勘探, 35(6): 40-43.

李强, 孙继东, 杨兴科. 2005.新疆石英滩金矿床原生晕特征与隐伏矿预测[J]. 地质与勘探, 41(4): 66-72.

李学虎, 高仁品, 王雁鹤. 2010a.青海省同德县加吾金多金属矿区微量元素在地质体中的分布特征[J].矿产与地质, 24(6): 513-519.

李学虎, 杨涛, 程真. 2013.青海加吾矿区地层控矿作用及赋矿特征[J]. 矿产勘查, 4(2): 146-153.

李学虎, 张东林, 高仁品. 2010b.青海加吾金矿地质特征及找矿标志[J].矿产勘查, 1(4): 354-359.

刘崇民, 马生明, 胡树起. 2010.金属矿床原生晕勘查指标[J]. 物探与化探, 34(6): 765-771.

刘崇民. 2006.金属矿床原生晕研究进展[J]. 地质学报, 80(10): 1529-1538.

潘桂堂, 李兴振, 王立全, 丁 俊, 陈智粱. 2002.青藏高原及邻区大地构造单元初步划分[J]. 地质通报, 21(11): 701-707.

邵跃. 1997.热液矿床岩石测量(原生晕法)找矿[M]. 北京: 地质出版社, 1-143.

孙华山, 孙玲, 曹新志, 王超, 谭俊, 刘福江, 杨开春. 2008.胶西北上庄金矿床原生晕轴(垂)向分带特征及深部矿体预测的勘查地球化学标志[J]. 矿床地质, 27(1):64-70.

王斌, 张栋, 路英川, 孙昊. 2015.青海加吾金矿床花岗斑岩脉锆石U-Pb年龄、Hf同位素特征及其地质意义[J]. 矿物岩石, 35(3): 52-60.

杨怀辉, 杨振京, 崔来运. 2007.豫西嵩县东湾金矿床地质及原生晕地球化学特征[J]. 四川地质学报, 27(2): 107-111.

袁士松, 葛良胜, 金宝义, 闫家盼, 喻万强. 2011.青海同德加吾金矿区印支期花岗岩类地球化学特征及其地质意义[J]. 矿物岩石地球化学通报, 30(增刊): 245.

张嘉升. 2015.青海加吾金矿地质地球化学特征研究[D]. 长安大学, 1-95.