上扬子西南缘下志留统细粒沉积岩有机质富集机理研究<br/>——以盐津牛寨剖面为例

上扬子西南缘下志留统细粒沉积岩有机质富集机理研究
——以盐津牛寨剖面为例

2021-06-17何佳伟谢渊刘建清何利

沉积学报 2021年3期

何佳伟,谢渊,刘建清,何利

1.成都理工大学沉积地质研究院，成都 610059

2.中国地质调查局成都中心，成都 610081

细粒沉积岩是指由粒径小于62μm细粒沉积物固结而成的沉积岩，这些细粒沉积物主要包括黏土类矿物、粉砂质、碳酸盐矿物和有机质等[1-2]。在我国细粒沉积岩主要分布于四川盆地、塔里木盆地和华北地区古生界的相关地层中[3]，随着非常规油气藏理论的不断完善，分布广、厚度大的细粒沉积岩进入了人们的视线[4]。位于上扬子板块的四川盆地内外发育多套细粒沉积岩，以五峰组—龙马溪组为代表的细粒沉积岩的研究，取得了巨大的成功，该套细粒沉积岩以页岩为主，地域分布稳定、沉积岩层厚度大、有机质含量高、脆性矿物发育，具有较大的页岩气资源潜力[5-6]，而以该套富有机质细粒沉积岩为代表的非常规油气作为国家的重要资源储备物资，对于非常规油气沉积学而言以往的研究主要集中岩石学(如：石英含量、黏土成分等)[7-10]、储层地质学(如：孔隙度、成熟度、渗透率等)[11-13]和展布特征[14]以及产量[15-16]等方面，而对其形成环境和富集模式的研究则略显不足，严重阻碍了非常规油气沉积学的进展。陆扬博等[17]对上扬子五峰组—龙马溪组的地质事件进行了论述并指出火山喷发是五峰组—龙马溪组甜点段形成的一个重要原因，Liet al.[18]对川南地区的四条剖面进行了详细的研究，并指出海平面变化控制下的沉积环境是该套细粒沉积岩有机质富集的主要原因，何龙等[19]则通过地球化学手段对南川地区的五峰组—龙马溪组的形成环境进行了探讨并指出沉积环境是该地区有机质富集的主控因素，严德天等[20]、李双建等[21]在对这套细粒沉积岩研究后认为较高的古生产力是造成有机质富集的根本原因；何利等[22]在对川南马边地区的五峰组—龙马溪组进行研究时，指出有机质的富集除了受控于氧化还原环境的改变还受控于较高的生产力，可见不同研究者对五峰组—龙马溪组有机质的富集模式不尽相同，并且，五峰组沉积厚度较薄，以往学者们常把五峰组和龙马溪组2套细粒沉积岩当做整体进行研究，忽视了在O-S过渡时期赫南特阶沉积的碳酸盐岩对这两套细粒沉积岩在沉积环境乃至富集模式的影响，也有研究者认为四川盆地观音桥组发育并不完全，特别是盆地西南缘代表冰期沉积的该套碳酸盐岩出现缺失[22-23]，故在探讨不同地区细粒沉积岩的有机质富集模式时，不能将两者混为一谈。此外，对于这两套细粒沉积岩有机质富集模式的探讨多集中在盆地内部的川南和川东两个沉降中心[24-28]，但对于靠近古隆起的盆缘细粒沉积岩中所蕴含的环境信息和页岩气等信息的挖掘还远不能与盆地腹部相同地层媲美，这严重制约了人们对早志留世非常规沉积学的客观评价，忽略了富有机质细粒沉积岩的沉积环境、富集模式在非常规沉积学中的重要地位，对非常规油气沉积学的认识和勘探造成了严重的失衡。基于此，本文在四川盆地西南缘盐津牛寨剖面系统采集了能揭示龙马溪组沉积环境和页岩气地质特征的系列样品，并与处于盆地腹部川南沉降中心长宁剖面和非沉降中心南川剖面的龙马溪组进行有机碳、岩石学特征和地球化学特征对比分析，以此阐述盆地西南缘早志留世细粒沉积岩的古环境特征与有机质富集模式，丰富了四川盆地边缘地区非常规油气沉积学的研究内容，为进一步拓展四川盆地边缘非常规资源勘探提供依据。

1 地质背景

研究区位于四川盆地外缘西南部，隶属于上扬子板块西南缘，其构造演化过程与上扬子板块息息相关，自震旦纪以来受多期构造作用的影响形成了现今的盆地格局(图1)。在奥陶世末期—志留世早期，由于华夏板块向扬子板块的强烈挤压作用，造成了川中—黔中—雪峰山三大古隆起的格局，在古陆隆升强制推进作用下形成了“三面环隆”的古地理格局，与此同时，发生了全球性的海侵事件，海平面迅速上升，形成了自川南至川东北为水体深度逐渐加深的扬子陆棚海环境[31-34]，并在川南、鄂西—渝东、川北—川东北3个深水陆棚相沉积区内[35-36]，形成了岩层厚度大、分布广、有机质高的富含笔石化石的龙马溪组细粒岩。

2 剖面概况

盐津牛寨剖面位于云南省盐津县牛寨乡，地层出露良好，地层总厚190.06 m，自下而上依次为上奥陶统五峰组、下志留统龙马溪组和中二叠统梁山组，各组之间均为平行不整合接触(图1)。五峰组主要为碳质页岩夹硅质页岩、生屑灰岩；龙马溪组主体为灰黑色—黑色碳质细粒沉积岩，总厚181.59 m，下部以黑色碳质页岩为主，发育水平层理富含笔石化石；上部以灰黑色碳质页岩为主夹薄层状泥岩，对比陈旭等[37]的笔石带划分标准，研究区龙马溪组笔石相带完整指示研究区早志留世鲁丹阶至埃隆阶沉积序列完整，具有良好的区域可对比性，梁山组则以微晶灰岩为主夹少许劣质薄煤层。

3 样品采集及测试方法

为了探讨上扬子局限海盆内的环境差异及其对有机质富集的影响程度，本文对盐津牛寨剖面志留系龙马溪组细粒沉积岩进行了系统采样，并根据邻区有机质特征和沉积环境的研究结果[22]，对龙马溪组下部进行了加密采样，选取了具有代表意义的样品19件，按照TOC=2.00%为界限，将牛寨剖面龙马溪组划分为上下两段，其中下段采样10件、上段采样9件(图1)，并将这19件样品送至国土资源部西南矿产资源监督检测中心进行岩矿分析和总有机碳(TOC)、主量元素以及微量元素的测试。其中岩石矿物相分析采用ZJ207 bruker D8 advance型X射线衍射仪(XRD)执行SY/T 5163—2010行业标准进行检测，应用High Score软件进行数据分析；主量元素采用X射线荧光光谱仪(XRF)进行检测，分析误差<1%，采用粉末压片法制备样品，先碎岩，然后称取经过高温烘干后的4.5 g样品倒入模具中，在35 t的压力下压制10 s进行制样，然后将样片取下，编号待测；TOC在总有机碳分析仪采用燃烧法进行测试；微量元素采ELEMENT XR等离子体质谱(ICP-MS)仪执行GB/T 14506.30—2010国标进项测定，首先将样品碎至200目以下，并采用加热法除去有机质的影响，然后称取0.10 g样品至于特定容器中并向容器内加入HNO3+HF+HClO4进行密封高温(140℃～150℃)24 h溶样，接着将密封容器敞开放置于加热板上高温加热除去溶样酸，最后向容器内加入HNO3复溶直至溶液完全澄清并将澄清溶液转移至容量瓶中定容，编号待测。

图1 四川盆地及其周缘晚奥陶世—早志留世古地理图及研究剖面柱状图(a)改自文献[29]；(b)改自文献[17，30]Fig.1 Paleogeographic map and research section histogram of Sichuan Basin and its surrounding area from Late Ordovician to Early Silurian

4 实验结果及分析

4.1 矿物岩石学特征

整体而言盐津牛寨剖面中石英含量在之间17.00%～53.00%，平均为33.11%；长石含量为2.00%～18.00%，平均为6.05%；碳酸盐矿物类型以方解石为主，平均含量可达11.58%；黏土矿物存在的主要类型为伊利石，其平均含量可达29.00%；含少量石膏和黄铁矿。垂向变化特征(图2)：龙马溪组从底到顶石英含量逐渐减少，而以绿泥石为代表的黏土矿物呈现逐渐增加的趋势，两者呈现负相关关系，碳酸盐矿物则呈现先减少后增加的两段式特征。

盐津牛寨剖面位于上扬子西南缘，并以长宁剖面和南川剖面为对比剖面，结合岩性特征及矿物在垂向上的变化，对在不同沉积环境下龙马溪组细粒沉积岩主要矿物特征进行分析。据石英＋长石＋黄铁矿、碳酸盐岩和黏土矿物三角图显示(图3)，三个地区的细粒沉积岩在矿物组成上具有明显的差异性，从分布模式上来看盆地西南缘矿物分布呈现陆缘斜坡式的分散状态，长宁地区呈现两点式而南川地区呈现集中式，从矿物成分上来看从川西南至川东北黏土矿物含量逐渐趋于稳定，以石英和长石为代表的陆源碎屑位置其含量总体呈现降低的特点。总体而言，盆地西南缘由于更靠近西边的隆起古陆，故其物质特点呈现杂乱式，在长宁和南川地区其沉积物经过海水的搬运和分选作用，物质成分趋于一致，呈现较为集中的特点，但南川地区较长宁地区而言，由于自川南至川东北为水体深度逐渐加深，导致物源的搬运过程中需要更长的时间和距离，在搬运过程中，黏土矿物发生络合作用和物理吸附作用并进行沉积，因此其黏土矿物含量明显低于川南沉降中心。

图2 四川盆地西南缘盐津牛寨剖面龙马溪组各参数垂向演化Fig.2 Vertical evolution of parameters of Longmaxi Formation in Niuzhai section of Yanjin,southwest margin of Sichuan Basin

图3 上扬子不同地区龙马溪组细粒沉积岩样品矿物组成三角图(a)为四川盆地西南缘，数据引自文献[23]；(b)为长宁剖面数据引自文献[7]；(c)为南川三泉子剖面数据引自文献[38]Fig.3 Mineral composition triangle plots of fine-grained sedimentary rock samples of Longmaxi Formation in different regions of Upper Yangtze

4.2 总有机碳

总有机碳含量(TOC)是衡量细粒沉积岩作为烃源岩的重要指标，国内现多以TOC含量>2.00%作为有效烃源岩的下限值。盐津牛寨剖面TOC含量分布在0.33%～5.87%之间，平均为2.08%，其中龙马溪组下段10件样品TOC含量介于0.87%～5.87%，平均为3.22%，远高于龙马溪组上段数值，为0.33%～1.40%，平均为0.82%。整条剖面在下段出现最大值5.87%，远高于国内有效烃源岩的下限值，整体表现出先快速增加后缓慢减弱的变化趋势(图2)，主要是由于龙马溪组下部为深水陆棚环境，水体安静易于浮游生物及菌藻类的生存，而上部地层水体变浅，陆源物质含量增多，因此在沉积相及地层岩性的影响下导致下段地层中TOC值高于上段。

对比长宁剖面和南川剖面发现，从川西南至川东北TOC含量分布逐渐趋于集中，TOC含量平均值和TOCmin均呈现增大的趋势，但TOCmax分布趋势呈现完全不同的特征(图4)。古地理位置所造成的水体循环差异已不能解释这种情况，此时，初级生产力和氧化还原条件的变化对有机质的富集更为关键。

图4 上扬子不同地区龙马溪组细粒沉积岩总有机碳含量盒型图长宁剖面数据引自文献[39]；南川剖面数据引自文献[19]Fig.4 Box diagram of total organic carbon(TOC)content in fine-grained sedimentary rocks of the Longmaxi Formation in different areas of the Upper Yangtze

4.3 主量元素

对盐津牛寨剖面的SiO2、Al2O3、Fe2O3、CaO、Na2O和TiO2等6种常见主量元素进行测试分析，测试结果表明盐津牛寨剖面细粒沉积岩主量元素分布特征存在一定的差异性(表1)。

自下而上，SiO2、Al2O3、TiO2和Na2O呈前期快速增加后期逐步减缓的变化特征(图2)；而CaO呈前期快速减小后期逐步增加的变化趋势，其规律与碳酸盐矿物变化规律相似，Fe2O3无明显变化趋势，表明海侵事件的发生对主量元素的构成影响明显，其中，SiO2在龙马溪组上、下段平均值分别49.90%、49.47%，整体表现出降低的趋势，Al2O3与TiO2在龙马溪组中的平均含量分别为12.22%、0.52%，且两者的变化趋势呈现良好的正向拟合关系。

4.4 微量元素

在评价海相细粒沉积岩的生产力时，常用生物钡(Babio)作为参考指标，一般用Ba的总量减去陆源碎屑Ba的估算值，可直接反应水体的生产力大小，计算公式为[40]：

式中：Ba样品和Al样品为样品中Ba和Al的含量，μg/g；(Ba/Al)PASS为后太古宙澳大利页岩中的Ba和Al含量比值[41]。

据公式(1)计算结算结果显示龙马溪组细粒沉积岩中Babio含量变化较大(图2)，其中龙马溪组下段为452.60～1 221.25μg/g，平均为911.34μg/g，上段为387.14～1 426.83μg/g，平均为868.39μg/g，明显低于底部黑色碳质页岩段。此外，采用元素比值法如Ba/Al或P/Al也能反映水体的生产力，研究区龙马溪组Ba/Al和P/Al呈现出先增后减的变化趋势，与Babio含量变化趋势相似(图2)。

细粒沉积物在沉积成岩过程中会比其他沉积岩更容易吸附微量元素，为了消除陆源成分的干扰，国际上常采用富集系数(EFs)来快速判别页岩中不同元素的富集程度，计算公式为[42]：

式中：X和Al为样品中某一元素X和Al的含量，μg/g；；PAAS为后太古宙澳大利亚页岩标准化[39]。若XEF>1则表明样品中X元素相对于PAAS页岩更富集，XEF<1则亏损[43]。

据公式(2)计算结算结果显MoEF、UEF在盐津牛寨剖面细粒沉积岩段上的分布存在明显差异，龙马溪组下段MoEF介于0.67～66.01，平均为41.80；UEF介于4.50～37.98，平均为19.82，明显大于龙马溪组上段灰黑色碳质页岩、泥岩段，其中MoEF介于1.15～9.20，平均为5.37；UEF介于3.84～5.58，平均为4.67。

Ni/Co、V/(V+Ni)、U/Th和V/Cr值自下而上在盐津牛寨剖面表现出先增后减的变化趋势(图2)，除氧化还原敏感元素外，Zr元素在古环境中也具有一定的指示意义，该元素在剖面中介于44.90～126.00μg/g之间，平均为84.51μg/g，与主量元素Al2O3和TiO2表现出一致的变化趋势。

5 讨论

5.1 地球化学特征

利用化学元素行为特征的差异性，可以定性—定量地描述海相细粒沉积岩的沉积环境特征，本次研究主要利用相关元素的含量或元素对比值来判定龙马溪组细粒沉积岩形成时的初级生产力、氧化还原条件、陆源供应的输入和盆地水体滞留程度。

5.1.1 初级生产力

在海洋环境中，Ba元素早在2009年就被证实生物来源的Ba可以定量的表征古海洋的生产力[44]，虽然P元素不仅是多数海洋生物组成骨骼的必备元素还是生物在生命活动过程中的营养元素之一[45]，但在缺氧环境下P元素会从有机质中释放出来，造成生产力指标失真[46]，采用Ba/Al或P/Al元素比值更能代表古海洋的初级生产力，因此本文采用Babio值、Ba/Al和P/Al综合判别古海洋的生产力。

表1 盐津牛寨剖面龙马溪组细粒沉积岩部分主量元素、微量元素特征Table1 Characteristics of major and trace elements of fine-grained sedimentary rocks from Longmaxi Formation in Niuzhai section of Yanjin

盐津牛寨剖面的Babio整体呈现与TOC一致的先增后减的趋势，但Babio高值段略微延迟于TOC高值段，Ba/Al和P/Al值除P/Al底部的一个样品外，其余样品的变化趋势与Babio一致(图2)，都呈现出先增后减的特征，但Ba/Al和P/Al值的减弱地起始点明显早于Babio值，与比陆源碎屑输入的判别指标相对比，Ba/Al和P/Al与其呈现出负向相关性。理论上这三者应具有同步变化的特征，但由于盐津牛寨剖面更靠近隆起古陆，受海平面变化的影响，大量陆源碎屑物质的涌入，致使盐津牛寨剖面的Ba/Al值、P/Al值与Babio值减弱出现不同步现象，但整体而言，研究区龙马溪组的生产力呈现出下段高于上段的特征，与有机质的富集层段相吻合，表明生产力是有机质富集的控制因素之一。

5.1.2 氧化还原条件

V、Mo、U等元素在氧化状态下，呈现为高价态离子易于溶解迁移而在还原状态下表现为低价态易于富集沉淀，当海水状态发生变化时这些元素能记录这一变化过程并保存在沉积物中，但也有极个别元素如Th等的溶解状态不受海水的氧化还原状态的控制，但整体而言Ni/Co、V(/V+Ni)、U/Th和V/Cr等参数是与缺氧还原程度呈正相关关系[47-51(]表2)。

表2 氧化还原条件地球化学判识指标[47-51]Table2 Geochemical identification indices of redox conditions

盐津牛寨剖面细粒沉积岩的Ni/Co、V(/V+Ni)和V/Cr平均值分别为6.34、0.73、2.81均指示盐津牛寨剖面处于缺氧—贫氧环境中，虽然U/Th平均值为0.62，但龙马溪组下部的13个样品平均值为0.79，同样处于缺氧—贫氧的环境中。对比四组数据分布模式，发现都具有从缺氧环境逐步过渡到氧化环境的特征(图5)，但在底部的一个样品值明显低于下段其他样品，其中Ni/Co、U/Th和V/Cr值为5.93、0.22、2.47指示氧化环境或接近于氧化环境，表明龙马溪组下段底部氧化还原条件的不稳定，之后由于海侵事件的发生其环境迅速演化为缺氧环境，造成V(/V+Ni)，U/Th和V/Cr值快速正偏，而到了龙马溪组上段由于海侵程度的减弱，相对海平面降低导致Ni/Co、U/Th和V/Cr值逐渐负偏指示海水氧化度逐渐趋于正常。这种变化趋势的与TOC变化趋势以及龙马溪组细粒沉积岩矿物分布的情况相吻合。

对比川西南至川东北的长宁剖面和南川剖面氧化还原参数在垂向上的变化趋势发现(图5)，龙马溪组细粒沉积岩的氧化还原条件都具有从缺氧—贫氧环境逐步过渡到氧化环境的特征，但其变化规律呈现截然不同的特征，盐津牛寨剖面和表现为底部从氧化环境快速过渡到缺氧然后缓慢转化为氧化环境，长宁剖面则从底部缺氧环境逐步过渡到氧化环境，而南川剖面表现为震荡式地从底部缺氧环境转换为氧化环境，造成上述不同现象的原因可能为海平面变化在不同古地理位置的差异响应。

5.1.3 陆源供应输入

陆源碎屑的输入主要用石英含量以及TiO2、Al2O3和Zr值的变化规律来表示，全岩矿物分析结果表明，盐津牛寨剖面石英含量整体自下而上呈现先升高后降低的趋势，与TOC的变化趋势呈现良好的正向拟合关系，而且石英含量超过33.11%的层段主要集中在龙马溪组下部(图2)。主微量元素分析表明，盐津牛寨剖面龙马溪组TiO2、Al2O3、Zr三者之间均表现为良好的一致性，且呈现从龙马溪组下段底部快速增加逐渐过渡到上部缓慢降低(图2)上述现象均反映了龙马溪组下段陆源供应量明显增加。

对比川西南至川东北方向的长宁剖面和南川剖面TiO2、Al2O3在垂向上的变化趋势发现(图6)，盐津牛寨剖面与长宁剖面都具有TiO2、Al2O3从底部到顶部呈现先增加后减小的趋势，表明陆源碎屑的供应量先增强后减弱，但从其增幅来看，盐津牛寨剖面明显大于长宁剖面，相比之下说明盆地西南缘在早志留世处于更靠近岸边、陆源输入物质更强的古地理位置。与上述变化规律相比南川剖面则表现为陆源碎屑的供应量呈现较弱的减缓趋势，说明南川地区在早志留世处于离岸较远、陆源物质输入较弱的古地理位置。综上这种现象与前文所论述的古地理位置状况相印证，盐津牛寨剖面更靠近剥蚀区故陆源物质的接受量明显较大，相对而言长宁剖面和南川剖面离剥蚀区最远，且自川南至川东北为水体深度逐渐加深，导致该地区的陆源碎屑注入量记录呈现微弱的减缓趋势。

图5 上扬子不同地区龙马溪组细粒沉积岩氧化还原条件判别图长宁剖面数据引自文献[39]；南川剖面数据引自文献[19]Fig.5 Identification of redox conditions for fine-grained sedimentary rocks of the Longmaxi Formation in different areas of the Upper Yangtze

图6 上扬子不同地区龙马溪组陆源碎屑输入变化图长宁剖面数据引自文献[39]；南川剖面数据引自文献[19]Fig.6 Variation of terrigenous debris input in the Longmaxi Formation in different areas of the Upper Yangtze

5.1.4 盆地水体滞留程度

U、Mo元素在上地壳的含量仅为2.7 ug/g、3.7 ug/g，相对含量较低[41]，在海洋生物中含量也很低，而海洋环境中的U、Mo元素主要通过陆地上的河流输入补充，因此富集在沉积物中的U和Mo元素一般来源于海水[52]。沉积过程中松散沉积对U和Mo元素的摄取开始时间不同，对U的摄取发生在氧化还原反应界线处，而对Mo的摄取则需要H2S的存在，故沉积物对U的摄取早于Mo元素，并且深海环境下的铁锰氢氧化物的吸附作用会吸附海水中的Mo元素而进入沉积物中，而U元素则不受吸附作用的影响[51,53]，基于两者地球化学行为的差异，现多采用Mo/TOC值和U-Mo协变模式来判别水体的滞留程度。

盐津牛寨剖面的Mo/TOC值介于0.66～10.08，平均为6.35，除一个样品可能因还原环境造成Mo的富集而升高外[19]，其余样品均说明了在早志留世盆地西南缘处于中等烈滞留环境中(图7)。从盐津牛寨剖面的U-Mo协变模式图可以看出UEF、MoEF值较小(图8)，Mo/U比值也小，主体落在0.3-1×SW区间内，从贫氧—缺氧环境UEF/MoEF值由0.3×SW逐步过渡到1×SW，随着水体逐还原程度的加剧UEF、MoEF值逐渐增大UEF/MoEF值却出现降低的趋势，表示海水中的Mo元素随着沉积作用的进行，得不到补充，使得UEF/MoEF值降低，这种变化模式介于东太平洋和现代黑海的UEF-MoEF协变模式之间(图8)，表明当时盆地西南缘处于一个中等滞留海盆环境中，符合Mo/TOC值得出的结论，也验证了上述龙马溪组从底部到顶部还原程度逐步减弱的趋势。

图7 上扬子不同地区龙马溪组细粒沉积岩TOC-Mo相关图黄线为盐津牛寨剖面线性回归线，长宁剖面数据引自文献[39]；南川剖面数据引自文献[19]Fig.7 TOC-Mo correlation of fine-grained sedimentary rocks of Longmaxi Formation in different areas of Upper Yangtze

图8 上扬子不同地区龙马溪组细粒沉积岩UEF-MoEF协变模式图长宁剖面数据引自文献[39]；南川剖面数据引自文献[19]Fig.8 UEF-MoEF covariation model of fine-grained sedimentary rocks of Longmaxi Formation in different areas of Upper Yangtze

盆地水体的滞留程度在川西南至川东北方向并没有表现出明显的差异性，从TOC-Mo相关性图中(图7)可以看出，当TOC含量<2.00%时，除牛寨剖面一个样品落在Mo/TOC＝4.5界线之下、其余样品均分布在Mo/TOC=9～4.5区间内，表明水体滞留程度呈现出相似的状态；当TOC含量>2.00%时，盐津牛寨剖面、南川剖面和长宁剖面Mo/TOC值基本绕Mo/TOC=9界线分布，仍然表明这三个地区的水体滞留程度相似，都为中等滞留环境。其次，对比U-Mo协变模式图(图8)，也发现盐津牛寨剖面、南川剖面和长宁剖面都具有UEF/MoEF值在高值段出现降低的趋势，表示Mo元素补充量不同于非滞留海盆那样持续补充也不同于强滞留环境下海水中的Mo沉积后得不到补充，再次说明这三个地区所在的扬子陆棚海在早志留世处于一个中等滞留的沉积环境中。综上表明，在早志留世的“三面环隆”的地质背景下，海平面的升降并没有明显影响水体的滞留程度，并且总体属于相对开放的海域环境即中等滞留的水体环境。

5.2 有机质富集机理

目前国际上在有机质富集的模式主要有三种模式：生产力模式、保存条件模式、生产力—保存条件协同控制模式，生产力模式认为生产力的高低是有机质富集的决定因素；保存条件模式强调在缺氧环境下不利于有机碳的分解，从而沉积物中会保存较多的有机质；而协同模式则认为有机碳的埋藏受上述条件的共同控制[54-56]。何龙等[19]认为氧化还原条件是影响南川非沉降中心细粒沉积岩有机质富集的主要因素，Liet al.[18]则认为海平面升降造成的氧化还原条件变化控制了川南沉降中心粒沉积岩有机质富集，选取Babio、Ni/Co、Al2O3、Mo/TOC代表初级生产力、氧化还原条件、陆源碎屑物质的注入和盆地水体的滞留程度分别与TOC进行线性拟合(表3)，发现无论是盆地西南缘还是川南沉降中心和南川非沉降中心，TOC与盆地的水体滞留程度相关性都极低，表明水体的滞留程度不是有机质的富集的主控因素，综合对比发现，川南沉降中心显示TOC与氧化还原条件呈现显著的正向相关关系，与初级生产力和陆源碎屑物质的注入呈负向相关关系，符合李艳芳提出的海平面变化—保存条件模式；南川非沉降中心与上述四个因素的相关系数均为正值，但TOC-Ni/Co相关系数远高于其余三者，印证了何龙提出的保存条件模式。

表3 不同研究区相关参数与TOC拟合度Table 3 Fitness of related parameters and TOC in different study areas

盐津牛寨地区TOC与不同参数之间的相关性投点图显示(图9)：TOC与氧化还原判别参数Ni/Co、V/Cr以及U/Th都呈现较好的正向相关，相关系数R2分别为0.61、0.49和0.61；与陆源碎屑输入判别参数Al2O3、TiO2和Zr的相关系数R2分别为0.12、0.06和0.00表明陆源碎屑的输入对该地区有机质的富集没有显著的影响；与初级生产力判别参数Ba/Al、Babio和P/Al的相关系数分别为0.27、0.02和0.04，其相关系数都较低，表明古海洋初级生产力不是有机质富集的主控因素。对比长宁剖面和南川剖面的TOC与各控制参数的相关性(表3)，发现盆地西南缘的表现形式和南川非沉降中心较为相似，都表现为TOC与氧化还原条件呈现显著的正向相关关系，与初级生产力、陆源碎屑的输入无明显相关性(图9)，综上，缺氧还原的沉积环境是影响研究区龙马溪组黑色页岩有机质富集的主要因素

图9 四川盆地西南缘盐津牛寨剖面龙马溪组细粒沉积岩TOC质量分数与不同参数相关性Fig.9 TOC of fine-grained sedimentary rocks of Longmaxi Formation in Niuzhai section of Yanjin,southwest margin of Sichuan Basin

在早志留世全球气候变暖，形成于晚奥陶世的冰川逐渐消融，在上扬子地区发生大面积海侵[57]，导致早龙马溪早期水体分层，底层水迅速转变为中等滞留—缺氧状态(图10a)，由于海平面的升高，处于盆地边缘的盐津牛寨地区水体滞留程度增加，同时由于大面积的海侵发生带来较多的笔石等浮游生物和菌藻类，随着生物代谢活动的进行这些浮游生物和菌藻类进入沉积物中进行埋藏，随着埋藏作用的进行，不断消耗本就处于缺氧环境下的溶解氧，最终形成缺氧环境，故此龙马溪组下段富集了大量的有机质。到了龙马溪组中—晚期随着海平面的下降(图10b)，原本处于缺氧环境的底层水随着海平面的降低含氧量增加，同时海水裹挟大量的陆源碎屑物进入盐津牛寨地区，一方面直接稀释了有机质的浓度，另一方面破坏了原有浮游生物及菌藻类的生存环境，导致生物丰度下降，不再消耗底层水的溶解氧，导致细粒沉积岩中上段有机碳含量的降低[58-59]。