长蛇山地处西藏羌塘地块中部，行政区划隶属尼玛县绒马乡。新中国成立后，大批科研机构和地质工作者对羌塘地区展开地质研究，在造山带地质理论研究等基础地质方面取得重大进展。例如，邓万明等[1]曾对茶布—双湖一线的超基性和基性岩浆岩进行了初步的岩石地球化学分析，认为其不同于一般的洋中脊基性岩，可能形成于板内裂谷环境，并认为双湖—龙木错一线并不存在古特提斯的缝合带，而是一条陆内裂谷；李才等[2-6]在1/25万玛依岗日幅区调的基础上提出了龙木错—双湖为冈瓦纳大陆与欧亚大陆间的一条古特提斯缝合带；Paul Kapp、尹安等[7-8]认为，羌塘在晚古生代为统一地块，古特提斯洋位于其北部的金沙江带，羌塘中部大面积分布的变质岩具有变质核杂岩的性质。而后通过许多学者（常承法[9]；邓希光[10-11]；李才等[12]；翟庆国等[13]；黄继钧[14]）的研究和地球物理资料的证实，现在普遍的观点认为龙木错-双湖缝合带是一条晚古生代古特提斯洋缝合带，同时在该带也发现了大量的蛇绿岩套和各类脉岩出露。其中，翟庆国等对角木日地区的基性岩，进行了系统的地球化学分析（李才等[15-16]；翟庆国等[17]），认为角木日地区的基性岩形成于洋中脊或者准洋脊环境。
前人对龙木错—双湖一线出露脉岩的研究，基本确定了角木日地区的基性岩形成于洋中脊或者准洋脊环境，但针对双湖冈玛日地区的地质研究依然较少，笔者通过对该地区长蛇山辉长辉绿岩的研究和区域对比，认为北西-南东向的辉长辉绿岩形成于大陆裂谷环境，与角木日东西向的基性岩所反映的洋中脊环境有较大差别。同时，对于古特提斯洋的宽度以及闭合时限等问题，地质学家仍有较大争论，本文通过冈玛日长蛇山研究区的辉长辉绿岩的地球化学和年代学的研究，也可为该问题的解决提供依据。
1 地质背景
研究区地处青藏高原羌塘盆地的腹地，为欧亚大陆与冈瓦纳大陆的结合地带，李才[18]等通过区地质调查和研究，认同了刘本培[19]对于该地区区域地质构造的认识，即双湖缝合带以及昌宁—孟连缝合带两者之间是相互连接的关系；并在此基础上，提出羌南和羌北处于不同的地质环境。
王根厚等[20～25]通过对野外地质构造现象以及岩性认识的划分，确定了龙木错—双湖缝合带的存在，并通过填图确定了其具体分布，由此划分出羌南与羌北两处不同的地质构造系统，羌南属于一套增生杂岩体系，而羌北属于史密斯地层体系。据[26]地球物理数据显示，羌塘盆地的构造格局为“两坳夹一隆”，即冈玛错-戈木日-玛依岗日-双湖一系列条带为羌塘盆地的“一隆”，区分隔了“两坳”，即北羌塘和南羌塘坳陷。
区域内基性岩墙群发育，玛依岗日以西岩脉多为北西-南东向，以东为东西向或近南北向。基性岩墙的岩石类型与岩脉（体）的规模有关，规模大者多为辉长岩，小型脉体以辉绿岩、石英辉绿岩为主。岩墙侵入于上石炭统含冰海杂砾岩的砂板岩中，岩墙对围岩没有明显的热蚀作用。
本文重点研究区域为长蛇山地区（图2），区内基性侵入岩多为辉长岩，主要呈北西-南东向小岩脉产出，笔者在4个地质点（Y001-Y004）采集了25件辉长岩样品。
图2 长蛇山地区地质简图（改自王根厚等[23～25]）
Fig.2 Sketch geological map of Changsheshan area（after Wang Genhou.et al.[23～25]）
1.辉长岩；2.细砾岩；3.千枚岩；4.石英砂岩；5.超基性岩；6.龙格岩组；7.推俯断层；8.平移断层
2 分析测试方法
为保证数据质量，样品采集尽量选取新鲜岩块，摒除一些表面风化较为严重的岩块。样品处理和测试由河北省廊坊地质区调所实验室完成，主量、微量及稀土分析使用的检测仪器主要有AxiosmasX射线荧光光谱仪，P1245电子分析天平，50 ml滴定管，X Serise2等离子体质谱仪，光栅摄谱仪等，其中主量元素是由X射线荧光光谱法（XRF）测得，而微量元素则根据元素的不同分别使用X射线荧光谱法和等离子体质谱法（ICP-MS）测定。
3 岩石学特征
所采样品经河北廊坊区调所实验室鉴定，均为蚀变辉长岩，辉长结构，局部嵌晶含长结构，块状构造，主要矿物成分为斜长石、辉石、角闪石（图3）。斜长石含量为60%左右，呈半自形板状、板条状，杂乱分布，局部似格架状，大小为0.2～1.8 mm，帘石化明显，局部绢云母化、高岭土化。辉石含量为35%左右，呈自形-半自形柱状，为单斜辉石，分布于斜长石间，大小0.05～2 mm，少绿泥石化，局部粒内及边缘嵌晶斜长石颗粒。角闪石含量大于5%，褐绿色，多分布于辉石边缘，大小为0.05～1 mm，绿泥石化，交代辉石。
图3 长蛇山辉长岩Y001～Y004镜下显微照片
Fig.3 Microscopy of gabbro Y001～Y004 in Changsheshan
4 岩石地球化学特征
4.1 主量元素地球化学特征
岩石化学分析结果见表1。辉长岩SiO2含量在45.29%～50.32%之间，平均为47.19%；MgO含量变化范围相对较大，介于6.22%～15.33%之间，平均为11.02%；Mg#范围在53.89～72.24之间，平均64.79，少部分（9组，总数25组）符合原生岩浆范围（Mg#=68～75）[27]，说明本区域内侵入岩结晶分异作用不强烈。TiO2含量在1.41%～2.27%之间，平均为1.84%，高于中国区域内辉长岩（1.67%）；Al2O3含量较高，在8.62%～14.45%之间，平均为11.24%，CaO含量较高，在9.02%～13.11%之间，平均为11.02%。岩石中全铁含量较高，平均值TFeO为10.35%，样品Fe2O3/FeO比值均小于1，说明岩石在弱还原环境下形成[28]；Na2O含量在1.10%～3.36%之间，平均为1.72%；岩石中K 2O含量为0.54%～1.68%之间，平均值为0.99%,；全碱Na2O+K 2O，含量为1.80%～4.16%，平均值为2.71；P2O5的含量较高介于0.11%～0.21%之间。MgO/(FeO+MgO)的比值介于0.44～0.64之间，相对于原始地幔的MgO/(FeO+MgO)比值（0.68～0.72）要低，说明其原生岩浆与原始地幔有一定差别[29]，推测其在岩浆演化过程中存在地壳物质的混染熔融。
表1 长蛇山辉长岩化学成分
Table1 Chemical composition of gabbro from changsheshan mountains
数据测试单位：河北廊坊区调所实验室
Sample Y001-1 Y001-2 Y001-3 Y001-4 Y001-5 Y001-6 Y002-1 Y002-2 Y002-3 Y002-4 Y002-5 Y002-6 Y002-7 Y002-8 Y002-9 Y002-10 Y003-1 Y003-2 Y003-3 Y004-1 Y004-2 Y004-3 Y004-4 Y004-5 Y004-6 SiO2%47.49 47.15 47.63 47.48 47.64 47.62 47.91 46.96 48.11 47.19 46.25 46.36 47.88 46.24 45.29 47.11 47.11 47.00 48.04 50.32 46.61 45.83 46.55 46.87 46.99 Al2O3%10.10 9.53 10.15 8.76 8.62 10.86 9.70 9.52 10.42 10.84 11.23 9.86 12.34 10.04 11.89 12.61 12.40 12.79 13.41 13.62 11.40 13.29 11.53 14.45 11.74 TiO2%1.70 1.68 1.73 1.47 1.42 1.76 1.46 1.41 1.52 1.69 2.27 1.51 1.78 1.70 1.98 1.90 2.18 2.14 2.18 2.02 2.00 2.15 2.01 2.23 2.01 Fe2O3%1.24 1.44 1.39 1.25 1.75 1.51 1.41 1.86 1.69 1.49 1.84 1.82 1.57 1.96 1.85 2.32 2.34 1.86 1.74 2.24 1.29 1.43 1.55 2.06 1.74 FeO%8.54 8.47 8.69 8.72 8.07 8.46 9.01 8.97 8.35 8.86 9.55 8.46 8.85 8.42 9.50 8.25 8.89 9.07 8.63 7.51 9.70 9.77 9.47 8.97 9.13 CaO%12.75 13.11 11.92 12.23 12.38 12.55 9.02 9.45 10.70 10.11 11.78 12.33 10.53 12.51 11.76 11.77 10.43 10.48 11.57 9.66 11.06 10.61 11.06 10.40 10.86 MgO%11.84 12.36 11.80 14.07 14.08 10.65 14.70 15.33 12.37 12.49 10.44 13.47 9.61 12.84 10.80 8.59 9.29 9.08 7.13 6.62 10.85 9.10 10.53 7.10 10.31 K 2O%0.79 0.70 0.98 0.61 0.69 0.98 0.98 1.36 0.97 1.33 0.94 0.88 1.35 0.54 0.96 1.12 1.16 1.13 0.75 0.80 1.05 0.94 1.00 1.68 1.09 Na2O%1.49 1.35 1.45 1.19 1.16 1.55 1.57 1.18 1.77 1.41 1.57 1.10 2.10 1.37 1.49 1.98 2.02 2.21 2.39 3.36 1.52 2.28 1.68 2.09 1.69 MnO%0.191 0.185 0.198 0.191 0.184 0.194 0.264 0.202 0.209 0.236 0.226 0.205 0.249 0.194 0.196 0.295 0.200 0.204 0.193 0.175 0.194 0.193 0.197 0.191 0.194 P2O5%0.151 0.147 0.163 0.135 0.132 0.151 0.145 0.145 0.144 0.170 0.118 0.108 0.146 0.125 0.113 0.164 0.209 0.202 0.194 0.186 0.194 0.200 0.178 0.179 0.187 H2O+%2.88 2.84 2.96 2.79 2.94 2.69 2.56 2.39 2.92 3.20 2.84 0.19 2.59 3.27 3.28 3.06 2.79 2.81 2.98 2.68 3.00 3.44 3.27 2.79 3.23 H2O-%0.20 0.22 0.29 0.20 0.27 0.28 0.26 0.18 0.18 0.30 0.31 0.30 0.27 0.20 0.21 0.29 0.30 0.26 0.20 0.17 0.16 0.23 0.24 0.19 0.23 LOI%3.50 3.66 3.69 3.64 3.64 3.49 3.55 3.31 3.51 3.91 3.55 3.65 3.29 3.84 3.98 3.69 3.51 3.58 3.63 3.33 3.91 4.02 4.04 3.62 3.85 Total%99.79 99.78 99.79 99.77 99.76 99.78 99.73 99.69 99.76 99.73 99.77 99.74 99.71 99.77 99.79 99.80 99.74 99.75 99.86 99.85 99.78 99.81 99.79 99.84 99.79
在TAS图中（图4），本区样品均投在辉长岩区，基本都属于拉斑玄武岩（亚碱性玄武岩）系列。在FAM图解中，具有明显的拉斑系列演化趋势（图5）。
图4 辉长岩样品的TAS图[30-31]
Fig.4 TAS diagram of basic intrusive rock[30-31]
图5 辉长岩FAM图
Fig.5 FAM diagram of gabbro
4.2 微量元素地球化学特征
4.2.1 稀土元素地球化学特征
稀土分析数据显示，辉长岩稀土元素质量分数∑REE=77.50×10-6～144.14×10-6，平均为105.34×10-6，LREE/HREE=4.16～5.85，平均值为4.99，轻重稀土分馏较为明显，LREE富集，(La/Yb)N比值在3.75～6.48之间，平均为4.99，(Ce/Yb)N比值在3.34～5.32之间，平均为4.39。与原生地幔岩浆相近的Mg#以及较低的(Ce/Yb)N值，都一定程度上说明，该研究区内的结晶分异作用不甚强烈。
样品的δEu为0.92～1.11，平均为1.00，基本无铕异常，说明基本无斜长石的分离结晶，基本偏向于原生岩浆。
表2 稀土地球化学数据
Tab.2 Rare earth geochemical data
数据测试单位：河北廊坊区调所实验室
Sample Y001-1 Y001-2 Y001-3 Y001-4 Y001-5 Y001-6 Y002-1 Y002-2 Y002-3 Y002-4 Y002-5 Y002-6 Y002-7 Y002-8 Y002-9 Y002-10 Y003-1 Y003-2 Y003-3 Y004-1 Y004-2 Y004-3 Y004-4 Y004-5 Y004-6 La μg/g 15.94 15.56 17.05 16.69 15.88 16.27 17.36 16.58 15.99 20.38 13.56 12.73 14.22 11.25 11.91 16.30 20.39 19.06 18.54 25.67 17.97 18.69 16.17 15.46 19.41 Ce μg/g 36.58 35.20 38.66 37.21 35.84 37.15 38.17 36.41 36.32 45.03 33.52 31.81 33.81 25.90 27.40 36.19 46.66 43.42 41.24 54.46 40.30 41.52 37.59 35.67 44.30 Pr μg/g 4.56 4.43 4.86 4.64 4.47 4.73 4.62 4.40 4.51 5.53 3.89 3.78 4.34 3.43 3.59 4.55 5.88 5.38 5.21 6.56 5.20 5.31 4.83 4.52 5.65 Nd μg/g 20.78 19.95 22.26 20.73 19.75 21.75 20.28 19.66 20.17 24.44 17.61 17.72 20.14 16.27 17.43 20.76 26.53 25.01 23.54 28.20 23.73 23.93 22.05 20.91 25.44 Sm μg/g 4.94 4.67 5.20 4.89 4.72 5.09 4.57 4.36 4.69 5.56 4.38 4.33 4.92 4.23 4.32 4.98 6.04 5.68 5.54 6.28 5.47 5.78 5.31 4.96 6.01 Eu μg/g 1.60 1.51 1.65 1.48 1.41 1.67 1.37 1.43 1.40 1.65 1.48 1.43 1.75 1.40 1.49 1.70 1.92 1.88 1.83 1.94 1.81 1.78 1.66 1.73 1.95 Gd μg/g 4.63 4.45 4.88 4.63 4.51 4.80 4.36 4.28 4.56 5.13 4.19 4.13 4.58 3.90 4.18 4.82 5.83 5.35 5.29 5.94 5.20 5.43 5.00 4.75 5.65 Tb μg/g 0.78 0.74 0.81 0.76 0.73 0.83 0.74 0.73 0.78 0.86 0.71 0.71 0.77 0.67 0.72 0.80 0.99 0.91 0.90 0.96 0.85 0.90 0.85 0.80 0.95 Dy μg/g 4.95 4.75 5.10 4.92 4.82 5.04 4.62 4.49 4.87 5.32 4.66 4.36 4.92 4.49 4.54 5.04 6.01 5.66 5.53 6.04 5.43 5.66 5.30 4.90 6.01 Ho μg/g 0.93 0.89 0.98 0.94 0.90 1.00 0.88 0.86 0.91 1.05 0.88 0.82 0.94 0.85 0.88 0.95 1.16 1.08 1.04 1.16 1.03 1.11 1.03 0.95 1.12 Er μg/g 2.48 2.40 2.62 2.49 2.33 2.59 2.39 2.37 2.51 2.78 2.33 2.16 2.49 2.23 2.31 2.50 3.09 2.79 2.82 3.07 2.71 2.82 2.70 2.51 3.03 Tm μg/g 0.39 0.37 0.41 0.41 0.38 0.41 0.38 0.36 0.40 0.43 0.37 0.36 0.40 0.35 0.36 0.40 0.48 0.45 0.44 0.48 0.43 0.46 0.43 0.40 0.45 Yb μg/g 2.23 2.12 2.42 2.30 2.26 2.43 2.23 2.18 2.31 2.63 2.16 2.00 2.31 2.15 2.13 2.40 2.89 2.62 2.63 2.84 2.45 2.64 2.54 2.28 2.76 Lu μg/g 0.42 0.42 0.45 0.43 0.42 0.45 0.46 0.44 0.44 0.48 0.39 0.40 0.41 0.37 0.41 0.45 0.51 0.47 0.46 0.54 0.48 0.50 0.49 0.45 0.54
在稀土配分模式图[47]上（图6），稀土曲线呈明显的右倾，表明轻稀土相对重稀土呈较强富集。
图6 稀土元素配分模式图解[47]
Fig.6 Diagram of rare earth elements distribution model[47]
4.2.2 不相容元素地球化学特征
前人研究表明，大多数基性熔岩由于受到地壳混染等作用，具有Nb、Ta、Ti的负异常，这些可以作为它们的判别因子来确定岩石产生的构造环境[32]。
研究区的辉长岩不相容元素（Cs、Ba、Hf、Ta、W、Tl、Pb、Bi、Th、Ti、U）的分析结果见表3，结合微量元素原始地幔标准化蛛网图[47]（图7），可以看出该区域辉长岩Nb明显具有负异常，负的Nb异常是大陆地壳的特征，同时也可能是地壳岩石对岩浆的混染作用，这点结论与之前的结论相一致。
表3 微量元素地球化学数据
Table 3 Geochemical data of trace elements
数据测试单位：河北廊坊区调所实验室
Sample Y001-1 Y001-2 Y001-3 Y001-4 Y001-5 Y001-6 Y002-1 Y002-2 Y002-3 Y002-4 Y002-5 Y002-6 Y002-7 Y002-8 Y002-9 Y002-10 Y003-1 Y003-2 Y003-3 Y004-1 Y004-2 Y004-3 Y004-4 Y004-5 Y004-WL6 Sample Y001-1 Y001-2 Y001-3 Y001-4 Y001-5 Y001-6 Y002-1 Y002-2 Y002-3 Y002-4 Y002-5 Y002-6 Y002-7 Y002-8 Y002-9 Y002-10 Y003-1 Y003-2 Y003-3 Y004-1 Y004-2 Y004-3 Y004-4 Y004-5 Y004-6 Li μg/g 27.45 30.40 33.94 33.80 36.36 25.53 34.12 37.16 35.80 36.50 36.03 35.55 38.51 35.44 31.87 25.68 39.78 38.43 31.22 29.42 43.07 42.81 37.09 29.17 39.72 In μg/g 0.07 0.07 0.07 0.06 0.06 0.06 0.06 0.06 0.06 0.07 0.07 0.05 0.07 0.07 0.06 0.07 0.08 0.08 0.08 0.07 0.07 0.10 0.07 0.08 0.07 Be μg/g 0.64 0.79 0.63 0.71 0.89 0.66 1.43 0.86 0.91 1.03 0.79 0.55 0.88 0.57 0.53 0.62 0.99 1.21 1.15 1.33 0.77 0.99 0.79 0.87 0.88 Cs μg/g 1.81 2.57 4.53 3.68 3.76 4.31 2.71 5.17 1.95 3.38 2.05 2.51 3.48 1.58 1.88 1.53 1.94 1.84 1.75 1.79 3.10 2.10 2.07 3.21 2.09 Sc μg/g 30.41 35.90 32.40 33.36 35.01 32.25 25.62 27.78 32.39 25.96 32.54 30.42 31.05 33.66 32.25 28.86 28.15 27.19 26.14 23.48 27.45 30.69 26.92 24.90 29.08 Ba μg/g 261.9 277.4 364.6 198.6 286.0 421.1 380.7 467.8 480.0 424.7 626.4 452.4 1054 384.1 375.0 474.4 692.4 701.6 360.4 333.3 447.4 354.8 417.9 448.5 503.7 V μg/g 210.3 243.9 226.5 227.0 231.8 230.9 191.1 197.2 218.6 204.9 291.7 193.5 241.3 229.2 259.6 224.0 238.5 238.5 255.2 234.4 231.8 267.6 216.3 257.5 234.4 Hf μg/g 3.73 3.65 3.75 3.20 3.38 4.37 3.57 2.89 3.54 3.92 3.27 3.30 5.12 3.33 3.10 3.86 5.39 5.01 4.53 6.77 4.72 4.56 4.72 4.15 5.13 Cr μg/g 916.1 1102 769.6 1097 1295 731.6 1310 1544 1136 1082 701.7 1024 539.4 1063 725.7 397.6 551.1 538.2 177.0 257.4 543.3 499.1 548.7 71.0 604.1 Ta μg/g 0.95 1.18 1.09 1.00 1.05 1.23 1.13 1.11 1.11 1.66 1.41 0.82 1.24 0.98 0.94 1.11 1.89 1.53 1.51 1.92 1.80 1.54 1.27 1.84 1.41 Co μg/g 45.4 52.4 51.2 57.4 56.4 46.6 62.6 70.5 55.8 51.9 54.7 52.7 50.6 58.2 53.3 56.8 50.4 46.5 40.4 41.7 53.8 52.3 49.3 42.5 52.7 W μg/g 0.35 0.42 0.37 0.39 0.41 0.41 0.53 0.66 0.49 0.54 0.45 0.42 0.64 0.38 0.36 0.56 0.59 0.57 0.58 0.77 0.54 0.64 0.47 0.75 0.47 Ni μg/g 167.9 205.8 190.0 275.1 266.7 160.8 383.5 500.2 218.7 238.6 162.9 222.4 137.9 231.9 155.0 317.7 208.8 172.0 94.6 104.0 300.2 184.1 238.7 79.5 242.3 Tl μg/g 0.24 0.25 0.37 0.23 0.29 0.24 0.41 0.59 0.33 0.36 0.28 0.32 0.33 0.17 0.22 0.41 0.34 0.25 0.27 0.13 0.51 0.20 0.19 0.28 0.21 Cu μg/g 36.2 37.0 33.6 34.0 39.0 47.9 83.0 105.5 52.7 54.2 50.6 44.1 56.0 45.3 28.6 94.5 52.1 49.2 44.4 54.8 65.2 32.3 52.4 30.3 54.2 Pb μg/g 3.30 5.05 2.97 6.60 4.44 4.37 5.64 9.36 5.55 3.95 6.04 4.72 18.74 4.02 3.72 8.43 6.87 7.20 6.81 6.13 5.82 6.90 4.05 3.93 5.16 Zn μg/g 70.7 85.5 80.8 85.0 79.6 84.3 84.5 97.7 87.3 92.4 104.4 82.1 167.5 90.7 99.2 83.9 110.9 105.2 109.5 98.9 105.1 132.1 102.3 119.2 108.1 Bi μg/g 0.02 0.04 0.04 0.06 0.04 0.06 0.10 0.14 0.06 0.04 0.08 0.08 0.06 0.04 0.03 0.07 0.05 0.06 0.05 0.04 0.04 0.05 0.04 0.03 0.04 Ga μg/g 13.71 15.09 14.67 13.41 13.65 15.45 15.02 14.77 15.23 15.20 16.19 12.12 17.59 14.39 15.17 17.83 18.76 18.06 20.26 20.05 17.34 19.18 15.79 19.43 17.86 Th μg/g 2.89 3.14 3.26 3.47 3.61 3.27 4.50 4.58 4.64 4.92 3.21 2.37 2.84 2.53 2.15 3.02 4.48 4.52 4.10 7.17 3.73 3.83 3.34 3.27 3.99 Rb μg/g 21.0 22.1 28.6 21.0 24.9 31.7 32.6 53.8 35.1 49.7 31.0 25.7 41.7 16.8 32.5 42.3 33.7 33.9 23.3 27.0 36.7 33.7 30.1 54.6 36.8 U μg/g 0.59 0.70 0.66 0.70 0.77 0.64 0.85 0.96 0.85 1.00 0.56 0.59 0.59 0.57 0.48 0.69 0.91 0.90 0.79 1.44 0.75 0.80 0.72 0.72 0.81 Sr μg/g 209.7 277.5 236.9 292.4 263.9 276.0 200.7 196.0 203.3 210.4 312.1 180.5 360.9 236.8 215.9 209.5 385.6 284.4 221.2 251.0 226.6 258.3 194.5 334.7 235.2 Ag μg/g 0.03 0.03 0.04 0.03 0.04 0.06 0.08 0.09 0.08 0.06 0.06 0.07 0.06 0.05 0.02 0.12 0.05 0.04 0.06 0.04 0.03 0.03 0.04 0.02 0.05 Zr μg/g 123.3 117.9 118.1 102.1 103.9 129.0 119.5 97.7 110.9 126.6 101.5 103.0 133.0 104.3 94.7 120.9 172.7 154.4 147.4 220.3 148.1 154.9 151.0 138.7 152.1 Sn μg/g 4.5 2.9 3.1 4.9 3.3 4.8 3.8 4.0 2.4 6.6 1.6 2.4 1.8 2.7 2.2 1.5 2.1 2.0 1.0 2.3 2.2 2.5 5.5 1.4 2.2 Nb μg/g 10.66 12.70 12.06 11.16 11.05 14.08 12.77 11.93 11.83 14.56 13.24 8.10 13.37 10.43 10.49 12.56 18.47 17.18 17.15 19.76 16.30 18.05 14.24 16.11 14.39 B μg/g 14.3 8.1 8.3 8.4 13.7 14.3 35.2 22.0 8.5 33.1 16.9 21.5 23.9 19.6 13.5 29.2 14.3 14.6 13.2 14.2 11.3 14.4 17.2 21.2 10.7 Mo μg/g 0.84 0.56 0.41 0.80 0.53 0.69 0.28 0.47 0.20 0.39 0.37 0.41 0.49 0.43 0.25 0.60 0.65 0.67 0.74 0.59 0.54 0.46 0.48 0.58 0.46 Ti μg/g 10220 10071 10378 8837 8514 10567 8781 8425 9095 10140 13640 9055 10663 10190 11892 11380 13048 12817 13042 12100 11988 12901 12054 13356 12050
图7 辉长岩微量元素蛛网图[47]
Fig.7 Spidergram of gabbro trace elements[47]
大离子亲石元素（LILE）Rb、Ba等较富集，高场强元素（HFSE）Hf、Th、Ta相对较为富集，不相容元素Nb亏损。由表中微量元素计算可知，研究区辉长岩的Nb/La<1，表明其基性熔岩遭受到了比较强烈的地壳混染作用[32]。
将本区辉长岩与典型环境不相容元素比值特征做对比（表4），本区的辉长岩与大陆壳成分各不相容元素的比值最为相近，由此可以推断本区辉长岩产自陆壳环境。
表4 本区辉长岩与典型环境不相容元素比值特征对比[33]
Table 4 Comparison of ratio characteristics of incompatible elements between gabbro and typical environment in this area(According to Weaver B L[33]）
元素比原始地幔N-MORB大陆壳HIMUOIB EMIOIB EMIIOIB OIB本区辉长岩Zr/Nb 14.8 30 16.2 3.2～5.0 4.2～11.5 4.5～7.3 5.83 9.54 La/Nb 0.94 1.07 2.2 0.66～0.77 0.86～1.19 0.97～1.09 0.77 1.24 Ba/Nb 9 4.3 54 4.9～6.5 11.4～17.8 7.3～11.0-33.47 Ba/Th 77 60 124 49～77 103～154 67～84-128.73 Rb/Nb 0.91 0.36 4.7 0.35～0.38 0.88～1.17 0.59～0.85 0.65 2.45 K/Nb 323 296 134 77～179 213～432 248～378--Th/Nb 0.117 0.071 10.44 0.080～0.101 0.105～0.122 0.111～0.157 0.08 0.27 Th/La 0.125 0.067 0.204 0.107～0.133 0.107～0.128 0.122～0.163 0.11 0.22 Ba/La 9.6 4 25 6.8～8.7 13.2～16.9 8.3～11.3 9.46 27.64
5 讨论
5.1 构造环境分析
由于地壳混染作用对于岩石中Zr和Y的含量影响不大，所以通过Zr/Y-Zr图解，判断出研究区的大致构造环境[34-35]。根据研究区所采样的微量元素，作图如下，从中可看到Y001-Y004地质点样品均为板内玄武岩环境（图8）。
图8 Zr/Y-Zr构造环境判别图[36]
Fig.8 Zr/Y-Zr tectonic environment discrimination diagram[36]
在辉长岩Zr-Nb-Y图解（图9）上，其点也基本上在C区和B区，即板内拉斑玄武岩区和E-MORB区，结合Zr/Y-Zr图解，判定其为板内拉斑玄武岩区域。
图9 Nb＊2-Zr/4-Y构造环境判别图[37]
Fig.9 Nb＊2-Zr/4-Y tectonic environment discrimination diagram[37]
AⅠ.板内碱性玄武岩；AⅡ.板内碱性玄武岩和板内拉斑玄武岩；B.E型MORB；C.板内拉斑玄武岩和火山弧玄武岩；D.N型MORB和火山弧玄武岩
在Ti-Zr-Y图解（图10）中，可以看出所投环境地质点基本落于D区，即板内玄武岩环境，此环境与之前的投图所指示的环境相一致。
图1 研究区构造单元划分图（据王根厚（2012年）[20～22]）
Fig.1 Division of tectonic units in the study area[20～22]
二个一级构造单元：I.羌北-昌都板块，属欧亚大陆构造域；II.羌南-保山板块，属冈瓦纳构造域；二个二级构造单元：I-1.羌北地块；II-2.羌南增生杂岩带；六个三级构造单元：I-1-1.俄久卖糜棱岩化片麻岩；I-1-2.羌北晚二叠世-中三叠世、晚侏罗世沉积地层；II-2-1.增生杂岩基质，片理化石英-石英、千枚岩；II-2-2.刚性洋岛岩块；II-2-3.高压变质带；II-2-4.中二叠世碳酸盐岩岩层；N1-2—中新世—上新世陆相沉积。
图10 Ti-Zr-Y图解[35]
Fig.10 Ti-Zr-Y diagram[35]
A.岛弧拉斑玄武岩；B.MORB、岛弧拉斑玄武岩和钙碱性玄武岩C.钙碱性玄武岩；D.板内玄武岩
根据以上图解以及微量元素分析，可知研究区的辉长岩基本属于陆壳内构造环境，其间可能伴有围岩混染，致使微量元素有轻微异常。
在研究大地构造背景时，需要分析区内辉长岩的源区性质，并研究地壳物质对岩浆岩的影响。岩石地球化学研究证实，如果铁镁质岩浆受来自陆壳的混染，其大离子亲石元素及一些主量元素含量相对较高，同时Zr/Nb、La/Nb比值较高，Ce/Pb、Nb/U比值相对较低[38]。典型地幔：Ce/Pb=25±5，地壳：Ce/Pb＜15[39]，本研究区辉长岩的Ce/Pb比值为1.80～13.03之间，平均值为7.27，从而得知其为典型的地壳环境，与前文所得结果相同。根据前人研究成果，当岩浆源区遭受来自俯冲板片来源的流体作用或地壳物质的混染时都可能使得Ce/Pb比值较低，受来自俯冲消减板片脱水或熔融流体流体改造的地幔，它们的Ta、Zr、Hf和Nb将发生不同程度的分异作用[40]。
由赵崇贺[41]的中基性火山岩成分ATK图解（图11），可以看出本研究区的所有辉长岩均产出于大陆玄武岩环境，结合之前所得结论，可得出研究区内辉长岩形成环境为大陆板内玄武岩。板内玄武岩广泛定义为产出于板块内部，远离板块边缘的地幔源火山岩，其主要分为拉斑玄武岩及碱性玄武岩。本研究区的辉长岩经微量地球化学投图，显示其主要成分为拉斑玄武岩系列。大洋内的板内玄武岩，主要构成洋岛、海底高原及海山。而本区的大陆板内玄武岩，首先为大陆溢流玄武岩，其主要组成为拉斑玄武岩，大陆溢流型玄武岩为陆上大火成岩省；其次是大陆碱性玄武岩带，其主要成分以碱性玄武岩为主，分布较为分散，基本每个大陆都有分布，但出露面积不大[42]。
图11 TiO2-Al2O 3-K2O图解[41]
Fig.11 TiO2-Al2O3-K2O diagram[41]
I.大洋玄武岩；II.大陆玄武岩、安山岩；III.岛弧、造山带玄武岩、安山岩
研究表明板块内部构造，主要分为大洋壳以及大陆壳两类，大洋壳为大洋板块，大陆壳为大陆板块，板内构造，顾名思义，为构造上构造相对稳定地区。本区辉长岩属于大陆板块，主要由拉斑玄武岩系列组成。大陆板内玄武岩主要有大陆裂谷和陆内热柱两种成因，而陆内热柱属于碱性玄武岩，结合研究区内辉长岩地球化学数据，Th/Ta为1.77～4.19，平均值为2.89，Ta/Hf为0.24～0.44，平均值为0.32，而大陆板内玄武岩Th/Ta>1.6，Ta/Hf>0.1,且典型裂谷玄武岩的Th/Ta值为1.6～4[43]，说明研究区内辉长岩属于典型的裂谷玄武岩（图12）。
图12 Th/Hf-Ta/Hf构造环境判别图解[43]
Fig.12 Th/Hf-Ta/Hf diagram[43]
通过对印度洋洋中脊玄武岩的研究，Le Roex等[44]认为Nb、Zr和Y的丰度可以直接反应地幔源类型，富集地幔的Zr/Y比值低于18，而亏损地幔具有大于18。在本研究区中，Zr/Y为4.09～7.02，平均值为4.94，Zr/Nb为7.66～12.72，平均值为9.54，由此可知该岩浆岩为富集型地幔源。
5.2 形成时代
本研究区共从Y001-Y004四个点位各采集大块岩体，在河北廊坊地质调查研究所制取人工重砂，由于辉长岩中锆石获取较难，所以所获得锆石数量有限，通过阴极发光、透射光、反射光图像，发现本次采集的锆石环带不清晰，且多具有核幔边结构，为捕掳锆石，仅有几粒锆石具有环带结构，但较为破碎，无法完成激光测年实验。因李才、翟庆国等在本区或相邻区域采集了相同岩性的岩体进行年龄分析实验，故本文使用其年龄分析结果作为参考。
本区域内所采集的辉长岩，宏观上看，其属于一系列大的基性岩墙群的部分组成，该岩墙群整体延伸方向与龙木错-双湖缝合带相同，且分布较为广泛，延伸亦较远，故推测岩墙群相同延伸方向在同一构造环境下产生。
李才[4]等在玛依岗日地区进行区调工作时，对区域内的基性岩墙群有一定的认识了解，并采取样品分析该区域内基性岩墙群的年龄，其中锆石U-Pb年龄为312±4 Ma，全岩Sm-Nd年龄为299±13 Ma、314±5 Ma，并认定此结果为整个基性岩墙群的年龄。
翟庆国[45]对羌塘这一套基性岩墙群进行进一步的研究分析，通过锆石激光测年得出羌塘中部的基性岩墙群年代为283 Ma。相比于其之前[46]采集戈木日地区得出的加权平均年龄284±3 Ma，数据较为接近可靠。
综上所述，由于长蛇山地区辉长岩属于基性岩墙群的一小部分，故结合前人的研究成果以及野外地质观察，推测其形成时代应与李才等人的研究结果相同，年龄300 Ma左右，即晚石炭-早二叠之间。
5.3 构造意义
对于该区域的野外特征、室内研究以及其地球化学特征在前文中都有了较为详细的论述，进一步探讨其构造环境背景以及在羌塘构造体中的意义。
冈玛日图幅内的基性侵入岩墙，主要分为东西向和北西-南东向两个方向，野外工作时，笔者发现两者的野外特征有所区别，根据所采样品进行的地化分析显示，它们的化学组分也有所差异，玛依岗日地区内东西向基性侵入岩构造投图为洋中脊环境，而本研究区基性侵入岩（辉长岩）投图为板内大陆裂谷环境，洋中脊环境说明其处于成熟的大洋时期，可能为威尔逊旋回的成熟或衰退时期，而大陆裂谷环境为大陆形成的威尔逊旋回的萌芽阶段，两个方向的基性侵入岩应为不同时期的产物，它们的岩浆源亦为不同阶段的岩浆源。
长蛇山地区的辉长岩岩体位于玛依岗日雪山峰南部，增生杂岩体系内，该处的基性侵入岩规模较大，且分布广泛，研究区的辉长岩岩体只是整个青藏高原中的小部分，通过地化填图得知，它的形成环境为板内环境，且为大陆板内环境，亦板内裂谷。
本人通过阅读相关文献，结合整个区域地质演化史，认为长蛇山地区的辉长岩，其所反应的大陆裂谷环境，有如下可能：
（1）其发育环境为中特提斯洋的萌芽阶段，在该阶段，古特提斯洋已经完成闭合，大陆重新处于拉伸环境下，地壳减薄，地幔基性岩浆岩随之喷发而出，其形成时代厘定了中特提斯洋的开启。
（2）由于与其相近的东西向辉长岩墙，其所反映的构造环境为洋中脊环境，由此可以推断，处于大陆裂谷时期的基性岩喷出后，经过后期地质改造，地壳抬升等一系列过程，使得两处不同构造环境的辉长岩在现阶段基本处于同一区域。
（3）关于特提斯洋系列发育的时代及演变期次，在地质界一直都有较大争议，本文认为长蛇山地区的基性侵入岩岩墙与大洋形成的初始大陆拉张减薄环境相关。
6 结论
（1）长蛇山地区北西-南东向的辉长岩MgO/(FeO+MgO)的比值介于0.44～0.64（53.9～72.2）之间，该比值相对于原始地幔的MgO/(FeO+MgO)比值（0.68～0.72）要低，熔浆成份与原生岩浆有差异，推测在岩浆演化过程中有一定程度的结晶分异或地壳混染。
（2）研究区辉长岩轻重稀土总量分数比值：LREE/HREE=4.16～5.85，平均为4.99，其分馏作用较为明显，LREE富集。本区的辉长岩与大陆壳成分（富集Pb，亏损Nb、Ta、P）最为接近，本文推断长蛇山地区辉长岩产自陆壳环境。
（3）长蛇山地区辉长岩构造环境判别形成于大陆裂谷环境，岩浆源来自富集型地幔。结合前人研究基础，该区域岩墙群形成于威尔逊旋回的萌芽阶段，属于大陆伸展背景下的初始裂谷环境，对龙木错-双湖缝合带的构造演化具有重要的制约作用，有助于深化对特提斯洋演化的认识。
[1]邓万明，尹集祥，呙中平.羌塘茶布一双湖地区基性超基性岩和火山岩研究[J].中国科学,1996,26(4)：296-301.
[2]李才.龙木错一双湖一澜沧江板块缝合带与石炭二叠系冈瓦纳北界[J]..长春地质学院学报.1987，17（2）：155-166
[3]李才，程立人，胡克，等.西藏龙木错-双湖古特提斯缝合线研究[M].北京：地质出版社.1995，23（3）：268-274
[4]李才，程立人，于介江，等.1：250000玛依岗日幅（I45C003002）区域地质调查报告[R].吉林大学地质调查研究院,2006,1-205.
[5]李才，董永胜，翟庆国，等.青藏高原羌塘早古生代蛇绿岩——堆晶辉长岩的锆石SHRIMP定年及其意义[J].岩石学报,2008,24(1)：31-36.
[6]李才，翟刚毅，王立全，等.认识青藏高原的重要窗口-羌塘地区近年来研究进展评述（代序）[J].地质通报，2009,28（9）：1170-1173.
[7]KAPP P,YIN A,MANNING C E,et al.Blueschist-bearing metamorphic core complexes in the Qiangtang block reveal deep crustal structure of northern Tibet[J].Geology,2000,28(1)：19-22.
[8]KAPPP,YIN A,MANNINGCE,et al.Tectonic evolution of the early Mesozoic blueschist-bearing Qiangtang metamorphic belt,central Tibet[J].Tectonics,2003,22(4).
[9]常承法，青藏高原的板块构造[J].矿物岩石地球化学通讯，1985,（2）：1007-2802
[10]邓希光，丁林，刘小汉.藏北羌塘中部冈玛日-桃形错蓝片岩的发现[J].地质科学,2000,35(2)：227-232.
[11]邓希光，丁林，刘小汉.青藏高原羌塘中部蓝片岩的地球化学特征及其构造意义[J].岩石学报,2002,18(4)：517-525.
[12]李才，黄小鹏，翟庆国，等.龙木错—双湖—吉塘板块缝合带与青藏高原冈瓦纳北界[J].地学前缘,2006,04：136-147.
[13]翟庆国，李才，王军，等.藏北羌塘地区基性岩墙群锆石SHRIMP定年及Hf同位素特征[J].科学通报,2009(021)：3331-3337.
[14]黄继钧.藏北羌塘盆地构造特征及演化[J].中国区域地质,2001,20(2)：178-186.
[15]李才，翟庆国，董永胜，等.青藏高原龙木错-双湖板块缝合带与羌塘古特提斯洋演化记录[J].地质通报，2007，26(1)：13-21.
[16]李才，翟庆国，陈文，等.青藏高原龙木错-双湖板块缝合带闭合的年代学依据—来自果干加年山蛇绿岩与流纹岩Ar-Ar和SHRIMP年龄制约[J].岩石学报,2007，26(5)：911-918.
[17]翟庆国，李才，黄小鹏.西藏羌塘中部角木日地区二叠系玄武岩的地球化学特征及其构造意义[J].地质通报,2006,25(12)：1419-1427.
[18]李才，翟庆国，程立人.青藏高原羌塘地区几个关键地质问题的思考.地质通报，2005,24（4）：295-301
[19]刘本培，冯庆来，C.chonglakmani,等.滇西古特提斯多岛洋的结构及其北西-南东延伸[J].地学前缘，2002，9(3)：163-164.
[20]王根厚，韩芳林，杨运军，等.藏北羌塘中部晚古生代基性超基性岩的发现及其地质意义[J].地质通报，2009，28(9)：1181-1187.
[21]王根厚，梁定益，张维杰，等.藏东北构造古地理特征及冈瓦纳北界的时空转换[J].地质通报.2007，26（8）：921-928.
[22]王根厚，周详.喜马拉雅造山带变质杂岩表露机制[J].地质力学学报，1996，2（3）：27-28.
[23]王根厚，周详，曾庆高，等.西藏东部嘉玉桥变质杂岩中片理化岩组主期面理形成及构造意义[J].中国区域地质，1998，17（2）：177-180.
[24]王根厚，贾建称，万永平，等.藏东巴青县北部酉西岩组构造片理形成及构造意义[J].地学前缘，2006,13（4）：181-185.
[25]王根厚，张维杰，周详，等.西藏东部嘉玉桥变质杂岩内中侏罗世高压剪切作用：来自多硅白云母的证据[J].2008,24（2）：396-399.
[26]卢占武，高锐，李永铁，等.青藏高原羌塘盆地基底结构与南北向变化——基于一条270km反射地震剖面的认识[J].岩石学报2011（11）：618-631.
[27]朱弟成，潘桂棠，莫宣学，等.特提斯喜马拉雅带中段桑秀组玄武岩的地球化学和岩石成因[J].地球化学，2005（1）：588-594.
[28]LASSITERJCANDDEPAOLODJ.Plume/lithosphere interaction in the generation of continental and oceanic flood basalts：chemical and isotope constraints.In：Mahoney J,ed.Large Igneous Provinces：Continental,Oceanic,and Planetary Flood Volcanism.Geophysical Monography 100[J].American Geophysical Union,1997,335～355.
[29]钟辉，傅俊彧.塔源地区变辉长岩岩石化学地球化学特征及成因[J].地质与资源，2006（1）：42-47.
[30]Wilson M.Igneous petrogenesis(Global tectonic approach)[M].London：Oxford University Press.1989.
[31]COX KG,BELL JD,PARKHURSTRJ.Theinterpretation of igneousrocks[M].London：George Allen and unwin ltd.1979.
[32]夏林圻，夏祖春，徐学义，等利用地球化学方法判别大陆玄武岩和岛弧玄武岩[J].岩石矿物学杂志，2007（1）：77-89.
[33]WEAVER B L,Trace element evidence for the origin of ocean-island basalts,Geology[J],1991(19),123-126.
[34]PEARCE J.A.,CANN R.Ophiolite origin investigated by discriminant analysis using Ti,Zr and Y[J].Earth and Planetary Science Letters,1971(12)：339-349.
[35]PEARCE J.A.,CANN R.Tectonic setting of basic volcanic rocks determined using trace element analyses[J].Earth and Planetary Science Letters,1973(19)：290-300.
[36]PEARCE J.A.and NORRY M.J.Petrogenetic implications of Ti,Zr,Y and Nb variations in volcanic rocks,Contrib[J].Mineral.Petrol.,1979(69)：33-47.
[37]MESCHEDEM.A method of discriminating between different types of mid-ocean ridge basalts and continental tholeiites with the Nb-Zr-Y diagram[J].Chemical Geology,1986(56);207-218.
[38]CAMPBELL IH and GRIFFITHS RW.The evolution of mantles chemical structure[J].Lithos,1993(30)：389-399.
[39]HOFMANN A W.Chemical differentiation of the Earth：the relationship between mantle,continental crust,and oceanic crust[J].Earth and Planet.1988
[40]STOLZA Z,JOCHUM K P,SPETTEL B,Hofmann A W.Fluid-and melt-related enrichment in the subarc mantle：evidence from Nb/Ta variations in island-arc basalts[J].The Journal of Geology.1996
[41]赵崇贺.中基性火山岩成分的ATK图解与构造环境[J].地质科技情报，1989（4）：1-5.
[42]FARMER G L.Continental basaltic rocks[G].//Treatise on Geochemistry,Elsevier Pergamon,Oxford：2003.85-121.
[43]汪云亮，张成江，修淑芝.玄武岩类形成的大地构造环境的Th/Hf-Ta/Hf图解判别[J].岩石学报，2001（3）：413-421.
[44]LEROEX,A PDICK,BRLANK H JB,et al.Geochemistry,mineralogy and petrogensis of lavas erupts along the southwest Indian Ridge between the Bouvet tripe junction and 11 degree east[J].Journal of Petrology,1983,24：267-318.
[45]翟庆国.青藏高原羌塘中部三叠纪高压变质带[A].中国地质学会.第七届世界华人地质科学研讨会摘要集[C].中国地质学会：中国地质学会地质学报编辑部,2013：1.
[46]翟庆国，李才，王军，等.藏北羌塘地区基性岩墙群锆石SHRIMP定年及Hf同位素特征[J].科学通报,2009,54(21)：3331-3337.
[47]SUN SS,MCDONOUGH W F.Chemical and isotopic systematics of oceanic basalts：implications for mantle composition and processes[J].Geological Society,London,Special Publications,1989,42(1)：313-345.