2890 ActaPetrologicaSinica 岩石学报 2009,25(11) 1 引言东沟钼矿是我国近年来新发现的超大型钼矿床, 位于世界最大的钼矿带东秦岭钼成矿带之中 ( 图 1) 对矿石辉钼矿 Re Os 同位素定年表明, 东沟钼矿成矿年龄为 116± 2Ma( 叶会寿等,2006

Size: px

Start display at page:

Download "2890 ActaPetrologicaSinica 岩石学报 2009,25(11) 1 引言东沟钼矿是我国近年来新发现的超大型钼矿床, 位于世界最大的钼矿带东秦岭钼成矿带之中 ( 图 1) 对矿石辉钼矿 Re Os 同位素定年表明, 东沟钼矿成矿年龄为 116± 2Ma( 叶会寿等,2006"

盐毓冀
5 years ago
Views:

1 /2009/025(11) ActaPetrologicaSinica 岩石学报河南东沟钼矿花岗斑岩成因 : 岩石地球化学锆石 U Pb 年代学及 Sr Nd Hf 同位素制约戴宝章蒋少涌王孝磊 DAIBaoZhang,JIANGShaoYong,WANGXiaoLei 南京大学地球科学系内生金属矿床成矿机制研究国家重点实验室, 南京 StateKeyLaboratoryforMineralDepositsResearch,DepartmentofEarthSciences,NanjingUniversity,Nanjing210093,China 收稿, 改回. DaiBZ,JiangSY,WangXL.2009.Petrogenesisofthegraniticporphyryrelatedtothegiantmolybdenum depositin Donggou,Henanprovince,China:Constraintsfrom petrogeochemistry,zirconu PbchronologyandSr Nd Hfisotopes. ActaPetrologicaSinica,25(11): Abstract Inthisstudy,wecariedoutacomprehensiveresearchincludingzirconU Pbdating,petrogeochemistryandSr Nd Hf isotopeanalysisonthegraniticporphyryrelatedtothegiantmolybdenum depositindonggou,henanprovince,china.la ICP MS zirconu PbchronologyindicatedthattheDonggouporphyryemplacedat114~117Ma,identicalwiththetimewhenthemolybdenum depositwereformed(116±2ma,re Osage).BoththeDonggouporphyryandthelargeTaishanmiaobatholithtoitssouthare enrichedinsi,k,lilessuchasrb,thandu,andhfsessuchasnb,ta,zrandhf,whiledepletedinfe,mg,ca,srandba. TheyareweaklyperaluminousandshowhighGa/Alratios,andthuscanbeclasifiedasaluminousA typegranite.theprimitive magmaexperiencedveryhighextentofcrystalfractionationandthenemplacedtoformthedonggouporphyryinanextensionseting. TheNdisotopiccompositionsoftheDonggousamplesrangefrom to ,withcorespondingε Nd (t)valuesof-17.3to Theε Hf (t)valuesalsoshowalargerange(-3.4to-18.7)withaxenocrystyieldingε Hf (t)valueof-2.4at1715ma. Basedonaboveisotopicdata,wearguethatthemagmasourcesoftheDonggouporphyrywerecomposedofmainlyoldcrustalrocks,but mayalsoinvolvesomeportionsofmantlematerials.wesuggestedthatinputofmantle derivedmeltandhighlyevolvednatureofthe Donggouporphyryareimportantfortheformationofthegiantmolybdenumdeposit. Keywords Donggougiantporphyry typemolybdenumdeposit;aluminousa typegranite;zirconu Pbdating;Sr Nd Hfisotopes; extensionseting;mantlematerials 摘要本研究对东沟超大型钼矿床的成矿母岩东沟花岗斑岩开展了系统的年代学岩石地球化学及 Sr Nd Hf 同位素分析工作 LA ICP MS 锆石 U Pb 定年结果表明, 东沟花岗斑岩成岩年龄为 114~117Ma, 与已有的成矿年龄 (116±2Ma,Re Os 法 ) 一致, 证实了东沟钼矿为一斑岩型矿床详细的岩石地球化学分析显示东沟花岗斑岩岩体与区域上太山庙大型花岗岩基为同源演化关系, 它们均为弱过铝质, 具有富 Si 富 K 富 Rb Th U 等大离子亲石元素富 Nb Ta Zr Hf 等高场强元素, 贫 Fe Mg Ca, 贫 Sr Ba,Ga/Al 比值较高等地球化学特征, 属铝质 A 型花岗岩, 形成于伸展构造体制, 东沟岩体是母岩浆经历了强烈结晶分异高度演化的产物 ; 东沟岩体 Nd 同位素组成为 ~ ,ε Nd (t)) 值在 -17.3~-14.3 之间, 锆石的 ε Hf (t) 值变化较大, 由 -3.4 至 -18.7, 另有一颗年龄为 1715Ma 的捕获锆石的 ε Hf (t) 值为 -2.4,Nd Hf 模式年龄分别为 1.5~1.8Ga 与 1.3~1.7Ga 我们认为东沟岩体的岩浆源区以古老地壳物质为主, 但也有少量幔源组分参入, 并且幔源物质的加入及很高的岩浆演化程度可能对东沟钼矿的成岩成矿过程具有重要作用关键词中图法分类号东沟超大型斑岩钼矿 ; 铝质 A 型花岗岩 ; 锆石 U Pb 定年 ;Sr Nd Hf 同位素 ; 伸展体制 ; 幔源物质 P ;P597.3;P 科技部 973 计划课题 (2006CB403506) 资助. 第一作者简介 : 戴宝章, 男,1977 年生, 博士, 研究方向 : 矿床地球化学,E mail:daibz@nju.edu.cn 通讯作者 : 蒋少涌,E mail:shyjiang@nju.edu.cn

2 2890 ActaPetrologicaSinica 岩石学报 2009,25(11) 1 引言东沟钼矿是我国近年来新发现的超大型钼矿床, 位于世界最大的钼矿带东秦岭钼成矿带之中 ( 图 1) 对矿石辉钼矿 Re Os 同位素定年表明, 东沟钼矿成矿年龄为 116± 2Ma( 叶会寿等,2006), 与东秦岭钼矿带中其他主要钼矿床如金堆城南泥湖上房沟雷门沟等存在 20Ma 以上的成矿时差 ( 东秦岭钼矿带主成矿期为 130~150Ma, 如黄典豪等, 1994; 李永峰等,2003,2006; 郭保健等,2006; 叶会寿等, 2006; 李诺等,2007,2009), 因此东沟超大型钼矿床的发现不仅极大地增加了东秦岭钼矿带的矿产储量, 拓宽了中生代钼成矿的时代, 具有非常重大的经济价值, 而且这一新进展也极大地推动了对东秦岭地区钼成矿规律更深入的认识 ( 叶会寿等,2006; 李诺等,2007), 具有非常重要的地质意义斑岩型钼矿及斑岩夕卡岩型钼矿是东秦岭钼矿带的主要矿床类型 ( 如李永峰等,2003,2005; 郭保健等,2006; 李诺等,2007,2009), 东沟钼矿属斑岩型矿床 ( 叶会寿等,2006) 东沟钼矿成矿母岩位于下铺村东, 出露面积不足 0.01km 2 ( 图 1), 岩体 SHRIMP 锆石 U Pb 年龄为 112±1Ma( 叶会寿等,2006), 辉钼矿 Re Os 年龄为 116±2Ma, 前者略晚为更好的理解东沟钼矿成因及成岩成矿环境与动力学, 已有研究者对东沟花岗斑岩进行了部分岩石地球化学工作, 得到了该岩体具有铝质 A 型花岗岩属性的认识, 并指出东沟钼矿成岩成矿物质来源明显不同于与 I 型花岗岩相关的金堆城雷门沟南泥湖等钼矿床 ( 叶会寿等,2006) 已有研究还表明东沟花岗斑岩可能是其南侧出露的大规模太山庙岩体 (~290km 2 ) 的分枝岩脉或岩浆分异晚期的产物, 后者同样具有铝质 A 型花岗岩地球化学特征, 且成岩年龄 (115±2Ma) 与东沟花岗斑岩基本一致 ( 叶会寿等,2006,2008) 东沟超大型钼矿成岩成矿年龄虽然大体一致, 但前人报道的成岩年龄晚于成矿年龄 3~5Ma, 因此还需要进一步的年代学工作进行限定 ; 东沟花岗斑岩虽然规模很小, 但由于与成矿过程密切相关, 需要进行精细解剖, 已有的岩石地球化学数据尚无法满足这种需要, 尤其是缺乏同位素数据严重影响了对岩浆源区等问题的探讨 ; 对东秦岭钼矿带的深入认识有赖于对含矿岩体与不含矿岩体及不同时代不同成因含矿岩体的成岩过程构造环境及地球动力学的对比研究, 这一方面的工作目前仍比较薄弱本文针对以上几方面问题, 对系统采集的东沟花岗斑岩样品进行了岩石地球化学 Sr Nd 同位素 LA ICP MS 锆石 U Pb 定年及锆石原位 Hf 同位素组成研究, 进而对东沟钼矿成矿斑岩的岩石成因及相关问题进行探讨图 1 东秦岭钼矿带主要钼矿床分布图及东沟钼矿地质简图 ( 据叶会寿等,2006,2008 修编 ) Fig.1 Geologicalsketchmapshowingdistributionofmajormolybdenum depositsineasternqinlingarea,withinset showingsimplifiedgeologicalmapofthedonggoumolybdenumdeposit

3 戴宝章等 : 河南东沟钼矿花岗斑岩成因 : 岩石地球化学锆石 U Pb 年代学及 Sr Nd Hf 同位素制约地质概况与样品采集东沟钼矿位于河南省汝阳县, 地理坐标为东经 ~ , 北纬 ~ , 大地构造位置位于华北克拉通南缘的东秦岭钼矿带东段, 属外方山地块, 东沟钼矿区具体位置在区域太山庙钾长黑云母花岗岩体北东外接触带约 10km 处 ( 图 1) 目前该矿床已控制钼资源量超过 62 万吨, 属超大型规模东沟矿区出露的地层主要为熊耳群火山岩系, 为鸡蛋坪组二三四段部分层位, 主要岩性为安山岩英安岩, 及少量火山碎屑岩 ( 凝灰质粉砂岩 ); 矿区主要侵入岩有两期, 早期中酸性石英二长岩为成矿前侵入岩, 晚期花岗斑岩为钼矿床的直接成矿母岩, 它是燕山晚期太山庙岩体向北东方向侵入形成的小岩枝 ; 矿区构造以北东向为主, 是主要的导岩导矿构造东沟钼矿容矿围岩除花岗斑岩外, 主要为中元古界熊耳群火山岩, 矿体范围远大于斑岩体, 地球物理及钻探揭示的深部隐伏岩体面积可达 1.32km 2, 略呈北西南东展布 ( 付治国等,2005) 地表出露的东沟花岗斑岩呈肉红色, 块状构造, 斑状结构为主该岩体虽然规模很小, 但仔细观察后发现岩体可分为两相, 大部分岩石具有斑状结构, 为主体花岗斑岩, 而西南部局部岩石不显斑状结构, 为细粒花岗岩另岩体内局部被石英脉等穿插本研究中, 我们沿矿区内南西北东向公路所揭示的岩体新鲜露头 ( 东经 '47, 北纬 ) 系统采集了 8 件样品,DG 01~DG 06 号样品为东沟岩体样品, 其中 DG 03 样品中见很多石英脉穿插,DG 07 号样品为岩体与围岩接触部位,DG 08 为外接触带熊耳群火山岩岩相观察表明东沟岩体样品多具有斑状结构 ( 图 2A 2C), 斑晶主要为石英斜长石钾长石或条纹长石, 石英斑晶边部常见熔蚀结构 ; 基质具显微花岗结构, 主要由钾长石或条纹长石石英斜长石及少量黑云母构成, 副矿物有磁铁矿榍石金红石和锆石等样品 DG 03 中可见成岩后期石英脉 ( 图 2D), 该样品中无或少见斑晶, 主要呈细粒花岗结构, 局部见石英与钾长石构成文象结构 3 分析方法本研究主要测试工作除锆石原位 Hf 同位素分析在中科院地质与地球物理所完成, 其余 ( 主量元素微量元素全岩 Sr Nd 同位素与 LA ICP MS 锆石 U Pb 定年 ) 均在南京大学内生金属矿床成矿机制研究国家重点实验室进行新鲜全岩样品磨碎至 200 目供后续测试分析使用主要氧化物含量分析采用等离子发射光谱 (ICP AES) 方法, 将约 0.1g 全岩粉末碱熔 ( 四硼酸锂 ) 后用 10%HCl 溶样, 以溶液状态在 JY38S 型发射光谱仪上进行测试, 测试中相对标准偏差 (relativestandarddeviation) 一般 2% 二价铁含量用湿化学方法单独测定后可计算出样品中三价铁的含量微量元素测定采用酸溶法 (HF+HNO 3 ), 用 Element2 型高分辨等离子质谱 (HR ICP MS) 对样品溶液进行测定, 测试过程中用标准岩石样品对仪器状态和数据质量进行监控, 结果显示对所测的大部分元素分析精度优于 10%, 详细测试流程见高剑峰等 (2003) 为进行全岩 Sr Nd 同位素分析, 将约 50mg 全岩粉末用 HF+HNO 3 完全溶解后用 BioRadAG50W 8 阳离子树脂分离纯化 Sr Nd 元素, 详细化学分离流程参阅濮巍等 (2005) Sr Nd 同位素分析在 FinniganTritonTI 型表面热 86 电离质谱仪 (TIMS) 上进行 Sr 同位素测试过程中采用 Sr/ 88 Sr= 对质量分馏进行校正, 对 Sr 同位素国际标准物质 NISTSRM987 长期的测试值 87 Sr/ 86 Sr= ± (2σ,n=65), 与文献报道值一致 ( ± ;Weisetal.,2006);Nd 同位素测试过程中采用 Nd/ 144 Nd= 校正质量分馏,Nd 同位素国际标准物质 JNdi 1 长期的测试值 143 Nd/ 144 Nd= ± (2σ,n=67), 在误差范围内与文献报道值 ( ± ;Tanakaetal.,2000) 一致本研究中对 DG 02 和 DG 05 两个样品进行了锆石分选, 将锆石制靶及进行透射光反射光及阴极发光 (CL) 研究之后, 对选定的锆石颗粒先后进行了 LA ICP MSU Pb 定年和 Hf 同位素分析锆石 LA ICP MS 定年采用的仪器型号为 Agilent7500a, 配备 UP213 型固体激光剥蚀系统, 分析过程激光束斑直径为 32um, 频率为 5Hz 实验原理和详细测试方法可参照 Jacksonetal.(2004) 数据处理使用 GLITTER4.0 程序, 计算获得同位素比值年龄和误差, 普通铅校正采用 Andersen(2002) 的方法进行锆石 Hf 同位素使用连接 193nm 激光取样系统的 Neptune 多接收电感耦合等离子质谱仪 (MC LA ICP MS) 进行测定分析时激光束斑直径为 40 或 50um, 相应频率为 6Hz 或 8Hz, 仪器的运行条件及详细的分析流程见 Wuetal.(2006) 测试过程中标准锆石的 176 Hf/ 177 Hf 分析结果为 ± , 与推荐值 ( ±8;Goolaertsetal.,2004) 在误差范围内一致 4 分析结果 4.1 主要氧化物含量东沟斑岩体主要氧化物含量列于表 1 由主量元素含量及计算出的相关比值或参数可以看出, 该岩体具有高 Si 富 K, 而贫 Fe Mg Ca 的特征, 各样品的 A/CNK 值在 0.98~ 1.06 之间, 因此应属富钾的铝质花岗岩标准矿物计算结果表明, 该岩体岩石的主要组成矿物为石英钾长石与斜长石, 其中,DG 03 号样品中因晚期石英细脉的存在石英的比例明显升高在 TAS 图解中, 东沟岩体样品均落入花岗岩区域 ( 图 3A), 大多属弱过铝花岗岩 ( 图 3B)

4 2892 ActaPetrologicaSinica 岩石学报 2009,25(11) 图 2 东沟岩体代表性显微照片 A DG 06 号样品, 斑晶为石英, 边部略融蚀 ;B DG 01 号样品, 斑晶为斜长石 ;C DG 04 号样品, 石英集合体及斜长石构成斑晶 ; D DG 03 号样品, 见晚期石英脉, 石英与钾长石构成微文象结构 Fig.2 PhotomicrographsshowingtypicalmineralsandtexturesoftheDonggougraniticporphyry A quartzphenocrystwithslightlyerosededgesinsampledg 06;B plagioclasephenocrystinsampledg 01;C phenocrystofquartzaggregate andplagioclaseinsampledg 04;D quartzveininsampledg 03,quartzandK felsparintergrowthcanbefoundinthissample. 4.2 稀土元素及其它微量元素东沟花岗斑岩稀土元素及其它微量元素含量见表 2 可见该岩体富集 Rb K U Th 等大离子亲石元素及 Zr Hf Nb Ta 等高场强元素, 而明显亏损 Sr Ba 及 P 等元素东沟花岗斑岩样品稀土元素总量为 47~90?g/g, 低于 Whalen(1987) 计算出的 A 型花岗岩的平均值, 在球粒陨石标准化的稀土元素配分曲线 ( 图 4A) 中, 这些样品总体呈现轻重稀土分异的右倾型, 所有样品均具有强烈的 Eu 的负异常, 并且中稀土元素 (MREE) 相对轻稀土 (LREE) 及重稀土 (HREE) 均更为亏损, 因此它们的配分曲线呈现出燕翼型中间下凹两侧扬起的特征在原始地幔标准化的微量元素蛛网图 ( 图 4B) 中, 可见明显的 Sr Ba 及 P 的凹谷以及 Sm Gd Tb Ti Dy Ho Y 等元素构成的连续的低谷 4.3 全岩 Sr Nd 同位素组成本次研究中对 4 件东沟花岗斑岩样品进行了 Sr Nd 同位素分析 Sr 同位素测试值在 ~ 之间, 各样品 87 Rb/ 86 Sr 均很高, 其中 3 件样品在 100 左右 (99.6~110.0), 而另一件样品明显较低 (58.3) 根据锆石 U Pb 定年所得到的年龄 (115Ma, 见后文详述 ) 进行放射性成因 Sr 扣除后得到的初始 Sr 同位素组成均小于 ( 表 3),Wuetal(2000,2002) 曾指出, 高 Rb/Sr 比值岩石样品往往会得到明显偏低的初始 Sr 同位素组成, 这样的结果不能用以进行岩石成因的讨论此外, 由于不同样品间的 Rb/Sr 比值存在明显差异, 因此不能排除其 Rb Sr 体系在成岩后发生扰动的可能 4 件样品的 Nd 同位素组成及 147 Sm/ 144 Nd 比值则比较均一, 初始 Nd 同位素组成在 ~ 之间, 计算得到的 ε Nd (t)) 值在 -17.3~-14.3 之间

5 戴宝章等 : 河南东沟钼矿花岗斑岩成因 : 岩石地球化学锆石 U Pb 年代学及 Sr Nd Hf 同位素制约 2893 表 1 东沟岩体主要氧化物含量 (wt%) Table1 Majoroxidecomposition ofthedonggou granitic porphyry(inwt%) 样品号 07DG 07DG 07DG 07DG 07DG 07DG A ave SiO TiO Al 2 O Fe 2 O FeO MnO MgO CaO Na 2 O K 2 O P 2 O LOI SUM K 2 O/Na 2 O K 2 O+Na 2 O /CaO A/CNK FeO /MgO A.I 表 3 东沟花岗斑岩 Sr Nd 同位素组成 Table3 Sr Ndisotopiccompositionofrepresentativesamplesof thedonggouporphyry 样品号 07DG 02 07DG 04 07DG 05 07DG 06 SiO 2 (wt%) Rb(μg/g) Sr(μg/g) Rb/ 86 Sr Sr/ 86 Sr σ Sr(i) Sm(μg/g) Nd(μg/g) Sm/ 144 Nd Nd/ 144 Nd σ Nd(i) ε Nd (t)) σ T(DM) T 2 (DM) 表 2 东沟花岗斑岩微量元素含量 (μg/g) Table2 TraceelementcompositionoftheDonggougranitic porphyry(μg/g) 样品号 07DG 07DG 07DG 07DG 07DG 07DG A ave Sc Ti V Cr Zn Ga Rb Sr Y Zr Nb Mo Cs Ba La Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb Lu Hf Ta Pb Th U REE Eu/Eu (La/Yb) N Zr+Nb+ Ce+Y

6 2894 ActaPetrologicaSinica 岩石学报 2009,25(11) 图 3 A 侵入岩岩石类型 TAS 判别图解 ( 据 Middlemost, 1994);B A/NKvs.A/CNK 图解 Fig.3 A TASplotforplutonicrocks(afterMiddlemost, 1994);B A/NKvs.A/CNKplot 图 4 A 球粒陨石标准化稀土元素配分曲线 ;B 原始地幔标准化微量元素蛛网图, 标准化数据据 McDonough& Sun(1995) Fig.4 A Chondrite normalizedrareearthelementpaterns of representative samples of the Donggou porphyry; B Primitve mantle normalized trace elementpaterns of representativesamplesofthedonggouporphyry 4.4 锆石 U Pb 定年对 DG 02 和 DG 05 两件样品中锆石 U Pb 定年的结果列于表 4, 可见大部分锆石年龄较谐和, 主要集中 115Ma 左右, 仅 DG 05 号样品中一个锆石颗粒得到 1709Ma 的年龄图 5 显示了锆石的阴极发光 (CL) 照片及 U Pb 定年与 Hf 同位素分析的点位和结果, 可见大部分锆石具有较好的岩浆振荡环带在 U Pb 谐和图 ( 图 6) 中可以看到, 两件样品均有一组锆石年龄集中于 115Ma 附近, 其 206 Pb/ 238 U 年龄的加权平均结果分别是 114±1Ma 和 117±1Ma, 它们应该代表了东沟花岗斑岩体的成岩年龄 4.5 锆石 Hf 同位素组成本研究中对 DG 02 和 DG 05 两件样品中进行过 U Pb 定年的锆石相应地进行了 Hf 同位素组成测定, 结果如表 5 所示进行数据处理时, 176 Lu 的衰变常数采用 a -1 (Schereretal.,2001),ε Hf (t) 的计算利用 Blichertetal. (1997) 推荐的球粒陨石 Hf 同位素值 ( ) 及 Lu/Hf 比值 (0.0332),Hf 模式年龄计算时采用当前亏损地幔的 176 Hf/ 177 Hf 比值 ( ;Noweletal.,1998) 和 176 Lu/ 177 Lu 比值 (0.0384;Grifinetal.,2000) 计算结果表明所有锆石均具有负的 ε Hf (t) 值, 主要集中在 -10.1~ 之间 ( 图 7), 仅有两个点具有较高的 Hf 同位素组成,ε Hf (t) 值分别为 -3.4(t=117Ma) 和 -2.4(t=1715Ma) 5 讨论 5.1 成岩年龄及其意义根据东沟钼矿体与花岗斑岩小岩体的空间位置关系以及成岩成矿的同时性, 前人将东沟超大型钼矿厘定为斑岩型 ( 叶会寿等,2006,2008) 但与东秦岭钼矿带中其他斑岩型矿床以及全国数十个大中型钼矿床相比, 东沟钼矿成矿母岩出露面积最小, 虽然勘探工程已揭示深部有隐伏岩体, 但规

7 戴宝章等 : 河南东沟钼矿花岗斑岩成因 : 岩石地球化学锆石 U Pb 年代学及 Sr Nd Hf 同位素制约 2895 表 4 LA ICP MS 锆石 U Pb 定年数据 Table4 ResultofzirconU PbdatingusingLA ICP MS 点号 Th 232 (μg/g) U 238 (μg/g) Th/U 207 Pb/ 206 Pb σ 207 Pb/ 235 U σ 年龄 (Ma) 206 Pb/ 238 U σ 208 Pb/ 232 Th σ dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg 模仍较小 (1.32km 2, 付治国等,2005), 因此成岩成矿年龄的精确限定成为确认其为斑岩钼矿的一个关键依据叶会寿等 (2006) 的研究显示东沟钼矿成矿年龄为 116±2Ma, 东沟花岗斑岩成岩年龄为 112±1Ma, 二者基本同时, 但仍存在 3~5Ma 的时差本研究中, 用 LA ICP MS 方法对两件东沟花岗斑岩样品进行锆石 U Pb 定年获得的结果分别为 114± 1Ma 与 117±1Ma, 锆石内部结构特征 ( 图 5) 等表明它们代表的是岩浆结晶年龄, 即东沟花岗斑岩体成岩时代为 114~ 117Ma, 这与矿石辉钼矿 Re Os 定年结果完全一致, 从而进一步证实了成岩成矿过程的同时性 5.2 岩石属性铝质 A 型花岗岩如前所述, 东沟花岗斑岩属富钾的铝质花岗岩就微量元素地球化学特征来看, 该岩体富集 Rb K U Th 等大离子亲石元素及 Zr Hf Nb Ta 等高场强元素, 而明显亏损 Sr Ba 和 P 等元素, 其稀土元素总量为 47~90μg/g, 总体呈现轻重

8 2896 ActaPetrologicaSinica 岩石学报 2009,25(11) 表 5 锆石 Hf 同位素组成 Table5 HfisotopecompositionofzirconsfromtheDonggouporphyry 点号 t(ma) 176 Yb/ 177 Hf 176Lu/177Hf 176 Hf/ 177 Hf 2σ 176 Hf/ 177 Hf(t) ε Hf (t) T DM1 T DM2 T DM3 f Lu/Hf DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG DG 稀土分异的右倾型, 所有样品均具有强烈的 Eu 的负异常, 并且中稀土元素相对轻稀土及重稀土均更为亏损在 SiO 2 vs.ar 图解 ( 图 8A) 中, 东沟花岗斑岩样品均落入碱性岩区域 ; 在各种花岗岩构造判别图解 ( 图 8B-8D) 中, 东沟花岗斑岩均落入板内花岗岩区或其与同碰撞花岗岩交界区域东沟花岗斑岩与其南侧太山庙花岗岩基具有密切的成因联系, 可能是后者的小岩枝或分异产物 ( 付治国等, 2005; 叶会寿等,2006,2008) 在标准化稀土元素配分图及微量元素蛛网图 ( 图 4A 4B) 中均可看到东沟花岗斑岩与太山庙花岗岩的化学组成呈现出逐渐过渡演化的关系, 特别是与太山庙岩体中的碱长花岗斑岩不论在化学组成还是岩相特征上都非常相似, 因此东沟花岗斑岩小岩体与太山庙大型

9 戴宝章等 : 河南东沟钼矿花岗斑岩成因 : 岩石地球化学锆石 U Pb 年代学及 Sr Nd Hf 同位素制约 2897 图 6 东沟花岗斑岩锆石 U Pb 谐和图 Fig.6 U PbConcordiaplotsofzirconsfrom thedonggou graniticporphyry 图 5 东沟花岗斑岩锆石阴极发光 (CL) 照片图中标出了 U Pb 定年及 Hf 同位素测试的点位及相应结果 Fig.5 Cathodoluminescence(CL) imagesofselected zircongrainsfromthedonggougraniticporphyry 花岗岩基属同一侵入系列, 它们均为板内环境下形成的花岗岩东沟花岗斑岩相对于太山庙岩体偏离板内花岗岩区域是由于其具有相对更高的 Rb 含量及相对较低的 Yb 及 Y 含量, 这些特征可能是高度分异演化的结果东沟花岗斑岩强烈亏损 Eu Sr Ba, 指示着极高程度的斜长石的结晶分异 ; P 的亏损及轻稀土中稀土相对太山庙岩基的明显降低表明发生了磷灰石及其他富稀土 ( 特别是富中稀土 ) 矿物的分离如 Sr 同位素结果所显示的, 东沟花岗斑岩中 Rb 含量可能还在成岩后的地质过程中发生了变化图 7 东沟花岗斑岩锆石 Hf 同位素组成 U Pb 年龄相关关系图解 Fig.7 Hfisotopiccompositionsvs.U Pbagesofzircons fromthedonggougraniticporphyry 前述讨论表明东沟花岗斑岩具有板内碱性花岗岩的属性, 其 Ga/Al 比值均大于 2.0, 在 Ga/Al 相关的判别图解中

10 2898 ActaPetrologicaSinica 岩石学报 2009,25(11) 图 8 A SiO 2 vs.alkalinityratio(ar) 图解 ( 据 Wright,1969);B C D 花岗岩构造环境判别图解 ( 据 Pearceetal., 1984;Pearce,1996) Fig.8 A SiO 2 vs.ar(alkalinityratio)plot(afterwright,1969);b C D tectonicdiscriminantplotsforgranitoids (afterpearceetal,1984;pearce,1996) ( 图 9A 9B) 东沟花岗斑岩样品及太山庙花岗岩样品均落入 A 型花岗岩区域而在 Zr+Nb+Ce+Y 相关的判别图解 ( 图 9C 9D) 中, 东沟花岗斑岩样品及太山庙岩体碱长花岗斑岩样品落入长英质或高分异的花岗岩区域, 而太山庙岩基主体仍落入 A 型花岗岩区域 Whalenetal.(1987) 指出, 在进行 A 型花岗岩判别时, 应对时空关系密切有成因联系的岩石系列进行综合考虑, 并以最主要的岩石作为依据因此, 我们认为东沟花岗斑岩确属铝质 A 型花岗岩 ( 如叶会寿等, 2006) 利用 Eby(1992) 的图解 ( 图 10) 进一步判别可看出, 东沟花岗斑岩及太山庙花岗岩均属 A1 型图 10 中, 太山庙大型花岗岩基三个期次的花岗岩 ( 中粗粒正长花岗岩细中粒正长花岗岩和碱长花岗斑岩 ) 及东沟花岗斑岩之间呈现出明显的演化关系, 东沟岩体的演化程度明显高于太山庙晚期碱长花岗斑岩,Nb 明显富集而 Y Ce 及 Ga 含量则不同程度的降低, 可见主要富 Nb Ta 的矿物可能为金红石, 它在岩浆演化过程中没有分离出去东沟斑岩极高演化程度与钼的成矿可能存在重要联系, 对此还需进一步研究 5.3 岩浆源区及成岩成矿的联系东沟超大型钼矿与东秦岭钼矿带其他主要斑岩型钼矿存在着约 20Ma 的时差, 而与成矿密切相关的东沟花岗斑岩及太山庙大型花岗岩基具有铝质 A 型花岗岩的特点, 明显不同于其他成矿斑岩及区域上的大型花岗岩基 ( 如合峪花山等 ), 可见东沟钼矿的成岩成矿具有很强的特异性对东沟花岗斑岩而言, 其特殊性可能来源于独特的岩浆源区, 也可能是特殊的形成条件或演化过程的产物, 或者是其叠加 A 型花岗岩的概念自其提出之初, 所强调的主要是地球化学特征熔融条件和构造背景等 ( 如碱性贫水非造山或伸展高温等,Colinsetal.,1982;Whalenetal.,1987; Creaseretal.,1991;Eby,1992) 而事实上,A 型花岗岩的源区则可能非常复杂 : 可以是造山型花岗岩熔体被提取后残留长英质麻粒岩相下地壳 (Colinsetal.,1982;Whalenet al.,1987) 或下地壳脱水熔融残留的紫苏花岗岩 (Landenberger&Colins,1996), 或地壳中的英云闪长岩或花岗闪长岩等岩浆岩 (Creaseretal.,1991;Skjerlieetal.,1992,

11 戴宝章等 : 河南东沟钼矿花岗斑岩成因 : 岩石地球化学锆石 U Pb 年代学及 Sr Nd Hf 同位素制约 2899 图 9 A 型花岗岩判别图解 ( 据 Whalenetal.,1987) Fig.9 DiscriminantplotsforA typegranites(afterwhalenetal.,1987) 图 10 A1 A2 型花岗岩判别图解 ( 据 Eby,1992, 图例同图 9) Fig.10 PlotsforsubdivisionofA typegranites(afterebbyetal.,1992) 1993), 因此可以归为 I 型花岗岩的一个亚类 (Creaseretal., 1991); 可以是类似 OIB 地幔源区的陆下地幔 (Eby,1992) 或亏损地幔源区 (Hanetal.,1997), 可作为 M 型花岗岩的一类 ; 还可能是碰撞造山过程中形成的混杂的岩石圈地幔 (Whalenetal.,1996) 或通过底侵形成的新生的基性地壳物质 (Frostetal.,2001),A 型花岗岩还可以是幔源岩浆与上地

12 2900 ActaPetrologicaSinica 岩石学报 2009,25(11) 壳熔体共同分异演化的产物 (Steinitzetal.,2008) Kinget al.(1997) 指出应将铝质 A 型花岗岩单独作为一类以区别于过碱性的 A 型花岗岩, 它是壳内长英质岩石高温熔融的产物, 如实验岩石学表明钙碱性花岗岩类岩石 ( 英云闪长岩和花岗闪长岩 ) 高温低压熔融产物的特征与准铝质 A 型花岗岩非常相似 (Douce,1997) 因此,A 型花岗岩岩浆源区的确定需要综合考虑地质地球化学等各方面情况, 而同位素对准确限定 A 型花岗岩岩石成因具有至关重要的作用 ( 如 Hanet al.,1997;frostetal.,2001;mushkinetal.,2003;kempet al.,2005;steinitzetal.,2008) 东沟花岗斑岩的 Nd 同位素组成为 ε Nd (t)=-17.3~-14.3, 相应的模式年龄为 1.5Ga ~1.8Ga; 锆石的 ε Hf (t) 值主要集中在 -10.4~-18.7 之间 ( 图 7), 另有两个点 Hf 同位素组成比较特殊, 一颗为古老的继承锆石,t=1715Ma, 相应的 ε Hf (t) 值为 -2.4, 另一颗锆石具有最高的 Hf 同位素组成,t=117Ma 时的 ε Hf (t) 值为 -3.4 ( 图 7), 除开以上两个点, 东沟花岗斑岩中锆石的 Hf 模式年龄在 1.3Ga~1.7Ga 之间, 与全岩 Nd 同位素分析得到的结果基本一致或略显年轻全岩 Nd 同位素锆石 U Pb 年代学及 Hf 同位素均表明, 古老的地壳物质是东沟花岗斑岩的一个最重要岩浆源区郭波等 (2009) 的研究表明, 区域燕山期大型岩基如合峪岩体华山岩体是挤压体制向伸展体制转换过程中增厚的下地壳部分熔融的产物, 并且推断区域上存在 1.7 ~2.1Ga 的增生地壳合峪与华山岩体 ( 成岩年龄为 130~ 135Ma) 锆石 ε Hf (t) 值主要为 -16~-20, 最低可达 -30, 明显低于东沟花岗斑岩, 因此东沟岩体的岩浆源区除了与合峪华山岩体类似的古老地壳物质或其熔融残留体外, 还需要一个具有较高 Hf 同位素组成的端元, 或者说含有一定的幔源或新生地壳组分由于 130Ma 后, 东秦岭地区转变为伸展体制, 我们推断东沟岩体岩浆源区中高 Hf 同位素组成的组分可能直接来源于地幔, 幔源岩浆与地壳熔体相互作用导致了锆石中 Hf 同位素组成较大范围的变化 (ε Hf (t)=-3.4 ~-18.7), 这也可以解释东沟岩体属于与陆内幔源岩浆有关的 A1 型花岗岩的判别结果从图 7 中也可以看出, 东沟斑岩中锆石的 Hf 同位素组成偏离古老地壳物质产生的熔体 ( 以合峪华山花岗岩基为代表 ) 不太大, 表明幔源岩浆并不占据主导作用另一种可能的情况是在东沟斑岩形成之前发生过幔源岩浆的底侵作用而有部分新生地壳形成, 后来在伸展条件下新生地壳与古老地壳物质共同熔融产生了东沟斑岩的母岩浆综上所述,115Ma 左右时由于区域伸展构造体制的增强, 地壳岩石在增温减压的双重效用下也发生熔融, 其中包含一定量的新生地壳或者幔源物质, 生成的熔体最终侵位形成了太山庙大型花岗岩基, 高度分异演化的部分熔体稍晚侵位形成岩枝状产出的具有铝质 A 型花岗岩地球化学特征的东沟花岗斑岩幔源组分的加入及充分的分异演化对东沟钼矿的成岩成矿具有重要作用 5 结论 LA ICP MS 锆石 U Pb 定年结果表明东沟花岗斑岩成岩年龄为 114~117Ma, 与东沟超大型钼矿成矿年龄 (116± 2Ma) 完全一致, 成岩成矿的同步性有力的支持了东沟钼矿属斑岩型矿床的认识东沟斑岩与太山庙花岗岩为同源演化的岩石系列, 它们均为铝质 A 型花岗岩, 东沟花岗斑岩是母岩浆高度分异演化的产物全岩 Nd 同位素与锆石 Hf 同位素研究的结果表明, 东沟斑岩的岩浆源区以古老地壳物质为主外, 还可能有一定量的幔源物质参入, 岩浆源区的独特性及非常高的演化程度可能对东沟钼矿的形成具有重要作用 References AndersenT.2002.CorectionsofcommonleadinU Pbanalysesthatdo notreport 204 Pb.Chem.Geol.,192:59-79 BlichertTJ,ChauvelC,AlbaredeF.1997.SeparationofHfandLufor high precisionisotopeanalysisofrocksamplesbymagneticsector multiplecolectoricp MS.Contrib.Mineral.Petrol.,127: ColinsWJ,BeamsSD,WhiteAJRetal.1982.NatureandoriginofA typegraniteswith particularreferencetosoutheastern Australia. ContribMineralPetrol,80: CreaserRA, Price RC, and Wormald RJ A type granites revisited:asesementofaresidual sourcemodel.geology,19: DouceAEP.1997.GenerationofmetaluminousA typegranitesbylow presuremeltingofcalc alkalinegranitoids.geology,25(8): Eby GN Chemical subdivision of the A type granitoids: Petrogeneticandtectonicimplications.Geology,20: FrostCD,BelJM,FrostBRetal.2001.Crustalgrowthbymagmatic underplating:isotopicevidencefromthenorthernshermanbatholith. Geology,29(6): FuZG,LvQW,TianXQetal.2005.Studyongeologicalcharacteristics andprospectingfactorsofdonggoumodeposit.chinamolybdenum Industry,29(2):8-16 GaoJF,Lu J,LaiMYetal.2003.Analysisoftraceelementsinrock samplesusinghr ICPMS.Jour.NanjingUniv.(NaturalSci.),39 (6): (inChinesewithEnglishabstract) GoolaertsA,MatieliN,DejongJetal.2004.HfandLuisotopic referencevaluesforzirconstandard91500bymc ICP MS.Chem. Geol.,206:1-9 GrifinWL,PearsonNJ,BelousovaE etal.2000.thehfisotope compositionofcratonicmantle:lam MC ICPMSanalysisofzircon megacrystsinkimberlites.geochim.cosmochim.acta,64: GuoB,ZhuLM,LiBetal.2009.ZirconU PbageandHfisotope compositionofthehuashanandheyugraniteplutonsatthesouthern marginofnorthchinacraton:implicationsforgeodynamicseting. ActaPetrologicaSinica,25(2): (inchinesewith Englishabstract) GuoBJ,MaoJW,LiHW etal.2006.re Osdatingofthemolybdenite fromtheqiushuwancu ModepositintheeastQinlingandits geologicalsignificance.actapetrologicasinica,22(9): (inchinesewithenglishabstract)

13 戴宝章等 : 河南东沟钼矿花岗斑岩成因 : 岩石地球化学锆石 U Pb 年代学及 Sr Nd Hf 同位素制约 2901 HanBF,WangSG,JahnBMetal.1997.Depleted mantlesourceforthe UlungurRiver A type granites from North Xinjiang, China: geochemistryand Nd Srisotopicevidence, and implicationsfor Phanerozoiccrustalgrowth.Chem.Geol.,138: HuangDH, W CY, Du AD etal Re Osisotopicagesof molybdenumdepositsineastqinlingandtheirsignificance.mineral Deposits,13(3): (inChinesewithEnglishabstract) JacksonSE,PearsonNJ,GrifinWLetal.2004.Theapplicationoflaser ablationmicroprobe inductivelycoupled plasma masspectrometry (LAM ICP MS) toin situ U Pb zircon geochronology. Chem. Geol.,211:47-69 KempAIS,WormaldRJ,WhitehouseMJetal.2005.Hfisotopesin zirconrevealcontrasting sourcesand crystalization historiesfor alkalinetoperalkalinegranitesoftemora,southeasternaustralia. Geology,33(10): KingPL,WhiteAJR,ChappelAW etal.1997.characterizationand originofaluminousa typegranitesfrom thelachlanfoldbelt, SoutheasternAustrilia.J.Petrol.,38(3): LiN,ChenYJ,ZhangH etal.2007.molybdenum depositsineast Qinling.EarthScienceFrontier,14(5): (inChinese withenglishabstract) LiN,ChenYJ,SunYLetal.2009.MolybdeniteRe Osisochronageof theyuchiling porphyry Mo deposit, Henan Province and its geologicalimplications.actapetrologicasinica,25(2): (inchinesewithenglishabstract) LiYF, Mao JW, BaiFJetal Re Osisotopic dating of molybdenitesinthenannihumolybdenum(tungsten)orefieldinthe EasternQinlinganditsgeologicalsignificance.GeologicalReview, 49(6): (inChinesewithEnglishabstract) LiYF,MaoJW,LiuDY etal SHRIMP zircon U Pb and molybdenitere OsdatingsfortheLeimengouprophyrymolybdenum deposit,westernhenananditsgeologicalimplication.geological Review,52(1): (inChinesewithEnglishabstract) McDonoughWF,andSunSS,1995.ThecompositionoftheEarth. Chem.Geol.,120: MiddlemostEAK,1994.Namingmaterialsinthemagma/igneousrock system.earthsci.rev.,37(3 4): MushkinA,NavonO,HaliczLetal.2003.ThepetrogenesisofA type magmasfrom theamdam Masif,southernIsrael.J.Petrol.,44 (5): NowelGM,KemptonPD,NobleSR etal.1998.highprecisionhf isotopemeasurementsofmorbandoibbythermalionizationmas spectrometry:insightsintothedepletedmantle.chem.geol.,149: PearceJA, HarisN BW, and Tindle AG Trace element discriminationdiagramsforthetectonicinterpretation ofgranitic rocks.j.petrol.,25(4): PearceJA.1996.Sourcesandsetingsofgraniticrocks.Episodes,19 (4): PuW,GaoJF,ZhaoKDetal.2005.SeperationmethodofRb Sr,Sm NdusingDCTAandHIBA.Jour.NanjingUniv.(NaturalSci.), 41(4): (inChinesewithEnglishabstract) SchererE,MunkerC,MezgerK.,2001.Calibrationofthelutetium hafniumclock.science,293: SkjerlieKP,andJohnstonAD,1992.Vapor absentmeltingat10kbarof abiotite andamphibole bearingtonaliticgneis:implicationsforthe generationofa typegranites.geology,20: SkjerlieKP,andJohnstonAD.1993.Fluid absentmeltingbehaviorofan F richtonaliticgneisatmid crustalpresures:implicationsforthe generationofanorogenicgranites.j.petrol.,34(4): SteinitzA,KatzirY,ValeyJW etal.2008.theorigin,coolingand alteration ofa type granites in southern Israel(norhternmost Arabian Nubian shield): amulti mineraloxygen isotopestudy. Geol.Mag.,146(2): TanakaT,TogashiS,KamiokaH etal.2000.jndi 1:aneodymium isotopicreferenceinconsistencywithlajolaneodymium.chem. Geol.,168: WeisD,KieferB,MaerschalkCetal.2006.High precisionisotopic characterizationofusgsreferencematerialsbytimsandmc ICP MS.Geochem.Geophys.Geosyst.,7,Q08006,doi: / 2006GC WhalenJB, CurieKL, and ChappelBW,1987. Atypegranites: geochemical characteristics, discrimination, and petrogenesis. ContributionstoMineralogyandPetrology,95: WrightJB,1969.Asimplealkalinityratioanditsapplicationtoquestions ofnon orogenicgranitegenesis.geologicalmagazine,106: WhalenJB,JennerGA,LongstafeFJetal.1996.Geochemicaland isotopic(c, Nd, Pb and Sr) constraints on A type granite petrogenesisbasedonthetopsailsigneoussuite, Newfoundland Appalachians.J.Petrol.,37(6): WuFY,JahnBM,WildeSetal.2000.Phanerozoiccontinentalcrustal growth:sr Ndisotopicevidencefromthegranitesinnortheastern China.Tectonophysics,328, WuFY,SunDY,LiHM etal.2002.a typegranitesinnortheastern China: age and geochemicalconstraintson theirpetrogenesis. Chem.Geol.,187, WuFY,YangJH,XieLW etal.2006.hfisotopiccompositionsofthe standardzirconsand baddeleyitesused in U Pb geochronology. Chem.Geol.,234: YeHS,MaoJW, LiYF etal SHRIMP Zircon U Pb and molybdenitere OsdatingforthesuperlargeDonggouporphyryMo depositineastqinling,china,anditsgeologicalimplication. ActaGeologicaSinica,80(7): (inChinesewith Englishabstract) YeHS,MaoJW,XuLGetal.2008.SHRIMPzirconU Pbdatingand geochemistryofthetaishanmiaoaluminousa typegraniteinwestern Henanprovince.GeologicalReview,54(5): 附中文参考文献付治国, 吕伟庆, 田修启等东沟钼矿矿床地质特征及找矿因素研究. 中国钼业,29(2):8-16 高剑锋, 陆建军, 赖鸣远等岩石样品中微量元素的高分辨等离子质谱分析. 南京大学学报 ( 自然科学 ),39(6): 郭保健, 毛景文, 李厚文等秦岭造山带秋树湾铜钼矿床辉钼矿 Re Os 定年及其地质意义. 岩石学报,22(9): 郭波, 朱赖民, 李等华北陆块南缘华山和合峪花岗岩岩体锆石 U Pb 年龄 Hf 同位素组成与成岩动力学背景. 岩石学报, 25(2): 黄典豪, 吴澄宇, 杜安道等东秦岭地区钼矿床的铼锇同位素年龄及其意义. 矿床地质,13(3): 李诺, 陈衍景, 孙亚莉等河南鱼池岭钼矿床辉钼矿铼锇同位素年龄及其地质意义. 岩石学报,25(2): 李诺, 陈衍景, 张辉等东秦岭斑岩钼矿带的地质特征和成矿构造背景. 地学前缘,14(5): 李永峰, 毛景文, 白凤军等东秦岭南泥湖钼 ( 钨 ) 矿田 Re Os 同位素年龄及其地质意义. 地质论评,49(6): 李永峰, 毛景文, 刘敦一等豫西雷门沟斑岩钼矿 SHRIMP 锆石 U Pb 和辉钼矿 Re Os 测年及其地质意义. 地质论评,52(1): 濮巍, 高剑锋, 赵葵东等利用 DCTA 和 HIBA 快速有效分离 Rb Sr Sm Nd 的方法. 南京大学学报 ( 自然科学 ),41(4): 叶会寿, 毛景文, 李永峰等东秦岭东沟超大型斑岩钼矿 SHRIMP 锆石 U Pb 和辉钼矿 Re Os 年龄及其地质意义. 地质学报,80(7): 叶会寿, 毛景文, 徐林刚等豫西太山庙铝质 A 型花岗岩 SHRIMP 锆石 U Pb 年龄及其地球化学特征. 地质论评,54(5):

, 3.1 PM303, 10(2)( 5, 5 ), 1 w (SiO 2 ) 68% ~89%, 2( 1, 1) (Al 2O 3,Fe 2O 3,MnO) (MgO,K 2O, Na 2O),w(TiO 2 ) w (SiO 2 ) ;,w (TiO 2 ) 0.31%

, 3.1 PM303, 10(2)( 5, 5 ), 1 w (SiO 2 ) 68% ~89%, 2( 1, 1) (Al 2O 3,Fe 2O 3,MnO) (MgO,K 2O, Na 2O),w(TiO 2 ) w (SiO 2 ) ;,w (TiO 2 ) 0.31% 36 3 2017 5 GeologicalScienceandTechnologyInformation Vol.36 No.3 May 2017 doi:10.19509/j.cnki.dzkq.2017.0302,,,. [J].,2017,36(3): 7-17. 1, 1, 1, 2, 1 3 (1., 264004;2. ( ), 430074; 3., 710000) :, (Al 2O

2890 ActaPetrologicaSinica 岩石学报 2009,25(11) 1 引言 东沟钼矿是我国近年来新发现的超大型钼矿床, 位于世界最大的钼矿带 东秦岭钼成矿带之中 ( 图 1) 对矿石辉钼矿 Re Os 同位素定年表明, 东沟钼矿成矿年龄为 116± 2Ma( 叶会寿等,2006

2890 ActaPetrologicaSinica 岩石学报 2009,25(11) 1 引言东沟钼矿是我国近年来新发现的超大型钼矿床, 位于世界最大的钼矿带东秦岭钼成矿带之中 ( 图 1) 对矿石辉钼矿 Re Os 同位素定年表明, 东沟钼矿成矿年龄为 116± 2Ma( 叶会寿等,2006