河南栾川鱼库Mo矿床黄铁矿热电性特征及成矿规律

贾文娟，王功文，韩江伟，燕长海，宋要武 ，刘传权，冯 源，王宏卫，杨 帆，孔 亮，郭娜娜，张旭晃，李军军

(1.中国地质大学(北京)地球科学与资源学院，北京 100083；2.河南省地质调查院，河南 郑州 450007；3.西北有色地质研究院，陕西 西安 710000；4.河南省栾川县地质矿产局，河南 栾川 471500)

河南栾川鱼库Mo矿床黄铁矿热电性特征及成矿规律

贾文娟1，王功文1，韩江伟2，燕长海2，宋要武2，刘传权2，冯 源3，王宏卫4，杨 帆1，孔 亮1，郭娜娜4，张旭晃4，李军军3

(1.中国地质大学(北京)地球科学与资源学院，北京 100083；2.河南省地质调查院，河南 郑州 450007；3.西北有色地质研究院，陕西 西安 710000；4.河南省栾川县地质矿产局，河南 栾川 471500)

栾川鱼库钼矿床位于华北克拉通南缘，矿床的形成与燕山期岩浆-矽卡岩热液成矿作用有关，是近年来发现的大型隐伏钼多金属矿床。对鱼库矿区骨干剖面的深部钻孔(>1 000 m)中黄铁矿开展热电性标型特征分析，发现样品中黄铁矿以N型为主，总含量达96%。根据黄铁矿导型及热电性参数变化特征推测出矿体向深部有一定规模的延伸。根据热电性实验数据计算出成矿温度范围在270～332 ℃之间，结合矿床矿石特征判断该矿床属于浅成中温热液成矿。并利用地质统计学知识，运用ArcGIS对计算出来的成矿温度进行克里格插值分析，探讨矿区深部成矿规律。

黄铁矿； 热电性特征； 成矿预测；鱼库Mo矿床；河南

0 引 言

黄铁矿是最常见的金属硫化物，同时也是多金属硫化物矿石典型矿物组分之一[1]。黄铁矿标型特征不仅可以用于判断成矿物源，也是重要的找矿标志，因此其成为矿床成因和找矿信息的重要研究内容[2-3]。黄铁矿标型特征在判断成矿温度、剥蚀程度、矿床规模和寻找隐伏矿体方面均起到重要作用[4-6]。

近年来，根据“十二五”科技支撑课题和整装勘查项目，河南省地质调查院与中国地质大学(北京)联合，通过三维多元地学信息集成找矿研究[7]，开展了栾川矿集区Mo、Pb-Zn矿增储的矿产勘查工作。截至目前，深度达1 000 m钻孔超过了25个，为该区深部钼多金属矿增储奠定了基础。本文在矿集区新发现的鱼库矿区优选典型钻孔，开展深部找矿规律研究。

图1 栾川钼钨铅锌银矿床地质简图[10]Fig.1 Regional geological map of the Luanchuan Mo-W-Pb-Zn-Ag deposit1.下古生界陶湾群；2.上古生界栾川群；3.中元古界官道口群；4.中元古界宽坪群； 5.太华群； 6.燕山期花岗岩；7.村镇； 8.压性断裂带； 9.张性断裂带； 10.张扭性断裂带； 11.压扭性断裂带； 12.性质不明断裂带； 13.背斜轴； 14.向斜轴； 15.倒转背斜轴； 16.地质界线； 17.平行不整合界线； 18.角度不整合界线； 19.钼矿床 ；20.铅锌矿床； 21.隐伏钼(钨)矿

本文借鉴以往黄铁矿在金矿成矿预测中运用的技术方法，首次利用钻孔样中黄铁矿较系统地研究栾川鱼库钼 Mo多金属矿床的成矿规律。结合地质、地球化学信息，综合分析黄铁矿形态、热电性标型在垂向空间上不同深度、不同矿化程度、不同岩性样品中的变化规律，探讨黄铁矿标型特征与成矿之间的关系，并对矿床的成矿温度、剥蚀程度进行了计算。在此基础上，利用地质统计学分析，运用ArcGIS地质统计学空间分析模块，对研究区黄铁矿样品估算的成矿温度进行了克里格插值，分析研究区钼多金属矿床成矿环境及其深部找矿地段。

1 地质概况

1.1 区域地质背景

栾川鱼库隐伏斑岩钼(钨)矿床位于南泥湖—上房特大型钼矿田的西南侧，是近年来新发现的斑岩—矽卡岩矿田(图1)。区域位于中朝准地台南缘的豫西断隆、三川—栾川陷褶断带中，陷褶断带北部为熊耳山隆褶带，南部为伏牛山隆褶带，是典型的活动性大陆边缘[8-9]。区内主要发育下古生界奥陶系陶湾群、新元古界青白口系栾川群、中元古界蓟县系官道口群和长城系熊耳群地层。本区构造复杂，最显著的构造是NNW向展布的两条区域性大断裂，即北部的马超营大断裂和南部的黑沟—栾川—维摩寺大断裂及一系列近于平行向南逆冲的推覆断层，对卢氏—栾川台缘褶皱带的地质发展和成岩成矿有着重要的控制作用[9]。

区域岩浆岩活动剧烈，整个南泥湖岩体是由燕山期斑状钾长花岗岩与斑状黑云母花岗岩组成的复式岩体。呈不规则椭圆形小岩株产出，岩体的延伸方向与构造线方向一致，主要为富硅、高钾的花岗岩体，与研究区内的 Mo-Pb-Zn矿关系密切[8-10]。

1.2 矿床地质特征

图2 栾川钼铅锌银矿带地质简图(据Cao等[9]修编)Fig.2 Sketch geologic map of the Luanchuan Mo-W-Pb-Zn-Ag deposit1.下古生界陶湾群三岔口组；2.新元古界栾川群鱼库组；3.栾川群大红口组；4.栾川群煤窑沟组；5.栾川群南泥湖组；6.栾川群三川组；7.中元古界官道口群白术沟组；8.燕山早期斑状二长花岗岩；9.前加里东期变辉长岩；10.前加里东期变正长斑岩；11.断层；12.层控型矿化带；13.矽卡岩型矿化带；14.矿化带编号

鱼库矿床属于石宝沟矿田，与其西北方向的南泥湖矿田属于同一成矿系统(图1)。鱼库钼铅锌矿床含矿地层主要为新元古界栾川群三川组(Pt3s)碳酸盐岩夹碎屑岩地层(图2)。三川组下段(Pt3s1)以灰白色变石英砂岩为主，夹钙质粉砂岩。岩石致密坚硬，内有细粒星点状黄铁矿，风化后呈“火烧皮”状，底部含燧石，厚度326.55 m。三川组上段(Pt3s2)上部主要为灰白色中厚层大理岩，岩石质纯，下部为青灰色条纹条带状大理岩夹钙质粉砂岩薄层，局部岩石矽卡岩化，可见透闪石及绿帘石化斑点，厚度289.86 m(图2)。

鱼库隐伏钼(钨)矿体周围构造发育。矿体位于一个倒转背斜的北西翼。西北侧和东南侧各发育两个北东—南西向压性断裂带，西南部发育多个北西西向压性断裂带[6]。沿断裂带热液活动强烈，控制着区域钼、铅锌、银的矿化[7]。钼主矿体主要分布于二长花岗斑岩体内，矿体在空间分布上与花岗斑岩顶面起伏形态相“整合”，主要赋存于斑岩体外接触带，总体向北西方向侧伏，侧伏角20°左右。目前勘探钻孔见矿最大厚度265.74 m，最小厚度4.38 m，矿体品位0.03%～0.10%。钼矿体岩控作用明显。矿石类型主要为斑岩型，其次为矽卡岩型。钼钨矿石和钼矿石主要赋存于矽卡岩中。矿石结构为鳞片状结构、镶嵌结构、自形—它形晶粒状结构、交代环状结构及放射状结构，构造为浸染状、细脉浸染状和细脉状，矽卡岩型及大理岩型矿石呈细脉状、网脉状，矽卡岩矿石以浸染状为主，细脉状次之。金属矿物有辉钼矿、白钨矿、黄铁矿、磁黄铁矿、磁铁矿、闪锌矿、方铅矿、黄铜矿；脉石矿物有钾长石、石英、斜长石、黑云母、透辉石、阳起石、次闪石、电气石、锆石、萤石、绿帘石、绿泥石、绢云母等[10]。

1.3 黄铁矿产出特征

黄铁矿的结晶习性对其形成环境的变化反映较为敏感，在低的过饱和度、低硫逸度以及过高或过低于其最佳形成温度的条件下，{100}形较为发育，且晶形相对简单; 在适中的温度、缓慢的冷却速度和有充分物质来源的高过饱和度、高硫逸度条件下，{210} 和{111}较为发育，且晶形相对复杂[2-4,11-12]。

图3 鱼库矿床黄铁矿微形貌特征Fig.3 The characteristics of pyrites from the Yuku deposita.花岗岩中自形黄铁矿；b.花岗岩中侵蚀结构黄铁矿；c.钾化花岗岩中碎裂结构黄铁矿脉；d.矽卡岩中骸晶结构黄铁矿；e.矽卡岩中侵蚀结构黄铁矿；f.辉长岩中黄铁矿脉；Py.黄铁矿；Mo.辉钼矿；Pl.斜长石；Kfs.钾长石

研究区样品中黄铁矿大部分被交代在矿石中呈浸染状分布，或呈颗粒状及脉状生长，其中花岗岩中含少数自形-半自形立方体黄铁矿，颗粒较小(图3a)，一般为0.1～0.5 mm，侵蚀结构黄铁矿颗粒较大(图3b)，可达1 mm以上，钾化花岗岩中还生长碎裂结构的黄铁矿脉，脉宽在0.7～1.5 mm之间(图3c)；矽卡岩中黄铁矿也多为交代结构(图3d)，且颗粒粒径一般较大，为1～2 mm，图中黄铁矿沿大颗粒石榴子石环带生长。矽卡岩中还生长被辉钼矿交代的黄铁矿(图3e)。辉长岩中黄铁矿脉(图3f)，脉宽2.2～4 mm之间。少数辉长岩中黄铁矿为它形粒状，粒径为0.5 mm左右。

2 黄铁矿热电性标型

2.1 测试原理及应用

矿物的热电性是指矿物晶体在温差条件下产生电势差的物理性质，包括热电系数和导电类型(简称导型)两层含义。利用热电仪测量矿物在冷端与热端的电位差，可以实现矿物热电性的测量。热电性系数α是指单位温差的热电势[13]，表达式为α=ΔE/ΔT。式中：α代表热电系数；E代表热电动势；ΔT代表温差。

黄铁矿是半导体矿物，其导电类型有两种，一种是电子型(N型)，另一种是空穴型(P型)。当热电动势E为负值时，矿物表现为N型导电；E为正值时，矿物表现为P型导电。黄铁矿主元素含量比为亏Fe(即S过剩)时显示空穴导型(P型)，所含的阴电荷过剩显示空穴型导型(N型)[2-4,14-15]。一般情况下，黄铁矿形成深度越大，其成矿温度越高，压力越大，热电动势趋于负值，黄铁矿表现为电子型(N型)；反之，黄铁矿形成深度越浅，其成矿温度越低，压力越小，热电动势趋于正值，黄铁矿表现为空穴型(P型)[16-18]。通过对黄铁矿热电系数和热电导型的测量与分析，可以有助于研究矿床成矿温度和剥蚀率及矿体深度等问题，进而推测矿体深部成矿前景[2-4]。

2.2 测试方法及结果

笔者在鱼库矿床对矿体进行了坑道与钻孔的样品采集工作。采样钻孔分别为分布于03线勘探线上的ZK0301、ZK0307与ZK0315，样品标高垂向分布范围200～1 330 m。

样品中的黄铁矿多呈浸染状与脉状分布，粒度与晶形均差别较大，笔者共选取了37个岩石样品中的黄铁矿单矿物。在中国地质大学(北京)成因矿物实验室，利用BHTE—06型热电系数测量仪进行实验。冷端温度设置为20 ℃左右，热端温度设置为80 ℃左右，活化温度设置为60 ℃[19]。在每个样品中精选50粒黄铁矿作为实验对象，实验结果详见表1。

研究结果表明，研究区的矿体黄铁矿热电系数α变化范围为-308～321.9 μV·℃-1。以N型为主，热电性系数介于-308～-12 μV·℃-1之间；少量P型黄铁矿，其热电系数介于36.8～321.9 μV·℃-1。

3 分析与讨论

3.1 成矿温度

大量经验证明黄铁矿的热电系数可能与其形成温度有一定关联，1964年，A.戈尔巴乔夫通过对大量数据的统计分析，做出了黄铁矿的热电性-温度图，并从中获得线性方程，即t=(704.51-|α|)/1.818 (N型)与t=3×(122.22+α)/5.0 (P型)。根据上述两式，用黄铁矿热电性数据计算出矿床的成矿温度。鱼库矿区矿化样品中黄铁矿成矿温度变化范围为270～332 ℃。这与流体包裹体测温结果(280～380 ℃)基本吻合*河南省地质调查院.栾川铅锌矿区深部资源勘查技术集成.2015.，矿区无矿化或少量矿化岩石样品中黄铁矿成矿温度变化基本分布于257.9～300 ℃之间。总体来看，矿化样品中黄铁矿的成矿温度大于无矿化或少量矿化岩石样品，由此也能指示出黄铁矿具有不同期次和成因。综合前面对黄铁矿晶形统计分析认为该矿床为中温矿床。

3.2 矿体剥蚀率

根据黄铁矿热电系数可以计算黄铁矿的热电性参数XNP,利用XNP可以计算出矿体剥蚀率γ (矿体相对矿化总长度的剥蚀百分比)[11]。黄铁矿的热电性参数XNP = (2fⅠ+fⅡ)-(fⅣ+ 2fⅤ)，其中：fⅠ表示黄铁矿样品中α>400 μV·℃-1的百分比；fⅡ表示α为200～400 μV·℃-1的百分比；fⅣ表示α为0～-200 μV·℃-1的百分比；fⅤ为α<-200 μV·℃-1的百分比。

最后利用γ=50-XNP /4计算鱼库矿体剥蚀率，显示标高1 283.86 m矿体剥蚀率为81%，由标高1 283.86 m至1 090 m，矿体的剥蚀率波动上升，而由标高964 m到225 m处波动下降至54%。以上数据说明矿体被剥蚀到下部，在钻孔垂向上，由浅至深呈现低→高→低的轻微韵律式变化表明矿体向深部有可能有一定延伸。这与黄铁矿热电导型特征所反映的情况一致。

3.3 黄铁矿热电性的垂向空间分布特征

鱼库矿区钻孔ZK0301黄铁矿的导型多为N型，热电性系数在垂向上呈波动式上升，在浅部具有一定的稳定性(图4)。在接触带和蚀变较严重的964 m和224 m处，黄铁矿热电性系数出现特高值。取样深度由地下1 300 m至200 m，热电系数峰值由-200 μV·℃-1波动上升至-100 μV·℃-1附近。孔深854 m和孔深224 m处出现P型黄铁矿，854 m深处黄铁矿热电性系数数值分散，基本平均分布在0～300 μV·℃-1之间，224 m处黄铁矿热电性系数集中在20 μV·℃-1。

表1 鱼库矿区黄铁矿热电性特征

图4 钻孔ZK0301热电性系数与深度变化关系Fig.4 Line chart of thermoelectric coefficients of pyrites from different elevations in ZK0301

鱼库矿区钻孔ZK0307中大理岩和花岗岩样品中，黄铁矿导型都为N型，热电性系数分散分布于-100～-300 μV·℃-1之间(图7)。本钻孔样品中未出现P型黄铁矿。在垂向上黄铁矿热电性系数呈波动式递增并具有一定的稳定性。取样深度从孔深840 m至孔深710 m处的黄铁矿热电系数峰值由-200 μV·℃-1增长至-150 μV·℃-1。

图5 钻孔ZK0307热电性系数与深度变化关系Fig.5 Line chart of thermoelectric coefficients of pyrites from different elevations in ZK0307

鱼库矿区ZK0315钻孔N型黄铁矿的热电性系数峰值出现在-100～-200 μV·℃-1之间(图8)，P型黄铁矿热电性系数峰值范围分布于80～100 μV·℃-1之间， P型黄铁矿占样品总数的2%。透辉石大理岩中N型黄铁矿的热电性系数峰值出现在-240～-260 μV·℃-1之间， P型黄铁矿热电性系数分散分布于80～180 μV·℃-1之间，P型黄铁矿占样品总数12%。在垂向上黄铁矿热电性系数呈波动式上升并具有一定的稳定性。取样深度由深部1 066 m至浅部565 m，热电系数峰值由-200 μV·℃-1增长至-100 μV·℃-1附近，有时峰值低至-250 μV·℃-1。

图6 钻孔ZK0315热电性系数与深度变化关系Fig.6 Line chart of thermoelectric coefficients of pyrites from different elevations in ZK0315

总体上看，3个钻孔从地表至深部，热电性系数变化趋势一致，在孔深0～1 000 m 内都呈波动式变化，矿体剥蚀率呈现高—低—高的韵律变化，所以推测矿体向深部可能会有剥蚀百分比低的富矿段重复出现。此外，黄铁矿热电系数变化梯度很小，且出现P型黄铁矿，这表明矿体具有稳定延伸和产出的可能。

矿化较少或无矿化岩石中的黄铁矿热电系数则显示出方差相对较大、分布更为分散的特点。这指示在鱼库矿区，黄铁矿热电系数波动程度与矿化存在一定关系，在黄铁矿热点系数变化较快的区域易出现矿化。

3.4 黄铁矿热电系数离散性特征

黄铁矿热电系数离散度(σα′)可以较准确地反映出热电系数相差较大的不同样品其热电系数数值相对集中与分散的情况。离散度的计算采用σα′=|σα/α|× 100%[20-22]，其中α代表黄铁矿样品热电系数的平均值；σα为热电系数的标准差。σα'数值越大，指示黄铁矿热电系数的分散程度越大；σα'数值越小，指示黄铁矿热电系数的分散程度越小。σα'值可反映成矿条件的稳定性。在稳定性较好的成矿条件下形成的黄铁矿，其晶体性质较为接近，热电性系数离散程度较小。σα′值也可以反映出不同阶段黄铁矿的叠加程度。形成以N 型黄铁矿为主的成矿阶段，与形成以P型黄铁矿为主的成矿阶段比较，两阶段的成矿条件均相对稳定，具有较小的σα′值; 而在两成矿作用阶段相叠加时，会造成成矿条件的波动，σα′值也会相应增高。

图7 黄铁矿深度与离散度关系图Fig.7 The characteristics of thermoelectric coefficient dispersions of pyrites

根据研究区黄铁矿的产出深度与离散度关系图(图7)，深度1 100～1 300 m 之间的离散度投点相对分散，600～900 m之间的离散度投点相对集中。总体上呈现由浅至深，离散度高值数据所占比例逐渐增多的现象。因此，推测1 100～1 300 m 之间的成矿阶段热液活动有多次脉动，造成不同区域的叠加程度有所差异； 而600～900 m，低值数据所占比例较高，指示不同阶段叠加影响较小，成矿环境相对较为稳定。

3.5 黄铁矿热电性对钼品位的标识

通过分析黄铁矿热电系数离散度以及钼品位数据的变化情况，结合地质特征，对鱼库矿体内黄铁矿热电性分布特征和钼矿化的关系进行了研究。ZK0315 共采集12件样品，样品岩性为大理岩、花岗岩、辉长岩，垂向分布范围为孔深565～1 143 m，通过对比离散度以及钼品位的关系，由于钼品位的高值区域与离散度的高值区域相对应，可知钼的富集区域位于该矿体中部且经历了多次成矿作用，进行了多次叠加。

3.6 克里格插值

克里格插值(Kriging)又称空间局部插值法，是以变异函数理论和结构分析为基础，在有限区域内对区域化变量进行无偏最优估计的一种方法，是地质统计学的主要内容之一。由于观测数据满足内蕴假设,且其区域化变量的平均值是未知的常数，作者选择了普通克里格插值，由于数据不符合正态分布，在插值之前对数据进行了对数处理[4]，其插值采用的变差函数及参数如图8所示。

图8 鱼库矿区03勘探剖面成矿温度插值参数及变差函数Fig.8 Variation function and parameter of mineralization temperature of exploration line 03 in Yuku using Kriging interpolation

根据样品统计分析，样品品位空间服从球形分布，块金值为0.002 6，垂向连续性较好，垂向样品的相似性与连续性范围达112.49 m，最大连续范围达774.26 m。上述统计特征表明研究区矿体在深部具有较好的连续性，样品品位在垂向变化不大，具有较好的继承性特征，这与研究区岩浆-矽卡岩-脉状矿体隶属于一个成矿系统的理论成矿模式认知一致[4，23-24]。

根据图9可以看出，研究区域地层为三川组及白术沟组地层，由于这两个地层岩性均主要为有利成矿的大理岩，因此，浅部脉状矿体的发育不受地层控制，横向切穿倾斜的地层，深部块状矿体与岩体形状基本一致。用黄铁矿计算出来的成矿温度进行克里格插值，结合薄片鉴定以及地质分析，得出该地区矿床成矿温度为270～332 ℃之间，根据克里格插值结果(图10)显示，中部低温区域与细层状矿体相对应。左下角成矿温度呈环带状逐渐升高，此处成矿应与下部岩体有关。结合地质背景发现，浅部矿体生长于有利成矿的三川组及白术沟组大理岩中，受构造控制影响呈脉状分布切穿地层，深部矿体受构造-岩体控制，与褶皱核部上侵的岩体形状重合度较高。这与研究区岩浆岩成矿能力与成矿特征趋于一致[23]。

图9 鱼库03勘探剖面图Fig.9 Map of the section of exploration line 03 in Yuku

图10 鱼库钼矿床03勘探剖面成矿温度克里格插值图Fig.10 Mineralization temperature plot of exploration line 03 using Kriging interpolation in the Yuku Mo deposit

4 结 论

通过对研究区黄铁矿的热电性标识特征的研究，得出以下结论:

(1)样品中黄铁矿的导型以N型为主，热电性系数介于-308～-12 μV·℃-1之间。含少量P型黄铁矿，热电系数介于36.8～321.9 μV·℃-1之间。由浅至深离散度值逐渐升高，表明成矿环境由浅至深逐渐复杂。

(2)根据黄铁矿热电性数据，计算出样品中黄铁矿成矿温度变化范围为257.9～332 ℃，矿体中黄铁矿成矿温度变化范围为270～332 ℃，结合地质背景和矿床地质特征综合判断，该矿床属于中温岩浆-矽卡岩型热液矿床。

(3)钼品位的高值区域与黄铁矿热电系数离散度的高值区域相互对应，黄铁矿导型及热电性参数变化特征反映鱼库矿体向深部应有一定规模的延伸。

(4)根据黄铁矿成矿温度克里格插值与实测勘探剖面对比分析，结合研究区 Mo多金属成矿模式，认为鱼库钼矿体主要受构造-岩体约束。

[1] ABRAITIS P K, PATTRICK R A D, VAUGHAN D J. Variations in the compositional textural and electrical properties of natural pyrite: a review [J]. International Journal of Mineral Processing, 2004, 74:41-59.

[2] XUE J L, LI S R, SUN W Y, ZHANG Y Q, et al. Characteristics of the genetic mineralogy of pyrite and its significance for prospecting in the Denggezhuang Gold Deposit, Jiaodong Peninsula, China [J]. Science China, Earth Sciences, 2014, 57(4): 644-661.

[3] 于洪军, 申俊峰, 王磊, 等. 玲珑金矿黄铁矿热电性与微量元素标型及深部金矿化评价[J]. 地质与勘探,2011,47(4):615-623.

[4] WANG G W, FENG Y, CARRANZA E J M, et al. Typo morphic characteristics of pyrite: criteria for 3D exploration targeting in the Xishan Gold Deposit, China[J]. Journal of Geochemical Exploration,2016,64:136-163.

[5] 陈光远, 孙岱生, 邵伟, 等. 胶东金矿成因矿物学与找矿[M]. 重庆: 重庆科技出版社,1989：1-20.

[6] 申俊峰,李胜荣,马广钢, 等. 玲珑金矿黄铁矿标型特征及其大纵深变化规律与找矿意义[J]. 地学前缘,2013,20(3):55-75.

[7] WANG Gongwen, LI Ruixi,CARRANZA E M J, et al. 3D geological modeling for prediction of subsurface Mo targets in the Luanchuan District, China[J]. Ore Geology Reviews, 2015,71:592-610.

[8] 罗铬玖. 中国钼矿床[M]. 郑州: 河南科学技术出版社, 1991:106-119.

[9] CAO Huawen, ZHANG Shouting, SANTOSH M,et al. The Luanchuan Mo-W-Pb-Zn-Ag magmatic-hydrothermal system in the East Qinling metallogenic belt, China: Constraints on metallogenesis from C-H-O-S-Pb isotope compositions and Rb-Sr isochron ages[J]. Journal of Asian Earth Science,2015,111:751-780.

[10] 韩江伟, 郭波, 王宏卫, 等. 栾川西鱼库隐伏斑岩型 Mo-W矿床地球化学及其意义[J]. 岩石学报,2015, 31(6):1789-1796.

[11] 刘坤,刘家军,吴杰,等. 甘肃马坞金矿床8号矿体黄铁矿热电性特征及其地质意义[J]. 现代地质, 2014,28(4): 711-720.

[12] 姜绍飞,杜振国,刘之洋. 黄铁矿的晶形与热电性的标型及其实际意义[J]. 东北大学学报(自然科学版), 1995(5):458-463.

[13] 高来之. 金属矿床标型矿物热电性研究状况[J]. 矿物岩石地球化学通报, 1990(3):200-201.

[14] 张胜男,黄菲,黄涛, 等. 红透山黄铁矿的热电性研究[J]. 沈阳师范大学学报(自然科学版), 2011(2):181-184.

[15] 孙国胜, 初凤友, 胡瑞忠, 等. 我国主要金矿类型中黄铁矿“电子-空穴心特征及影响因素[J]. 矿物学报, 2004(3): 211-217.

[16] 李胜荣,陈光远,邵伟, 等．胶东乳山金矿田成因矿物学[M]．北京: 地质出版社,1993: 98-101.

[17] 曾祥涛,王建平,刘必政, 等. 陕西省双王金矿床五号矿体黄铁矿标型特征研究[J]． 地质与勘探,2012,48(1): 76-84.

[18] 王鹏,董国臣,李志国, 等. 辽西北票二道沟金矿的成矿特点和黄铁矿热电性特征[J]． 现代地质,2013,27(2):314-323.

[19] 刘云龙,申俊峰,高海军, 等. BHTE-06型热电性测试仪活化温度的标定及其意义[J]. 西部探矿工程,2011(8):146-149.

[20] 邵伟,陈光远,孙岱生．黄铁矿热电性研究方法及其在胶东金矿的应用[J]．现代地质,1990,4(1): 46-57.

[21] 刘平,邱朝霞．某些金矿床中黄铁矿热电系数值的分布特点及其意义[J]．矿物学报,1991,11(1): 60-69.

[22] 杨竹森,李红阳,高振敏, 等．胶东北部超高品位金矿黄铁矿热电性研究[J]．矿床地质,2000,19(4): 307-314.

[23] 武宗林,张寿庭,许腾, 等. 栾川矿集区黄背岭中鱼库花岗岩中黑云母成分特征及成岩成矿意义[J]. 矿物岩石, 2015, 35(3): 11-19.

[24] 段士刚,薛春纪,刘国印, 等. 河南栾川地区铅锌银矿床地质、流体包裹体、同位素地球化学与年代学[J]. 矿床地质, 2010(s1):429-430.

Thermoelectric Characteristics and Metallogenic Regularity of Pyrite in Yuku Molybdenum Deposit, Luanchuan,Henan

JIA Wenjuan1，WANG Gongwen1，HAN Jiangwei2，YAN Changhai2，SONG Yaowu2，LIU Chuanquan2，FENG Yuan3，WANG Hongwei4，YANG Fan1，KONG Liang1，GUO Nana4，ZHANG Xuhuang4，LI Junjun3

(1.School of Earth Sciences and Resources, China University of Geosciences, Beijing 100083, China；2.HenanInstituteofGeologicalSurvey,Zhengzhou,Henan450007,China；3.NorthwestInstituteofNonferrousMetalGeological,Xi’an,Shaanxi710000,China；4.LuanchuanBureauofGeologyandMineralResources,Luanchuan,Henan471500,China)

Luanchuan Yuku molybdenum deposit is located in the southern margin of the North China Craton, whose Yanshanian magma-skarn type mineralization is related to the large concealed Mo polymetallic deposit discovered in recent years. The samples from deep drilling mining (>1,000 m) were used to analyze thermoelectric characteristics of pyrite and the aim was to provide an important scientific basis for the exploration of deep prospecting in this area. The research contents and methods are as followed:the thermoelectric properties of pyrites are mainly of the N type with a total content of 96%.The thermoelectric parameters characteristics reflect the orebody denudation to the middle and lower position, but the statistical data shows the depth potential extension. According to the pyroelectricity of pyrite ore in the study area to calculate the forming temperature of fan between 270-332 ℃, the deposit belongs to shallow medium temperature hydrothermal mineralization. On the basis of geological statistics, ArcGIS Geostatistical Analyst modular was used to analyze the trend of metallogenic temperature, and the ore-forming regularity and potential targets were discussed.

pyrite; thermoelectric characteristic; metallogenic prediction； Yuku Mo deposit;Henan

2015-12-25；改回日期：2016-05-12；责任编辑：戚开静。

河南省栾川县整装勘查项目(12120114035001)；国家自然科学基金面上项目(41572318)。

贾文娟，女，硕士研究生，1992年出生，矿物学、岩石学、矿床学专业，主要从事黄铁矿与成矿规律关系研究。

Email: 1479662844@qq.com。

王功文，男，副教授，1972年出生，博士生导师，矿物学、岩石学、矿床学专业，主要从事三维成矿预测研究。

Email:gwwang@cugb.edu.cn。

P618.6；P571

A

1000-8527(2016)06-1209-10