Suitability assessment of geological sequestration of CO2 in Songliao Basin based on gray relational analysis method
-
摘要:
松辽盆地是中国主要的能源基地，随着中国振兴东北老工业基地政策的实施，如何应对工业和经济发展带来的CO2排放量的增加，将成为亟需解决的问题。在盆地内进行CO2地质储存将成为盆地内解决该问题的主要技术选择之一。以松辽盆地为研究对象，在深入分析松辽盆地CO2地质储存条件的基础上，从储存规模、储存安全性和社会环境风险及经济适宜性4个方面综合考虑，构建了由16个指标组成的松辽盆地CO2地质储存适宜性评价指标体系，利用灰色关联分析法对松辽盆地33个二级构造单元进行了盆地级适宜性评价工作，评价结果为进一步场地级目标靶区的筛选提供了科技支撑。
Abstract:Songliao Basin is one of China's major energy bases. With the implementation of the policy of revitalizing the old industrial bases in Northeast region, how to deal with the increase of CO2 emissions arising from industrial and economic development has become an urgent problem. Performing CO2 geological sequestration in this basin will be one of the major techniques to solve the problem. Therefore, with Songliao Basin as the study object and in consideration of the geological conditions of Songliao Basin, the authors selected indexes from the four aspects which include storage's size, storage's security, social environmental risk, economic suitability, and established an index system that includes 16 indicators by applying the gray relational analysis method to the basin level suitability assessment of the 33 Ⅱ tectonic units of Songliao Basin. The results provide a basis for target area selection of ground level.
-
新生代印度-欧亚大陆的碰撞汇聚作用，导致晚始新世—中新世喜马拉雅地壳物质发生大规模的熔融作用，并形成巨量的淡色花岗岩、花岗岩体[1]。这些岩体沿碰撞造山作用有关的伸展构造分布[2-11]，并与造山带演化过程密切相关，是探索构造-岩浆相互作用关系和机制的重要“岩石探针”[12]。
在喜马拉雅山脉的脊部，自西向东出露一系列不连续的淡色花岗岩侵入体，岩体的分布受区域性伸展构造藏南拆离系（STDS）控制，并构成STDS韧性剪切带的一部分[13-14]，因此，又称其为STDS淡色花岗岩。前人对造山带西部的淡色花岗岩进行了较多的年代学和成因机制方面的研究，取得的主要认识和存在的争议包括：①淡色花岗岩的时代，基本厘定STDS淡色花岗岩的年代格架，时代集中于36~13Ma[2, 4-5, 7, 9, 15-22]，暗示长时间的地壳深熔作用；但多数年龄来自喜马拉雅中西部地区的岩体，缺少东部岩体年龄的补充和限定。②淡色花岗岩的源区和源岩，大多数学者认为，STDS淡色花岗岩源区和源岩为大喜马拉雅变质沉积岩[23-26]；而杨晓松等[27]根据部分熔融实验研究，认为源岩为大喜马拉雅中的黑云斜长片麻岩；Guo等[28]通过系统的地球化学和实验岩石学模拟研究，认为源区为大喜马拉雅变质沉积岩和来自小喜马拉雅流体的两端元混合区。③淡色花岗岩的形成机制，淡色花岗岩源自白云母脱水熔融这一结论已基本达成共识[29-30]，但对触发部分熔融的因素仍有争议；Le Fort等[31]认为，流体（如水）是导致部分熔融的因素，然而，实验岩石学研究表明，在无流体参与时，部分熔融行为仍可进行[24]；Har⁃ris等[25]、Davidson等[32]认为，构造减压是部分熔融的直接触发因素。后期，又有学者提出断裂的剪切摩擦生热和放射性同位素生热可以导致部分熔融[17, 30-33]，但理论计算表明，单纯的剪切摩擦热或地壳放射性元素生热难以触发大规模的岩浆作用[34]。由此可见，STDS淡色花岗岩的研究仍存在不足，已取得的部分认识多基于对造山带中西部岩体的分析，而缺少来自东部岩体的年代学、地球化学等的约束。
本次研究采集喜马拉雅造山带东部错那地区淡色花岗岩样品，进行锆石U-Pb年龄、全岩地球化学和Sr-Nd同位素研究，为STDS淡色花岗岩研究提供了新的数据，以期厘定深熔作用时代，查明其源区和源岩，并探讨成因机制和构造意义。
1. 区域地质背景
狭义的喜马拉雅造山带指雅鲁藏布江缝合带与主前锋逆冲断裂（MFT）之间的、由新生代印度-欧亚大陆碰撞形成的强烈变形、变质带。造山带自北向南发育的构造-岩石单元依次为特提斯喜马拉雅沉积系（THS）、STDS、大喜马拉雅结晶杂岩（GHC）、主中央逆冲断裂（MCT）、小喜马拉雅沉积岩系（LHS）、主边界逆冲断裂（MBT）、西瓦里克前陆盆地沉积（SS）和MFT（图 1）[35]。THS位于最北部，南北分别以STDS和雅鲁藏布江缝合带为界，主要发育一套古生代—始新世、经历极低级变质的硅质碎屑岩和碳酸盐岩沉积，在中部自西向东分布一系列片麻岩穹窿——北喜马拉雅片麻岩穹窿（NH⁃GD），是北喜马拉雅伸展构造的重要组成部分[36]。GHC位于MCT和STDS之间，由一套源岩时代为古元古代—奥陶纪的中高级变质结晶杂岩系组成[37]，上部紧邻STDS处出露大量淡色花岗岩。LHS位于MBT和MCT之间，主要由碎屑沉积岩和低级变质岩组成。最南部为SS，由古近纪—中新世的海相和陆相沉积组成。
![]()
图 1 喜马拉雅造山带中东段地质简图（据参考文献[35]修改）GCT—大反冲断层；STDS—藏南拆离系；MCT—主中央逆冲断裂；MBT—主边界逆冲断裂；MFT—主前锋逆冲断裂；NHGD—北喜马拉雅片麻岩穹窿Figure 1. Sketch geological map of the middle and east Himalayan orogen研究区位于喜马拉雅造山带东部的错那地区（图 1）。研究区以STDS为界，北部为THS，由三叠纪—白垩纪沉积岩组成，主要岩性为砂岩、粉砂岩、泥页岩，以及低级变质的板岩和千枚岩（图 2）。南部为GHC中高级变质的肉切村群和聂拉木群，其中，肉切村群主要岩石包括钙质石英片岩、泥质片岩、大理岩等（图 2），内部发育少量切层的淡色花岗岩脉（图 3-a、3-b）。聂拉木群主要由黑云斜长片麻岩、透辉大理岩、含石榴子石矽线黑云斜长变粒岩、含矽线黑云石英片岩、石英岩等组成。在STDS附近出露淡色花岗岩体（图 2、图 3-c），岩石整体呈灰白色，显微镜下观察主要矿物成分为石英（40%）、斜长石（25%）、钾长石（20%）、白云母（10%），以及少量的黑云母和石榴子石（图 3-d）。
![]()
图 3 藏南错那地区主要岩石野外及显微照片a—泥质片岩及侵入其中的淡色花岗岩脉；b—大理岩及侵入其中的淡色花岗岩脉；c—淡色花岗岩野外露头；d—淡色花岗岩矿物组合，正交偏光。Qtz—石英；Pl—斜长石；Kfs—钾长石；Mus—白云母；Bi—黑云母；Grt—石榴子石Figure 3. Photographs and photomicrograph of the rocks in Cuona area, southern Tibet2. 实验方法
2.1 LA-MC-ICP-MS锆石U-Pb定年
采用常规方法分选出锆石单矿物并制作成环氧树脂靶，打磨靶面使锆石中心暴露，拍摄锆石的可见光和阴极发光显微照片，然后进行LA-MCICP-MS锆石U-Pb同位素测试。锆石阴极发光图像（CL）在北京大学造山带与地壳演化教育部重点实验室的扫描电镜上完成。LA-MC-ICP-MS锆石U-Pb同位素测试在天津地质矿产研究所同位素实验室进行。所用仪器为Thermo Fisher公司生产的Neptune型多接受器等离子体质谱仪（MC-ICPMS），并结合美国ESL公司生产的UP193-FX ArF准分子激光器，激光剥蚀束斑直径为35μm，激光能量密度为2.5J/cm2，频率为10Hz，激光剥蚀物质以氦为载气送入Neptune的电感耦合等离子体，GJ-1作为外部锆石年龄标准进行U-Pb同位素分馏校正。在测试过程中，每测定7个样品点之前、之后各测定2次锆石标样GJ-1。数据的离线处理采用ICPMS⁃DataCal程序[38]完成，锆石U-Pb谐和图用Isoplot/ Ex3.00程序绘制[39]。
2.2 主量和微量元素
野外采集的样品挑选较新鲜、无蚀变的部分碎至200目，在加拿大ACME实验室完成全岩主量和微量元素测定。主量元素采用电感耦合等离子体发射光谱仪（ICP-ES）测定，测试精度优于5%；微量元素采用Agilent 7500a型电感耦合等离子体质谱仪（ICP-MS）测定，测试精度优于5%~10%。
2.3 Sr-Nd同位素
Rb-Sr、Sm-Nd同位素测试样品的化学分离在北京大学造山带与地壳演化教育部重点实验室完成，Sr、Nd同位素比值测定在天津地质矿产研究所热电离质谱仪TRITON上完成。同位素比值测试过程中，扇形磁分析器的有效半径为81cm，加速电压为10kV时分析质量数范围为3~320amu，分辨率大于450（10%峰谷定义），灵敏度大于3ion/100μmol，丰度灵敏度≤10×10-9。在样品测试过程中，测定的JN-DI Nd同位素标样和NBS-987 Sr同位素标样的同位素值分别为143Nd/144Nd=0.512107±0.000004 （2σ）和87Sr/86Sr=0.710265±0.000004 （2σ）。
3. 测试结果
3.1 LA-MC-ICP-MS锆石U-Pb年龄
错那淡色花岗岩样品CN-06的锆石阴极发光图像显示，锆石颗粒自形，长柱状，棱角清晰，颗粒长度为200~400μm，长宽比介于2∶1~4∶1之间。大多数锆石发育核-幔-边结构，核部和幔部为继承锆石，核部色调较暗，幔部形状不规则，色调灰白；边部宽度为40~60μm，色调黑暗并显示密集的韵律生长环带，表明为岩浆成因。少数锆石仅发育核-边结构，核部色调较亮，围绕核部发育具有韵律环带的暗色生长边（图 4-a）。
![]()
图 4 藏南错那淡色花岗岩代表性锆石CL图像(a)和U-Pb谐和图(b)Figure 4. CL images of representative zircons (a) and zircon U-Pb concordia diagram (b) of Cuona leucogranite, southern Tibet对样品CN-06共选取23个点位进行U-Pb同位素测试，全部分析点均位于锆石边部的暗色韵律环带上，测试数据列于表 1。23个分析点的Th和U含量分别为187×10-6~509×10-6和8186×10-6~8970×10-6，Th/U值极低，为0.02~0.06，均不大于0.10，与深熔成因的锆石特征吻合[40]。23个分析点的206Pb/ 238U年龄介于18.9~16.4Ma之间，给出的年龄加权平均值为17.7±0.3Ma（MSWD=2.0）（图 4-b），代表错那淡色花岗岩的最终结晶年龄。
表 1 藏南错那淡色花岗岩样品（CN-06）LA-MC-ICP-MS锆石U-Th-Pb同位素测试数据Table 1. LA-MC-ICP-MS zircon U-Th-Pb isotopic data of the leucogranite sample CN-06 from Cuona, southern Tibet分析点号 Pb Th U Th/U 同位素比值 表观年龄/Ma 10-6 207Pb/235U ±1σ 206Pb/238U ±1σ 208Pb/232Th ±1σ 207Pb/235U ±1σ 206Pb/238U ±1σ 208Pb/232Th ±1σ 01 2598 339 8357 0.04 0.0170 0.0009 0.00255 0.00013 0.00108 0.00006 17.1 0.9 16.4 0.8 21.9 1.3 02 2687 208 8481 0.02 0.0177 0.0009 0.00275 0.00014 0.00087 0.00010 17.8 0.9 17.7 0.9 17.6 2.0 03 2716 213 8475 0.03 0.0181 0.0010 0.00279 0.00014 0.00086 0.00009 18.2 0.9 17.9 0.9 17.5 1.9 04 2889 226 8970 0.03 0.0175 0.0009 0.00267 0.00013 0.00085 0.00007 17.7 0.9 17.2 0.9 17.3 1.3 05 2612 275 8418 0.03 0.0174 0.0009 0.00270 0.00014 0.00086 0.00007 17.5 0.9 17.4 0.9 17.5 1.4 06 2944 309 8375 0.04 0.0201 0.0010 0.00284 0.00014 0.00177 0.00011 20.2 1.0 18.3 0.9 35.7 2.3 07 2826 187 8499 0.02 0.0185 0.0011 0.00287 0.00015 0.00166 0.00027 18.6 1.1 18.5 1.0 33.6 5.4 08 2630 298 8395 0.04 0.0177 0.0009 0.00275 0.00014 0.00083 0.00006 17.9 0.9 17.7 0.9 16.7 1.1 09 2634 350 8344 0.04 0.0180 0.0009 0.00267 0.00014 0.00103 0.00007 18.1 0.9 17.2 0.9 20.9 1.3 10 2747 207 8481 0.02 0.0182 0.0010 0.00283 0.00014 0.00111 0.00010 18.3 1.0 18.2 0.9 22.4 2.0 11 3001 250 8433 0.03 0.0202 0.0011 0.00277 0.00014 0.00284 0.00022 20.3 1.1 17.9 0.9 57.3 4.5 12 2707 439 8254 0.05 0.0181 0.0010 0.00264 0.00014 0.00136 0.00010 18.2 1.0 17.0 0.9 27.4 2.1 13 2834 371 8321 0.04 0.0195 0.0010 0.00257 0.00013 0.00232 0.00015 19.6 1.0 16.5 0.9 46.8 2.9 14 2721 354 8336 0.04 0.0181 0.0009 0.00274 0.00014 0.00102 0.00007 18.2 0.9 17.6 0.9 20.7 1.4 15 2719 193 8496 0.02 0.0182 0.0009 0.00277 0.00014 0.00117 0.00012 18.3 0.9 17.8 0.9 23.7 2.4 16 3128 225 8456 0.03 0.0192 0.0019 0.00275 0.00014 0.00430 0.00048 19.3 1.9 17.7 0.9 86.7 9.6 17 2612 192 8501 0.02 0.0174 0.0010 0.00271 0.00014 0.00090 0.00025 17.5 1.0 17.4 0.9 18.3 5.0 18 2699 337 8354 0.04 0.0185 0.0010 0.00276 0.00014 0.00106 0.00008 18.6 1.0 17.8 0.9 21.4 1.6 19 2704 278 8412 0.03 0.0187 0.0010 0.00282 0.00014 0.00112 0.00009 18.8 1.0 18.1 0.9 22.7 1.7 20 2832 305 8379 0.04 0.0188 0.0010 0.00294 0.00015 0.00100 0.00009 18.9 1.0 18.9 1.0 20.2 1.8 21 3020 198 8483 0.02 0.0181 0.0010 0.00292 0.00015 0.00300 0.00024 18.2 1.0 18.8 1.0 60.6 4.8 22 2760 400 8290 0.05 0.0178 0.0009 0.00279 0.00014 0.00126 0.00011 18.0 0.9 18.0 0.9 25.4 2.3 23 2677 509 8186 0.06 0.0182 0.0011 0.00265 0.00014 0.00158 0.00022 18.3 1.1 17.1 0.9 32.0 4.4 3.2 主量和微量元素
样品全岩主量和微量元素数据见表 2。错那淡色花岗岩具有较高的SiO2（74.46%~75.57%）、Al2O3 （14.07%~14.64%）和K2O（4.19%~4.85%），以及较低的CaO（0.51%~0.67%）、MgO（0.07%~0.10%）和TiO2（0.02% ~0.04%）含量。K2O/Na2O值较高，介于1.09~1.31之间。在SiO2-K2O图解中，样品点位于高钾钙碱性系列（图 5-a）。A/CNK值为1.15~1.25，均大于1.1，在A/CNK-A/NK图解中，样品点全部落入过铝质区域（图 5-b）中。总体来看，错那淡色花岗岩属于高钾钙碱性过铝质S型花岗岩。
表 2 藏南错那淡色花岗岩全岩主量、微量和稀土元素测试数据Table 2. Whole-rock major, trace and rare earth element analyses of the Cuona leucogranites in southern Tibet样品号 CN-01 CN-02 CN-03 CN-04 CN-05 SiO2 74.59 74.58 74.46 74.87 75.57 Al2O3 14.48 14.64 14.60 14.31 14.07 Fe2O3 0.70 0.83 0.90 0.77 0.69 MgO 0.08 0.09 0.10 0.07 0.08 CaO 0.51 0.53 0.59 0.64 0.67 Na2O 3.69 3.64 3.51 3.90 3.84 K2O 4.85 4.68 4.41 4.48 4.19 TiO2 0.02 0.03 0.04 0.03 0.03 P2O5 0.20 0.19 0.18 0.18 0.15 MnOz 0.04 0.04 0.06 0.05 0.04 烧失量 0.8 0.7 1.1 0.7 0.7 合计 99.16 99.25 98.85 99.30 99.33 A/NK 1.28 1.32 1.38 1.27 1.30 A/CNK 1.18 1.22 1.25 1.15 1.16 Na2O+K2O 8.54 8.32 7.92 8.38 8.03 K2O/Na2O 1.31 1.29 1.26 1.15 1.09 CaO/Na2O 0.14 0.15 0.17 0.16 0.17 Be 7 12 10 16 14 Sc 3 3 4 3 3 Co 0.4 0.2 0.4 0.5 1.1 Ni 0.1 0.1 0.3 0.2 0.2 Zn 5 9 13 11 9 Ga 16.9 18.5 20.8 18.9 19.5 Rb 377.6 413.5 381.7 353.8 353.0 Sr 12.8 23.3 18.2 15.3 15.2 Y 9.7 12.8 20.1 13.3 11.2 Zr 32.6 39.9 52.9 41.0 25.2 Nb 14.2 17.7 20.5 18.3 18.0 Cs 25.3 30.9 31.2 25.7 23.7 Ba 17 47 37 24 23 La 3.4 5.5 6.8 5.2 4.0 Ce 7.2 11.2 14.9 10.0 8.9 Pr 0.86 1.37 1.79 1.17 0.99 Nd 3.6 3.5 6.4 3.3 3.3 Sm 1.03 1.43 2.15 1.55 1.28 Eu 0.09 0.14 0.14 0.13 0.10 Gd 1.29 1.83 2.57 1.49 1.42 Tb 0.27 0.39 0.56 0.35 0.32 Dy 1.39 2.07 3.24 1.95 1.93 Ho 0.26 0.36 0.63 0.41 0.36 Er 0.78 1.11 1.78 1.09 0.92 Tm 0.12 0.19 0.30 0.20 0.15 Yb 0.86 1.28 1.94 1.18 1.11 Lu 0.13 0.19 0.28 0.19 0.16 Hf 1.0 1.3 2.3 1.7 1.3 Ta 4.0 4.4 5.4 6.1 4.8 Pb 6.3 6.7 7.2 6.3 6.7 Th 2.4 3.4 5.0 2.9 2.4 U 5.4 5.7 8.1 4.1 4.2 Sn 14 17 20 14 12 Rb/Sr 29.50 17.75 20.97 23.12 23.22 LREE 16.18 23.14 32.18 21.35 18.57 HREE 5.10 7.42 11.30 6.86 6.37 TREE 21.28 30.56 43.48 28.21 24.94 (La/Sm)w 2.13 2.49 2.04 2.17 2.02 (Gd/Yb)N 1.24 1.18 1.10 1.04 1.06 (La/Yb)N 2.84 3.08 2.52 3.16 2.59 δEu 0.24 0.26 0.18 0.26 0.23 注：主量元素含量单位为%，微量和稀土元素含量单位为10-6；A/NK=摩尔Al2O3/ （Na2O+K2O），A/ CNK=摩尔Al2O3/ （CaO+Na2O+K2O）；δEu=2EuN/ （SmN+GdN），其中N为球粒陨石标准化值[41] ![]()
图 5 藏南错那淡色花岗岩的SiO2-K2O (a)、A/CNK-A/NK分类图解(b)及原始地幔标准化蛛网图(c)和球粒陨石标准化稀土元素配分模式图(d) (原始地幔和球粒陨石数值据参考文献[41])Figure 5. TAS diagram (a), A/CNK-A/NK diagram (b), primitive mantle-normalized trace element spider diagram (c) and chondrite-normalized REE patterns (d) of Cuona leucogranites, southern Tibet在原始地幔标准化微量元素蛛网图（图 5-c）中，显示Rb、Th、U和Pb的正异常，以及Ba、Nb、Sr和Zr的负异常。稀土元素总量较低，为21.28×10-6~ 43.48×10-6，轻、重稀土元素分馏较弱，（La/Yb）N为2.59~3.16。在球粒陨石标准化的稀土元素配分模式图（图 5-d）中，具有强烈的负Eu异常，δEu值为0.18~0.26。
3.3 Rb-Sr、Sm-Nd同位素
样品全岩Rb-Sr、Sm-Nd同位素测试结果列于表 3。经过t=18Ma校正后，4个样品的Isr值较高，为0.78982~0.79276，εNd(t)值较低，为-19.5~-18.2，SrNd同位素的初始比值变化不大。Nd同位素二阶段亏损地幔模式年龄为2469~2372Ma，表明淡色花岗岩可能源自古老地壳的重熔。
表 3 藏南错那淡色花岗岩全岩Rb-Sr和Sm-Nd同位素测试数据Table 3. Whole-rock Rb-Sr and Sm-Nd isotopic compositions of the Cuona leucogranites in southern Tibet样品号 t/Ma Rb/10-6 Sr/10-6 87Rb/86Sr 87Sr/86Sr 2σ Isr Sm/10-6 Nd/10-6 147Sm/144Nd 143Nd/144Nd 2σm TDM2/Ma εNd(t) CN-01 18.0 377.6 12.8 83.3240 0.81406 0.00020 0.79276 1.03 3.60 0.181524 0.511636 0.000004 2372 -19.5 CN-02 18.0 413.5 23.3 50.1266 0.81431 0.00006 0.80150 1.43 3.50 0.259220 0.511652 0.000004 2469 -19.4 CN-04 18.0 353.8 15.3 65.3153 0.80910 0.00001 0.79241 1.55 3.30 0.298001 0.511654 0.000004 2452 -19.4 CN-05 18.0 353.0 15.2 65.5963 0.80659 0.00001 0.78982 1.28 3.30 0.246091 0.511709 0.000006 2390 -18.2 注：87Rb/86Sr和147Sm/144Nd通过ICP-MS测试的微量元素Rb、Sr、Sm和Nd计算所得，计算公式为87Rb/86Sr=Rb/Sr×2.981，147Sm/144Nd=Sm/Nd×[0.531497+0.142521×（143Nd/144Nd）s]。ISr=（87Sr/86Sr）s+87Rb/86Sr（eλt–1），143Nd/144Nd（t）=（143Nd/144Nd)s+147Sm/144Nd （eλt-1）；εNd（t）=[ （143Nd/144Nd）s/ （143Nd/144Nd)CHUR-1]×104。（143Nd/144Nd）CHUR=0.512638，（147Sm/144Nd）CHUR=0.1967，（143Nd/144Nd)DM=0.51315，（147Sm/144Nd)DM=0.2137；λRb=1.42×10-12/年[42]，λSm=6.54×10-12/年[43]；二阶段模式年龄TDM2的计算见参考文献[41] 4. 讨论
4.1 STDS淡色花岗岩的时代
喜马拉雅淡色花岗岩的形成时代一直是地学界研究的重点，准确获得各岩体的年龄，不仅有助于限定相关构造的活动时代，而且有利于重建造山过程中岩浆作用的时空格局，进而反演造山带的岩浆演化历史。近期，吴福元等[1]根据目前已发表的年龄数据，重新勾画了喜马拉雅淡色花岗岩的年代格架，大致划分为始喜马拉雅阶段（Eo-Himalayan：44~26Ma）、新喜马拉雅阶段（Neo-Himalayan：26~ 13Ma）和后喜马拉雅阶段（Post-Himalayan：13~ 7Ma）。以上3个时期的花岗岩最显著的区别在于受造山带内不同构造的控制：始喜马拉雅阶段的花岗岩多分布于THS中部的NHGD核部，新喜马拉雅阶段的淡色花岗岩主要沿STDS分布，而后喜马拉雅阶段的花岗岩受控于南北向裂谷（NSTR）。错那淡色花岗岩锆石U-Pb年龄为17.7±0.3Ma （MSWD=2.0），且岩体沿STDS分布，应属于上述划分方案中的新喜马拉雅阶段淡色花岗岩。在造山带的中西部，相同峰期时代的STDS淡色花岗岩也广泛分布（图 1）。①在西部，玛纳斯鲁淡色花岗岩年龄为19.3~18.1Ma[30, 33, 45]。②在中部，吉隆STDS淡色花岗岩脉年龄为18.6~16.6Ma[14]，聂拉木淡色花岗岩年龄为17.1~16.8Ma[9]，希夏邦马淡色花岗岩年龄为17.3Ma[46]，达扎卡淡色花岗岩年龄为16.7Ma[4]，绒布淡色花岗岩年龄为16.8~16.2Ma[18, 47]，珠峰淡色花岗岩年龄为16.9~16.4Ma[2, 30]，定结淡色花岗岩年龄为16.9~15.1Ma[7, 48-49]。③在东部，东不丹地区淡色花岗岩年龄为17.6~16.8Ma[50]。这些淡色花岗岩均沿区域性伸展构造STDS分布，且具有近乎一致的峰期结晶年龄，暗示其形成可能与STDS的活动密切相关。
4.2 错那淡色花岗岩的源区和源岩
错那淡色花岗岩具有高SiO2和Na2O含量，铝过饱和，K2O＞Na2O，为高钾钙碱性过铝质花岗岩。全岩的Isr值较高（0.78982~0.79276），εNd(t)值较低（-19.5~-18.2），Rb/Sr值（17.75~29.50）远高于大陆地壳的平均值（0.32），与典型壳源成因花岗岩特征吻合。而淡色花岗岩体内部及附近均未发现幔源包体和相关的岩石，因此，笔者认为错那淡色花岗岩应为壳源成因。
目前，关于STDS淡色花岗岩的源岩问题仍存在分歧。杨晓松等[51]的高温高压岩石学模拟研究发现，在温压分别为770~980℃和1.0~1.4GPa，且无外来流体的条件下，GHC中黑云斜长片麻岩发生部分熔融，形成的熔体与STDS淡色花岗岩成分类似，且熔融后的残留物质与GHC下地壳麻粒岩成分相当，据此认为，淡色花岗岩源岩为GHC黑云斜长片麻岩[27]。Knesel等[52]通过实验岩石学研究发现，GHC变泥质岩早期发生白云母脱水熔融，形成的花岗质熔体成分与前人获得的STDS淡色花岗岩的地球化学特征一致，因此，认为源岩应该为GHC变质沉积岩。考虑到淡色花岗岩的壳源成因特性，将错那淡色花岗岩样品的Sr、Nd同位素成分与潜在的可能源岩，包括GHC变泥质岩、LHS变质沉积岩和北喜马拉雅片麻岩[53-56]，统一以t=18Ma重新进行计算并综合对比（图 6）。结果显示，4件样品均位于GHC变泥质岩区域。此外，错那淡色花岗岩的二阶段亏损地幔模式年龄TDM2值（2469~2372Ma）与GHC变泥质岩的相同，且与GHC变泥质岩中碎屑锆石获得的约2500Ma的峰值年龄重叠[57]，暗示其可能来自GHC变泥质岩的部分熔融。综上所述，错那淡色花岗岩的源岩为GHC变泥质岩。
![]()
图 6 藏南错那淡色花岗岩的Isr-εNd(t)图Figure 6. Isr-εNd(t diagram of the Cuona leucogranites, southern Tibet4.3 错那淡色花岗岩的形成机制
STDS淡色花岗岩是GHC变泥质岩部分熔融的产物，但部分熔融的机制仍有争议，争议点包括熔融的矿物是白云母还是黑云母？熔融过程是否有外来流体参与？
Knesel等[52]对GHC变泥质岩进行部分熔融实验研究发现，单一源区的变泥质岩在熔融的不同阶段具有不同的熔融机制，并产生化学成分迥异的熔体。早期在无水条件下，白云母发生脱水熔融，形成低Sr、高Isr的花岗质熔体。而后期伴随早期脱水熔融产生的含水流体的加入，主要发生水致部分熔融，形成高Sr、低Isr的熔体，分布于不同地区的花岗岩的Sr和Isr数据也证实了上述实验结论。本文搜集已发表的喜马拉雅花岗岩的Sr和Isr数据，统一以t=18Ma重新计算并归纳于图 7-a中。从该图可以看出，花岗岩主要集中分布于a1和a2两个截然不同的区域，其中a1区域内花岗岩表现出低Sr（10×10-6~ 80×10-6）、高Isr（0.75~0.80）的特征，与部分熔融实验中无水条件下白云母脱水熔融形成的熔体化学特征吻合。而a2区域内花岗岩显示高Sr（50×10-6~ 200×10-6）、低Isr（0.72~0.75）的特征，与水致部分熔融形成的熔体化学特征相符（图 7-a）。错那淡色花岗岩具有低Sr（12.8×10-6~23.3×10-6）和高Isr （0.78982~0.79276）的特征，在图 7-a中，与其他地区的STDS淡色花岗岩一样位于a1区域中，暗示其为无水条件下白云母脱水熔融的产物。
![]()
图 7 藏南错那淡色花岗岩的Sr-Isr (a)和Ba-Rb/Sr(b)图Mus(VP)—饱和蒸汽相的白云母熔融蚀变；Mus(VA)—缺乏蒸汽相的白云母熔融蚀变；Bi(VA)—缺乏蒸汽相的黑云母熔融蚀变Figure 7. Sr-Isr(a) and Ba-Rb/Sr(b) diagrams of the Cuona leucogranites, southern Tibet实验岩石学[52, 60]和理论计算结果[61-62]表明，随温压条件和含水量的变化，变泥质岩发生从白云母到黑云母的递进部分熔融，形成具有复杂RbSr组分的花岗质熔体[24, 52, 61]。在变泥质岩中，Rb主要赋存于云母类矿物中，而Sr主要赋存于长石矿物中，在有水参与的地壳熔融过程中，白云母和斜长石都发生部分熔融，且长石熔融的量远大于云母，导致熔体的Rb/Sr值偏低。在无水条件下，白云母大量熔融，形成的熔体Rb/Sr值较高[63]。因此，可以通过Rb/Sr值及其与Ba浓度的关系反推部分熔融机制。错那淡色花岗岩具有较高的Rb/Sr值（17.75~29.50），在Ba-Rb/Sr图解（图 7-b）中，Rb/ Sr值随着Ba浓度的增加而减小，二者呈负相关，与缺乏蒸汽相的白云母脱水熔融的趋势一致[63]，也与上述Sr-Isr系统分析结果一致。同时，由于无水条件下参与熔融的矿物主要为白云母，斜长石未进入熔体而残留在源区，导致错那淡色花岗岩具有较强的负Eu异常。综上所述，错那淡色花岗岩是GHC变泥质岩在无水条件下白云母脱水熔融的产物。
4.4 构造动力学意义
喜马拉雅造山带内发育多种类型的花岗岩，其形成过程对应不同的造山阶段，反映不同的构造背景和动力学过程。如晚始新世—渐新世（45~28Ma）的花岗岩主要出露于NHGD，如然巴穹窿（44Ma）[64]、麻布迦穹窿（35Ma）[65]、萨迦穹窿（28.1Ma）[66]、佩枯错穹窿（28.2Ma）[67]等。这些花岗岩的形成年龄与喜马拉雅变质岩的进变质时代相同[68-69]，是早期地壳逆冲增厚条件下部分熔融的产物[37]，熔体的产生使地壳弱化，并直接影响与STDS相关的中下地壳的伸展拆离，是触发STDS启动的关键因素。另一类为晚渐新世—早中新世（26~13Ma）[1]的淡色花岗岩，以STDS淡色花岗岩为代表，其地球化学成分与始新世—渐新世花岗岩不同，以高Rb/Sr值和高Isr为特征，结晶时代与喜马拉雅峰期变质作用时代相同，其形成可能与地壳伸展背景下GHC的快速折返/隆升导致的减压熔融有关[25]。最后一类为晚中新世（13~7Ma）花岗岩，这类花岗岩的矿物和化学成分与STDS淡色花岗岩类似，主要的区别在于前者主要沿NSTR分布，其形成与造山后期的东西向伸展有关。
在GHC变泥质岩部分熔融的过程中，温度和压力是主控因素。错那淡色花岗岩的锆石饱和温度为657~715℃，这一温度低于白云母脱水熔融的温度（720~770℃）。在缺少深部热源证据和剪切摩擦生热，地壳中放射性同位素衰变生热无法产生大规模岩浆作用的情况下[34]，减压成为部分熔融的必要条件，而STDS的启动不仅可以起到减压的作用，还可为岩浆上侵提供空间和通道，可能是触发熔融的重要因素。已有研究表明，STDS强烈活动的时间集中在28~13Ma[14]，与STDS淡色花岗岩的时代基本一致，且淡色花岗岩主要沿STDS分布，佐证了以上结论。错那淡色花岗岩的年龄（17.7±0.3Ma）与STDS活动时代吻合，其地球化学特征也与晚渐新世—早中新世淡色花岗岩相同。从构造-岩浆作用角度综合考虑，笔者认为，错那淡色花岗岩可能形成于STDS启动之后的伸展减薄背景下，反映了构造减压导致的中下地壳变泥质岩中含水矿物（白云母）脱水熔融的动力学过程。
5. 结论
（1）藏南错那淡色花岗岩的结晶年龄为17.7±0.3Ma，岩石属于高钾钙碱性过铝质花岗岩，富集Rb、Th和U，亏损Ba、Nb、Sr和Zr，具有强烈的负Eu异常；Rb/Sr值较高，Isr和εNd(t)值分别为0.78982~0.79276和-19.5~-18.2。
（2）错那淡色花岗岩为壳源成因，源岩是GHC变泥质岩，是无水条件下白云母脱水熔融的产物。熔体的形成与STDS的启动有关，反映了地壳伸展减薄背景下，构造减压导致中下地壳变泥质岩中含水矿物脱水熔融并沿STDS上升侵位的动力学过程。
-
![]()
图 2 松辽盆地储盖层组合示意图
Figure 2. The sketch map of reservoir-cap assemblages in Songliao Basin
![]()
图 3 松辽盆地CO2地质储存适宜性区划
Figure 3. Suitability evaluation zoning map of CO2 geological sequestration in Songliao Basin
表 1 松辽盆地CO2地质储存适宜性评价指标体系
Table 1 Suitability assessment index system of CO2 geological sequestration in Songliao Basin
准则层 指标层 含义 储存安全性B1 地壳动峰值加速度（g）C1 地震时地面运动的加速度 地震安全性C2 地震的活跃程度 盖层埋深（m）C3 盖层距地表面的埋深 盖层岩性C4 盖层的岩石类型 盖层累计厚度（m）C5 盖层中泥岩的总厚度 断裂密度（km·km-2）C6 二级构造单元内所有断裂总长度与二级构造单元面积的比值，原则上二级构造单元越完整，越利于CO2地质储存 水文地质条件C7 地下水的水动力条件 储存规模B2 储层厚度（m）C8 储层中砂岩的总厚度 储层岩性C9 储层的岩石类型 储层压力系数C10 地层压力与静水压力的比值 地热流值（MW·m-2）C11 单位时间通过地球表面单位面积的热流值 地温梯度（℃·100m-1）C12 地层温度随着深度增加的增长率 有效孔隙度（%）C13 砂岩的有效孔隙度 有效储存量（104t）C14 CO2在埋存过程中受到储层的非均质性、CO2浮力等因素的影响，理论储存量是无法达到的，有效储存量可以利用有效储存系数换算得到 社会环境风险B3 人口密度（人·km-2）C15 单位面积土地上居住的人口数 经济适宜性B4 碳源规模（104t·a-1）C16 碳气源的规模 
下载: 导出CSV
表 2 CO2地质储存适宜性评价指标分级
Table 2 Suitability index classification of CO2 geological sequestration
准则层 指标层 好 一般 差 储存安全性B1 地壳动峰值加速度（g）C1 < 0.05 ≤0.15 ＞0.2 地震安全性C2 安全（8） 中等（5） 较危险（2） 盖层埋深（m）C3 1100~1900 1900~2900
700~1100＞2900
< 700盖层岩性C4 蒸发岩类（8） 泥质岩类（5） 页岩和致密灰岩（2） 盖层累计厚度（m）C5 ＞300 150~300 < 150 断裂密度（km· km-2）C6 < 0.3 0.3~0.4 ＞0.4 水文地质条件C7 流动规模长、深度深、更新慢（8） 流动规模中等、深度中等、更新一般（5） 流动规模短、深度浅、更新快（2） 储存规模B2 储层厚度（m）C8 ＞80 30~80 < 30 储层岩性C9 碎屑岩（8） 碎屑岩与碳酸盐混合（5） 碳酸盐（2） 储层压力系数C10 < 0.9 0.9~1.1 ＞1.1 地热流值（MW·m-2）C11 < 50 50~70 ＞70 地温梯度（℃· 100m-1）C12 < 2 2~3 ＞3 有效孔隙度（%）C13 ＞15 10~15 < 10 有效储存量（104t）C14 ＞9000 3000~9000 < 3000 社会环境风险B3 人口密度（人·km-2）C15 < 25 26~50 ＞51 经济适宜性B4 碳源规模（104t· a-1）C16 ＞10000 5000~10000 < 5000
下载: 导出CSV
表 3 松辽盆地CO2地质储存适宜性评价结果
Table 3 Suitability assessment results of geological sequestration of CO2 in Songliao Basin
一级构造 二级构造 λi1 (好） λi2(—般） λi3（差） 适宜性 北部倾没区 三兴背斜带 0.538 0.694 0.676 一般 克山依龙背斜带 0.559 0.656 0.653 一般 乾元背斜带 0.603 0.645 0.641 一般 乌裕尔凹陷 0.556 0.661 0.658 一般 西部斜坡区 富裕构造带 0.775 0.843 0.833 一般 泰康隆起带 0.598 0.606 0.674 差 西部斜坡带 0.570 0.649 0.661 差 东北隆起区 海伦隆起带 0.599 0.649 0.641 一般 绥棱背斜带 0.578 0.650 0.641 一般 绥化凹陷 0.569 0.661 0.649 一般 呼兰隆起带 0.589 0.638 0.669 差 东南隆起区 长春岭背斜带 0.506 0.658 0.684 差 登娄库背斜带 0.515 0.657 0.658 差 宾县王府凹陷 0.615 0.719 0.576 一般 钓鱼台凸起 0.509 0.735 0.659 一般 九台阶地 0.497 0.718 0.682 一般 榆树德惠凹陷 0.535 0.710 0.625 一般 青山口隆起 0.562 0.717 0.600 一般 中央凹陷区 黑鱼沟凹陷 0.715 0.636 0.542 好 明水阶地 0.729 0.611 0.531 好 龙虎泡大安阶地 0.696 0.596 0.582 好 齐家古龙凹陷 0.713 0.616 0.568 好 大庆长垣 0.707 0.597 0.568 好 三肇凹陷 0.723 0.592 0.557 好 朝阳沟阶地 0.638 0.659 0.575 一般 扶余隆起 0.689 0.601 0.584 好 双坨子阶地 0.641 0.653 0.587 一般 长岭凹陷 0.714 0.603 0.558 好
下载: 导出CSV
-
张菊. 华北南部盆地二氧化碳地质储存潜力与适宜性评价[D]. 石家庄经济学院硕士学位论文, 2013. http://cdmd.cnki.com.cn/Article/CDMD-10077-1014016120.htm Bachu S. Screening and ranking of sedimentary basins for sequestra-tion of CO2, in geological media in response to climate change[J]. Environmental Geology, 2003, 44(3):277-289. doi: 10.1007/s00254-003-0762-9
李甫成, 张杨, 张晓娟, 等.深部咸水层CO2地质储存适宜性评价方法研究[J].冰川冻土, 2014, 36(3):649-660. 沈平平, 廖新维.二氧化碳地质埋存与提高石油采收率技术[M].北京:石油工业出版社, 2009. 张森琦.中国二氧化碳地质储存地质基础及场地地质评价[M].北京:地质出版社, 2011. 巫润建, 李国敏, 黎明, 等.松辽盆地咸含水层埋存CO2储存容量初步估算[J].工程地质学报, 2009, 17(1):100-104. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=gcdz200901016&dbname=CJFD&dbcode=CJFQ 杨国强, 苏小四, 杜尚海, 等.松辽盆地CO2地质储存适宜性评价[J].地球学报, 2011, 32(5):570-580. Su X, Xu W, Du S. Basin-scale CO2, storage capacity assessment of deep saline aquifers in the Songliao Basin, northeast China[J]. Greenhouse Gases Science and Technology, 2013, 3(4):266-280. doi: 10.1002/ghg.1354
刁玉杰, 张森琦, 郭建强, 等.深部咸水层CO2地质储存地质安全性评价方法研究[J].中国地质, 2011, 38(3):786-792. 罗毅, 周创立, 刘向杰.多层次灰色关联分析法在火电机组运行评价中的应用[J].中国电机工程学报, 2012, 32(17):97-103. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=zgdc201217016&dbname=CJFD&dbcode=CJFQ 侯启军.松辽盆地陆相石油地质学[M].北京:石油工业出版社, 2009. 李国玉.中国含油气盆地图集[M].北京:石油工业出版社, 2002. 张洪涛, 文冬光, 李义连, 等.中国CO2地质埋存条件分析及有关建议[J].地质通报, 2005, 24(12):1107-1110. doi: 10.3969/j.issn.1671-2552.2005.12.004 范基姣, 贾小丰, 张森琦, 等. CO2地质储存潜力与适宜性评价方法及初步评价[J].水文地质工程地质, 2011, 38(6):108-112. 邓聚龙.灰色系统基本方法[M].武汉:华中工学院出版社, 1987. 施萍.利用灰色关联法评价苏州河水质状况[J].环境科技, 2008, 21(a01):60-63. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=jshj2008s1022&dbname=CJFD&dbcode=CJFQ 王平, 王云峰.综合权重的灰色关联分析法在河流水质评价中的应用[J].水资源保护, 2013, 29(5):52-54. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=szyb201305011&dbname=CJFD&dbcode=CJFQ Zou Z H, Yi Y, Sun J N. Entropy method for determination of weight of evaluating indicators in fuzzy synthetic evaluation for water quality assessment[J]. Journal of Environmental Sciences, 2006, 18(5):1020-1023. doi: 10.1016/S1001-0742(06)60032-6
马雪倩. 黑龙江省水资源可持续利用评价研究[D]. 东北农业大学硕士学位论文. 2013. 国家地震局.1:4000000中国地震烈度区划图[M].北京:地震出版社, 1990. 楼章华, 程军蕊, 金爱民.沉积盆地地下水动力场特征研究——以松辽盆地为例[J].沉积学报, 2006, 24(2):193-201. 楼章华, 金爱民, 朱蓉, 等.论松辽盆地地下水动力场的形成与演化[J].地质学报, 2001, 75(1):111-120. http://d.wanfangdata.com.cn/Periodical/dizhixb200101013 杨绪充.论含油气盆地的地下水动力环境[J].石油学报, 1989, 10(4):30-37. 吉林油田石油地质志编写组.中国石油地质志卷二——大庆、吉林油田(上册), 吉林油田(下册)[M].北京:石油工业出版社, 1993. 赵志魁, 张金亮, 赵占银, 等.松辽盆地南部坳陷湖盆沉积相和储层研究[M].北京:石油工业出版社, 2009. 侯启军.深盆油藏:松辽盆地扶杨油层油藏形成与分布[M].北京:石油工业出版社, 2010. 吴乾蕃.松辽盆地地热场[J].地震研究, 1991, (1):31-40. 吴乾蕃, 谢毅真.松辽盆地大地热流[J].地震地质, 1985, 7(2):61-66. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzdz198502009&dbname=CJFD&dbcode=CJFQ 迟元林, 云金表, 蒙启安, 等.松辽盆地深部结构及成盆动力学与油气聚集[M].北京:石油工业出版社, 2002. Bachu S. Comparison between methodologies recommended for estimation of CO2 storage capacity in geological media by the CSLF task force on CO2 storage capacity estimation and the USDOE capacity and fairways subgroup of the regional carbon sequestration partnerships program[J]. Carbon Sequestration Leadership Forum, 2008:1-21.
涂华, 刘翠杰.标准煤二氧化碳排放的计算[J].煤质技术, 2014, (2):57-60. 李小春, 魏宁, 刘延锋, 等. 中国二氧化碳地质封存潜力评估和早期实施机遇. 中国科学院武汉岩土力学研究所, 2008. 中国地质调查局水文地质环境地质调查中心.重点地区二氧化碳地质储存适宜性与选址调查评价技术要求（征求意见稿）. 中国地质调查局水文地质环境地质调查中心, 2013.
计量
- 文章访问数: 2650
- HTML全文浏览量: 398
- PDF下载量: 2068
