盆地地层热历史主要受埋藏和热体制控制，岩浆活动作为非稳态的热源，其在地球内部积聚、上涌以及在近地表喷发的过程也是能量释放的过程，会对盆地内有机质的演化产生影响. 研究表明塔里木盆地古地温梯度自震旦纪以来整体呈下降趋势，其中志留纪至晚古生代从拉张盆地时期过渡到稳定地壳热演化时期，地温梯度逐渐降低；中生代地温梯度较稳定；新生代则叠加了盆地边缘强烈挠曲热演化过程(Chang et al.，2014; Li et al.，2010; Qiu et al.，2012; 李成等，2000; 刘绍文等，2006; 邱楠生等，2004; 左银辉等，2015). 二叠纪岩浆岩喷发于早、中二叠世，约300—275 Ma，其中大规模玄武岩发育在290—288 Ma(Li et al.，2012b; Polyakov et al.，2008; Wei et al.，2014; Xu et al.，2014; Yu et al.，2011; Zhang et al.，2010，2013; 陈汉林等，1997，2009; 杨树锋等，2014)，然而目前针对该期岩浆活动对古地温的影响范围和时间存在明显争议. 由于有机质成熟演化的不可逆性，其与非稳态高温热事件的可能关系值得关注. 为此，本文通过分析塔里木盆地R_o值随深度的演化规律，恢复典型单井热历史，探讨了石炭-二叠纪异常热演化及其成因，这将为进一步认识盆地深部构造-岩浆活动及其石油地质意义提供重要基础.

塔里木盆地是中国最大的沉积盆地，位于天山和昆仑山褶皱带之间，是在古老基底之上震旦纪-古生代克拉通叠覆中新生代前陆盆地而成的大型叠合盆地(贾承造，2004; 贾承造和魏国齐，2002)，经历了震旦纪-中泥盆世、晚泥盆世-三叠纪和侏罗纪-第四纪3个伸展-聚敛演化旋回(何碧竹等，2011; 何登发等，2005)，构造环境发生了多次重大变革. 复杂的构造运动导致了多期岩浆活动，盆地内主要发育四期火山岩，分别在震旦-寒武纪、晚奥陶世-志留纪、二叠纪、白垩纪，其中晚古生代火山活动强度很大(陈汉林等，1997; 温声明等，2005; 徐汉林等，2006; 张水昌等，2004). 晚古生代是塔里木盆地演化过程的一个重要阶段，此时塔里木克拉通处于强烈的拉张状态，盆地内部发育了大量岩浆岩，以基性岩类为主，如玄武岩、辉绿岩、玄武安山岩和超镁质岩类等，在满加尔凹陷西段、塔北隆起西段、塔中隆起、巴楚隆起和塔西南坳陷中段均有分布，面积达2.5×10⁵ km²以上(图 1)(刘晓等，2011; 潘赘等，2013; 杨树锋等，2005，2014).

图 1

Fig. 1

图 1 塔里木盆地构造单元和研究井位分布图 二叠纪火山岩分布据杨树锋等(2005). Fig. 1 Sketch map of structural units in the Tarim basin and well locations The distribution of Permian igneous rock is drawn according to Yang et al.(2005).

石炭-二叠纪是塔里木盆地内部沉积演化的重要时期. 塔里木盆地经历了震旦纪-中泥盆世活跃的构造变迁之后，在晚泥盆世进入了相对稳定的演化阶段，发育了由西南向东北的海侵，并在石炭纪初急剧扩张(何登发等，2005). 晚泥盆世-石炭纪总体上以海平面相对上升的海侵为主，期间出现两次快速下降和两次缓慢小幅度下降(王毅，1998; 朱如凯等，2002)，海平面的频繁变化造就了碳酸盐岩与碎屑岩交互发育的特点，自下而上依次沉积了含砾砂岩段、下泥岩段、生屑灰岩段、上泥岩段、砂泥岩段、顶灰岩段. 前期相对稳定的构造环境使得盆地内含砾砂岩段-砂泥岩段保存较完整，直至顶灰岩段在盆地内部及周缘地区出现不同程度的缺失，且向北向西愈加严重(朱如凯等，2002)，上石炭统剥蚀发生于石炭纪末小海子组沉积后-二叠纪初碎屑岩段沉积前，剥蚀量可达到800 m以上(Li et al.，2014). 本文选取了岩浆岩分布区中8口井进行研究(图 1)，自北向南的地层连井剖面显示石炭系顶剥蚀程度具有北强南弱的趋势，残余顶灰岩段向南增厚(图 2). 二叠纪是塔里木板块构造-沉积演化一个重要的转折期，早二叠世结束了海相沉积，进入了陆相沉积阶段，发育了以陆相碎屑岩为主的沉积以及火山岩和火山碎屑岩.

图 2

Fig. 2

图 2 石炭系-下二叠统连井剖面 Fig. 2 Stratigraphic profile of Carboniferous and Lower Permian

镜质体反射率是目前应用最广泛的古温标之一，是本文研究塔里木盆地晚古生代热事件的依据. R_o值受控于有机质受热温度和受热时间，以温度影响为主(Buchardt and Lewan，1990; Burnham and Sweeney，1989; 邱楠生等，2004)，稳定沉降的状态下R_o与深度(温度)在对数坐标中表现出线性关系(Baker and Pawlewicz，1986； Baker and Goldstein，1990; Sweeney and Burnham，1990)；有机质热演化具有不可逆性，能够指示如岩浆活动异常高温热事件，非稳态高热源侵入浅部地层引起R_o增大，若后期正常沉积下地层温度不超过非稳态热效应对有机质的影响，R_o演化不连续. 另外，R_o值演化可揭示地层的沉降-抬升信息，在相同的热状态下，地层抬升剥蚀是导致R_o在深度剖面上呈现不连续演化的主要原因.

综合分析塔里木盆地不同地区的大量实测R_o数据，多口钻井的R_o值在深度剖面上的分布具有规律性. 二叠系与前二叠系(SL1井为泥盆系，其余井为石炭系)中的R_o值不连续，不整合面附近多具零星的极高R_o异常值，记录了明显的构造-热事件，揭示出石炭纪末地层经历了较高的古地温(表 1、图 1和图 3).

表 1(Table 1)

表 1 塔里木盆地镜质体反射率 Table 1 Vitrinite reflectance in Tarim basin 井号 深度(m) 层位 Ro 岩性 SL1 4403 K 0.54 粉砂质泥岩 SL1 4451 T 0.56 泥岩 SL1 4452 T 0.59 泥岩 SL1 4515 T 0.60 粉砂质泥岩 SL1 4518 T 0.62 粉砂质泥岩 SL1 4550 T 0.63 粉砂质泥岩 SL1 4911 T 0.72 泥岩 SL1 4940 T 0.72 泥岩 SL1 5236 P 0.78 粉砂质泥岩 SL1 5470 P 2.64 泥质粉砂岩 SL1 6290 S 1.61 泥岩 SL1 6364 S 1.63 泥质粉砂岩 SL1 6395 S 1.64 泥质粉砂岩 Sh1 2995 T 0.38 泥岩 Sh1 3270 T 0.50 泥岩 Sh1 3574 T 0.51 泥岩 Sh1 3742 P 0.60 泥岩 Sh1 3889 P 0.65 粉砂质泥岩 Sh1 4913 P 1.70 灰质泥岩 Sh1 4222 C 0.85 泥岩 Sh1 4222 C 0.89 泥岩 Sh1 4224 C 0.90 泥岩 Sh1 4836 D 1.32 粉砂质泥岩 H4 954 P 0.41 泥岩 H4 1058 P 0.48 泥岩 H4 1061 P 0.51 泥岩 H4 1070 P 0.52 泥岩 H4 1240 C 0.96 泥岩 H4 1243 C 0.96 泥岩 H4 1244 C 0.98 粉砂质泥岩 H4 1260 C 1.00 泥岩 H4 1264 C 0.97 泥质粉砂岩 H4 1285 C 0.76 泥岩 H4 1312 C 0.74 泥岩 H4 1340 C 0.75 泥岩 H4 1542 C 0.80 泥质粉砂岩 H4 1544 C 0.80 泥质粉砂岩 H4 1570 C 0.81 泥岩 H4 1653 C 0.80 灰质泥岩 H4 1815 D 0.83 粉砂质泥岩 H4 1818 D 0.83 粉砂质泥岩 H4 1820 D 0.84 粉砂质泥岩 HT1 1068 N 0.3 泥岩 HT1 1173 P 0.31 泥岩 HT1 1898 P 0.51 泥岩 HT1 1900 P 0.50 泥岩 HT1 2588 P 1.06 含砾泥岩 HT1 2762 C 0.72 泥岩 HT1 2806 C 0.74 灰质泥岩 HT1 2849 C 0.76 灰质泥岩 HT1 2890 C 0.76 泥岩 HT1 2920 C 0.77 泥岩 HT1 2958 C 0.79 泥岩 HT1 2985 C 0.80 粉砂质泥岩 BD2 1537 T 0.32 泥岩 BD2 1771 P 0.42 泥岩 BD2 2313 P 0.45 泥岩 BD2 2766 P 0.51 泥岩 BD2 2932 P 0.52 泥岩 BD2 2976 P 0.80 泥岩 BD2 3504 C 0.62 灰质泥岩 BD2 3520 C 0.59 灰质泥岩 BD2 3622 C 0.60 粉砂质泥岩 TZ1 2495 T 0.46 粉砂岩 TZ1 2497 T 0.46 泥岩 TZ1 2705 P 0.50 泥岩 TZ1 2711 P 0.51 泥岩 TZ1 2882 P 0.53 泥岩 TZ1 2915 P 0.58 泥岩 TZ1 3145 C 0.74 泥岩 TZ1 3463 C 0.81 泥岩 TZ1 3478 C 0.89 灰质泥岩 TZ1 3482 C 1.00 泥岩 TZ1 3579 C 0.95 泥岩 TZ1 3581 C 1.04 泥岩 TB2 2206 R 0.42 泥岩 TB2 2502 T 0.48 泥岩 TB2 3411 P 0.60 泥岩 TB2 3711 P 1.65 泥岩 TB2 4043 C 0.79 泥岩 TB2 4162 C 0.90 泥岩 TB2 4242 C 0.81 泥岩 TB2 4549 C 0.85 泥岩 TB2 4706 C 1.00 泥岩 TC1 2680 T 0.56 泥岩 TC1 2844 T 0.57 泥岩 TC1 3128 P 0.61 泥岩 TC1 3283 P 0.62 泥岩 TC1 3556.5 P 0.63 泥岩 TC1 3665 P 0.87 泥岩 TC1 4049 P 1.03 泥岩 TC1 4135 C 1.05 泥岩 TC1 4242 C 0.97 泥岩 TC1 4339 C 0.83 泥岩 TC1 4356 C 0.83 泥岩 TC1 4442 C 0.89 泥岩 TC1 4816 C 1.00 泥岩 TC1 5031 S 1.08 泥岩 注：Sh1井、BD2井和TC1井来源于中石化西北分公司，TZ1井数据来源于王铁冠等(2011)和肖贤明等(1995).

表 1 塔里木盆地镜质体反射率 Table 1 Vitrinite reflectance in Tarim basin

图 3

Fig. 3

图 3 塔里木盆地典型井镜质体反射率随深度的演化关系 Fig. 3 Ro versus depth of typical wells in Tarim basin

本文根据EasyR_o%模型(Burnham and Sweeney，1989; Sweeney and Burnham，1990)，利用R_o值进行热史反演. 为探讨二叠系与石炭系之间R_o值不连续与古地温和地层剥蚀量的关系，选择不同构造单元中4口井为基础进行详细模拟. 如巴楚隆起H4井二叠系与石炭系间R_o的差值若仅由地层抬升剥蚀造成的情况下，模拟结果显示石炭系顶的剥蚀量需超过3500 m(图 4a)，综合前人研究成果其剥蚀量为~1000 m(Li et al.，2014；Lin et al.，2012; 齐永安和刘国臣，1999; 张一伟等，2000; 李丽贤等，2011; 丛培泓等，2015). 若存在高温热事件的条件下，考虑岩浆上涌释放充沛热量烘烤地层，限定剥蚀量为1000 m，模拟结果显示石炭纪末地层温度显著升高，达到峰值，促进了有机质的热演化；而后地层温度逐渐下降，到中二叠世恢复到正常热状态(图 4b). Sh1井、HT1井和TB2井模拟亦揭示出相同的结果：石炭-二叠系间R_o值的不连续演化是高温热事件与抬升剥蚀共同作用的结果，石炭纪末地温出现短暂的高峰；二叠纪初地层温度快速降低，之后地层温度主要受埋藏深度的控制.

图 4

Fig. 4

图 4 H4井、Sh1井、TZ1井和TB2井埋藏史和热史模拟结果 (a，c，e，g)分别为H4井、Sh1井、TZ1井和TB2井仅有剥蚀的假设情况;(b，d，f，h) 分别为H4井、Sh1井、TZ1井和TB2井在剥蚀和热事件控制下的实际情况. Fig. 4 Burial and thermal history of wells H4，Sh1，TZ1 and TB2 (a，c，e，g)are simulated results only with erosion for Well H4，Sh1，TZ1 and TB2，respectively.(b，d，f，h)are simulated results controlled by erosion and magmatic activity for wells H4，Sh1，TZ1 and TB2，respectively.

塔里木盆地自震旦纪至今总体上逐渐变冷，前石炭纪稳定演化，邱楠生等(2004)研究成果显示出地温梯度逐渐降低的规律. 研究井热史模拟显示石炭纪末盆地地温梯度短暂地快速提升，至300 Ma左右达到高峰值，分布在4.8~5.6 ℃/100 m之间. 距今约285 Ma时高地温梯度快速回落至3 ℃/100 m左右，进入持续缓慢降低的阶段(图 4). 地温梯度高峰期在区域上呈现出西北高东南低的规律：在塔北地区(SL1井)最高，巴楚隆起北部(H4井)次之，塔中地区(TZ1井)最低(图 5). 由于塔里木盆地地质历史复杂，多期构造运动叠加使得盆地内地层遭受了不同程度的剥蚀，东部大部分地区二叠系缺失，南部二叠系则俯冲至昆仑山下，因此盆地东南部对石炭-二叠纪热事件的记录缺失，有待进一步分析确认.

图 5

Fig. 5

图 5 塔里木盆地石炭纪末地温梯度等值线图 剥蚀区据Li et al.(2014). Fig. 5 Thermal gradient contours at end of Carboniferous in Tarim basin The erosion zone is from Li et al.(2014).

镜质体是由高等植物演化而来，陆生植物大发展始于晚古生代. 泥盆纪晚期植物界完成脱离水体的变革，乔木型植物已占相当的优势，并可形成小型森林；石炭纪陆生植物进一步繁荣，地球上首次出现大规模的森林；晚二叠世开始，全球许多地区裸子植物兴起，占主导地位(杜远生和童金南，1998; 刘本培和全秋琦，1996; 肖传桃，2007). 因此晚志留世之后的沉积地层中包含镜质体，可以直接测量；而晚古生代之前沉积的地层中只含有沥青、镜状体等，这些古温标必须通过公式转换成等效镜质体反射率，准确性相对较差. 在文中重点比较的二叠系、石炭系和上泥盆统中的R_o值是直接测量的结果，采用这些测量值进行研究具有较高的可靠性.

盆地内岩浆活动对地温场的影响范围和程度一直有争议，本文中-下二叠统和上石炭统中零星R_o极高值点和石炭-二叠系R_o值错断的现象为这一问题提供了解释的切入点. 侵入体对围岩热影响的范围与岩体性质和规模有关. 酸性岩浆的温度约为700~900 ℃，中性岩浆的温度约为900~1000 ℃，基性岩浆的温度约为1000~1300 ℃，高温熔融物质在浅部地层运移过程中冷却结晶成岩释放的热量作用于围岩有机质热成熟演化，导致R_o值急剧上升，可达6%以上，远远超出有机质正常热演化所能达到的程度(Annen and Sparks，2002; Fjeldskaar et al.，2008; Galushkin，1997; 王民，2010). 岩浆沿构造薄弱带上涌时以不同形式侵入浅部地层或喷出地表，研究表明受火山岩体热影响的围岩厚度约为其厚度的0.3~4倍，通常火山侵入体规模较小，对多个地区火山侵入体周边围岩内有机质成熟程度统计结果显示距侵入体数米至百米范围内的围岩被严重烘烤，有机质演化程度明显增高(Fjeldskaar et al.，2008; Galushkin，1997; Stewart et al.，2005; 陈荣书等，1989; 高福红等，2009; 朱传庆等，2010a). 本文研究井中除TZ1井外均钻遇了岩浆岩，以玄武岩为主，安山岩、英安岩和辉绿岩也可见，厚度范围在63~429 m之间，岩浆岩附近出现了R_o高值“漂点”(如Sh1井4913 m测量点为1.70%)，离岩浆岩距离越近R_o值越高，但高温作用范围有限，距岩浆岩较远的有机质演化正常，R_o值不受影响(图 3). 塔里木盆地中邻近岩浆岩的R_o演化与前人研究规律一致：围岩温度随与岩体距离的增加急剧降低，导致围岩内R_o值随着快速下降，呈陡倾状连续演化(图 6)(Clayton and Bostick，1986; Galushkin，1997; Simoneit et al.，2004; Simonet et al.，1981). 岩浆活动难以大范围提高地层温度，尤其是在盆地尺度上，仅能作为局部热源. 塔里木盆地石炭-二叠系之间R_o值不连续，热史模拟显示塔里木盆地北部及西部地区石炭纪末地层温度升高，意味着在岩浆喷发前盆地大范围地区处于高地热背景下. 地球内部活动性大的岩浆易向低压区域聚集形成岩浆房，其分布范围远远超出岩浆侵入带的范围，该阶段岩浆熔融形成过程中积聚的热量通过热传导作用于上部地层，对区域热演化有一定的影响，且影响时间超过岩浆侵入的影响时间(范桃园等，1999; 朱传庆等，2010a). 塔里木盆地石炭纪末岩浆自地球深部携带热量上涌的过程导致了区域性地温升高，而二叠纪岩浆喷发活动对地温场的影响局限在岩体周围但效果显著.

图 6

Fig. 6

图 6 受岩浆岩影响的Ro值演化剖面 图中实线为Ro变化趋势. Fig. 6 Vitrinite reflectance profile influenced by igneous rocks Solid line shows changing trend of Ro.

塔里木早-中二叠世大火成岩省是在我国继峨眉山大火成岩省之后相对比较确认的又一个大火成岩省，随着研究的深入越来越多的证据指出塔里木大火成岩省与地幔柱关系密切. 朱传庆等(2010a，2010b)在恢复四川盆地热流演化史的基础上，证实古热流的时间演化规律、空间分布特征与峨眉山岩浆活动密切相关，为峨眉山地幔柱在上扬子地区的存在提供地热学依据.塔里木盆地构造-热演化亦具有相似的特征. 通常认为地幔柱起源于核幔边界，地幔柱柱头直径最大可达千米以上，具有异常高的温度(Campbell，2001)，由于上升地幔柱对岩石圈的动力和能量冲击，可造成大规模的地壳抬升和地层温度升高(Griffiths and Campbell，1990)，有别于一般的岩浆活动和地层抬升.塔里木盆地近半区域内异常的构造-热历史从地热学角度一定程度地指示出大火成岩省与地幔柱的关系. 塔里木盆地大火成岩省的岩浆喷出在300—275 Ma(Li et al.，2012b; Polyakov et al.，2008; Wei et al.，2014; Xu et al.，2014; Yu et al.，2011; Zhang et al.，2010，2013; 陈汉林等，1997，2009; 杨树锋等，2014)，大规模玄武质岩浆岩形成于290—288 Ma. 通过沉积地层学研究确定塔里木岩浆喷发前地壳快速大幅度抬升，强烈隆升发生于~300 Ma(Li et al.，2014)，造成石炭-二叠系之间的不整合. 盆地地温梯度值自石炭纪末期开始迅速升高，最高值出现在~300 Ma，与地层强烈隆升的时间一致，呼应了高温地幔柱柱头的上升过程(图 7b). 热史模拟显示石炭纪末期塔里木盆地地温受热效应的影响程度自北向南方向减弱(图 5)，该时期地层剥蚀程度亦表现出北强南弱的特点，热史演化和地层剥蚀在空间上相呼应. 随后地温梯度开始迅速降低(图 7c)，这是由于地幔柱的热量在与塔里木岩石圈地幔发生持续作用和岩浆喷发的短暂过程中大量损耗. 由于地幔柱与岩石圈地幔的相互作用，大火成岩省从早期的玄武岩到晚期的侵入岩体现出从富集地幔组分转变为亏损地幔组分的演化趋势，侵入岩均具有源于亏损地幔源区的特征(Li et al.，2012a; 励音骐，2013; 余星，2009). 前人通过沉积学、岩石学和地球化学方法得出，地幔柱中心柱头位于盆地北部地区(Tian et al.，2010; Zhang et al.，2013; 陈汉林等，2006)；值得注意的是天山二叠系玄武岩与塔里木玄武岩具有相似的地球化学特征，推断是同一地幔柱的产物，但天山玄武岩低钛高镁、高ε_Nd(t)值的特点表明其可能是来自高程度部分熔融的岩石圈地幔的产物，与地幔柱头更为邻近，塔里木盆地玄武岩高钛低镁、低ε_Nd(t)值的特点意味着其可能来自于柱头边缘低程度部分熔融的区域(Qin et al.，2011; 夏林圻等，2004). 塔里木盆地石炭纪末地温梯度自北向南降低，反映出异常热作用北强南弱的特点，这种空间上的分布规律应是对石炭纪末地幔柱柱头烘烤浅部地层的响应，意味着盆地西北部更靠近热源中心.

图 7

Fig. 7

图 7 ~315—285 Ma塔里木盆地古地温和地幔柱演化示意图 Fig. 7 Schematic illustration for ~315—285 Ma evolution of thermal history and mantle plume in Tarim basin

(1) 塔里木盆地西北部多口钻井显示二叠系与前二叠系内镜质体反射率(R_o)值随埋深出现不连续演化，R_o分析和热史模拟表明，盆地内自晚石炭世发育一期高温热事件，岩浆喷出前岩浆房的烘烤对盆地地温场产生了重要影响，导致区域地温升高，并在~300 Ma达到峰值，随后至早二叠世温度迅速下降；而二叠纪大规模岩浆喷发过程(290—288 Ma)的热效应仅能作用有限区域.

(2) 塔里木盆地石炭纪末-二叠纪初的构造-热事件，从地热学角度提供了塔里木地幔柱存在的证据，盆地地温梯度具有北高南低的特点，异常热事件对盆地的影响向南部地区减弱，指示塔里木盆地北部更靠近热源中心区.