高小卫,吴秀荣,杨振强
GAOXiao-Wei,WUXiu-Rong,YANGZhen-Qiang
(中国地质调查局武汉地质调查中心,武汉430205)
(Wuhan Center of China Geological Survey,Wuhan 430205,China)
苏门答腊岛是印度尼西亚盛产金、银、铜和锡等矿产的火山岛。岛上贵金属和有色金属矿分布于苏门答腊大断裂的两侧;锡矿产于东苏门答腊、廖内岛、班加岛和勿里洞岛(锡岛),以河流和滨外砂锡矿为主[1]。海西期岩浆旋回以来的金属矿产类型有:斑岩型Cu-Mo矿、密西西比型(MVT)铅-锌矿、花岗岩型锡矿(“锡岛”)、火山喷气型Au-Ag矿等。区域金属矿产的分布受控于不同地体中的岩浆-构造带。
苏门答腊岛在构造上位于巽他古陆西南边缘。根据不同的地层系统、沉积古地理、古生物地理区系、岩浆旋回和构造运动特征可将该岛晚海西-印支期的板块构造划分为东苏门答腊地体(亲冈瓦纳)和西苏门答腊地体(亲华夏地体),以及燕山期的沃伊拉群洋壳-洋岛推覆体。两者之间的界限为梅迪亚(中)苏门答腊构造带(MSTZ)[2-6]。西苏门答腊地体侵入岩的岩浆活动特征与东马来半岛I-型相似,具有强烈的石炭-二叠纪中-酸性岩浆侵入和基性火山活动,对贵金属和有色金属成矿和分布起着重要的控制作用,而东苏门答腊地体的S-型花岗岩与冈瓦纳体系的“暹缅马苏”地体的岩浆活动特征相似,对锡矿的成矿起着控制作用[7-8]。
本文是中国/印度尼西亚地学国际合作研究项目的初步成果。项目的重点研究区域为西苏门答腊地体巴东-明古鲁地区。笔者根据2010年以来在该研究区野外实地地质工作期间所采集的火成岩样品的岩石化学分析结果,并参考研究区周围火山-侵入岩的岩石化学分析数据(表1),对苏门答腊岛上不同地体的含矿性及其岩浆作用进行分析研究,探讨其含金(银)、铜和锡多金属矿产构造环境和源区。目的在于揭示本区含铜的岛弧型火成岩与相邻地体的含锡花岗岩的地球化学特征,并对其岩浆成因进行对比。
苏门答腊岛在构造上被划分为东苏门答腊地体和西苏门答腊地体。根据地体划分的标志和原则,大体上可划分为两条成矿带:西苏门答腊地体铜-金矿成矿带、东苏门答腊地体锡矿带成矿带、燕山期的沃伊拉群推覆体铜-金成矿带以及新生代金-银成矿带。各成矿带的分布特征受岩浆岩分布规律的控制[1]。位于东、西苏门答腊地体两者之间的印支期梅迪亚(中间)苏门答腊构造带(MSTZ)的含矿性特征更加接近东苏门答腊地体锡矿成矿带。东苏门答腊地体锡矿成矿带与西苏门答腊地体的含矿差异很大,而与“暹缅马苏”地体及其相邻地体的含矿性可以对比。新生代金-银矿沿苏门答腊-巴厘散大断裂两侧成带分布,受控于大陆边缘火山弧的岩浆活动,为印度洋俯冲的结果[1,6]。
铜-金矿成矿带呈NW-SE向延伸,在构造上隶属于西苏门答腊地体(断块)和沃伊拉群推覆体。该矿带的走向分布与海西-燕山期以来的岩浆活动带的走向一致,其中包括有石炭纪-中三叠世火山-侵入岩(Rb-Sr等时线年龄和白云母K-Ar法年龄287~246 Ma)带、晚三叠-早侏罗世火山-侵入岩带(黑云母K-Ar法、全岩Rb-Sr法和长石Rb-Sr法年龄224~181 Ma)、中侏罗世-白垩纪(沃伊拉群)火山岩带(黑云母、白云母K-Ar法年龄169~75 Ma)和新生代火山岩带(全岩K-Ar法同位素年龄65~1.6 Ma)[2-6]。海西-印支期的成矿带受印支期的MSTZ构造带的控制。燕山期的沃伊拉群及其加积复合体成矿带受中侏罗世-白垩纪火山弧的控制(见参考文献[4]中图1),与燕山早期巽他古陆边缘古裂谷的海底扩张、闭合和洋岛逆冲碰撞作用有关。
在本矿带内的矿化富集地带和地球化学异常带,自北西方向的穆阿拉西邦基铅-锌矿起,向东南方向依次分布的有:Mudik斑岩型铜金矿床,马塔比(Martabe)金矿群(Purnama金-银矿床等)、Padang铜(锌)矿、Lubuksulasih大型铜锌矽卡岩型矿床、苏利特河(Sulit Air)矽卡岩型铜-金多金属矿、苏里安矽卡岩型铜-铅-锌矿点,再向南东至邦科花岗岩南部一带,长度达250 km以上[9]。
铜-金矿成矿带的形成是海西-印支期以来多期次的火山-侵入活动的结果。岩浆活动的范围与成矿带的走向一致,起始于北苏门答腊省的实武牙花岗岩(Rb-Sr法年龄264~257 Ma;219~144 Ma),向东南方向依次分布的有:穆阿拉西邦基闪长花岗岩(黑云母K-Ar年龄197~158 Ma)[8]、辛卡拉克(Ombilin,翁比林)花岗岩(白云母K-Ar法、全岩Rb-Sr法等时线和白云母Rb-Sr法年龄246~287 Ma和206 Ma)、苏利特河(Sulite Air)花岗岩(白云母K-Ar法和角闪石40Ar/39Ar法年龄209~138 Ma)、苏里安(Surian)花岗岩和燕山期至喜山期的班科南(Bango South)花岗岩(白云母和角闪石K-Ar法年龄169~131 Ma)和拉西(Lassi)花岗岩(白云母、角闪石K-Ar法和白云母Rb-Sr法素年龄122 Ma;52~57 Ma)等。岩浆岩的岩石类型为岛弧型的基性至中-酸性火山-侵入岩套,包括闪长岩、二长花岗岩、闪长花岗岩、花岗岩、石英斑岩,以及辉石粗安岩、安山岩和辉石玻基安山岩、安山质火山角砾岩、晶屑玻屑凝灰岩。岩石产状为中-酸性侵入岩、脉岩和层状基-中性火山岩。据报道,分布于北苏门答腊的马迪纳勒根西地区(Madina Regency)的玄武岩、安山岩和I-型闪长花岗岩是属于岛弧型-弧后钙碱性中-基性火山岩和活动大陆边缘火山弧侵入岩[8]。本项目的最新研究结果也表明,该成矿带出露于辛卡拉克湖东南的诸多花岗侵入岩体和火山喷发岩均为埃达克质岩[9]。
白垩-新近纪的铜-金(银)矿成矿带沿苏门答腊-巴厘散大断裂两侧成带分布,受控于大陆边缘火山弧的岩浆活动,为印度洋俯冲的结果。侵入岩为I-型花岗岩、花岗闪长岩和闪长岩,火山岩为玄武岩、安山岩和高钾橄榄玄粗岩系列,构成典型的陆缘火山弧型岩浆类型[3-6]。
据文献记载,在中苏门答腊省钻孔中曾获得花岗岩的Rb-Sr法年龄(427~335 Ma,相当于志留纪~早石炭世)[5-7],证明了东苏门答腊地体是一个具有前石炭纪结晶基底(加里东期造山褶皱带)的古大陆壳的一部分。而与锡有关的海西-印支期花岗岩广泛分布于东苏门答腊地体,个别分布于西苏门答腊地体(例如二叠纪实武牙含锡花岗岩和加巴含锡花岗岩等等)。
1.2.1 “锡岛”锡矿成矿带
东苏门答腊地体是印度尼西亚重要的锡矿成矿带。地理位置位于北苏门答腊省、中苏门答腊省和南苏门答腊省的范围内,在构造上夹于梅迪亚(中间)苏门答腊构造带(MSTZ)和著名的文冬-劳勿缝合线之间的广大区域,可以划分为两条分支:一条分支包括宾坦(Bintan)岛、廖内群岛、新及岛、班加岛至勿里洞岛(“锡岛”)等一系列滨外岛屿,是位于西马来半岛上的“暹缅马苏”(Sibumasu)地体的主要山脉(Main Range)锡矿带向东南延伸的部分(图1),成矿时代为印支期(252~193 Ma),与主要山脉锡矿带比较,其成矿时代(247~143 Ma)要早得多。典型矿产地为:位于“锡岛”北西方向的丹绒潘当(Tanjungpandan)锡矿和克拉帕坎皮特(Kelapa Kampit)锡矿。
丹绒潘当锡矿为云英岩型和砂锡矿。含矿母岩为中三叠世(215 Ma)岩基。含锡岩浆岩为钛铁矿系列,由石英二长岩套和花岗岩套以及微粒花岗岩和细晶岩组成。黑云母花岗岩套包括花岗岩主相(19个平均样品)和伟晶微粒花岗岩(19个平均样品);石英二长岩套包括了辉长堆晶岩(3个平均样品)、含黑云母角闪石石英二长岩(8个平均样品)和碱性长石角闪花岗质伟晶岩[11]。
克拉帕坎皮特锡矿是东南亚最具有重要经济价值的层控型锡矿之一。它是同沉积物质来源的含锡沉积-变质岩,后期被花岗岩岩浆热液运移和叠加的例子,也是解释陆壳元素矿化与镁铁火山物质共生概念的重要范例[11]。含矿层(Nam Salu)是石炭-二叠纪低变质岩石系列(碎屑沉积与灰岩、燧石和铁矿互层,伴随海相镁铁火山物质),包含有变玄武岩(K K7)和绿泥石片岩(K 3),类似板块碰撞缝合线所常见的岩石类型。整个层系被三叠纪含锡花岗岩侵入。
1.2.2 梅迪亚(中)苏门答腊构造带锡矿成矿带
东苏门答腊地体锡矿成矿带的另一条分支位于北苏门答腊省和中苏门答腊省,包括梅迪亚(中)苏门答腊构造带(MSTZ)的两侧。主要矿产地有:哈塔邦(Hatapang)锡矿、洛干(Rokan)、双溪伊萨汉(Sungeilsahen)锡矿、Tigapuluh、加巴(Garba)等,成矿时代为印支晚期至燕山期,属于S-型含锡花岗岩成矿带(K-Ar法年龄203~186 Ma)[12],比西马来半岛上主要山脉锡矿带的时代要晚。成矿时代与古特提斯消亡和东、西苏门答腊地体拼合、碰撞,以及碰撞后地壳伸展时期完全一致。
1.2.3 “暹缅马苏”地体及其相邻的锡成矿带
东南亚三条重要锡矿带由主要山脉锡矿带西花岗岩省锡矿带和东花岗岩省锡矿带所组成。西马来半岛的主要山脉花岗岩带与东苏门答腊地体属于同一个岩浆岩省(图1),构成“暹缅马苏”地体锡矿带的主体,是马来西亚最具经济价值的锡矿带。
(1)主要山脉花岗岩带锡成矿带
主要山脉花岗岩带(MR)是由过铝质、S-型黑云母花岗岩组成(Rb-Sr法年龄220~200)Ma[11],为陆-陆碰撞和后碰撞的产物。该带向北方向延伸,在泰国境内演变为混合花岗岩带,以罗勇(Rayong)岩基为例[13]。
坚打谷锡矿区主要的含锡花岗岩相为三叠纪的斑状黑云母花岗岩(Rb-Sr等时年龄207 Ma)。细粒的含电气石花岗岩年龄为203~206 Ma(K-Ar法)。白云母花岗岩为强烈过铝质S-型花岗岩。吉隆坡周围和坚打谷锡矿区是马来西亚最大的砂锡矿产地。
泰国的罗勇岩基的侵入时代也为三叠纪。岩性为K-长石伟晶黑云母和黑-白云母花岗岩过铝质(S-型),(87Sr/86Sr)i为0.7193,未见锡矿化。
(2)西花岗岩省锡成矿带
西花岗岩省是与主要山脉花岗岩带相邻的另一条锡矿带,位于中南半岛的中西部(泰国的西北部)和缅甸。含锡花岗侵入体由准铝质角闪石-黑云母花岗岩(I-型)和过铝质黑云母花岗岩(S-型)组成,主要成矿期为白垩-新近纪(Rb-Sr法年龄95~55 Ma),锶同位素初始值为0.708~0.735,与印度洋板块俯冲成因有关。代表性的分界山脉(Border Range)岩体有:皮洛克(Pilok)、赫敏之(Hermiyingi,Rb-Sr年龄95.5 Ma)浅色花岗岩[10]和邦克拉辛(Pongkrathing)花岗岩。主要锡矿区为:位于缅甸土瓦东北的赫敏之钨锡矿和泰国南部的普吉(Phuket)和潘纳(Phangna)锡矿田,其中Takua Pa岩基的Rb-Sr年龄为78Ma[11]。
(3)东花岗岩省地体成矿带
东花岗岩省位于东马来半岛(即主要山脉花岗岩省之东),在构造上属于印支地体的海西-印支造山带的南延部分,此带锡矿不发育。该花岗岩省由二叠-三叠纪(白云母、角闪石K-Ar和Rb-Sr法年龄265~230 Ma)的角闪石-黑云母花岗岩和黑云母花岗岩组成;含角闪石花岗岩为I-型花岗岩,而庄他武里黑云母花岗岩更加靠近S-型花岗岩,为火山弧花岗岩。但是,它是否可能延伸入印度尼西亚境内值得怀疑,因为上述“锡岛”花岗岩也是以S-型花岗岩为主,叠加有东花岗岩省特有的I-型花岗岩。
地球化学分析资料证明,分布于北苏门答腊省马迪纳勒根西的钙碱性中-基性火山岩和活动大陆边缘陆缘火山弧花岗侵入岩的主量元素SiO2变化于46.28%~68.57%,MgO为1%~5.93%。一般K2O低于Na2O,只有穆阿拉西邦基闪长花岗岩K2O(5.24%)高于Na2O(2.75%)。安山岩和花岗岩类具有较高的不相容元素(Ba、Rb、K)、Th和P,含有较高的高场强元素Zr、Nb、Hf、Ta,而Sr和Ti较低(表1)。除了闪长花岗岩以外,在蛛网图上所有样品都显示Nb和Ta的低谷,为岛弧和主动大陆边缘环境的特征。
从辛卡拉克(亦称为翁比林,Ombilin)花岗岩、苏利特河闪长花岗岩至班科南(Bango)花岗岩一带是本项目的重点研究岩体。该花岗岩带中出露的埃达克质岩有花岗岩、闪长花岗岩、含铜石英斑岩、含矿石英脉岩和凝灰质岩等等[9]。其中,埃达克质花岗岩以高硅、高钾和低钛、低镁、低磷为特征(表1)。
本成矿带的埃达克质岩地球化学特征与典型埃达克岩的地球化学特征接近,其主量元素SiO2(63.81%~95.25%)和MgO(0.09%~1.88%)变化很大,且钛低(0.01%~0.73%)。在SiO2-K2O判别图上,大多数样品落在拉斑玄武岩系列范围,少数落在钙碱性岩石系列中。印支期的苏里安至班科南(Bango)的埃达克质花岗岩为高钾钙碱性系列(图2a、2c),在A/NK-A/CNK图解上,落在过铝质花岗岩类的范围内(图2b)。
本区埃达克质岩的微量元素表现为Th、Hf、Sr较低,亲石元素Rb不富集,也不太富集高场强元素,Sr比阿留申群岛的C-型埃达克岩低(表1)。Th变化于0.09×10-6~10.5×10-6之间。在微量元素蛛网图上显示为Ba、Ta、Nb、Sr、P、Ti低谷(图2c)。岩石的Sr/Y较低(图4a),大多数埃达克岩的Sr/Y小于C-型埃达克质岩(≤20)。在Harker图解上,埃达克质岩的Sr和TiO2与SiO2略成反比关系(图3a、b);Y和Nb与SiO2不成比例关系(图3c、d)。在Rb-Y+Nb)图解上(图6b)大多数埃达克岩样品落在火山弧花岗岩范围内。
表1 苏门答腊岛及其相邻地体岩浆岩主量元素(wt%)和微量元素(×10-6)组成Table 1 Major and trace element composition for igneous rocks in Sumatra and terrains nearby
续表1 表1 苏门答腊岛及其相邻地体岩浆岩主量元素(wt%)和微量元素(×10-6)组成
埃达克质岩的Y(1.16×10-6~16.7×10-6)变化很大,Yb(0.11×10-6~2.02×10-6)也很低(表1)。在La/Yb-Yb图解上,埃达克质岩分别落在C-型和O-型埃达克质岩的范围内(图4b)。其中,苏里安埃达克质花岗岩La/Yb(=10.2~14.58)比阿留申群岛的埃达克岩La/Yb(>20)低。根据埃达克质花岗岩REE配分曲线类型为右倾铲状型确定其为C-型埃达克质岩。
东苏门答腊地体岩浆岩的分布比较零星,主要岩体有哈塔邦花岗岩和双溪伊萨汉及其邻近地区花岗岩。前者(Rb-Sr等时线年龄为80 Ma),位于托巴(Toba)湖以南,Sr的原始比值为0.7151,属于大陆板内S-型花岗岩的性质;后者及其邻近地区为云英岩化和钾长石巨晶黑云母花岗岩(白云母石K-Ar法年龄198~193 Ma),SiO2值皆大于71%[14]。
位于巴东NE方向巴里散山脉以东的锡朱琼(Sijunjung)花岗岩岩基(白云母、角闪石K-Ar法年龄246~206 Ma)与西苏门答腊地体的岩浆岩带十分接近,具有高SiO2(72.71%)的地球化学特征,但其含矿性却类似于西马来半岛和“锡岛”的S-型花岗岩。
图1 东南亚锡矿带的花岗岩省的地理分布图(据Cobbing等,1986[11])Fig.1 Geographic distribution of granitic provincein south-eastern Asian tin belt(after Cobbinget al.,1986[11])
东苏门答腊地体岩浆岩的地球化学分析数据主要来源为“锡岛”丹绒潘丹含锡花岗岩套和克拉帕坎皮特含锡变质岩套[11]。
2.2.1 丹绒潘丹岩浆岩套
丹绒潘丹锡矿的含锡岩浆岩的主量元素变化很大。花岗岩套的主量元素以高硅、高钾、低钙、低镁铁、低钛、低锰、低磷为特征。SiO2(75.549%~76.10%)变化不大,MgO低(0.13%~0.02%)、K2O(5.10%~6.72%)>Na2O(2.15%~3.01%),和低的TiO2(0.07%~0.18%)。在A/NK-A/CNK图解上,落在准铝质-和过铝质花岗岩类过渡的范围内(图2b)。微量元素Ba、Nb、Th、Y有富集的趋势。Sr(26×10-6~102×10-6)和Zr(53×10-6~165×10-6)较低。F(715×10-6~1856×10-6)相当富集。在Rb-(Y+Nb)图解上(图5b)样品落在同碰撞带花岗岩和板内花岗岩的分界线上。
图2 火成岩的地球化学图解Fig.2 Geochemical discrimination diagramsfor igneousrocks
石英二长岩套的主量元素SiO2变化于48.65%~65.14%之间。石英二长岩则以低钙、低镁铁、低钛和低磷为标志,K2O(6.66%)>Na2O(3.62%),微量元素富集Zr(653×10-6),稀土元素La(212×10-6)、Ce(322×10-6) 和Y(43×10-6) 相当富集;在Rb-Y+Nb)图解上(图5b),样品落在同碰撞带花岗岩的范围内;辉长岩以高钙、高镁铁和富集钛为标志,Na2O(2.56%)>K2O(1.68%),微量元素富集Ba(403×10-6)和Sr(330×10-6),稀土元素La(83×10-6)、Ce(120×10-6)和Y(38×10-6)有富集的趋势,富集F(707×10-6);在Rb-(Y+Nb)图解上(图5b)样品落在同火山弧花岗岩(VAG)的范围内丹绒潘丹含锡花岗岩套的Sr同位素初始比值为0.7140,表明物质来源为地壳,而石英二长岩套Sr同位素初始比值为0.7049~0.7153,显示壳-幔物质之间的混合过程。
2.2.2 克拉帕坎皮特含锡花岗岩和变质岩套
克拉帕坎皮特锡矿含矿岩系显示同碰撞带岩浆侵入和热液交代的的地球化学特征,其主量元素SiO2(32.49%~90.33%)、MgO(0.81%~3.64%)、K2O(0.31%~4.16%)和TiO2(0.14%~1.463%)变化很大。微量元素Rb、Ba、Nb、Th、Y有富集的趋势。Ta(小于5×10-6)、Sr(4×10-6~201×10-6)和Zr(7×10-6~89×10-6)较低。在Rb-Y+Nb)图解上(图5b)落在同碰撞带花岗岩的范围内。
图3 火山岩的Harker图解Fig.3 Harker diagramsfor igneousrocks
图4 火成岩的Sr/Y-Y(a)和La/Yb-Yb(b)图解Fig.4 Sr/Y-Y(a)and La/Yb-Yb diagrams(b)for igneous rocks
2.2.3 “暹缅马苏”地体及其相邻地体的岩浆岩(1)主要山脉花岗岩带
在主要山脉花岗岩带内坚打谷锡矿区花岗岩的锶同位素初始值(87Sr/86Sr)i为0.715~0.751。其中,斑状黑云母花岗岩的(87Sr/86Sr)i为0.7193,过铝质的白云母花岗岩的微量元素模式为高度的结晶分离和强烈的热水叠加作用[11]。
“暹缅马苏”地体花岗岩的地球化学分析数据来自罗勇(Rayong)岩基。其主量元素以高硅、高钾、低钙、低镁铁、低钛、低锰、低磷为特征。大多数样品(25个样品平均)SiO2为72.60%,K2O(5.00%)>Na2O(2.74%),MgO为0.66%,钛0.313%(表1,样品号为R-3)。在SiO2-K2O判别图上,样品落在高钾钙碱性岩石系列和钾玄岩交界处(图2a)。微量元素表现为Sr含量较低,亲石元素Rb(351×10-6)、Ba(476×10-6)和Th(24×10-6)比较富集,高场强元素不太富集。在微量元素蛛网图上显示为Sr、P、Ti低谷(图2b)。Nb(15×10-6)趋向富集。Y(42×10-6)也比较富集,岩石的Sr/Y比值较低(图4a),表明岩浆源的长石结晶分离作用非常发育。
图5 火成岩的地球化学-构造环境判别图解Fig.5 Discrimination diagramsof geochemic-tectonic settingfor igneousrocks
(2)西花岗岩省含锡花岗岩
西花岗岩省含锡花岗岩的代表性样品的地球化学分析数据见于表1(样品号为B-4、P-5和H-6)。花岗岩的主量元素特征与泰国南部的罗勇岩基相似,都是以高硅、高钾、低钙、低镁铁、低钛、低锰、低磷为特征。大多数样品(平均)SiO2为73.42%~76.18%,K2O(4.38%~5.11%)>Na2O(2.92%~3.92%)。MgO为0.01%~0.35%,钛0.03%~0.25%(表1,样品号为R-4、P-5和H-6)。微量元素表现为Sr较低(4×10-6~65×10-6),亲石元素Rb(442×10-6~979×10-6)比较富集,Ba(13×10-6~255×10-6)变化较大。Nb(24×10-6~47×10-6)、Y(63×10-6~190×10-6)、Th(21×10-6~43×10-6)、U(19×10-6~31×10-6)比较富集。在Rb-(Y+Nb)图解上(图5b)样品落在同碰撞带花岗岩和板内花岗岩的分界线上。
该花岗岩省锶同位素初始值为0.708~0.735。普吉-潘纳锡矿田的Takua Pa岩基的锶同位素初始值为0.7346。
(3)东花岗岩省花岗岩
东花岗岩省花岗岩代表性样品的地球化学分析数据见于表1,分别为庄他武里花岗岩(样品号为C-2) 和黎府花岗闪长岩(样品号为L-1和C-2)。
黎府花岗闪长岩SiO2(平均)为63.32%,Na2O(3.43%)>K2O(2.99%),微量元素富集Ba(538×10-6)和Sr(381×10-6),而Rb(99×10-6)、Nb(7×10-6)、Y(22×10-6)和Th(10×10-6)不富集。在ANK-ACNK图解上,落在准铝质花岗岩类的范围内(图2b)。在Rb-(Y+Nb)图解上(图5b)样品落在火山弧花岗岩的范围内。
庄他武里花岗岩以高硅、高碱、低钙、低镁铁、低钛、低锰和低磷为特征。SiO2(平均)为70.26%,K2O(3.96%)>Na2O(3.76%),微量元素富集Ba(368×10-6)、Rb(368×10-6)、Y(51×10-6)和Th(23×10-6),而不富集Sr(118×10-6)、Nb(9×10-6)。在Rb-(Y+Nb)图解上(图5b)样品落在火山弧花岗岩和板内花岗岩的分界线上,而在ANK-ACNK图解上,落在过铝质花岗岩类的范围内(图2b)。
研究苏门答腊岛晚古生代不同地体的岩浆作用及其含矿性和成因具有重要的理论意义。它涉及到古特提斯构造域中与板块碰撞有关的岛弧和主动大陆边缘岩浆岩(包括埃达克质岩)岩浆-构造环境的判别和岩浆源区地球化学特征的对比研究。
众所周知,微量元素的比值及其图解是解释大陆边缘岩浆-环境重要手段。据文献记载,微量元素Th-Nb-Zr、Th-Ta-Hf和La-Nb-Zr地球化学图解[14,17-18]对板块汇集边缘(或大陆边缘)玄武-英安岩的大地构造环境有很好的判别效果,对与其共生的大陆边缘中-酸性岩浆系列也能起到很好的判别作用。笔者实践证明,这些地球化学-构造环境判别图(图5 c、d)以及Th/Yb-Ta/Yb判别图[20]不但适合于中-基性岩浆系列,也适合于主动大陆边缘长英质岩浆岩的环境判别和解释[9,19,21]。
3.1.1 铜-金矿带岩浆岩
本岛铜-金矿成矿带的大多数岩浆岩样品在Rb-(Y+Nb)图解上(图5e)落在火山弧花岗岩类的范围内。在La/Zr-Nb/Zr和Nb/Zr-Th/Zr判别图(图5c、d)上,多数岩浆岩样品分布于岛弧及其相邻的海、陆分界线(Nb/Zr=0.04)附近,落在岛弧和陆-陆碰撞带的范围内,少数落在大陆环境的拉张带(初始裂谷)、板内裂谷内,个别样品落在大洋板块离散边界N-MORB区的范围内,形成于洋壳环境(图5c)。这些图解表明,本矿带的岩浆岩主要形成于陆缘火山弧构造环境,与北苏门答腊省的马迪纳勒根西地区的岩浆岩套和南苏门答腊省楠榜第四纪火山岩的构造环境(图5c、d)几乎一致。
3.1.2 锡矿带岩浆岩
锡矿带岩浆岩的地球化学-构造环境判别数据主要以“暹缅马苏”地体花岗岩套为依据。其中,包括了取自东苏门答腊地体“锡岛”的丹绒潘丹含锡花岗岩套、克拉帕坎皮特含锡变质岩套和泰国南部罗勇花岗岩岩基,以及西花岗岩省的分界山脉岩体、皮洛克岩体和赫敏之岩体的岩石化学分析样品的平均值。该带的地球化学分析数据的缺点是缺乏高场强元素Hf、Ta和Yb,因而在地球化学-构造环境判别图的作图方面受到了极大的限制,没办法完成Th/Yb-Ta/Yb的图解的判别。
主要山脉罗勇花岗岩岩基和西花岗岩省3个岩体的样品(平均值),在Rb-(Y+Nb)图解上,都落在同碰撞带花岗岩带与大陆板内花岗岩带分界线上(图5e)。在La/Zr-Nb/Zr和Nb/Zr-Th/Zr判别图(图5c、d)上,这些岩浆岩样品分布于陆内拉张带(初始裂谷)和陆-陆碰撞区内(图5d)。其中,有一个西花岗岩省样品的Nb(平均值)高达到47×10-6,落在地幔热柱的范围内(图5c;表1,样品号P-5),表明“暹缅马苏”地体锡矿带的成矿的构造环境与地幔热柱隆起或陆内地壳减薄和拉伸引起的地幔热流活动有关。
东苏门答腊地体“锡岛”的的样品(平均值)在Rb-(Y+Nb)图解上(图5e)落在同碰撞带花岗岩带与大陆板内花岗岩带分界线附近的范围内(图5e)。在La/Zr-Nb/Zr和Nb/Zr-Th/Zr判别图(图5c、d)上,大多数样品分布于陆-陆碰撞区内,个别基性岩样品(TJPD-1)落在陆缘岛弧环境的范围内。其中,克拉帕坎皮特锡矿含矿岩系似乎富集微量元素Rb和Nb,而丹绒潘丹岩浆岩有更加富集微量元素Ba、Th、Y的趋势。前者的陆-陆同碰撞的特征更加明显。
3.1.3 东花岗岩省岩浆岩
东花岗岩省岩浆岩的地球化学-构造环境可与北苏门答腊省马迪纳勒根西的岛弧-陆缘火山弧岩浆岩带对比,为准铝质至过铝质花岗岩类的过渡类型(图2b)。在Rb-(Y+Nb)图解上,2个样品(平均值)落在火山弧花岗岩和与大陆板内花岗岩带分界线上(图5e)。在La/Zr-Nb/Zr和Nb/Zr-Th/Zr判别图(图5c、d)上,样品分布于陆缘岛弧及其相邻的大陆板内分界线(Nb/Zr=0.04)附近,显示岛弧和活动大陆边缘环境。
岛弧型钙-碱性岩浆岩(特别是埃达克岩)以及活动大陆边缘的火山弧花岗岩的成因和源区是近十年来化学地球动力学所探讨的重要内容。通过对与洋壳板块俯冲有关的埃达克岩的Sr/Y和La/Yb比值以及稀土元素图解的研究,可以认为岛弧型岩浆岩(包括埃达克岩)具有岩石类型和成因多样化的特点[22],其岩浆岩区可以归纳为:俯冲板块的部分熔融、地幔楔局部交代和部分熔融、下地壳的拆沉和熔融MASH(熔融-混染-储存-均一化)和上地壳AFC(混染-分异-结晶)作用有关[23]。
从以上Harker图解(图3)所示,本岛岩浆岩样品的Sr、MgO和TiO2与SiO2含量呈反比关系,暗示其岩浆源区,除了俯冲板块部分熔融作用以外,还有上地壳AFC作用(包括角闪石和长石的结晶分离作用)。但是,在Y-SiO2图解(图3c)上,随着SiO2含量达到60%以后,Y明显地显示出有两种不同的演化轨迹:西苏门答腊地体的Y出现与SiO2成反比关系的现象,而东苏门答腊地体(以及“暹缅马苏”地体)的Y却与SiO2成正比关系,暗示前者的岩浆作用以俯冲板块部分熔融作用为主,而后者则以长石结晶沉淀作用为主。这表明西苏门答腊地体(主要为埃达克质岩)的岩浆成因与洋壳板块俯冲部分熔融作用有关,其岩浆源区达到地幔楔之下的深度(相当于80~100 km),位于俯冲板片的残留相为高压矿物相(榴辉岩相)的位置。由于在A/NK-A/CNK图解上,有部分埃达克质岩落在过铝质花岗岩类的范围内(图2b),可以肯定其岩浆成因与上地壳的长石的结晶分离作用也有密切联系。与此相似,K2O、Sr和Sr/Y也与SiO2呈反比关系(图略),暗示上地壳岩浆源区钾长石结晶分离作用比较强烈,其源区的位置深度不可能超过70 km。
另一方面,岛弧型岩浆岩中MgO含量变化,也是岩浆岩物源区的重要标志。西苏门答腊地体岩浆岩的MgO不高于6%,其MgO的高值样品出现在中侏罗世-白垩纪(沃伊拉群)火山岩和新生代始新世火山岩中。本岛岩浆岩MgO含量随着SiO2含量增加呈现递减的趋势(图3)。Zr/Nb-MgO图解是揭示熔体中俯冲沉积物混入到地幔楔中的程度的可靠方法[24-25]。从Zr/Nb-MgO图解(图6a)上可以看出,本岛大多数岩浆岩样品的Zr/Nb比值小于40,落在洋壳MORB及以下的范围内(40~20),其物质源区受俯冲洋壳MORB板片的局部熔融作用影响十分强烈。而沃伊拉群始新世火山岩的样品普遍具有Zr/Nb>40比值,充分证明新生代火山岩的岩浆确实起源于被交代的地幔楔的部分熔融,极少有俯冲沉积物混入。在巴东地区,有一个变质火山岩样品(HQ-5)具有较高的Zr/Nb比值(=39),并且MgO含量较低(1.65%),表明俯冲洋壳之上的地幔楔受到一定的交代作用,其主要来源为俯冲的洋壳MORB的部分熔融和地幔楔局部熔融区。本区内大多数样品MgO含量小于2%,显然代表俯冲的洋壳板片和海沟深海沉积物局部熔融作用。对比之下,北苏门答腊马迪纳勒根西地区岛弧型火山-侵入岩的样品,其MgO(3.07%~3.86%)相对较高,暗示其物质源区是来自被交代作用的地幔楔,局部熔融的俯冲岩层熔体物质和热水对地幔楔的混染作用较强,与研究区内埃达克质岩的成因和岩浆岩源区略有不同,可能比西苏门答腊地体遭受过更强烈地幔楔局部熔融和交代作用。
除了俯冲洋壳板片和海沟深海沉积物局部熔融作用以外,上地壳AFC作用也很重要。本岛样品在Zr/Nb-Zr图解(图6b)上分别显示两种不同的岩浆演化趋势:强烈的俯冲洋壳板片局部熔融的演化趋势(左箭头)和相对较弱的来源于上地壳分离结晶(右箭头)演化趋势。本岛岩浆演化是两种不同岩浆作用结果:即西苏门答腊地体为洋壳板片和海沟深海沉积物局部熔融叠加弱的MASH作用和东苏门答腊地体强烈的AFC作用的结果。在Ba/Nb-Nb图解和Sr/Y-Y图解(图略)上,都显示与Zr/Nb-Zr图解上所表示的岩浆演化趋势。
图6 火成岩的Zr/Nb-MgO图解(a)[25]和Zr/Nb-Zr图解(b)[26]Fig.6 Zr/Nb-MgO(a)[25]and Zr/Nb-Zr(b)discrimination diagrams[26]for igneous rocks
图7 含锡花岗岩及变质岩石的Sn含量变化图解Fig.7 Variation diagramof Sn contant for tin-bearinggraniteas well as metamorphic rocks
但是,不管岩浆演化的途径不同,其岩浆演化机制和源区却相同。因此,铜-金矿成矿带的岩浆岩套以及丹绒潘丹岩浆岩表现在Ce-La图上(图5a)呈线性关系,与南苏门答腊省楠榜第四纪火山岩的岩浆演化规律一致,显示其岩浆成因都是与俯冲板片机制有关,属于来源于深部相似性质的岩浆源区[13],而克拉帕坎皮特锡矿含矿岩系和西花岗岩省浅色花岗岩则排列成另一条直线,表示其岩浆源区不同于西苏门答腊地体。
根据上述资料,笔者认为:西苏门答腊地体铜-金矿成矿带的成因与埃达克质岩的成因相同,岩浆演化过程以俯冲洋壳板片局部熔融叠加上弱的MASH为主,以上地壳AFC作用为次,而东苏门答腊地体锡成矿带以及与其相邻的“暹缅马苏”地体锡成矿带则都是AFC的结果。铜-金矿成矿带与壳-幔相互作用形成的I-型(岛弧型)花岗岩岩浆侵入有关[27],而锡成矿带与大陆碰撞带的后造山岩浆AFC作用和大陆边缘的扩张过渡带(弧后盆地)以及陆内裂谷的构造环境有关。实际上,锡矿成因机制十分复杂,包括成矿环境、岩浆系统、岩浆-热液交换系统和热液系统等因素。
本岛含锡母岩(花岗岩和变质岩)样品在Sn含量的图解上(图7),可以划分为热液变质区域和岩浆分异区域。在图7a、c和d上显示,母岩的Sn与Y、Rb/Sr和SiO2都呈正比关系,而与TiO2含量呈反比关系,其中巴东地区的埃达克质岩含有最低的Sn含量。
本岛丹绒潘丹锡矿岩浆岩主要相带Sn含量很低,平均为(2~4)×10-6。含锡母岩的Sn与Y(或Y+Nb)和SiO2和Rb/Sr比值成正比关系,而SiO2与TiO2成正比关系(图7)。含锡岩浆岩套的F含量(273~1856)×10-6富集,表明岩浆-热液交换系统对成矿的富集作用。
克拉帕坎皮特锡矿含矿岩系样品Sn含量落在热液变质的范围内。Sn与微量元素Ba、Rb以及稀土元素La和Ce皆成正比关系。由图7上所示,该锡矿含矿岩系Sn与Y(以及Sn与Rb/Sr比值)之间成正比关系,而Sn与SiO2(或与TiO2)呈分散分布(图7)。资料表明,含矿岩系的F(542~603)×10-6十分富集,岩浆热液对成矿起着富集作用。
相比之下,“暹缅马苏”地体主要山脉含Sn花岗岩的样品Sn的平均值为10×10-6,在TiO2-Sn图解上(图7b),总体上显示比较分散,其中的白云母花岗岩的Sn在该图解上比黑云母花岗岩的含量高些。表明其含锡模式受碰撞后岩浆的结晶分离作用所控制,并受到热水交代作用的影响。
西花岗岩省岩浆岩主要相带Sn较高,平均为(14~72)×10-6。Sn富集受热水交代作用的影响比主要山脉和丹绒潘丹锡矿岩浆岩更强。
(1)根据地体划分的标志和原则,本岛大体上可划分为西苏门答腊地体铜-金矿成矿带、东苏门答腊地体锡矿成矿带和燕山期的沃伊拉群推覆体铜-金成矿带以及新生代金-银成矿带。
(2)铜-金成矿带的岩浆岩系列为准铝质至过铝质岩浆岩类型,其中含矿母岩为高硅、高钾的埃达克质(过铝质)花岗岩(属SI-型花岗岩)。锡矿成矿带为准铝质(I-型)和过铝质(S-型)花岗岩的过渡岩浆岩系列,以S-型花岗岩为主。I-型为岛弧-活动大陆边缘(ACM)型火山弧花岗侵入岩,S-型为大陆碰撞带的后造山岩浆分异的产物。构造环境为碰撞带的弧后盆地和陆内裂谷。
(3)岛弧-活动大陆边缘(ACM)型火山弧花岗岩的源区为俯冲洋壳板片局部熔融和弱的MASH为主,上地壳源区的AFC作用比较次要;而东苏门答腊地体锡成矿带及其相邻的“暹缅马苏”地体锡成矿带的的岩浆源区来自地壳重熔和岩浆AFC作用,其中西花岗岩省锡矿成因与局部的地幔柱(或幔隆)构造环境引起的陆内地壳减薄和拉伸构造和地幔热流活动有关。
(4)西苏门答腊地体和东苏门答腊地体岩浆岩在Ce-La图上呈线性关系,显示都是来源于地壳深部与俯冲板片机制有关的岩浆源区和相似的岩浆演化机制。
[1]Crow M J,van Leeuwen.Metallic mineral deposits[M].//In:Barker A J,Crow MJand Milson JS(eds).Sumatra:Geology,Resources and Tectonics.Geological Society,London,Memoir,2005,31:147-174.
[2]Hutchison CS.Gondwanaand Cathaysian blocks,Palaeotethys sutures and Cenozoic tectonics in Southeast Asia [M].New York:Offord University Press,1994:388-405.
[3]高小卫,吴秀荣,杨振强.苏门答腊(印度尼西亚)火山岩的地层时代及其中、新生代地球化学-构造环境判别[J].华南地质与矿产,2012,28(2):107-113.
[4]高小卫,吴秀荣,杨振强.初论苏门答腊(印尼)的岩浆-构造旋回及其板块构造背景[J].华南地质与矿产,2013,29(4)259-270.
[5]高小卫,杨振强,吴秀荣.苏门答腊岛岩浆旋回中的成矿作用[J].华南地质与矿产,2013,29(4):299-307.
[6]高小卫,杨振强,吴秀荣.苏门答腊岛(印尼)两类地体划分及其区域成矿背景[J].地质通报,2015,待刊.
[7]Barley ME,Rak Pand Wyman D.Tectonic controls on magmatic-hydrothermal gold mineralization in the magmatic of SEAsia[M].//In:Blundell DJ,Neubauer Fand von Quadt A(eds).The Timing and Location of Major Deposits in a Evolving Orogen.Geological society,London,Special Publication,2002,204:39-47.
[8]Zulkarnain I.Geochemical Signature of Mesozoic Volcanic and Granitic Rocks in Madina Regency Area,North Sumatra,Indonesia,and its Tectonic Implication[J].Indonesian Journal on Geoscience,2009,14(2):117-131.
[9]朱章显,杨振强,胡鹏.印尼苏门答腊岛巴东地区埃达克质岩地球化学特征和构造环境意义[J].中国矿业,2014,23(8):92-99,109.
[10]王 宏,林方成,施美凤,石洪召,王疆丽.滇缅腾冲-毛淡棉构造岩浆岩带钨锡成矿系统与典型矿床研究[J].地质与勘探,2013,49(1):89-98.
[11]Lehmann B.Metallogeny of Tin.Berlin:Spring-Verlag,1990:1-211.
[12]Koning T,Darmono FX.Thegeology of the Beruk Northeast Field,Central Sumatra:Oil productin from Pre-Tertiary basement rocks[M].//In:Indonesian Petroleum Assosiation Proceeding of the 13th Annual convention,Jakarta,1984,1:385-406.
[13]Zulkarnain I.Geochemical evidence of island-arc origin for Sumatra Island:a new perspective based on volcanic rocks in Lampung Province,Indonesia [J].Jurnal Geologi Indonesia,2011,6(4):213-225.
[14]Cobbing E J.Granite[M].//In:Barker A J,Crow M Jand Milson JS(eds),Sumatra:Geology,Resourcesand Tectonics.Geological Society,London,Memoir,2005,31:54-62.
[15]武莉娜,王志畅,汪云亮.微量元素La、Nb、Zr在判别大地构造环境方面的应用[J].华东地质学院学报,2003,26(4):343-317.
[16]朱章显,杨振强.巴布亚新几内亚新生代两类埃达克岩的构造环境意义[J].华南地质与矿产,2007,(2):1-6.
[17]汪云亮,张成江,修淑芝.玄武岩类形成的大地构造环境的Th/Hf-Ta/Hf图解判别[J].岩石学报,2001,17(3):413-421.
[18]孙书勤,张成江,黄润秋.板块汇聚边缘玄武岩大地构造环境的Th、Nb、Zr判别[J].地球科学进展,2006,21(6):593-598.
[19]朱章显,赵财胜,杨振强.苏拉威西埃达克岩、类埃达克岩分布和特征[J].吉林大学学报(地球科学版),2009,39(1):80-88.
[20]Gorton M P,Schandl ES.From continents to island arcs:a geochemical index of tectonic setting for arc-related and within-plate felsic to intermediate volcanic rocks.Canadian Mineralogist,2000,38(5):1065-1073.
[21]胡鹏,朱章显,杨振强.爪哇和加里曼丹(印尼)新生代埃达克岩地球化学特征、成因及构造环境对比[J].华南地质与矿产,2015,31(2):125-135.
[22]张旗,许继峰,王焰,肖龙,刘红涛,王元龙.埃达克岩的多样性[J].地质通报,2004,23(9-10):959-965.
[23]Richards JP,Kerrich R.Special paper:adakite-like rocks:their diverse origins and questionable role in metallogenesis[J].Economic Geology,2007,102(4):537-576.
[24]Crow MJ,Tertiary volcanicity[M].//In:Barker A J,Crow M J,and Milson J S(eds),Sumatra-Geology,Resources and Tectonics.Geological Society,London,Memoir,2005,31:98-119.
[25]Macpherson CG,Hall R.Tectonic controls of geochemical evolution in arc magmatismin SEAsia(C).//In:Proceedings PACRIM 99 congress,Australian Institute of mining and metallurgy Publications,Series,1999:359-368.
[26]Maulana A.Petrology,Geochemistry and Metamorphic Evolution of South Sulawesi Basement Rock Complexes,Indonesia[D].Canberra:The Australian National University,2009:1-89.
[27]向文帅,高小卫,程 湘,印尼苏门答腊岛苏利特河铜矿地质特征[J].华南地质与矿产,2014,30(4):361-367.