小秦岭地区金矿床成矿构造地球化学动力学研究

1、页 ( Pa ge s) 155162 ,2000 ,6 (J u n . ,2000)Geo t ec to nica e t M e t allo ge ni a小秦岭地区金矿床成矿构造地球化学动力学研究方维萱 1 ,3 ,李亚林 2 ,黄转莹 1 ,31 . (中国科学院地球化学研究所矿床地球化学开放研究实验室 ,贵州 贵阳 550002)2 . (成都理工学院沉积学研究所 ,四川 成都 610059)3 . (西北有色金属地质勘查局 ,陕西 西安 7100541)摘 要 :小秦岭地区金矿床与构造有着密切的关系 ,不同方向的断裂中元素富集规律不同 ,构造带中流体是十分活跃的 。构造动力成

2、矿作用主要有动热退变质作用 、剪切 - 扭动裂隙成矿作用 、构造带内 岩石与流体交换作用 、构造扩容带成矿作用 , 这些成矿成晕作用均发生在特定的构造动力学系统 中 。关键词 :构造地球化学形迹 构造地球化学特征 成矿成晕作用 成矿作用动力学中图分类号 : P618 . 2/ 618 . 52 ; P595文献标识码 : A文章编号 :1001 - 1552 (2000) 02 - 08近年来随着构造地质学的进展 , 构造地质与矿床学 、地球化学之间不断相互渗透 , 形成了一些新兴交叉学科 , 如成矿构造学 、构造地球化学 、成矿作用动力学 , 这些交叉学科的研究成果 使构造与成矿之间的关系更

3、为密切 1985 , 1992 , 陈国达 ; 1982 , 1992 张湘炳 ; 1990 , 杨国清 ;1993 ,孙岩等 ; 1993 ,彭少海 ; 1994 ,李德威 。小秦岭地区是我国金 、钼的重要矿产地 ,在金成矿过程中 ,构造作用是十分重要的 ,它不仅控制了矿产形成 ,而且为矿体定位提供了储存空间 。构 造是导矿 - 储矿系统 , 作为成矿作用的统一整体组成部分 , 直接参与了金 、钼的成矿作用 。因 此 , 研究本区成矿作用地球化学动力学可以探索金富集成矿规律 , 指导成矿预测和隐伏矿床的 勘查 。成矿构造地球化学形迹1地球化学异常 、含金石英脉及含金构造带主要沿断裂破碎带分布

4、 , 地球化学异常所显示的构造地球化学行迹是十分显著地 ,其特征详见表 1 。收稿日期 :2000 - 01 - 13 ;改回日期 :2000 - 02 - 28基金项目 : 本文研究由国家重点黄金地质科技攻关项目 ( 90051 - 09 - 08) 及中科院地化所矿床地球化学开放实验室2 断裂构造地球化学特征2 . 1 断裂构造地球化学特征及金富集规律不同产状断裂中元素的富集成矿规律不同 (表 2) ,东西向断裂中富集 A u 、A g 、P b 、Mo 、Bi ,倾角较缓 (小于 30°) 的东西向断裂中 A u 最富集 , 而倾角较陡 (大于 50°) 的东西向断裂

5、中比前者更富集 A g 、P b 、Mo 、Bi 。近南北向断裂中富集 A u 、A g 、C u 、P b 、Zn 、Mo , A g 和 C u 含量最 高 。北东向断裂中富集 A u 、A g 、P b 、Mo 。北西向断裂中赋存独立 A g 矿体 。A u 主要富集于东 西向 、南北向 、北东向断裂中 , A u 在断裂中富集程度由东西向 南北向 北东向 北西向依次 降低 。A g 在断裂中富集程度由北西向 南北向 东西向 北东向依次降低 。Mo 富集程度由 东西向 北东向 北西向 南北向依次降低 。北东向断裂组是本区活动历史最长的断裂 ,在燕 山期活动性仍很强 ,其他方向断裂中含 H

6、 g 较低 ,唯有北东向断裂中有线状 H g 异常 ,含 H g 为176 ×10 - 9 1100 ×10 - 9 , H g 异常多出现于断裂下盘附近碎裂状岩石及矿石中 。 断裂构造的规模及发育情况与金的富集成矿关系密切 , 本区含金石英脉及含金蚀变体均产于断裂构造中 ,断裂构造规模及发育密度愈大 ,金的富集程度和规模愈大 。如潼峪金矿田内 ,层间滑动剪切构造带已发现有五层 ,东西向储矿断裂构造 505 、401 号长分别为 4630 m 、3000m ,已发现两个中型金矿 ( Q 505 、Q 401) 产于其中 。葫芦沟金矿床内 ,已发现七条含矿蚀变带 ,有两条长度

7、超过 2000 m ,葫芦沟中型钾质蚀变岩型金矿的主要储量分布于这两条蚀变带中 。 断裂构造活动期次愈多 , 金富集程度愈高 。石英脉充填于断裂构造中 , 未经破碎含 A u 性差 , A u 多在 100 ×10 - 9 以下 。石英脉多次破碎 , 裂隙发育 , 而裂隙又被硫化物网脉所胶结 、充填 ,金富集成矿 。只有构造作用的地段 ,常发育构造片岩 ,A u 不富集 ,或只发生轻微富集 ,A u< 100 ×10 - 9 。有流体参与常形成蚀变作用 ,充填有石英脉 ,A u 发生富集 。如绿泥石化构造片岩中 A u 可达 47 . 6 ×10 - 9 ,

8、构造片岩中夹有石英细脉时 ,A u 为 168 ×10 - 9 980 ×10 - 9 。形成蚀 变后或蚀变岩经再次断裂破碎 ,A u 的富集程度进一步升高 ,形成蚀变岩型金矿石 (见表 3) 。多次活动的断裂 , 其各期次动力学结构面不相同时 , 对金富集成矿有明显影响 。浅层次的 脆性断裂以张性 、张扭性的断裂为主 , 常是对深层次韧性剪切构造带的叠加变形 , 对含金石英脉形成有利 , 是矿体的定位提供构造扩容 - 充填的空间 。以充填 - 扩散交代成矿作用为主 , 成 矿流体常沿构造带有限空间活动 ,所形成的异常宽度小 ,异常长度大而强度高 。以压性 、压扭性 为主的

9、断裂 , 挤压片理化带和碎裂蚀变岩发育 , 有利于形成含金蚀变岩 , 常具有脆性构造带特 征 , 是对深层次韧性剪切带的叠加变形 。由于成矿流体在构造带内受到外界应力作用 , 沿片理化带和碎裂蚀变岩形成弥漫性粒间流体 , 成矿流体以渗滤交代作用为主 , 异常宽度可达 200 m 以上 ,异常浓集中心宽度大 ,但长度多在 500 m 以内 。当张性 、张扭性断裂与压性 、压扭性 - 韧 性带叠加或复合 ,形成蚀变岩型金矿和含金石英脉型金矿复合类型 。断裂构造带内及附近派生微小构造 (如挤压片理化带 、裂隙密集带 、碎裂岩带 、矿物线理)2 . 2 含金构造带的分带特征及金银富集规律含金构造带的水

10、平分带现象明显 , 如回马坪 - 朱家沟壳型断裂是本区高序次构造 , 属导矿 构造 。在孙家沟一带的发育宽度在 1000 余 m , 在垂直断裂带方向构造岩及元素有明显水平分带 (从外向中心随着应变的增强) , 碎裂蚀变混合花岗岩 初碎裂岩化混合花岗岩 ( Zn 、Hg)碎裂岩 ( A g) 糜棱岩 ( P b 、A g) 超糜棱岩 ( Si 、K) 。在含金石英脉岩金矿中 ,表现为脆性构造的叠加变形分带 , 在垂直断裂带方向构造岩及元素水平分带 (从外向中心) 为混合岩化片麻岩( N i 、Co 、C r 、V 、Ti ) 糜棱岩化片麻岩 ( A u 、A g) 碎裂岩 ( A u 、A g

11、 、K) 蚀变碎裂岩 ( A u 、A g 、 Mo 、P b) 破碎石英脉 ( A u 、A g 、Mo 、P b 、Si 、Fe ) 。葫芦沟钾质蚀变岩型金矿中 ,含金蚀变岩及 元素水平分带为外带黄铁绢云岩化带 ( K 、Fe 、Si 、A u ) ,中带钾长石化 - 黄铁矿化 - 硅化带 ( K 、 Si 、Fe 、A u ) ,内带黄铁矿化 - 硅化带 ( Fe 、Si 、A u 、Mo 、C u 、P b 、Zn ) 。铁族元素 (表 4) C r 、V 、Ti 、 N i 、Co 被活化带出 ,A u 、Mo 富集成矿 。西沟 - 大银硐沟含银韧性剪切带的水平分带为 (自外向内)

12、: 绢英岩化混合花岗岩 碎裂蚀变辉绿岩 ( A g 、Ti 、C r 、V 、N i 、Co ) 绿泥石线理带 ( A g 、C u 、P b 、Zn ) 含银蚀变岩带 ( A g 、C u 、P b 、Zn 、H g 、W 、F 、B a ) 。本区虽然由于不断抬升剥蚀 , 深部韧性剪切带不断遭受脆 - 韧性变形 , 但总体上仍显示有垂直分带的特点 ,从深部到浅部为 : 混合岩带 ( K 、N a 、Si ) 混合岩化带 ( K 、Si 带入 , M g 、Ca 带 出) 韧性剪切带 (浅色矿物集中 ,暗色矿物形成线理 ,塑性流变发育 。K 、Si 、Fe 、M g) 脆 - 韧 性断裂带

13、( K 、Si 、Mo 、Fe ) 脆性断裂带 ( C u 、P b 、Zn 、Bi 、A g 、A u 、Si ) , 成矿元素集中于脆 - 韧 性断裂带 、脆性断裂中富集成矿 , K 、Si 趋向于脆 - 韧性断裂带富集 , Si 趋向于脆性断裂带富 集 。这种宏观垂直分带中 , 产于脆性 、脆 - 韧性含金构造带中金矿类型也有垂直分带现象 , Q 04 、Q 2142 、Q 2144 、Q 539 等出露标高为 18521655 m , 属自然金 - 多金属硫化物 - 石英脉 型金矿 。岳王庙地区的标高为 1400 m ±,形成含金石英脉型 、蚀变千糜岩型及蚀变岩型三种复合的金

14、矿 。葫芦沟地区的标高为 1272 m ±, 形成钾质蚀变岩型金矿 , 呈现“上脉型 、中复合型 、下蚀变岩型”的规律 。岳王庙地区的含金构造带中 ,成矿元素 ( A u 、C u 、P b 、A g) 具有自地表向深部富集的趋势 (见表 5) , 它是隐伏金矿体的构造 - 矿物 - 地球化学垂直分带 , 是预测隐伏金矿 体的构造 - 矿物 - 地球化学的标志 。3构造带中流体动力学笔者 1991 ,方维萱 估计过混合岩化过程中矿物脱水为变质热流体形成提供的 H2 O 量 ,3从角闪斜长片岩 、片麻岩类到混合岩类 ,H2 O 降低了 0 . 11 % , 每 k m岩石可形成的变质水

15、达275 万 t (按岩石平均比重为 2 . 5t / m3 ) ; 当发展到中等程度的混合岩类时 , H O 降低了 1 . 1 % ,2每 k m3 岩石可形成 2750 万 t 变质水 , 假定古大气降水的参与忽略不计 ; 混合岩化作用形成的变质流体中含 A u 可达 6 . 5 ×10 - 6 ,形成了富 A u 的变质流体 。显然 ,混合岩化过程中变质流体 的规模宏大 。中期混合岩化大约在 69 亿年间发生 , 表现为长英质脉体沿片麻理 、裂隙注入充填 , 碱质升高 , 碱质 、硅质交代作用进一步增强 , 地幔柱热蚀变作用导致碱性流体的交代作用强烈 , 由于垂向隆升作用而使

16、金罗斑 大月坪复背斜形成初型 。变质流体与地幔热蚀作用的热传输集中 于金罗斑 大月坪背斜的核部 , 形成构造 - 热穹隆 , 但仍是比较封闭的构造流体系统 。因扬子 板块向北俯冲而形成的近南北向水平挤压作用 , 形成与金罗斑 大月坪背斜相配套的北东 、北 西 、东西 、南北向断裂构造 。这些配套的断裂构造成为热减阀构造 ,使封闭的构造流体系统有了 减压释放的构造通道 。沿这些构造通道迁移的变质流体有利于形成绿片岩相的动力变质岩 。绢 云母化是本区普遍发育的一种围岩蚀变和含金构造带中的蚀变现象 , 且与金矿化关系密切 , 在断裂带中有含 H2 O 流体活动时 ,钾长石 、钠长石在含水的碱性流体作

17、用下发生绢云母化 。假 设含金构造带长度 10 k m ,深度 5 k m ,蚀变宽度 10 m , 0 . 5 k m3 的含金构造带中长石经绢云母 化后加入了 453 万 t 水 , 显然在含金构造带中的成矿流体活动也是规模宏大的 。据候连珍 ( 1992) 对小秦岭金矿流体包裹体研究 , 选用 H2 O - CO2 体系 , 应用密度法查定金矿成矿压力 , 含金石英脉型金矿的成矿压力为 0 . 04220 . 0635 GPa . , 平均 0 . 056 GPa . , 按围岩静压力计- 6算属中浅成 ( 12 k m ) ,CO2 含量平均为 248 ×10。葫芦沟钾质蚀变

18、岩型金矿压力为 0 . 067- 60 . 0902 GPa . ,平均 0 . 0811GPa . ,属中深成 ( 2 . 73 . 6 k m ) ,CO2 含量为 1441 ×10。成矿流体中 N a Cl 为 1 . 8 %17 . 8 % ,平均 9 . 55 % ,属低盐度 。密度值平均 0 . 00087 k g/ m3 。成矿流体是富碱金属 N a + 及 K + ,含 CO 、SO 、Cl - 、F - 及多种金属复杂流体 。在中深成条件下 ,含金钾质2蚀变岩型金矿中 , 压性 - 压扭性剪切带中富 H2 O 、CO2 的成矿流体在应力驱动下 , 沿岩石粒间发生了弥

19、漫型的渗滤流动交代作用 。4构造成矿作用断裂构造成矿作用主要有断裂变质作用 、构造带内流体与岩石交换作用 、动力分异作用 、构造扩容作用 ,这些成矿作用是发生在特定的构造动力学系统中 。4 . 1 断裂变质作用强烈的构造应力作用引起动力变质 ,使断裂带内中 - 高级变质岩发生动力退变质作用 。伴 随着构造带内流体的产生 , 形成动力热退变质作用 , 以绿片岩相为主的退变质岩沿断裂带分布 。这种作用主要形成于压性 、压扭性 、剪性断裂带内 。主要表现为暗色矿物(角闪石 、辉石) 被 浅色矿物 (绿泥石 、绿帘石 、石英 、方解石) 替代 , 发生铁质析出 (磁铁矿 、赤铁矿 、镜铁矿 、钛铁 矿

20、) ,链状 、架状无水硅质酸盐矿物 (长石等) 被层状含水硅酸盐矿物 (绢云母 、高岭土) 所替代 。在小秦岭回马坪 - 朱家沟壳型导矿断裂带上 , 普遍发育绿泥石化 、绢云母化 、绿帘石化 、硅化 、高岭土化 、碳酸盐化 , 局部有磁铁矿化 、钛磁铁矿化 , 含星散状黄铁矿 、赤铁矿 、镜铁矿等 , 形成A u 、A g 、C u 、P b 的弱异常 。地球化学成晕作用实质是 K 、Si 、A u 、A g 、C u 、P b 富集 ,而 Ca 、M g 、Fe3 + 、Fe2 + 常被带出断裂带 ,在其附近的围岩中发生碳酸盐化 、含铁矿物沉淀 。在含金剪切带中 , 可见到剪切带或石英脉上

21、、下盘岩石中发育羽毛状裂隙和矿体分枝 , 岩 石发生绢云母化 、碳酸盐化 。沿羽毛状裂隙中充填有方铅矿细脉 、黄铁矿 、镜铁矿细脉 。这种现象在脆性较强的混合花岗岩中比片麻岩 、混合岩发育 , 它们是有剪切 - 扭动裂隙动力学系统中 ,超临界富 H2 O 、CO2 的成矿流体在剪切 - 扭动构造裂隙作用及扩容构造共同作用下产生储 矿空间 。石英脉常见到被剪切 - 扭动裂隙破碎 ,黄铁矿中发育“X”型裂理 ,沿其有方铅矿显微细 脉充填交代时 , 裂理宽度增大 。在方铅矿交代溶蚀黄铁矿后而形成港弯状结构中 , 分布有叶片 状 、枝状银金矿 ,使金的品位增高 。这种成矿成晕作用主要发生在剪切 - 扭

22、动力学系统中 ,而且在缓倾斜含金构造带比陡倾斜的含金构造带中明显发育 。在压扭性动力学系统中 ,金偏向于压力减小的部位富集 。如含金构造带 (石英脉) 走向转折处 、分枝处 , 这种动力学系统中处于相对 还原条件而引起络合物的分解 ,发生金的沉淀 。4 . 2 构造带内流体 - 岩石相互交换作用钾质蚀变岩型金矿中 , 在压性 - 压扭性动力学系统中 , 构造带内流体 - 岩石相互交换作用 十分强烈 。当成矿流体进入构造带中 ,二者成分差异很大 ,其交代作用强烈 ,迅速改变交代岩石的成分 , Si4 + 是带入量最多的离子 , 首先发生硅质交代作用 , 在构造带内形成了 0 . 20 . 8 m

23、 厚的含金石英脉 。成矿流体中高化学位的 Si4 + 、K + 、A u 等组分向外侧发生弥漫型渗滤扩散 ,因石 英脉的形成相对降低了化学位 , 故以钾质交代为主 , 形成钾长石化带 , 在蚀变交代作用最强烈的部位 ,钾长石达 80 % 。硅钾质交代作用的结果使成矿流体中与岩石的成分已基本接近 ,因此 ,成矿流体向外侧主要是活性组分 H2 O 、CO2 的交代作用 ,形成黄铁矿绢英岩化带和各类岩石的 蚀变 。富 Ca2 + 、M g2 + 、Fe2 + 、CO 、SO 的残余成矿流体最终形成方解石化 、铁白云石化 、重晶石 化 。钾长石化中碱性增强和黄铁矿化引起 A u 的沉淀 。N a +

24、进入矿区内霓霞正长岩中霓辉石 , F - 、Cl - 则进入氟碳铈矿 、磷灰石中 。在弱碱性 - 碱性 、还原条件下 ,压性 - 压扭性动学系统中 , 金的亲铁性大于亲铜性 , 流体中 C u2 + 、P b2 + 、Zn2 + 、Mo2 + 与硫结合 , 形成黄铁矿 、方铅矿 、闪锌 矿 、辉钼矿沉淀 ,并不引起 A u 富集 ,但因其硫化物与黄铁矿共生 ,所以 C u 、P b 、Zn 、Mo 异常可 以指示金的富集部位 。所形成的原生异常规模大 ,A u 异常宽度可达 50110 m ,比蚀变体规模 大 25 倍 ,比金矿体规模大 510 倍 ,元素组合为 A u 、A g 、C u 、

25、P b 、Zn 、Mo 、F 在金矿体 、含金 蚀变体上形成宽阔正异常 。4 . 3 动力分异作用地球化学S - C 组构中元素富集规律 : 对微构造 ( S - C 组构) 中的 C 面和 S 面系统采样分析发现1995 , 方维萱等 : C面中富集 P b 、A g 、Bi 、H g ,而 S 面中富集 A u 。说明 P b 、A g 、Bi 、H g 在压溶面富集 , 而 A u 则富集于张性面 ( S 面) 上 , 这种富集规律与宏观上近东西向缓倾斜的含金石英脉中易形成中型以上金矿的事实相吻 。缓倾斜的含金石英脉主要形成于层间滑剪切带( S 面) , 矿石类型主要有黄铁矿 - 石英脉

26、型 。而在陡倾斜构造中 , S - C 组构的交角变小 ,C 面发育程度 比 S 面好 ,矿石类型产要是黄铁矿 - 多金属硫化物 - 石英脉型矿 。中粗粒 、晶体弯曲 、挠曲的方铅矿含 A g 高达 549 . 0 ×10 - 6 ,中粒自形晶方铅矿含A g253 . 0 ×10 - 6 ,而细粒它形方铅矿中含 A g 167 ×10 - 6 1993 , 方维萱 。方铅矿是塑性极强的矿物 ,含 A g 高的方铅矿晶体弯曲 、挠曲显示经历应力作用 ,而细粒它形方铅矿呈脉状产于石英脉 边部 ,断口呈参差状 ,有折射劈理多组 ,反映方铅矿至少受两期应力影响 。方铅矿中

27、的银主要呈机械混入物分布于裂理 、晶粒之间 ,其解理和裂理愈发育 ,为 A g 的赋存提供空间愈多 。受应力 而细粒化后 ,导致赋存空间的减少 ,使 A g 被排出 。细粒化的方铅矿中包裹自然银 、黄铜矿 、银金 矿 、硫铜银矿 , 说明方铅矿受应力作用后而发生了 A g 的动力分异作用 , A g 形成独立矿物包含在方铅矿中 。虽然 A g 在方铅矿中发生发贫化 ,但银的独立矿物形成却使 A g 发生了再次富集 ,银金矿中含 A g 11 . 34 %25 . 90 % , 硫铜银矿中含 A g 51 . 82 % 。方铅矿的应变现象在构造带中 比较发育 ,在岳王庙地区可见到方铅矿围绕钾长石

28、 、石英碎斑形成了“”碎斑系 ,在垂直主应力方向方铅矿发生细粒化 , 形成细粒致密块状方铅矿 , 而在斜交主应力方向上形成拔丝状方铅矿 的面理效应 ,这是剪切带中受应力方铅矿的应变现象 。形成金富集和富矿体 , 而缓倾斜含金石英脉易形成层间“S”面的层滑张性扩容带 。由于成矿流体已发生了碱质交代作用 ,导致成矿流体的酸度增高 。在近地形表的构造扩容带中 ,压力降低 ,氧逸度升高 , Fe2 + 稳定性减小 , Fe3 + 增加 ,出现了镜铁矿 + 磁铁矿组组合 , 这种矿物组合常出现在陡倾斜的含金石英脉的外侧 , 不形 成金的富集 ,镜铁矿中含 A u 0 . 11 g/ t 。由于 Fe2

29、+ 浓度不足 ,C u2 + 、P b2 + 则代替 F2 + 与 S2 - 结合 ,形 成多金属型矿石 ,黄铜矿 、方铅矿的沉淀使金发生沉淀 。而在层滑张性扩容带中 ,由于顶板围岩封闭使 Fe2 + 不易发生氧化逃逸 ,黄铁矿沉淀使 A u 发生富集成矿 ,并含有较多的黄铜矿 。 笔者在工作过程中 , 得到了西北有色金属地质勘查局 712 队的支持与帮助 , 欧阳宗圻高工给予了热情指导 ,特此致谢 。参考文献 :张湘炳 . 成矿构造理论研究中的几个问题J .大地构造与成矿学 ,6 (2) :141152 .1982198519901991199219921993199319931994199

30、5陈国达 . 成矿构造研究法M .北京 :地质出版社 ( 第二版) ,166 .杨国清编著 . 构造地球化学M .桂林 :广西师范大学出版社 ,1304 .方维萱 . 陕西秦岭地区混合岩化成矿作用地球化学J , 地质与勘探 , (5) :47535 . 陈国达 . 历史 - 因果论大地构造学刍议J .大地构造与成矿学 ,16 (1) :1716 . 张湘炳 . 论构造成矿规律及其动力学机制J .大地构造与成矿学 ,16 (2) :1131227 . 方维萱 . 小秦岭含金石英脉的矿物地球化学研究J .地质与勘探 . 32 (3) :4050 .彭少海 . 粤北新洲推覆构造同推覆期构造地球化学的

31、动态学研究J .大地构造与成矿学 ,17 (1) :7181孙岩 ,戴春森 . 论构造地球化学研究J .地球化学进展 , (3) :1610 .李德威 . 构造动力成矿研究的进展J .地学前缘 , (1) :184190 .方维萱 ,樊子珲 ,兀鹏武 . 论小秦岭岳王庙地区成矿预测与找矿预测M .见 :傅水兴 ,祝新友主编 ,走向 地学新世纪 - 首届有色系统青年地质工作者学术讨论会论文集 . 北京 :冶金工业出版社 . 371376O N D Y NAM ICS O F M ETAL L O GE NIC GEO C HEM IS T RYFO R GOL D D EPOS I TSI N

32、T HE XIAO QI NL I N G A REAFA N G We i - x ua n1 ,3 , L I Ya - l in2 , HUA N G Zhua n g - y in g1 ,31 . (L O O G I n st i t u t e ofGeo c h e mi st r y , C hi n e se A c a de m y of Sci e nc e s , Gui ya n g , 5500022 ,C hi n a )2 . ( C h e n g do u U ni ve r si t y of Tec h no lo g y , C h e n g d u 610059 ,C hi n a )3 . ( No r t h w e