娄里金矿位于邯郸永年洪山碱性杂岩体北侧火山凝灰岩及凝灰熔岩中。目前已在4 条含金矿化蚀变带中圈出12 个矿体。李惠等通过对胶东、豫陕小秦岭和燕辽等中国主要金矿集中区典型金矿成矿作用的研究,将金矿床的形成归纳为4 个阶段,分别为黄铁矿-石英阶段、金-石英-黄铁矿阶段、金-石英-多金属硫化物阶段和石英-碳酸盐阶段^[1]。已有的研究结果显示,娄里金矿床形成温度为150℃左右,成矿深度在59~191 m,平均为92 m。金矿成因属浅成低温热液矿床,矿床为多期热液活动叠加的结果^❶。其中,金-石英-黄铁矿阶段为娄里金矿的主要成矿阶段,金-石英-多金属硫化物阶段不发育。矿床原生晕理论研究表明,热液中的成矿物质在溶解、迁移、沉淀过程中,由于内部和外界物理化学条件的变化,导致成矿指示元素及其组合的空间分布,形成规律性分带,每个矿体有其特定的前缘晕、近矿晕和尾晕元素组合。研究热液矿床的原生叠加晕特征,即研究识别每一期次形成的矿体-晕的轴向分带及不同期次形成的矿体-晕在空间上的叠加结构。矿化在空间上表现为成矿晕的叠加、再组合。可通过矿床成矿晕的空间组合变化规律,判断金矿体剥蚀程度、分析深部矿化富集情况,据此推测深部矿化特征^[2-7]。因此矿区元素轴向分带性的研究具有重要的理论和实际意义。本次分别挑选了娄里金矿区北侧ZK16-1 和南侧ZK11-2 两个钻孔,按10 m间隔采取岩心样品,分析有关指示元素的轴向变化规律,为矿区找矿勘查提供理论依据。

❶河北省地矿局地质一队. 河北省邯郸市洪山岩体及外围金矿成矿规律研究报告[R]. 2015.

1 区域地质

研究区位于华北地台山西断隆与华北断坳两个Ⅱ级构造单元结合部位的武安凹断束中,地壳运动频繁,岩浆活动十分强烈,形成一系列北东向褶皱和断裂及十分复杂且规模巨大的构造岩浆岩带（图 1）。研究区火山喷发作用强烈,火山机构受区域性构造所控制分布在邯郸深大断裂与鼓山—紫山大断裂之间的北北东向构造岩浆岩带内。

与金矿成因关系密切的洪山碱性杂岩体沿紫山—鼓山断裂侵入,且分布于断裂两侧。陈斌等^[8]对洪山正长岩中的锆石采用SHRIMP方法测年,确定其侵位年代为（135±2.7）Ma,黑云母Rb-Sr同位素等时年龄为（120.3±2.4）Ma。测年结果显示,洪山岩体的形成时代为白垩纪。杂岩体岩性较为复杂,以正长岩为主。岩体北侧火山凝灰岩发育,野外实测地质剖面显示,岩体北侧娄里组火山活动的性质表现为6个完整的火山爆发相-岩浆喷溢相的演化序列。含矿热液受特定旋回的岩性岩相控制,金矿化主要赋存于娄里组二段的火山凝灰岩中。受洪山岩体火山机构的影响,矿区范围内多数断裂带方向常在NW至NNE之间变化,平面上具弧形变化特点。

2 矿床地质特征

目前已在洪山火山机构北部娄里一带初步圈出4 条含金矿化蚀变带。蚀变带主要发育于白垩系娄里组二段（K₁L²）火山碎屑岩地层中,岩性以粗面质含角砾凝灰岩或粗面质角砾凝灰熔岩为主。蚀变带中岩石较破碎、片理发育,蚀变带受构造控制明显。4 条蚀变带总体上呈现出以洪山火山机构为南部交汇点,向北呈NWW、SN、NEE放射状分布态势（图 2）。矿体产在矿化蚀变带中,受蚀变带控制,以Ⅱ号矿化蚀变带矿体规模最大。4 条矿化蚀变带均呈陡倾斜,倾角多在70°以上（图 3）。4 个矿化蚀变带矿化类型均为破碎带蚀变岩型,矿脉主要由劈理、片理发育的碎裂蚀变凝灰熔岩组成。主要矿化蚀变为黄铁矿化,地表多已褐铁矿化,其次为硅化,多表现为微细石英脉呈条纹状分布。金矿化严格受限于构造热液蚀变范围,矿体多呈细脉或两向延长薄饼状产出,形态较为复杂。水平和剖面上均呈舒缓波状弯曲,与矿化蚀变带一致。4 条矿化蚀变带中12 个矿体统计结果显示长度在80~308 m,厚度在1.09~9.4 m,品位在1.02×10^-6~4.8×10^-6,平均品位2.49×10^-6。

矿石物质组成较为简单。金属矿物主要为自然金,含银自然金,银金矿、黄铁矿、磁铁矿等。脉石矿物主要为长石、石英,绢云母。其中长石以钾钠长石、透长石、微斜长石、钠长石、中长石等偏碱性长石为主。矿石中常见晶粒结构（包括自形晶结构、半自形晶结构、他形晶结构）、碎裂结构、交代结构、包含结构等。矿石中常见的构造有块状构造、细脉浸染状构造等。

3 原生晕垂向分带特征 3.1 元素垂向分带序列

元素地球化学活动性的差异使得不同指示元素的原生晕在空间上呈现序次分布。在20 世纪70年代前苏联学者格里戈良提出指示元素垂向分带指数计算法的基础上,国内化探界学者先后提出了广义衬值法、浓集重心法、浓集指数法、比重指数法、含量梯度法等多种计算方法^[9-12]。对上述计算元素垂向分带法的优缺点我国学者也进行探讨和对比分析^[13-15]。如对目前应用较多的格里戈良分带指数法,叶庆森等认为该方法存在某截面无异常时,其分带指数梯度值△G无法直接计算而影响排序,此外元素在分带序列中的位置会因参与计算的指示元素不同或多寡而改变等缺陷^[13]。为便于对比分析和验证结果,本次采用格里戈良分带指数法和广义衬值法两种方法计算了元素垂向分带序列。所谓广义衬值法分带指数法是从衬值或衬度（异常含量Ca与背景值Cb之比）的概念扩展而来。叶庆森认为广义衬值法不仅克服了上述衬值法的缺陷,而且可揭示各指示元素在各截面的富集程度,从而能更客观地反映指示元素分带现象^[13]。叶庆森等通过数个具体矿床实例证明,该法更能反映客观规律。两种计算方法均需要指示元素在矿区范围内的背景值。本次根据娄里矿区地表 386 件岩石样品的分析结果,在反复剔除特异值（大于平均值加3 倍标准差）后,使数据基本服从正态或对数正态分布情况下,确定了相关指示元素的背景值（表 1）。

从表中可以看出,矿区金元素背景值为25.99×10^-9,远高于中国碱性正长岩（0.9×10^-9）、华北地台正长岩（0.61×10^-9）和大陆地壳（1.21×10^-9）^[16]。洪山杂岩体及其附近各类岩石中Au元素本底含量测试结果表明,娄里矿区容矿围岩（粗面凝灰熔岩）及其他岩性段（尤其是黑云辉石粗面斑岩）Au 成矿元素含量较高,这种成矿元素高背景特征是本区金矿是成矿的物质基础。

3.1.1 ZK11-2钻孔元素垂向分带

ZK11-2 钻孔采用格里戈良分带指数法计算的元素分带序列为：

Pb-As-Cd-Zn-Mo-W-Sb-Co-Au-Ag-Mn-Bi-Ba-Ni-Sn-Cr-Cu

根据广义衬值计算法得到ZK11-2 钻孔元素垂向分带序列为：

Pb-As-Bi-Cd-Zn-Co-Mo-Cu-W-Sb-Mn-Au-Ag-Sn-Ni-Ba-Cr

对比两种计算方法所得分带序列可以看出,两种方法虽在所得的元素排列顺序上存在一定差异,但在反映热夜活动期次及不同温度元素在空间叠加模式上基本相同。均能显示出钻孔上部主要为低温元素（As、Cd）组成的前缘晕和中温元素（Pb、Zn）组成的近矿晕,钻孔中部为高温元素组成的尾晕元素（W、Co、Mo）。上述元素组合反映了一次较为完整的热液活动过程。钻孔中下部又出现一次以矿前晕（低温元素组合）开始到高温元素结束的热夜活动。两种分带序列计算结果均显示,ZK11-2 孔由2 次热液活动元素叠加而成。钻孔ZK11-2测试数据中Au、Ag 矿化与As 元素关系密切,相关系数分别为r_Au-As=0.748,r_Ag-As=0.88。这与已有的研究认为矿床属浅成低温热液成因相吻合。

3.1.2 ZK16-1钻孔元素垂向分带

ZK16-1 钻孔采用格里戈良分带指数法计算的元素分带序列为：

Ba-Cu-Zn-Mn-Bi-Sn-Mo-Cd-Pb-Co-Cr-Ni-W-Au-Hg-Sb-As-Ag

根据广义衬值计算法得到ZK16-1 钻孔元素垂向分带序列为：

Cu-Zn-Bi-Sn-Mo-Sb-Cd-Pb-Cr-Co-Ni-W-Au-Hg-As-Ag-Mn-Ba

两种方法计算所得的元素分带序列基本一致,钻孔中上部及中下部出现2 处高温元素（Bi-Sn-Mo和Cr-Co-Ni-W）组合,反映了至少存在两次热夜活动过程。与ZK11-2 钻孔不同的是钻孔顶部未出现与金矿化关系密切的低温元素,这是由于矿区含矿地层为娄里组二段（K₁l²）,而ZK16-1 钻孔上部岩层为娄里组一段（K₁l¹）,其上部的金矿体含矿层位已风化剥蚀。钻孔深部虽出现低温元素组合,根据成矿条件推测出现金矿化的可能性较小。

3.2 指示元素含量垂向变化规律 3.2.1 ZK11-2钻孔指示元素垂向变化

娄里金矿区ZK11-2 钻孔元素垂向变化规律图（图 4）显示,地表附近低温元素（As、Sb、Hg）含量较高,部分中温元素如Pb 发育,高温元素（W、Sn、Mo）含量较低。钻孔中金元素含量较高区域与低温元素关系密切,而与高温元素关系较差。这与矿床成矿温度较低的现象相吻合。

ZK11-2 全孔岩性以粗面质晶屑凝灰熔岩为主（图 3）,全孔采样分析显示,该钻孔分别在据地下27 m和100 m左右处有2 处较强的金矿化。其中距地表 27 m左右的金矿化,前缘低温元素含量较高,尾晕元素含量低。低-高温元素比值（HgAsSb/WSnMo）的地化参数突出反映了这一现象,该地段也是金矿化的主要富集地段。对比钻孔高-中-低温元素轴向分布可发现,从该处向下,中温元素（Cu-Pb-Zn）、高温元素（W-Sn-Mo）含量依次升高。反映了一个较为完整的热液活动过程。

在距地表 100 m（标高50 m处）左右地段,低温元素含量较高,中高温元素含量较低,高温元素不仅含量较高,范围也较宽。表明第二次热夜活动时,前缘低温元素与第一次热液活动时形成的尾晕（高温）元素叠加。该期热液活动强度较高,金矿化可达0.84 g/t。随后钻孔中显示第二次热液活动的前缘晕低温元素含量降低,而近矿晕中温元素含量相对升高。反映了第二次较为完整的热液活动过程。

第二次矿化后,钻孔向下近100 m范围（标高20 m到-90 m）内未出现高含量的前缘元素。表明钻孔已穿透含矿层,深部出现金矿化的可能性较小。这也与矿床属浅成低温热液的认识相吻合。

3.2.2 ZK16-1钻孔元素垂向分带

ZK16-1 钻孔岩性以硅化晶屑凝灰岩为主（见图 3）。该钻孔取样分析结果显示（图 5）,第一次热夜活动形成的前缘晕元素在地表附近（标高130 m附近）含量较低,而近矿晕的中温元素含量相对较高。暗示着第一次金矿化作用形成的矿体可能已受剥蚀。这与钻孔所处位置为娄里组一段地层,上覆金矿体赋存层位（娄里组二段）已风化剥蚀的事实相吻合。

在钻孔深部80 m处（标高50 m）,高-中-低温元素分布区间相吻合,结合元素轴向分带序列,说明第二次热液活动的前缘晕元素叠加在第一次热液活动的尾晕元素上,但热液含矿性较差。从标高-100 m开始向深部,矿前晕低温元素、近矿晕中温元素和尾晕高温元素含量依次降低,反映了完整的第三次热夜活动过程,该期矿化作用也相对较强。根据钻孔中元素轴向分布规律的变化趋势,标高-150 m以下,中-高温元素含量将会依次升高。而区内金矿化与低温元素关系密切,深部出现金矿化体的可能性不大。结合区域地质资料分析,深部应注意（正长）斑岩型铜钼矿化发生。

3.3 矿床原生叠加晕模型

用两种不同计算方法均得出元素分带序列存在明显的“前尾晕共存”现象,即由低温元素组成的前缘元素不仅出现在矿体前缘,而且在矿体中部或尾部也有出现。李惠等在研究原生晕垂向分带模型时发现很多矿床,尤其是金矿床原生晕垂向分带多具有反分带现象^[1-2]。若用热液床一次成矿作用形成的原生晕垂向分带理论是无法解释的,并根据热液矿床成矿晕严格受构造控制,热液矿床具有多期多阶段叠加成矿成晕特点,创造性地提出了原生叠加晕理论及其技术方法。该理论认为热液矿床每次成矿作用都有自己的前缘晕、近矿晕和尾晕。不同期次成矿形成的原生晕在空间上互相叠加。据此不仅要研究每次成矿形成的原生晕垂向分带特征,而且还要研究与识别不同成矿阶段形成的原生晕在空间上叠加呈现的结构形态,建立矿床垂向分带的原生叠加晕模型。

根据娄里金矿区ZK11-2 和ZK16-1 两个钻孔垂向元素分带序列及相关元素分布规律,绘制了娄里金矿体原生叠加晕模型示意图（图 6）。其中地化参数的计算采用前缘晕(HgSbAs)元素含量标准化累乘/尾晕(WSnMo)元素含量标准化累乘。

4 结论

（1）采用分带指数法与广义衬值法对钻孔数据的指示元素垂向分带序列计算结果显示,不同方法虽在元素排列顺序上存在一定差异,但在反映热夜活动期次及元素空间叠加模式上基本相同。

（2）钻孔元素分带序列存在明显的前尾晕共存现象,表明矿体的形成非一次热夜活动所致。两个钻孔的垂向分带都由两个完整的小循环（前缘晕、近矿晕和尾晕的元素都存在）组成,说明金矿体的形成至少经历了两次热夜活动叠加。

（3）原生晕轴向分带序列的确定及深部成矿潜力的预测应建立在全面正确理解和掌握区域成矿条件和控矿因素基础上。结合区域成矿条件,推断ZK11-2 钻孔已穿透含矿层,且深部出现金矿化的可能性较小。ZK16-1 钻孔顶部矿体已风化剥蚀,根据钻孔元素分布趋势,推测钻孔深部（标高-150 m以下）中-高温元素含量会依次升高,出现金矿化体可能性不大,但应注意（正长）斑岩型铜钼矿化的发生。