岩石学报  2018, Vol. 34 Issue (8): 2380-2392   PDF    
新疆东昆仑达拉库岸镁铁-超镁铁质岩体年代学、地球化学及成矿条件
夏明哲1,2,3 , 范亚洲4 , 夏昭德1,2 , 芮会超1 , 姜常义1,2     
1. 长安大学地球科学与国土资源学院, 西安 710054;
2. 西部矿产资源与地质工程教育部重点实验室, 西安 710054;
3. 国土资源部岩浆作用成矿与找矿重点实验室, 西安 71005;
4. 中国地质科学院地质研究所, 北京 100037
摘要:达拉库岸镁铁-超镁铁质岩体位于东昆仑造山带南带之喀拉米兰晚古生代沟弧系,主要由二辉橄榄岩、单辉橄榄岩、橄榄二辉岩、单辉辉石岩和辉长岩组成。辉长岩的锆石U-Pb谐和年龄为244.4±1.5Ma,属于中三叠世。岩石普遍弱富集稀土元素和大离子亲石元素(Rb、Ba、Sr),亏损不相容元素(Nb、Ta)。岩体原生岩浆为高镁拉斑玄武质岩浆(MgO 11.4%,FeO 10.8%)。元素地球化学和Nd、Sr同位素组成特征表明岩浆源区为富集岩石圈地幔,岩浆运移和侵位过程中遭受不同程度地壳物质的同化混染作用。依据橄榄石组分模拟获得橄榄石结晶过程中母岩浆达到硫饱和,并发生硫化物的熔离作用。从岩体特征、矿石结构、原生岩浆性质、深部硫化物熔离和物探信息等方面综合分析,岩体具有形成岩浆型铜镍硫化物矿床的良好条件。
关键词: 成矿条件     地球化学     年代学     达拉库岸岩体     东昆仑    
Geochronology, geochemical characteristics and ore-forming conditions of the Dalakuan mafic-ultramafic intrusion, East Kunlun, Xinjiang
XIA MingZhe1,2,3, FAN YaZhou4, XIA ZhaoDe1,2, RUI HuiChao1, JIANG ChangYi1,2     
1. College of Earth Science and Land Recourses, Chang'an University, Xi'an 710054, China;
2. Key Laboratory of western China's Mineral Resources and Geological Engineering, Ministry of Education, Xi'an 710054, China;
3. Key Laboratory for the Study of Focused Magmatism and Giant Ore Deposits, Ministry of Land and Resources, Xi'an 71005;
4. Institute of Geology, Chinese Academy of Geological Sciences, Beijing 100037, China
Abstract: The Dalaku'an mafic-ultramafic intrusion is located in the Late Paleozoic Karamilan arc zone, the southern active zone of the East Kunlun. It consists of lherzolite, wehrlite, olivine websterite, clinopyroxenite and gabbro. U-Pb zircon analysis yield an age of 244.4±1.5Ma from gabbro, indicates the intrusion emplaced in Anisian of Middle Triassic. Most of the samples are slightly enriched in LREE, moderately depleted in HFSE (Nb, Ta). The estimated primary magma for ultramafic rocks is tholeiitic and contains 11.4% MgO and 10.8% FeO. The characteristics of trace elements geochemical and Nd-Sr isotope indicate the magmatic source is enriched lithosphere mantle, and contamination at different levels is involved as magma migration and emplacement. The sulfide segregation and olivine fractional crystallization occurred almost at the same time during the magma evolution. Comprehensive analysis based on detailed fieldwork, ore texture/structures and sulfide segregation in the deep demonstrates the intrusion is conductive to form magmatic Cu-Ni sulfide deposits.
Key words: Ore-forming conditions     Geochemistry     Geochronology     Dalaku'an mafic-ultramafic intrusion     East Kunlun    

东昆仑造山带是中央造山带西段重要的组成部分,也是青藏高原北缘一条巨型岩浆岩带和我国重要的多金属成矿带(莫宣学等, 2007)。近年来,我国地质工作者在青海东昆仑发现了与镁铁-超镁铁质岩石密切相关的夏日哈木超大型镍矿床,探明金属镍储量超过100万吨,成为了我国岩浆矿床找矿工作的重大突破(李世金等, 2012)。此外,在该构造带陆续发现了石头坑德、冰沟南等较好的铜镍矿化岩体(周伟等, 2016; 董俊等, 2017; 张照伟等, 2017),使东昆仑地区成为了寻找岩浆型铜镍硫化物矿床的有利地区。

新疆东昆仑位于阿尔金断裂之南,北邻塔里木板块、柴达木陆块、南邻羌塘-巴颜喀拉构造带(姜春发, 1992; 新疆维吾尔自治区地质矿产局, 1993)。1:25万区域地质调查中发现该地区存在诸多镁铁-超镁铁质岩体,前人多将这些岩体归为蛇绿岩套的组成部分,认为这些岩体具有蛇绿岩属性(何国琦, 1994; 赖绍聪等, 1996; 王焰等, 1999)。然而,2012年至今,在该区域相继发现了一些与围岩呈显著侵入接触关系的镁铁-超镁铁质岩体,例如达拉库岸、几克里阔勒等,其中最具有代表性是达拉库岸岩体,岩石种类丰富,地表矿化特征显著,橄榄岩相中赋存岩浆型铜镍硫化物矿化体,单件样品的Ni最高品位为0.65%,Cu最高品位1.4%(新疆百川新源矿业开发有限责任公司, 2012)。本文以达拉库岸岩体为研究对象,通过详细岩石学、矿物晶体化学、年代学、地球化学的研究,明确岩体形成时限、探讨岩浆演化过程,综合分析岩体成矿潜力和找矿方向。

① 新疆百川新源矿业开发有限责任公司. 2012.新疆且末县达拉库岸铜多金属矿勘查工作报告

1 岩体地质概况

达拉库岸岩体位于东昆仑造山带南带之喀拉米兰晚古生代沟弧系(弓小平等, 2004; 韩红卫等, 2007)(图 1a)。区域内出露地层有:中泥盆世布拉克巴什组、早泥盆世托库孜达坂群、早-中石炭世喀拉米兰河群、二叠纪喀尔瓦组和第四系上更新统。区内岩浆岩较发育,主要为华力西中期的花岗岩、二长花岗岩;华力西晚期的黑云母花岗岩、二长花岗岩、花岗闪长岩、闪长岩、辉长岩、超基性岩。

图 1 达拉库岸岩体岩体构造位置(a据姜春发,1992)及地质略图(b, c) Fig. 1 Tectonic location of Dalakuan intrusion (a, after Jiang, 1992) and the geological sketch map (b, c)

达拉库岸岩体有东、西两个镁铁-超镁铁质小侵入体组成(图 1b)。其中,西侧侵入体出露长约90m,宽约8m,面积小于0.1km2,局部可见其侵入华力西期花岗闪长岩中,界线清晰;东侧侵入体呈不规则透镜状产出,长约96m,宽9~12m,出露面积约0.1km2,一侧侵入到中细粒结构花岗闪长岩中,另一侧与花岗闪长岩呈断层接触,断层带岩石破碎(图 1c)。此外,在岩体外围分布有二十余个辉长岩体,面积多小于0.1km2。辉长岩较新鲜,蚀变弱,以中粗粒结构、辉长结构为主,块状构造。

达拉库岸岩体地表出露岩石以超镁铁质岩石为主,东侧小侵入体主要岩石类型有橄榄二辉岩、单辉橄榄岩和辉长岩;西侧小侵入体主要为二辉橄榄岩、橄榄二辉岩和单辉辉石岩。其中,橄榄二辉岩和单辉辉石岩与二辉橄榄岩、单辉橄榄岩呈渐变接触,辉长岩多呈透镜状分布在超镁铁质岩石内部或在岩体边部不连续分布。岩石中常见结构有自形-半自形粒状结构、堆晶结构、包含结构、辉长结构等典型镁铁-超镁铁质岩石结构特征, 以块状构造为主(图 2)。岩石普遍遭受了较强蚀变, 有蛇纹石化、透闪石化、阳起石化、绿泥石化和钠黝帘石化。此外,各类岩石中含有5%~7%的普通角闪石,表明岩浆源区富水。

图 2 达拉库岸岩体典型岩石结构特征 (a)二辉橄榄岩中橄榄石堆晶,单斜辉石呈填隙状分布;(b)橄榄二辉岩中的大颗粒斜方辉石中包裹多个橄榄石,呈典型包橄结构. Ol-橄榄石;Opx-斜方辉石;Cpx-单斜辉石 Fig. 2 Microphotographs showing the textures of representative rock of the Dalakuan intrusion (a) all olivine and occur as cumulate phase, while clinopyroxene are intercumulate in the lherzolite; (b) olivine in websterite commonly occurs as inclusions in large orthopyroxene crystals. Ol-olivine; Opx-orthopyroxene; Cpx-clinopyroxene

达拉库岸岩体地表矿化特征明显,岩石表面强烈孔雀石化、局部有铜蓝,围岩花岗闪长岩裂隙中也发育孔雀石化(图 3a)。侵入体橄榄岩相中赋存铜镍硫化物矿化体,硫化物主要以磁黄铁矿、镍黄铁矿、黄铁矿和黄铜矿为主(图 3b, c),多呈不规则状集合体填充在橄榄石/辉石晶间空隙中,局部见有海绵陨铁结构。在橄榄石/辉石中可见包裹有乳滴状、浑圆状硫化物(图 3d-f)。

图 3 达拉库岸岩体地表矿化及原生硫化物特征 (a)橄榄二辉岩强烈孔雀石化;(b)二辉橄榄岩中硫化物局部呈稠密浸染状;(c)磁黄铁矿出溶链状镍黄铁矿;(d-f)橄榄石/辉石中呈乳滴状分布的硫化物.Sul-硫化物;Pn-镍黄铁矿;Po-磁黄铁矿;Cpy-黄铜矿 Fig. 3 The characteristics of sulfide in the intrusion (a) sulfide mineralization in olivine websterite; (b) interstitial sulfide assemblages in lherzolite; (c) pyrrhotite exsolute catenary pentlandite; (d-f) sulfide droplets in olivine or pyroxene minerals. Sul-sulfide; Pn-pentlandite; Po-pyrrhotite; Cpy-chalcopyrite
2 测试分析方法

锆石U-Pb年代学样品采自岩体内辉长岩,岩石新鲜,块状构造,辉长结构,无矿化。样品在河北省廊坊市区域地质调查研究所完成锆石挑选工作。在双目镜下选择颗粒大、晶型好、透明的锆石进行制靶。锆石测试工作在西北大学大陆动力学国家重点实验室完成。锆石LA-ICP-MS U-Pb年代学分析采用带有Shield Torch的Agilent 7500a的ICP-MS和德国Micro-Las公司生产的Geo-Las200M激光剥蚀系统。测试使用的标准锆石是91500和GJ-1,激光束斑直径为30μm,激光剥蚀样品的深度为20~40μm。

造岩矿物组成在长安大学西部矿产资源与地质工程教育部重点实验室采用JXI-8100型电子探针分析,加速电压15kV,束电流1.0×10-8 A,束斑直径1μm。主量元素和微量元素分析在长安大学西部矿产资源与地质工程教育部重点实验室采用3080E型X荧光光谱仪和美国X-7型ICP-MS测定。XRF熔片法按国家标准GB/T 14506.28—1993。Nd、Sr同位素比值测试在中国科学院广州地球化学研究所完成。Nd、Sr同位素分析仪器为Micromass Isoprobe型多接收器等离子体质谱仪。Sr、Nd的全流程空白分别为20pg、10pg,对样品的贡献可以忽略。该仪器测量标样La Jolla的143Nd/144Nd比值为0.511861±10,测量标样SRM987的87Sr/86Sr比值是0.710263±10。

3 锆石U-Pb年代学

选择岩体中新鲜辉长岩进行锆石U-Pb年代学研究,共测数据点15个,数据见表 1。辉长岩中锆石宽板状,内部较纯洁,无裂纹和包裹体,岩浆振荡环带较发育,U含量介于80×10-6~2408×10-6,Th含量介于127×10-6~5677×10-6,Th/U比值在1.6~2.4,具有典型岩浆成因特点(吴元保和郑永飞, 2004)。岩体辉长岩的锆石206Pb/238U-207Pb/235U谐和年龄为244.4±1.5Ma,206Pb/238U加权平均年龄为244.4±1.5Ma(图 4),由此可见,达拉库岸岩体的形成时代为244.4±1.5Ma,属于中三叠世。

表 1 达拉库岸岩体辉长岩锆石定年分析结果 Table 1 U-Pb zircon data for Dalakuan intrusion

图 4 达拉库岸岩体辉长岩锆石U-Pb谐和年龄 Fig. 4 The zircon U-Pb concordia diagrams of gabbro in Dalakuan intrusion
4 岩相学及矿物晶体化学特征

岩体中不同岩相橄榄石和辉石电子探针分析数据见表 2表 3。橄榄石成分变化范围较小,Fo值介于82.17~88.86,为贵橄榄石。斜方辉石的En为81~81.99;Fs为16.14~16.88;Wo为1.83~2.12,均为古铜辉石。单斜辉石的En为45.94~48.75;Fs为2.42~8.2;Wo为44.16~50.44,均属Ca-Mg-Fe辉石族,在Wo-En-Fs图中,主要位于透辉石区(图略)。

表 2 达拉库岸岩体橄榄石的电子探针分析数据(wt%) Table 2 The crystal chemical data of olivine of the Dalakuan intrusion (wt%)

表 3 达拉库岸岩体辉石的电子探针分析数据(wt%) Table 3 The crystal chemical data of pyroxenes (wt%) of the Dalakuan intrusion
5 地球化学特征 5.1 主量元素地球化学

达拉库岸岩体主量元素分析数据见表 4。各类岩石组成与岩相学特征变化相对应,Al2O3、Fe2O3T、MgO、CaO、Na2O含量变化大,MgO含量普遍高于Fe2O3T含量。岩体中所有岩石TiO2含量均 < 1.0%,平均0.27%。所有岩石的m/f值介于3.20~5.14之间,为铁质超基性岩。在Na2O+K2O-SiO2图解中,全部样品都投在了亚碱性区域(图 5a)。在FeOT/MgO-SiO2图解中,除有4件岩石样品投影在了钙碱性系列区,大部分样品均投影于拉斑玄武岩区(图 5b)。

图 5 Na2O+K2O-SiO2图解(a)和SiO2-FeOT/MgO图解(b)(据Miyashiro, 1974) Fig. 5 Na2O+K2O vs. SiO2 diagram (a) and SiO2 vs. FeOT/MgO diagram (b) (after Miyashiro, 1974)

表 4 达拉库岸岩体主量元素(wt%)和微量元素(×10-6)分析数据 Table 4 Major elements (wt%) and trace elements (×10-6) of the Dalakuan intrusion
5.2 微量元素地球化学

岩体中超镁铁质岩石稀土元素总量变化较大,∑REE=9.89×10-6~62.57×10-6,平均为16.10×10-6,2件辉长岩稀土元素总量相对高,∑REE=35.60×10-6~36.90×10-6。岩石的(La/Sm)N值1.08~4.26,(Gd/Yb)N值1.39~2.67,(La/Yb)N值1.85~5.42之间,表明轻、重稀土元素内部的分馏程度较小,而重轻稀土元素之间分馏程度较大。在稀土元素的配分图上,各类岩石具有相似的球粒陨石标准化的配分型式,均表现为属轻稀土弱富集型(图 6a)。

图 6 达拉库岸岩体球粒陨石标准化稀土元素配分曲线图(a, 标准化值据Taylor and Mclennan, 1985)和原始地幔标准化多元素配分曲线图(b, 标准化值据Sun and McDonough, 1989) Fig. 6 Chondrite-normalized REE patterns (a, normalizing values after Taylor and Mclennan, 1985) and PM-normalized trace elements spider diagram (b, normalizing values after Sun and McDonough, 1989) of the Dalakuan intrusion

微量元素蛛网图型式上,岩体中除部分样品的Sr和U外,大部分样品的不相容元素、适度不相容元素与原始地幔相应元素比值均不超过10。总体来看,岩石弱富集大离子亲石元素,有不同程度的Nb、Ta亏损和弱Ba亏损。部分样品大离子亲石元素(Rb、Ba、Sr、Th、U)的不规则波动,可能主要是受到热液蚀变作用的影响(图 6b)。

5.3 Nd、Sr同位素地球化学

岩体Nd、Sr同位素分析数据见表 5。岩体Nd同位素组成相对均一,143Nd/144Nd为0.512408~0.512603,而Sr同位素变化稍大,87Sr/86Sr为0.7056~0.7076,εSr(t)为15.8~36.4。采用岩体形成时代244Ma计算,获得二辉橄榄岩的εNd(t)为+0.15,具有弱亏损型地幔特征;其余样品εNd(t)在-4.02~-0.23,显示弱富集型特征。在εNd(t)-(87Sr/86Sr)i相关图上,样品位于第一象限和第四象限中,具有明显的EMⅡ变化趋势(图 7)。

图 7 达拉库岸岩体εNd(t)-(87Sr/86Sr)i相关图(据Zindler and Hart, 1986) Fig. 7 εNd(t) versus initial 87Sr/86Sr of the Dalakuan intrusion (after Zindler and Hart, 1986)

表 5 达拉库岸岩体Nd-Sr同位素组成 Table 5 The isotopic compositions of Nd and Sr of the Dalakuan intrusion
6 讨论 6.1 原生岩浆

与镁铁-超镁铁质岩浆有关的铜镍硫化物矿床中原生岩浆是成矿作用的关键(Naldrett, 1999)。由于在野外没有找到岩体的冷凝边,所以无法直接获得岩体原生岩浆的成分。单斜辉石的主量元素成分可以很好地反映原生岩浆特征(邱家骧和廖群安, 1987)。达拉库岸岩体中单斜辉石具有低Al2O3、TiO2、Na2O和K2O的特点,单斜辉石的SiO2-Al2O3和SiO2-Na2O-TiO2图解显示原生岩浆为大陆拉斑玄武质岩浆(图 8)。依据橄榄石-熔体平衡原理可以定量估算岩体原生岩浆的组成。达拉库岸岩体中二辉橄榄岩中橄榄石具有最高的Fo=88.86,选择不含硫化物且MgO含量高、CaO含量相对低的样品,利用Li and Ripley (2010)介绍的原生岩浆模拟计算方法,选取Mg-Fe在橄榄石-熔体之间的分配系数KD=X(FeO/MgO)Ol/X(FeO/MgO)melt=0.3±0.03(Roeder and Emslie, 1970)。获得岩体原生岩浆MgO为11.4%、FeO为10.8%。由此可见,达拉库岸镁铁-超镁铁质岩体的原生岩浆为高镁拉斑玄武质岩浆。

图 8 单斜辉石Al2O3-SiO2图解(a)和SiO2-Na2O-TiO2图解(b)(邱家骧和廖群安, 1987) Fig. 8 Al2O3 vs. SiO2 diagram (a) and SiO2-Na2O-TiO2 diagram (b) of clinopyroxene from the Dalakuan intrusion (after Qiu and Liao, 1987)
6.2 同化混染作用

达拉库岸岩体地表覆盖严重,出露面积小,很难观察到围岩捕虏体等显示岩浆与围岩之间可能有物质交换的直接证据。通常,幔源岩浆上升侵位过程中会不同程度受到陆壳物质的混染。从元素地球化学角度来看,总分配系数相同或很相近的元素比值在岩浆结晶过程中不会改变。因此,根据总分配系数相同或很相近、对同化混染作用又敏感的元素比值(例如Ce/Pb、Th/Yb、Nb/Ta、Ta/Yb、K2O/P2O5、Ti/Yb、Zr/Nb等)之间的协变关系,可以检验是否存在同化混染作用,并判断混染程度(Mecdonald et al., 2001; Barker et al., 1997; Campbell and Griffiths, 1992, 1993)。在图 9中Zr-Th、La/Yb-Ce/Yb和Ta/Yb-Th/Zr较好的相关性证实岩浆演化过程中存在陆壳物质的混染作用。此外,达拉库岸岩体中Sr同位素均显示了富集型特征,而橄榄岩相Nd同位素组成有弱亏损型特征,并且从橄榄岩相-辉石岩相-辉长岩相εNd(t)逐渐变小(-4.02~+0.15),εSr(t)逐渐变大(15.8~36.4),表明岩浆演化的早期阶段同化混染作用相对较弱,而演化到晚期阶段同化混染作用逐渐增强。大部分样品岩石化学系列属于拉斑玄武岩,少量的辉石岩和辉长岩显示了钙碱性系列特征,这可能也是同化混染作用所致。因为岩浆演化过程中同化混染作用达到一定程度时,会促使岩浆由拉斑玄武质系列转化为钙碱性系列,并导致岩浆演化方向由Fenner趋势(铁富集)转化为Bowen趋势(硅碱富集)。

图 9 达拉库岸岩体同化混染作用判别图 Fig. 9 Plots of selected trace elements for checking contamination of the Dalakuan intrusion
6.3 岩浆源区

一般来说,软流圈地幔和岩石圈地幔由于演化历史的差异存在系统区别。大陆岩石圈地幔,尤其是古老的大陆岩石圈地幔往往具有低的Sm/Nd、U/Pb和高Rb/Sr比值。所以,随着时间的积累,大陆岩石圈地幔往往具有低的εNd(t)值和高的εSr(t)值(Farmer et al., 1989; Menzies, 1989; Rudnick and Fountain, 1995)。达拉库岸岩体遭受了一定程度的同化混染作用,通常同化混染过程中混染物往往具有较高的Sr丰度和87Sr/86Sr比值,混染作用对幔源岩浆Sr同位素组成的影响程度往往会明显大于对Nd同位素的影响程度。在这种情况下,只有εNd(t)最高值才能在某种程度上反映源区的特征,而Sr同位素基本不代表源区性质。达拉库岸岩体岩石普遍弱富集轻稀土元素,具有较低的εNd(t)值(-4.02~+0.15),表明岩浆源区应该为富集的岩石圈地幔。

东昆仑地区晚古生代-早中生代是区域构造体制转换的重要时期,结合区域地层接触关系及花岗质岩石时空分布特点和相邻区域已发现镁铁质岩体的年代学和地球化学特征,表明至少251Ma开始,区域上已经开始由挤压环境转变为后碰撞伸展环境(莫宣学等, 2007; 奥琮等, 2015; Huang et al., 2014; Ding et al., 2014; 王亚磊等, 2017)。达拉库岸岩体形成于244Ma,应该是后碰撞伸展阶段的产物。通常,后碰撞伸展体制下形成的镁铁质岩体具有岛弧岩浆岩的地球化学特征,被认为是岩浆源区受俯冲板片流体交代改造的结果(Elliott et al., 1997)。在Th/Yb-Nb/Yb相关图中,达拉库岸岩体样品均落在了全球岛弧火山岩范围内,远离MORB-OIB地幔演化趋势线(图 10a)。大部分样品富集大离子亲石元素、有明显的Nb、Ta亏损,Th/Ta(4.79~57.67)、Th/Nb(0.23~1.99)比值较高,表明源区含有交代的地幔楔物质(Wilson, 1989)。由于Th、Zr和Nb具有相近的分配系数,分离结晶作用不会影响岩浆中这些元素的比值,因此,它们比值的显著不同可以用来指示源区特征(Saunders et al., 1992; Woodhead et al., 2001)。在Nb/Zr-Th/Zr图解中显示达拉库岸岩体源区存在俯冲流体交代作用(图 10b)。由此可见,达拉库岸岩体岩浆源区可能为板片流体交代改造的富集型岩石圈地幔。

图 10 达拉库岸岩体Nb/Yb-Th/Yb图解(a)和Nb/Zr-Th/Zr图解(b) Fig. 10 Plots of whole-rock Th/Yb vs. Nb/Yb (a) and Nb/Zr vs. Th/Zr (b) for the Dalakuan intrusion
6.4 成矿条件

达拉库岸岩体形成于后碰撞伸展阶段,与新疆东天山地区大量产出的岩浆型铜镍硫化物矿床形成环境基本一致(王京彬等, 2008; 夏明哲等, 2010)。岩体岩浆分异充分,岩石种类较丰富,原生岩浆为高镁拉斑玄武质岩浆,岩浆运移和侵位过程中遭受了一定程度地壳物质的同化混染作用,这些因素对形成岩浆型铜镍硫化物矿床是有利的(汤中立等, 2015)。岩体地表矿化特征明显,橄榄岩相中赋存铜镍硫化物矿化体。在岩体中开展的电法和磁法物探工作已圈定了一条宽约130m,长约400m,呈东西向展布的电磁综合异常。其中,激电异常和磁异常形态基本重合,激电剖面显示引起异常的极化体为南倾的中阻高极化体,磁测剖面结果也显示岩体南倾且向深部有一定的延伸,综合物性资料和地质特征推断异常主要为矿化超镁铁质岩石引起。

达拉库岸岩体中二辉橄榄岩、单辉橄榄岩和橄榄二辉岩中的橄榄石,具半自形粒状-浑圆状结构、堆晶结构,颗粒间分布有填隙状辉石,橄榄石裂纹发育,裂隙及边缘常发生蛇纹石化。在橄榄石/辉石中常见包裹有乳滴状、浑圆状硫化物,显示橄榄石结晶之前,母岩浆可能已经发生了硫化物的熔离。我们仅考虑只有橄榄石分离结晶以及橄榄石分离结晶的同时有硫化物熔离的情况来模拟橄榄石结晶过程。假设Ni在橄榄石与玄武质岩浆中的分配系数为7,在硫化物熔体与硅酸盐熔浆中的分配系数为500,达拉库岸岩体原始岩浆中MgO含量为11.4%、FeO为10.8%,Ni的含量为500×10-6,在总压力为0.2GPa,氧逸度为QFM-1,利用MELTS软件对岩浆演化和硫化物熔离过程进行模拟。如图 11所示,结果显示岩体中橄榄石均落于正常分离结晶线下方,表明橄榄石主要是从硫饱和的岩浆中结晶出来,而且在橄榄石的结晶过程中一直伴随有硫化物的熔离作用发生,橄榄石与硫化物的晶出质量比(橄榄石:硫化物)为15:1。综上所述,从岩体特征、矿石结构/构造、原生岩浆性质、深部硫化物熔离特征和物探信息等分析岩体具有形成岩浆型铜镍硫化物矿床的良好条件。目前,岩体地表出露岩相主要以辉石岩相为主,橄榄岩相较少,而铜镍硫化物矿化体常赋存在橄榄岩相或辉长苏长岩相中,因此需要进一步评价岩体深部橄榄岩相,追踪早期硫化物熔离时形成的硫化物熔体。

图 11 达拉库岸岩体橄榄石分离结晶模拟图解 Fig. 11 Compositional variations of olivine in the Dalakuan intrusion
7 结论

(1) 达拉库岸岩体主要由二辉橄榄岩、单辉橄榄岩、橄榄二辉岩、单辉辉石岩和辉长岩组成。辉长岩的锆石U-Pb谐和年龄为244.4±1.5Ma,属于中三叠世。

(2) 岩石普遍弱富集稀土元素和大离子亲石元素(Rb、Ba、Sr),亏损不相容元素(Nb、Ta)。岩体地球化学和Nd、Sr同位素组成表明岩浆源区为交代改造的富集岩石圈地幔。

(3) 岩体原生岩浆为高镁拉斑玄武质岩浆(MgO 11.4%、FeO 10.8%),岩浆运移和侵位过程中遭受不同程度地壳物质的同化混染作用。

(4) 从岩体特征、矿石结构/构造、原生岩浆性质、深部硫化物熔离特征和物探等综合分析,岩体具有形成岩浆型铜镍硫化物矿床的良好条件。

致谢      本次研究野外工作得到了新疆地质矿产开发局第三地质大队同行们的热情帮助;两名评审专家和编辑部老师提出的宝贵修改意见使本文更加完善;在此一并表示衷心的感谢。

参考文献
Ao Z, Sun FY, Li BY, Wang G, Li L, Li SJ and Zhao JW. 2015. U-Pb Dating, geochemistry and tectonic implications of Xiaojianshan gabbro in Qimantage Mountain, Eastern Kunlun Orogenic Belt. Geotectonica et Metallogenia, 39(6): 1176-1184.
Baker JA, Menzies MA, Thirlwall MF and MacPherson CG. 1997. Petrogenesis of quaternary intraplate volcanism, Sana'a, Yemen:Implications for plume-lithosphere interaction and polybaric melt hybridization. Journal of Petrology, 38(10): 1359-1390. DOI:10.1093/petroj/38.10.1359
Bureau of Geology Mineral Resources of Xinjiang Uigur Autonomous Region. 1993. Regional Geology of Xinjiang Uigur Autonomous Region. Beijing: Geological Publishing House: 657-717.
Campbell IH and Griffiths RW. 1992. The changing nature of mantle hotspots through time:Implications for the chemical evolution of the mantle. The Journal of Geology, 100(5): 497-523. DOI:10.1086/629605
Campbell IH and Griffiths RW. 1993. The evolution of the mantle's chemical structure. Lithos, 30(3-4): 389-399. DOI:10.1016/0024-4937(93)90047-G
Ding QF, Jiang SY and Sun FY. 2014. Zircon U-Pb geochronology, geochemical and Sr-Nd-Hf isotopic compositions of the Triassic granite and diorite dikes from the Wulonggou mining area in the Eastern Kunlun Orogen, NW China:Petrogenesis and tectonic implications. Lithos, 205: 266-283. DOI:10.1016/j.lithos.2014.07.015
Dong J, Huang HL, Yin JH, Li SJ and Tian CS. 2017. Geological characteristics of the Shitoukengde mafic-ultramafic rocks in East Kunlun and related metallogenic conditions. Northwestern Geology, 50(2): 50-59.
Elliott T, Plank T, Zindler A, White W and Bourdon B. 1997. Element transport from slab to volcanic front at the Mariana arc. Journal of Geophysical Research, 102(B7): 14991-15019. DOI:10.1029/97JB00788
Farmer GL, Perry FV, Semken S, Crowe B, Curtis D and DePaolo DJ. 1989. Isotopic evidence on the structure and origin of subcontinental lithospheric mantle in Southern Nevada. Journal of Geophysical Research, 94(B6): 7885-7898. DOI:10.1029/JB094iB06p07885
Gong XP, Ma HD, Yang XK, Li GT and Wang QM. 2004. Meaning and evolution and characteristic of Muztag-Cetacean Lake fracture zone. Geotectonica et Metallogenia, 28(4): 418-427.
Han HW, Wei MY, Mu LX, Yu JH, Zhang YS, Zhou ZY, Hu JM, Feng JX and Zhang CJ. 2007. The genesis of Katelixi copper zinc deposit in East Kunlun Mountains. Geotectonica et Metallogenia, 31(1): 77-82.
He GQ. 1994. Paleozoic Crustal Evolution and Mineralization in Xinjiang of China. Urumqi: Xinjiang People's Publishing House: 43-47.
Huang H, Niu YL, Nowell G, Zhao ZD, Yu XH, Zhu DC, Mo XX and Ding S. 2014. Geochemical constraints on the petrogenesis of granitoids in the East Kunlun Orogenic belt, northern Tibetan Plateau:Implications for continental crust growth through syn-collisional felsic magmatism. Chemical Geology, 370: 1-18. DOI:10.1016/j.chemgeo.2014.01.010
Jiang CF. 1992. People's Republic of China Ministry of Geology and Mineral Resources. Beijing: Geological Publishing House: 1-224.
Lai SC, Deng JF and Zhao HL. 1996. The Volcanism and Tectonic Evolution of the Northern Margin of Qinghai-Tibet Plateau. Xi'an: Shaanxi Science and Technology Press: 74-96.
Li C and Ripley EM. 2010. The relative effects of composition and temperature on olivine-liquid Ni partitioning:Statistical deconvolution and implications for petrologic modeling. Chemical Geology, 275: 99-104. DOI:10.1016/j.chemgeo.2010.05.001
Li SJ, Sun FY, Gao YW, Zhao YW, Li LS and Yang QA. 2012. The theoretical guidance and the practice of small intrusions forming large deposits:The enlightenment and significance for searching breakthrough of Cu-Ni sulfide deposit in Xiarihamu, East Kunlun, Qinghai. Northwestern Geology, 45(5): 185-191.
Mecdonald R, Rogers NW, Fitton JG, Black S and Smith M. 2001. Plume-lithosphere interactions in the generation of the basalts of the Kenya Rift, East Africa. Journal of Petrology, 42(5): 877-900. DOI:10.1093/petrology/42.5.877
Menzies A. 1989. Cratonic, circumcratonic and oceanic mantle domains beneath the western United States. Journal of Geophysical Research, 94(B6): 7899-7915. DOI:10.1029/JB094iB06p07899
Miyashiro A. 1974. Volcanic rock series in island arcs and active continental margins. American Journal of Science, 274(4): 321-355. DOI:10.2475/ajs.274.4.321
Mo XX, Luo ZH, Deng JF, Yu XH, Liu CD, Shen HW, Yuan WM and Liu YH. 2007. Granitoids and crustal growth in the East-Kunlun orogenic belt. Geological Journal of China Universities, 13(3): 403-414.
Naldrett AJ. 1999. World-class Ni-Cu-PGE deposits:Key factors in their genesis. Mineralium Deposita, 34: 227-240. DOI:10.1007/s001260050200
Qiu JX and Liao QA. 1987. The main characteristics and genetic information of clinopyroxene megacrysts in cenozoic basalt of eastern China. Acta Petrologica et Mineralogica, 6(1): 56-64.
Roeder PL and Emslie RF. 1970. Olivine-liquid equilibrium. Contributions to Mineralogy and Petrology, 29: 275-289. DOI:10.1007/BF00371276
Rudnick RL and Fountain DM. 1995. Nature and composition of the continental crust:A lower crustal perspective. Reviews of Geophysics, 33(3): 267-309. DOI:10.1029/95RG01302
Saunders AD, Storey M, Kent RW and Norry MJ. 1992. Consequences of plume-lithosphere interactions. In:Storey BC, Alabaster T and Pankhurst RJ (eds.). Magmatism and the Cause of Continental Breakup. Geological Society, London, Special Publication, 68:41-60
Sun SS and McDonough WF. 1989. Chemical and isotopic systematics of oceanic basalts:Implications for mantle composition and processes. In:Saunders AD and Norry MJ (eds.). Magmatism in the Ocean Basins. Geological Society, London, Special Publications, 42(1):313-345
Tang ZL, Jiao JG, Yan HQ and Xu G. 2015. Theoretical system for (large) deposit formed by smaller intrusion. Engineering Sciences, 17(2): 4-17.
Taylor SR and Mclennan SM. 1985. The Continental Crust:Its Composition and Evolution. Oxford: Blackwell: 1-312.
Wang JB, Wang YW and Zhou TF. 2008. Metallogenic spectrum related to post-collisional mantle-derived magma in North Xinjiang. Acta Petrologica Sinica, 24(4): 743-752.
Wang Y, Liu L, Che ZC, Chen DL and Luo JH. 1999. Geochemical characteristics of Early Paleozoic ophiolite in Mangnai area, Altun Mountains. Geological Review, 45(Suppl): 1010-1014.
Wang YL, Zhang ZW, Zhang JW, Qian B, Liu YG and You MX. 2017. Early Mesozoic mantle-derived magmatic events and their geological significance in the East Kunlun orogenic belt. Geology and Exploration, 53(5): 855-866.
Wilson M. 1989. Igneous Petrogenesis. London: Springer: 1-466.
Woodhead JD, Hergt JM, Davidson JP and Eggins SM. 2001. Hafnium isotope evidence for 'conservative' element mobility during subduction zone processes. Earth and Planetary Science Letters, 192(3): 331-346. DOI:10.1016/S0012-821X(01)00453-8
Wu YB and Zheng YN. 2004. Genesis of zircon and its constraints on interpretation of U-Pb age. Chinese Science Bulletin, 49(16): 1589-1640.
Xia MZ, Jiang CY, Qian ZZ, Xia ZD, Wang BY and Sun T. 2010. Geochemistry and petrogenesis of Huangshandong intrusion, East Tianshan, Xinjiang. Acta Petrologica Sinica, 26(8): 2413-2430.
Zhang ZW, Wang YL, Qian B, Liu YG, Qi CW, Dong J and Yin JH. 2017. Mineralogical characteristics of the Shitoukengde mafic-ultramafic intrusions in East Kunlun orogenic belt and the ore-forming indication. Geology and Exploration, 53(5): 825-837.
Zhou W, Wang BY, Xia MZ, Xia ZD, Jiang CY, Dong J and Xie ES. 2016. Mineralogical characteristics of Shitoukengde mafic-ultramafic intrusion and analysis of its metallogenic potential, East Kunlun. Acta Petrologica et Mineralogica, 35(1): 81-96.
Zindler A and Hart S. 1986. Chemical geodynamics. Annual Review of Earth and Planetary Sciences, 14: 493-571
奥琮, 孙丰月, 李碧乐, 王冠, 李良, 李世金, 赵俊伟. 2015. 东昆仑祁漫塔格地区小尖山辉长岩地球化学特征、U-Pb年代学及其构造意义. 大地构造与成矿学, 39(6): 1176-1184
董俊, 黄华良, 尹建华, 李世金, 田承胜. 2017. 东昆仑石头坑德镁铁-超镁铁质岩体地质特征及成矿条件分析. 西北地质, 50(2): 50-59.
弓小平, 马华东, 杨兴科, 李国填, 王庆明. 2004. 木孜塔格-鲸鱼湖断裂带特征、演化及其意义. 大地构造与成矿学, 28(4): 418-427.
韩红卫, 魏梦元, 牟伦洵, 余建华, 张永生, 周忠宇, 胡建明, 冯金星, 张春江. 2007. 东昆仑卡特里西铜锌矿成因. 大地构造与成矿学, 31(1): 77-82.
何国琦. 1994. 中国新疆古生代地壳演化及成矿. 乌鲁木齐: 新疆人民出版社: 43-47.
姜春发. 1992. 昆仑开合构造. 北京: 地质出版社: 1-224.
赖绍聪, 邓晋福, 赵海玲. 1996. 青藏高原北缘火山作用与构造演化. 西安: 陕西科学技术出版社: 74-96.
李世金, 孙丰月, 高永旺, 赵俊伟, 李连松, 杨启安. 2012. 小岩体成大矿理论指导与实践——青海东昆仑夏日哈木铜镍矿找矿突破的启示及意义. 西北地质, 45(5): 185-191.
莫宣学, 罗照华, 邓晋福, 喻学惠, 刘成东, 谌宏伟, 袁万明, 刘云华. 2007. 东昆仑造山带花岗岩及地壳生长. 高校地质学报, 13(3): 403-414.
邱家骧, 廖群安. 1987. 中国东部新生代玄武岩中单斜辉石巨晶的主要特征及成因信息. 岩石矿物学杂志, 6(1): 56-64.
汤中立, 焦建刚, 闫海卿, 徐刚. 2015. 小岩体成(大)矿理论体系. 中国工程科学, 17(2): 4-17.
王京彬, 王玉往, 周涛发. 2008. 新疆北部后碰撞与幔源岩浆有关的成矿谱系. 岩石学报, 24(4): 743-752.
王焰, 刘良, 车自成, 陈丹玲, 罗金海. 1999. 阿尔金茫崖地区早古生代蛇绿岩的地球化学特征. 地质论评, 45(增): 1010-1014.
王亚磊, 张照伟, 张江伟, 钱兵, 刘月高, 尤敏鑫. 2017. 东昆仑造山带早中生代幔源岩浆事件及其地质意义. 地质与勘探, 53(5): 855-866.
吴元保, 郑永年. 2004. 锆石成因矿物学研究及其对U-Pb年龄解释的制约. 科学通报, 49(16): 1589-1604. DOI:10.3321/j.issn:0023-074X.2004.16.002
夏明哲, 姜常义, 钱壮志, 夏昭德, 汪帮耀, 孙涛. 2010. 新疆东天山黄山东岩体岩石地球化学特征与岩石成因. 岩石学报, 26(8): 2413-2430.
新疆维吾尔自治区地质矿产局. 1993. 新疆维吾尔自治区区域地质志. 北京: 地质出版社: 657-717.
张照伟, 王亚磊, 钱兵, 刘月高, 祁昌炜, 董俊, 尹建华. 2017. 东昆仑石头坑德镁铁-超镁铁质岩体矿物学特征及成矿指示. 地质与勘探, 53(5): 825-837.
周伟, 汪帮耀, 夏明哲, 夏昭德, 姜常义, 董俊, 谢恩顺. 2016. 东昆仑石头坑德镁铁-超镁铁质岩体矿物学特征及成矿潜力分析. 岩石矿物学杂志, 35(1): 81-96.