江西万年银金矿床Ⅶ号脉原生晕特征与深部预测

2013-08-20吴建亮吕新彪杨永胜高荣臻王于键

地质与勘探 2013年5期

吴建亮，吕新彪，杨梧，杨永胜，高荣臻，刘洪，王于键

(1.中国地质大学，湖北武汉 430074;2.中国地质调查局成都地质调查中心，四川成都 610081)

赣东北乐华－德兴成矿带地处扬子板块江南地体东南缘，是我国东部矿集区的重要组成部分，其内发育有铜厂、朱砂红、富家邬铜矿床、银山铜铅锌金银多金属矿床、金山金矿床、万年银金矿床等一系列内生金属矿床。万年银金矿床地处乐德成矿带西南部，从20世纪90年代开始投入生产，目前开采到－235m中段。前人对该矿床已做过部分工作(徐磊明，1993a，1993b;董法先，1995)，但对该矿床原生晕空间分布特征及矿区深部成矿尚未有系统研究。因此，本文试图通过对万年银金矿床Ⅶ号主矿脉原生晕特征研究，建立元素地球化学模型，结合矿体定位规律分析，为矿区深部找矿提供有力依据。

1 矿床地质概况

万年银金矿床矿区范围内出露的地层较为简单(图1)，为中元古界双桥山群上亚群第四岩性组的一套浅变质的变质－沉积火山碎屑岩和第四系。赋矿围岩主要为凝灰质千枚岩和绢云母千枚岩。岩层产状稳定，常见保存完好的直线型及微波型薄层原生层理构造，部分地段可见变形层理。矿区构造以NE向断裂为主，受区域性构造应力场控制作用明显，形成在主干断裂(Ⅵ、Ⅶ号断裂)两侧不等间距分布的断续近平行断裂带(群)。断裂倾向SE，局部具双向倾斜特征，倾角陡立①。Ⅶ号断裂多期次继承性活动明显，成矿期多次张开闭合，构成矿区主要控矿断裂，控制了整个矿化过程。成矿期后断裂规模及影响较小。矿区范围内岩浆活动微弱，至今未见有岩体出现，仅发现少量中酸性岩脉，以石英闪长玢岩为主，少见石英安山玢岩及闪斜煌斑岩脉，均蚀变严重，局部可见硅化、碳酸盐化及矿化叠加，总体Au，Ag含量较低，为成矿前脉体①。矿区围岩蚀变种类简单，以硅化和碳酸盐化为主，少见绢云母化，总体沿矿脉呈带状线型分布。各蚀变带呈渐变过渡接触关系。硅化与成矿关系最为密切。矿区受区域构造作用表现明显，形成了较多碎裂岩化及糜棱岩化，在硅化叠加部位往往是矿化极佳场所。

矿体总体呈脉状，严格受构造破碎带控制，与围岩接触界线清晰。主矿脉为Ⅵ和Ⅶ号矿脉。Ⅶ号脉规模较大，总体走向北东，呈波状延展。倾向上具北西和南东双向倾斜的特征，倾角陡立，沿走向和倾向均膨胀收缩明显。Ⅵ号脉规模相对较小，矿脉受后期断裂影响较小。银矿物以银黝铜矿－黝锑银矿、深红银矿、辉银矿－螺状硫银矿等(徐磊明，1993)为主，粒径极小，主要以裂隙嵌布和粒间嵌布的方式分布在黄铁矿、毒砂、闪锌矿、方铅矿等金属矿物中①。各矿物表现出多世代及动力变形特征，反映成矿热液在构造应力环境下多期次活动特点。根据野外及镜下观察，矿区成矿期分为热液期及表生期，其中热液期又分为Ⅰ石英－碳酸盐阶段、Ⅱ石英－多金属硫化物阶段、Ⅲ细粒硫化物阶段和Ⅳ碳酸盐阶段，其中第Ⅱ、Ⅲ阶段为主要的成矿阶段。

2 矿床原生晕特征分析

图1 万年银金矿区1～64线区段地质略图(据资料②修改)Fig.1 Simplified geological map between 1st and 64th exploration lines of the Wannian silver-gold deposit area in Jiangxi Province(modified from data②)

原生晕方法又称为岩石地球化学方法，是利用在成岩成矿作用中形成、赋存于矿体或其他地质体中的地球化学异常而进行找矿的地球化学方法(蒋敬业等，2006)。从20世纪50年代起到现在，已发展成为地球化学找矿的最主要方法技术之一，尤其在寻找热液成因隐伏矿床方面具有明显直接的效果(刘崇明，2006;李惠等，2010)。

本次样品总体上是在Ⅶ号脉各中段沿脉方向进行系统布置和采集，共采集样品129件。样品由湖北省地矿局鄂东北实验室测试分析完成。选择与Au、Ag成矿有关的15种元素进行分析，包括B、Sn、Au、Ag、Bi、Hg、Mo、W、F、Cu、Cr、Mn、Ni、Pb、Zn。各元素采用的分析方法为:元素 Cu、Cr、Mn、Ni、Pb、Zn采用电感耦合等离子发射光谱法分析，B、Sn、Au采用发射光谱法分析，Bi、Hg采用原子荧光法分析，Mo、W采用极谱分析法分析，F采用电位分析法分析。测试结果精度满足分析要求。

2.1 成矿及伴生元素组合特征

2.1.1 相关分析

相关分析是利用元素间的相关系数来衡量各素间相关性和亲和性的一种数学方法。矿区各指示元素相关性如表1所示，以相关系数0.5为界，成矿组合 Ag、Au与 Mn、Hg相关性较好，表明 Mn、Hg可作为矿区深部预测的重要成矿指示元素，其分布具较大的指示意义。元素 Cu、Pb、Zn、W、Sn、Hg表现出较好的正相关性，可能是它们同为矿区紧密共生金属硫化物组合所致。元素B、F、Ni、Cr相关性较为密切，也是良好的伴生指示元素。

表1 万年银金矿床原生晕多金属元素相关系数表Table 1 Correlation coefficients of metal elements in the Wannian silver-gold deposit

2.1.2 聚类分析

R型聚类分析主要遵循“物以类聚”的原则来对变量归类。将万年银金矿原生晕各指示元素的聚类分析如图2所示。

图2 万年银金矿床原生晕指示元素R型聚类分析谱系图Fig.2 R type cluster hierarchical diagram of the indication elements in primary haloes of the Wannian silver-gold deposit

从上图可以看出，大约以相似系数0.29为界，可将元素分为 Ag、Mn、Pb、Zn、Cu、Hg、Sn、Au、W、Bi和 Ni、F、Cr、B、Mo 两个大群。当相似系数提高到0.51左右又可将第一个大群分为Ag、Mn亚群、Pb、Zn、Cu、Hg、Sn 亚群，Au 亚群，W 亚群以及 Bi亚群。Ag与Mn相关性较好，表明元素Mn可作为重要的成矿指示元素。Pb、Zn、Cu、Hg、Sn 具一定的相关性，该五种元素均具有较强的亲硫性，在矿物组合上表现为闪锌矿与方铅矿紧密共生，而Ⅰ世代黄铜矿、黝锡矿常以乳浊状出融于闪锌矿中。Au亚群与前两个亚群具一定相关性可能是因为Au以自然金形式主要赋存在毒砂和黄铁矿中，其次为方铅矿和闪锌矿中。Ni、F、Cr、B、Mo等元素为次要的伴生元素，虽然这些元素都是在成矿作用中发生了相对富集和运移，但它们沉淀时的物化条件和生成顺序是有差别的，R型聚类分群结果是这种差别的具体体现(张彦艳等，2006)。

2.1.3 因子分析

因子分析实际上是一种降维分析，是将存在复杂关系的较多变量依据某种内在联系生成几种新的变量，该变量提取了原来众多变量的主要信息，从而使标本或变量具有更明确的意义(Konstantinov et al，1995)。对于某一矿体来讲，所测元素含量实际上是在多个地质作用叠加条件下元素行为历史积累。R型因子分析的任务不仅要根据矿体中各种元素含量间的相互关系来识别一个过程，而且在多过程叠加的情况下来区分这些过程，并按不同的地质过程将元素总的历史成因进行分解(余金生，1985;赵鹏大等，1994)。对万年银金矿床因子分析结果中提取前3个主因子，累计百分数为76.2%，可以认为代表了原始变量中绝大部分信息。为使主因子地质意义更加明确，对因子载荷矩阵进行了方差极大旋转，其结果如表2所示。

从各指示元素的因子载荷可以发现:F1因子的元素组合为 Pb、Sn、Hg、Zn、Cu，成矿元素 Ag 有一定的贡献，而Au无载荷。因此，F1因子可能代表的是与贱金属硫化物相关的中低温热液银矿化，对应矿石中方铅矿、闪锌矿、黄铜矿、黝锡矿以及与方铅矿紧密共生的辉银矿、银黝铜矿等金属矿物组合，应为成矿期中的石英－多金属成矿阶段。F2因子的元素组合 Ag、Au、Mn，为成矿元素组合，Ag、Au 在成因上与元素Mn联系紧密，元素Hg、Sn也有一定载荷，这代表的是后期银矿化的叠加和金矿化的形成，应为主成矿阶段即细粒硫化物阶段，在矿石中表现为Ag、Au元素以含银黝铜矿、自然金等矿物的极细硫化物脉穿插交代毒砂、黄铁矿等早期硫化物裂隙，反映出该矿床具两期银矿化，一期金矿化的成矿特点，也是导致在元素相关性矩阵及聚类谱系图中Ag与Au仅具有一定相关性的原因，同时暗示出该矿床成矿物质具有多方面来源的可能性。F3因子的元素组合为 F、B、Cr、Ni，可能是代表围岩的元素组合。

表2 万年银金矿床R型分析因子旋转因子载荷矩阵Table 2 Factor loading matrix of R type factor analysis for the Wannian silver-gold deposit

2.2 矿床地球化学分带的空间分布特征

确定元素异常分带首先以确定异常下限为基础，合理确定原生晕异常下限是原生晕研究中的一个重要方面。在对原始数据进行统计分析的基础上，采取迭代法(平均值加上三倍标准离差)剔除高值，最终确定各元素背景值及异常下限，其具体结果见表3。

根据异常下限绘制的各元素原生晕剖面等值线图(图3)基本反映了万年银金矿床成矿及伴生元素的剖面发育特征。可以看出，各元素高值带普遍向北东侧伏，与矿体侧伏方向一致。元素Ag、Mo、Cr异常主要分布在矿脉东段，西段较低;Au、Mn元素异常在东、西两段上部均有分布;Hg、Pb、Zn、Cu、Sn、W、Bi元素异常在中段深部较发育呈现出中间高、两端低的分布特征;Ni元素异常主要在中段中上部发育，F、B元素异常在西段浅部和东段深部异常均有出现，根据指示元素异常的分布特征可初步划出三个原生晕带:

Ⅰ号原生晕带:位于矿体6号线以西地段，主要为近矿晕元素(Ag、Pb、Cu)，分布在矿体上部，向上未封闭。前缘晕元素异常分布相对较少，仅外带与中带发育，内带不发育。尾晕元素异常特征与前缘晕相似，但主要分布在矿体中下部。

Ⅱ号原生晕带:在矿区中段20～56线地区元素异常则表现为浅部主要分布为尾晕元素(Mn、Ni)和成矿元素(Au)，而在深部46线附近前缘晕元素(Hg)异常内、中、外带具有发育，面积较大且异常向下未封闭，在空间上与近矿元素(Pb、Zn)与尾晕元素(Sn、Mn、W)同位共存，有可能既是深部有盲矿体存在，也是两个矿化带或大矿体间有小矿体存在的反映(孙华山等，2004)。异常形态具北东向深部延伸趋势，显示出矿区深部46线以东较好的成矿远景。

表3 万年银金矿床原生晕多金属元素异常下限表Table 3 Anomaly threshold of metal elements in primary haloes for the Wannian Silver-Gold deposit

图3 万年银金矿床Ⅶ号脉指示元素品位纵向剖面等值线图Fig.3 Isolines of the primary haloes for indication elements in the Wannian Silver-gold deposit

表4 万年银金矿床原生晕元素轴向分带序列表(表示最大值)Table 4 Primary halos`axial zoning sequence in the Wannian sliver-gold deposit(* represents the maximum)

Ⅲ号原生晕带:位于矿体62线以东地段，近矿晕元素(Ag、Au、Bi)主要分布在浅部，异常较发育，具明显的分带特征。尾晕元素(Mo、Mn、Ni、Cr)在深部较发育。

原生晕的水平分带是指示元素在现代水平方向上异常发育的强度、范围的规律性变化。水平分带又可分为纵向分带和横向分带。纵向分带即指顺矿体走向反应的元素分带，横向分带则是垂直于矿体轴面方向上的元素分带，主要是扩散作用造成的分带。在同一次成矿热液活动中是不可能形成两个元素浓集中心的，假如成矿作用是在多期次矿化作用中形成的，各期矿化相互叠加，导致上述分带规律的破坏，形成“多建造晕”和“叠加晕”(蒋敬业等，2006)。采用格里戈良分带指数法(蒋敬业等，2006)计算万年银金矿床Ⅶ号脉纵向分带序列为由南西往北东依次为 Sn、W、Zn、Pb、Cu、Ag、Au、Mn、Bi、Cr、B、F、Hg、Ni、Mo。元素走向上显示矿体南西段以尾晕及近矿元素组合为主，往北东逐渐过渡到以成矿及低温元素组合为主，到矿体东部又出现高温元素和成矿元素组合，是多期次矿化的典型表现。以矿体为中心，其横向分带分为:边缘晕B、F、Hg、Ni;中间晕 Cr、Bi、Mo、Cu;内部晕 Ag、Au、Pb、Zn、Mn、Sn、W。

2.3 原生晕轴向分带序列

利用矿床原生晕轴向分带的规律在评价异常剥蚀水平，指导盲矿体的寻找和勘探工作上都具有重要意义。目前，原生晕轴向分带序列研究中常用的方法有浓集指数法、分带性衬度系数法、分带指数法、重心法、比重指数法、概率值法、金属量梯度法等(邱德同，1989;朴寿成等，1996，1994;解庆林，1992;赵琦，1989;张定源，1989)。本文采用解庆林(1992)提出的浓集指数法进行计算。与其他方法相比，该方法的原始数据标准化的结果是唯一的，且数据集中，能消除地壳中各元素丰度差可能给计算带来的影响;其数据标准化为无量纲化处理(解庆林，1992)，同时吸收了格里戈良分带指数法考虑不同截面数据的作用(蒋敬业等，2006)，在确定分带序列时具有较好的分辨性。

由表4可知，万年银金矿床原生晕轴向分带序列表现为Sn－Au－Cu－Mn－B－Mo－Ni－F－Ag－Cr－Hg－Pb－Zn－W －Bi。可以看出，矿体上部富集的元素为 Sn、Au、Cu、Mn、B、Mo，以中高温元素组合为主，表明矿体上部遭受了一定程度的剥蚀;矿体中部富集的元素为 Ni、F、Ag、Cr，而矿体下部富集的元素为 Hg、Pb、Zn、W、Bi，两者总体上表现为以中低温元素组合或近矿晕元素组合为主，是成矿阶段热液活动产物的主要反映。李惠等(1999，2006)等研究认为金矿的每一次成矿作用形成的矿体(晕)都有自己的前缘晕、矿体晕和尾晕，属正常分带序列，不同成矿阶段形成的矿体(晕)在空间上的叠加会使分带序列计算受到一定的影响，甚至出现“反(向)分带”或“逆(向)分带”，多期多阶段叠加形成的“反分带”是预测深部盲矿的重要依据。与李惠等(1999，2006)统计的中国金矿床原生晕正常垂向分带B－I－As－Hg－F－Sb－Ba－Pb－Ag－Au－Zn－Cu－W－Bi－Mo－Mn－Ni－Cr－Co－V－Bi相比，万年银金矿床其前缘晕元素B、F、Hg，近矿晕元素 Cu、Pb、Zn、Ag、Au，尾晕元素 Sn、Mn、Ni、Cr、W、Bi分布混杂，呈现出明显的逆向分带特征，在矿体的尾部出现了前缘晕和尾晕叠加的现象，表明该矿区存在多次成矿热液活动，反映了另一阶段形成的矿体头晕的叠加，预示矿体往深部还有较大的延伸。

2.4 原生晕轴向地球化学参数

原生晕轴向地球化学参数变化特征是原生晕各元素轴向分布的另一种表达方式，是预测隐伏矿体的重要依据，地球化学参数可以用前缘晕指示元素浓度标准化/尾晕指示元素浓度标准化、前缘晕指示元素浓度标准化累加/尾晕指示元素浓度标准化累加来表示(李惠等，2006)。矿床(体)在不同阶段形成矿体(晕)的叠加会破坏单一阶段成矿体(晕)参数在轴向上单一的升或降的规律，使其发生“转折”。大型、特大型金矿地化参数在轴向上往往有多个“转折”(李惠等，1998)。因此选择前缘特征元素和尾晕特征元素以及由此建立的地球化学参数可以对矿体(脉)的深部含矿性作出有效评价。

据此，得到江西万年银金矿床Ⅱ号原生晕带轴向地球化学参数 w(Hg)/w(Mn)、［w(Hg)×w(F)］/［w(Mn)× w(Ni)］、w(Hg)/w(Ni)、［w(Au)×w(Ag)］/［w(Mn)×w(Ni)］、［w(Hg)×w(F)×w(B)］/［w(Mn)×w(Ni)×w(W)］如图4所示(w表示元素标准化值的自然对数值)。地球化学参数可表示矿上晕元素相对于矿下晕元素的发育程度，其值越大，反映出矿体头部晕特征越明显，深部矿体存在的可能性越大;反之，则反映出矿体尾部晕特征发育，深部存在矿体的可能性越小(李惠等，1998)。

图4 万年银金矿床Ⅱ号原生晕带轴向地球化学参数变化曲线Fig.4 Curves showing geochemical parameters in theⅡprimary halo belt，Wannian sliver-gold deposit

矿区Ⅱ号原生晕带轴向地球化学参数显示由浅部到深部各地球化学参数具有相似的“降－升－降－升”变化特征，总体上呈振荡波动变化，波峰与波谷相间出现。结合该矿床成矿期具有多阶段的特点，其地球化学参数的强烈振荡波动也同时是成矿期热液多阶段活动的一种反映，是成矿多阶段反复叠加的结果。这种原生晕轴向化学参数叠加结构与李惠等(1999)总结的金矿床原生晕轴向地球化学参数理想模型中的C型相近，特别是由－185m中段到－235m中段各参数均表现为“升”的趋势，指示Ⅱ号原生晕带－235m中段以下矿体向深部还有很大延伸，具良好的成矿远景。

3 矿体侧伏特征

根据矿区钻孔资料绘制出9～72线区段Ⅶ号矿体厚度纵向等值投影图(图5)可见，厚度(＞2m)部分主要分布于矿体中、东部的0～－200 m区域地段。从西向东，矿体厚度总体由薄向厚增大的趋势。其中，极薄矿体主要分布在5～6线地段。薄矿体主要分布在10～38线地段，而较厚的矿体主要集中于34～60线的中深部，局部厚度在2m以上的集中段主要分布在14、26～30、46～54以及64～72线矿体的中下部。由图上可以看出高值区(＞2m)的厚矿体多呈不规则瘤状，往NE向深部侧伏延伸，且呈雁列状排布，侧伏角约在45°～55°左右，每个瘤状矿体厚度中心之间的距离大约在180～240m左右。可见矿体厚度在其侧伏方向上显示出周期性和规律性的膨胀和收缩。吴树仁(1993)、白万成(1995)和郭光裕(1999)等研究认为地质体内的断裂构造形迹呈波状分布是构造应力在地质体内作用的结果。不同作用方向和大小的力会形成不同的断裂裂隙性质和规模，这可以反映在受控于该断裂构造的矿脉的走向、倾向、侧伏向等产状要素上并在不同的空间上发生变异。同时，对于控矿断裂而言，某一成矿阶段的构造应力会造成矿石品位、矿脉厚度、矿体形态与规模等矿化要素在空间上有规律性和周期性的分布。不同成矿阶段的构造应力会使上述产状要素和矿化要素发生变异而在空间上周期性重复出现。因此，我们现在所观察到的这种矿脉形态的波形分布是成矿前及成矿期构造应力作用积累的结果，也是对矿区成矿期控矿断裂构造分布的直接反映。万年银金矿床矿体产出于断裂构造破碎带内部，断裂构造作为矿体的根本控制因素，对成矿元素、矿化类型和矿化带分布以及矿体形态、产出位置有着明显的控制作用，其在深部区域的发育和延伸对于盲矿体预测起着决定性作用。图5可见在46～54线、72线以东区段－235m以下深部区域具有厚矿体有明显的延伸趋势，显示出控矿断裂在该区域有良好的发育和延伸，是对深部盲矿体存在的构造标志和背景。

图5 Ⅶ号矿体厚度等值线图(m)Fig.5 Contours of thickness for the Ⅶ ore body(m)

4 矿体厚度与品位乘积变化规律

成矿元素品位与矿体厚度乘积是矿床勘探的重要指标之一，利用其变化规律可有效地对矿区深部成矿远景进行衡量。根据对矿区9线～72线各勘探线Ⅶ号脉Ag、Au品位与矿脉厚度统计，可得Ag品位与厚度乘积等值线图(图6)及Au品位与厚度乘积等值线图(图7)。可以看出，Ag、Au品位与矿体厚度乘积的高值区均呈近椭圆形或纺锤形雁列状排布，并向NE方向往深部侧伏延伸，其中以46线～54线、66线以东地区延伸趋势较明显，其中46线～54线地区与Ⅱ号原生晕带前缘晕异常套和较好，是深部找矿较有希望的地区。

5 讨论与结论

通过对万年银金矿床Ⅶ号矿体原生晕特征及矿体变化规律的系统研究，结合地质资料分析，可得出以下认识:

(1)万年银金矿床原生晕的轴向分带序列Sn－Au－Cu－Mn－B－Mo－Ni－F－Ag－Cr－Hg－Pb－Zn－W－Bi，与正常原生晕序列相比，其前缘晕、近矿晕、尾晕元素分布混杂，呈现出明显的逆向分带特征。在矿区中段Ⅱ号原生晕带46线地区深部Hg元素异常发育完全，面积较大且向下未封闭，出现前缘晕与尾晕叠加。同时矿区Ⅱ号原生晕带轴向地球化学参数［w(Hg)×w(F)×w(B)］/［w(Mn)×w(Ni)×w(W)］、w(Hg)/w(Mn)、［w(Au)×w(Ag)］/［w(Mn)×w(Ni)］、w(Hg)/w(Ni)、［w(Hg)×w(F)］/［w(Mn)×w(Ni)］呈“降－升－降－升”振荡波动变化特征，既是成矿期热液多阶段活动的一种反映，也是成矿多阶段反复叠加的结果。在矿体末端各地球化学参数均表现为“升”的趋势，反映出前缘晕元素富集程度增大，并且空间上与Pb、Zn、Sn、Mn、W 同位共存，与 F1、F2 因子代表的元素组合相近，表明该区域不仅发育有石英－多金属硫化物阶段的银矿化，还有主成矿阶段－细粒硫化物阶段的银、金矿化叠加，表明Ⅱ号原生晕带深部区域矿体还有较大的延伸，显示出良好的成矿远景。

(2)矿区厚矿体多呈不规则瘤状近雁列式近等距排布，往NE向深部侧伏延伸，侧伏角在45°～50°之间，每个瘤状矿体厚度中心之间的距离大约在180～240m之间。矿体厚度在其侧伏方向上显示出周期性和规律性的膨胀和收缩的波形变化，这是成矿前及成矿期构造应力作用积累的结果，也是对矿区成矿断裂构造分布的直接反映。矿区矿体厚度等值线图显示厚矿体在46～54线、72线以东区段－235m以下深部区域表现出明显的延伸发育，显示出以上区域发育有良好的成矿期断裂构造背景。Ag、Au品位与矿体厚度乘积的高值区均呈近椭圆形或纺锤形雁列状排布，并向NE方向往深部侧伏延伸，其中以46线～54线、66线以东深部地区延伸趋势较明显，可作为深部成矿远景的判别指示标志之一。

综合以上各成矿指示信息，矿区66线以东深部地段虽具有较好的成矿构造背景和矿体厚度与Ag、Au元素乘积高值延伸的指示，但无异常原生晕特征的表现，该区域可能为成矿期热液活动的导矿构造，不是主要的矿质沉淀堆积场所。相比之下，46～54线－235m中段以下深部地区则具有明显的原生晕异常表现、良好成矿构造背景以及矿体厚度与Ag、Au元素乘积高值延伸直接的指示等各种有利成矿信息的耦合，据此建立如下立体矿体预测模型(图8)。该区域在后续地质工作中取得了较好的效果。

图8 万年银金矿床Ⅶ号脉深部预测立体模式图Fig.8 Schematic diagram of three-dimensional deep orebody prediction model for the veinⅦin the Wannian silver-gold deposit

感谢两位审稿专家提出的宝贵建议!野外工作和资料收集得到了万年矿业公司领导的热情关心和大力支持，野外采样过程中得到了中国地质大学(武汉)硕士研究生刘阁、阮班晓的帮助，在此一并表示感谢。

［注释］

① 江西地质勘查局赣东北地质队.1987.江西省万年县虎家尖矿区Ⅵ、Ⅶ号矿带13－88线区段银、金矿勘探报告［R］

② 江西地矿资源勘查开发有限公司.2006.江西万年县虎家尖银、金矿区13－88线区段地形地质图［R］

Bai Wan-cheng，Ren Lin-zi，Qin Min，Wang Chun-hong.1995.The method of simulating facture waveform and application on deep prediction of vein gold deposits［J］.Journal of Precious Metallic Geology，4(3):214－221(in Chinese with English abstract)

Dong Fa-xian.1995.Geological structures and their control of silver-gold deposits in Hujiajian area，Jiangxi［J］.Journal of Geomechanics，2(1):50－58(in Chinese with English abstract)

Guo Guang-yu，Lin Zhuo-hong，Guo Wan-chao.1996.The application of waveform’s distribution in deep part prognostication of geo-environment［J］.Metallurgical Geological Dynamic，6:3 － 7(in Chinese with English abstract)

Jiang Jing-ye，Cheng Jian-ping，Qi Shi-hua，Xiang Wu-fu.2006.Applied Geochemistry［M］.Wuhan:China University of Geosciences Press:1－180(in Chinese with English abstract)

Konstantinov M M，Strujkov S F.1995.Application of indicator halos(signs of ore remobilization)in exploration for blind gold and silver deposits［J］.Journal of Geochemical Exploration，54(1):1217

Li Hui，Yu Bin，Li De-liang，Ma Jiu-ju，Xu Li-yun，Wei Jiang，Zhao Jiaxiang，Zhang Qiang，Si Shu-yun，Yang Ya-juan，He Rong-hua.2010.Summary of new methods on deep prediction of geochemical exploration［J］.Mineral Exploration，1(2):156 － 160(in Chinese with English abstract)

Li Hui，Zhang Wen-hua，Liu Bao-lin，Wang Jing-chen，Guo Rui-dong.1999.The study on axial zonality sequence of primary halo and some criteria for the application of this sequence for major types of gold deposits in China［J］.Geology and Prospecting，35(1):32 －35(in Chinese with English abstract)

Li Hui，Zhang Guo-yi，Yu Bin.2006.The structural superimposed halo model of deep blind ore body prediction and its prospecting effect［M］.Beijing:Geological Publishing House:1－146(in Chinese)

Li Hui，Zhang Wen-hua，Chang Feng-chi.1998.Models of overprint of primary halo for large，mega size blind gold ore deposits［M］.Beijing:China Metallurgical Industry Press:1－195(in Chinese with English abstract)

Li Hui，Zhang Wen-hua，Chang Feng-chi，Qu Yuan-fei，Zuo Hong-we.1999.The standards of the opposite direction zoning of primary halo and‘the turning points’of geochemical parameters in the direction of the axis in the blind gold deposits［J］.Journal of Guilin Institute of Technolgy，19(2):114 －117(in Chinese with English abstract)

Liu Chong-min.2006.Progress in studies on primary halos of ore deposits［J］.Acta Geologica since，80(10):1528 －1538

Piao Shou-cheng.Lian Chang-yun.1994.A new kind of method for determinating zoning sequence of primary halo:The method of Barycenter［J］.Geology and Prospecting，30(1):63 － 65(in Chinese with English abstract).

Piao Shou-cheng，Yang Yong-qiang，Lian Chang-yun.1996.A summary of studying methods in primary halo zoning sequence［J］.Global Geology，15(1):44 －48，53(in Chinese with English abstract)

Qiu De-tong.1989.A new method for zoning sequence determination of indicator elements within the primary halo of a deposit［J」.Geology and Prospecting，30(8):51－53(in Chinese with English abstract)

Sun Hua-shan，Sun Lin，Cao Xin-zhi，Wang Chao，Tan Jun，Liu Fujiang，Yang Kai-chun.2008.Axial/vertical zoning characteristics of primary halos and geochemical exploration indicators for deep ore body prognosis in Shangzhuang gold deposit，Northwest Jiaodong peninsula，Shandong Province［J］.Mineral Deposits，27(l1):64 － 70(in Chinese with English abstract)

Wu Shu-ren.1993.The research of ore-controlling fracture geometry and kinematic and its ore-controlling regularity［J］.Geology and Prospecting，(1):1－6(in Chinese with English abstract)

Xu lei-ming，Bai Ge，He Quan-bo.1993a.Geological characteristics and genesis of the hujianjian silver-gold deposit in Wannian county，Jiangxi Province［J］.Mineral Deposits，12(2):138 － 147(in Chinese with English abstract)

Xu lei-ming，Wu cheng-yu，Bai Ge，He Quan-bo.1993b.The characteristics of some major minerals and mineralization of the Hujiajian silvergold deposit in Wannian county，Jiangxi Province［J］.Acta Mineraloica Sinica，13(3):254－262(in Chinese with English abstract)

Xie Qing-lin.1992.Concentration index method to determine axial element sequence of primary halo［J］.Geology and Prospecting，33(6):55(in Chinese English abstract)

Yu Jin-sheng，Li Yu-wei.1985.Geological factor analysis［M］.Beijing:Geological Publishing House:1－418(in Chinese with)

Zhang Ding-yuan.1989.Study on methods calculating for element zonality of primary halo［J］.Geology and prospecting，25(5):45 － 49(in Chinese with English abstract)

Zhang Yan-yan，Wang Jian-xin，Zhao Zhi，Yin Zhan-jun，Xie Hai-dong，Zang Xing-yun.Application of R-type clustering analyses in subdivision of mineral Stages-an example from gold ore area in Damiaozi-Caiqiangzi of Huadian［J］.Global Geology，25(1)，29 －38(in Chinese with English abstract)

Zhao Peng-da，Hu Wang-liang，Li Zi-jin.1994.Statistics on deposit projecting［M］.2nd ed.Beijing:Geological Publishing House:260 －276(in Chinese)

Zhao Qi.1989.The index of element＇s proportion calculating method on the vertical primary halo［J］.Geophysical and Geochemical Exploration，13(2):157 －159(in Chinese with English abstract)

［附中文参考文献］

百万成，任林子，卿敏，王春宏.1995.断裂面波形模拟预测法在脉状金矿床深部定位预测中的应用［J］.贵金属地质，4(3):214－221

董法先.1995.江西虎家尖银、金矿区构造特征与矿化研究［J］.地质力学报，2(1):50－58

郭光裕，林卓虹.郭万超.1996.波形分布在深部预测区地质环境推测中的应用［J］.冶金地质动态，6:3－7

蒋敬业，程建萍，祁士华，向武副.2006.应用地球化学［M］.武汉:中国地质大学出版社:1－326

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

李惠，张文华，刘宝林，王敬臣，郭瑞栋.1999.中国主要类型金矿床的原生晕轴向分带序列研究及其应用准则［J］.地质与勘探，35(1):32－35

李惠，张国亿，禹斌.2006.金矿区深部盲矿预测的构造叠加晕模型及找矿效果［M］.北京:地质出版社:1－146

李惠，张文华，常凤池.1998.大型、特大型金矿盲矿预测的原生晕叠加模型［M］.北京:冶金工业出版社:1－195

李惠，张文华，常凤池，曲媛菲，左宏伟.1999.金矿盲矿预测的原生晕轴向“反(向)分带”和地化参数轴向“转折”准则［J］.桂林工学院院报，19(2):114－117

刘崇明.2006.金属矿床原生晕研究进展［J］.地质学报，80(10):1528－1538

朴寿成，杨永强，连长云.1996.原生晕研究方法综述［J］.世界地质，15(1)44－47