第一作者简介 刘招君,男,1951年生,教授,主要研究方向为油页岩成矿理论与资源评价。E-mail: liuzj@jlu.edu.cn。
油页岩作为一种重要的非常规油气资源,如何有效地利用油页岩成矿作用研究中的关键方法和技术实现油页岩的综合识别与评价极其重要。油页岩识别的关键技术主要包括基于岩石学的地质识别、精确的地球化学指标识别、测井技术的垂向高分辨率快速识别和地震技术的三维空间定量刻画。先进的测试技术、创新的理论方法以及高品质地球物理资料的二次开发利用将不断提升油页岩从定性到定量预测识别的精度。油页岩成矿作用研究的关键方法是从古构造条件、古气候条件、古地理条件以及古湖泊环境等地质因素入手,探讨其综合作用下有机质来源、保存条件及成岩作用的变化,由定性向半定量—定量地精细刻画油页岩中有机质聚集过程,从而有效地预测油页岩的展布及其富集规律。由于文中材料以中国陆相含油页岩盆地的研究成果为主,而陆相成因油页岩的识别特征及成因机制可能与海相等成因存在差异,因此,仍需油页岩地质学家不断深入研究与总结。希望此文对页岩油、页岩气及烃源岩研究也具有一定的借鉴意义。
About the first author Liu Zhao-Jun,born in 1951,is a professor of College of Earth Sciences,Jilin University. Now he is engaged in metallogenetic theory of oil shale and resource evaluation. E-mail: liuzj@jlu.edu.cn.
Oil shale is an important unconventional oil & gas resource. It is critical to comprehensively identify and evaluate oil shale recourse by utilizing key methods and technologies in the study of oil shale mineralization. The key technologies of oil shale identification include geology identification based on petrology,geochemistry identification,high-resolution logging,seismic technique for three-dimensional space description,advanced test technique,innovated theory & method and secondary exploitation of high-quality geophysics information,all of which will continuously improve the accuracy of oil shale identification and forecast from qualitative to the quantitative. The key methods in the study of oil shale mineralization are firstly to find out the geological factors,including palaeostructure,palaeogeographic and palaeoclimate conditions,and palaeolake environments,to discuss a combined effects on various organic matter sources,preservation conditions and diagenesis. Then,finely depict the process of organic matter enrichment of oil shale from qualitative to the quantitative,thus effectively forecast the patterns of oil shale formation and distribution. This paper was mainly based on the research results of continental oil shale-bearing basins in China. Due to the diversity of oil shale genetic type,there might exist differences of identifying characteristics between continental origin and marine origin oil shale,which still needs to be further studied and summarized. It is hoped that this paper could also be useful for researches an shale oil,shale gas and hydrocarbon source rocks.
油页岩是重要的非常规油气资源, 可以用来直接燃烧或者通过干馏得到油页岩油及相关产品, 具有极大的经济和社会效益。中国油页岩资源丰富, 2003— 2006年“ 全国油页岩资源评价” 结果表明, 埋深不大于1000 m油页岩油可采资源为120× 108 t, 相当于常规可采石油资源量的一半以上(刘招君等, 2009a, 2009b)。近年来, 国际油价的持续波动, 中国很多地方都掀起了油页岩勘探开发热潮, 促进了油页岩地质调查工作的发展, 地质找矿成果显著。但随着油页岩勘探难度的增大, 从以前的露头矿、浅层矿到目前的深层找矿, 以及原位开发技术的不断发展, 如何有效识别并准确地确定油页岩矿层位及展布规律, 是油页岩勘探与开发的关键地质任务, 为此, 笔者及其团队经过10余年来的理论与实践研究, 总结出以陆相盆地为主的油页岩成矿作用研究中的一些关键方法和技术, 希望能对今后的油页岩资源评价与成矿预测工作有所借鉴。
油页岩的识别是油页岩资源评价及成矿作用研究的一个关键问题。油页岩在岩石学特征上是一种富含有机质的细粒沉积岩, 其可靠识别主要依据野外露头与钻孔样品测试数据实现, 一般可以利用含油率、总有机碳含量(TOC)与产烃潜量(S1+S2)等指标来识别油页岩(贾建亮等, 2015)。然而, 鉴于样品数量、样品分布、样品代表性和测试分析效率等条件限制, 在垂向及横向上都难以实现客观全面地识别和评价盆地油页岩(贾建亮, 2012)。因此, 如何利用先进的地球物理技术实现油页岩的快速定量识别与评价也是油页岩综合研究的重要课题之一。
油页岩的综合地质与地球化学识别技术主要是基于油页岩的岩石学与地球化学特征, 是一种最直接可靠的油页岩识别方法, 而明确油页岩本身的含义是进行油页岩识别的前提。2006年“ 全国油页岩资源评价” 中定义了油页岩, 认为油页岩是一种高灰分的固体可燃有机沉积岩, 低温干馏可获得油页岩油, 含油率大于3.5%, 有机碳含量较高, 有机质类型主要为Ⅰ 型和Ⅱ 1型干酪根, 其发热量一般不小于4.18 mJ/Kg, 其工业指标在开发利用过程中可能随经济和技术条件的变化而变化(刘招君等, 2009a, 2009b, 2013; 孟庆涛, 2010; 贾建亮, 2012)。
1.1.1 油页岩的地质识别特征
油页岩的地质识别主要通过野外露头及钻孔岩心的最直观的岩石学特征进行。任何岩石分类都应该具备易被识别的关键岩石组成特征, 对于油页岩而言其特征应该是有机质, 主要包括有机质类型与丰度(Hutton, 1987; 贾建亮, 2012)。油页岩极其富含有机质使其在颜色、岩性、结构、构造和成分等方面区别于常规泥岩、低品质烃源岩、煤以及碳质泥岩等(刘招君等, 2009a; 孟庆涛, 2010; 贾建亮, 2012), 也使油页岩具有显著的地球化学与地球物理特征。
油页岩的岩石类型呈现多样性, Π о г р е б н о в 等(1983)认为油页岩是固体可燃矿产的一种岩石类型且是化学成因— 碎屑成因的沉积岩; Dyni(2003)认为油页岩是一种细粒沉积岩并包含大量有机物; 刘招君等(2009a)认为油页岩是一种高灰分的固体可燃有机沉积岩。由此可见, 世界上油页岩的岩石类型主要由富含有机质的页岩(泥岩)、泥灰岩和碳酸盐岩等岩石组成, 而中国油页岩主要以页岩或泥岩为主(刘招君等, 2009a)。
油页岩在野外露头或钻孔岩心上表现出质轻、具油腻感, 用指甲刻划时划痕呈暗褐色, 并且用小刀沿层面切削时, 常呈刨花状薄片, 用火柴点燃时冒烟且具油味等特征(图1)(刘招君等, 2009a, 2009b; 孟庆涛, 2010; 贾建亮, 2012)。油页岩手标本及岩心颜色主要为黑色— 灰黑色、褐黑色、黑褐色、灰褐色、棕褐色、褐黄色或深灰色(图1), 且其颜色随有机碳含量不同变化比较明显(刘招君等, 2009b)。油页岩多具泥质结构。油页岩的原生构造主要包括层理和块状构造, 层理主要以发育稳定水体环境中的水平层理为主(图1), 块状构造仅仅是用肉眼观察看不到任何内部构造, 但微观尺度可能存在一定的生物扰动构造(刘招君等, 2009a)。此外, 在油页岩上也可见到软体动物和植物等化石分布, 还可见到树皮结构、母质结构的镜质体。
中国典型陆相含油页岩盆地(如松辽盆地、桦甸盆地和抚顺盆地等)发育极其优质的油页岩, 颜色以深灰色、灰黑色、灰褐色以及暗褐色为主(图1-c, 1-e, 1-f), 含油率最高可达20.0%, 有机碳含量最高可达30.0%, 有机质类型以富含藻类与浮游生物的Ⅰ 型和II1型干酪根为主, 而成熟度普遍较低(RO< 0.5%)(柳蓉, 2007; 孟庆涛, 2010; 贾建亮, 2012; Jia et al., 2013a, 2013b; 孙平昌, 2013; Xu et al., 2015; 徐进军, 2015)。油页岩以泥质结构为主, 主要为黏土级别极细沉积物, 发育块状层理和水平层理, 有时含有介形虫类、叶肢介类等生物化石、碎屑条带及方解石、凝灰岩等夹层(图1-e, 1-f)。
1.1.2 油页岩的地球化学识别特征
含油率是定量识别与评价油页岩的最直接参数, 一般油页岩含油率边界品位大于3.5%, 并以5%与10%为界分为低、中、高3个品级(刘招君等, 2009a)。油页岩含油率的直接获取一般通过铝甑干馏测试方法, 而获取的油页岩油主要来自富烃有机显微组分, 如腐泥无定形体、藻类体、壳质组分和富烃镜质体等, 从而使油页岩油具有较高的转化效率, 而较低成熟度是保证其高效转化的基本保证。一般而言, 稳定水体环境中以浸油型干酪根为主要贡献的TOC含量往往与含油率的关系密切, 从而可以通过建立含油率与TOC含量之间的定量关系来识别与评价油页岩。此外, 以游离烃S1与裂解烃S2组成的产烃潜量的热升温获取机理无疑可以与含油率和TOC含量建立直接关系, 从而也使其成为一个判定油页岩地球化学识别的有效参数。因此, 含油率、TOC含量与产烃潜量是3个极其重要且具有内在成因联系的用地球化学方法识别与评价油页岩的关键参数。
要建立含油率、TOC含量与产烃潜量三者之间稳定可靠的定量关系模型, 需要做大量的测试数理统计。统计数据要满足以下条件: (1)不同地区具有不同成因类型的油页岩应分别统计, 也就是不同油页岩油转化效率会影响其对TOC含量与产烃潜量的贡献; (2)为避免构造埋深因素对油页岩形成的影响, 数据最好来自研究区不同位置的井; (3)为消除数据层段不全带来模型的不稳定性, 数据最好来自连续取心钻井(贾建亮, 2012)。
中国典型陆相大型湖盆— — 松辽盆地中发育的油页岩具有典型代表性, 稳定沉积背景下的巨大发育规模、高有机质富集、良好的有机质类型和较低有机质成熟, 使其能够在TOC含量与产烃潜量一定时具有较高的油页岩油转换效率, 从而为建立含油率、TOC含量与产烃潜量之间的稳定关系模型奠定了基础。为此, 选取松辽盆地不同构造单元4口全井段取心钻孔的青山口组层段(贾建亮, 2012; 孙平昌, 2013; 徐进军, 2015), 分别建立了TOC含量与含油率、TOC含量与产烃潜量之间的定量关系(图2)。需要指出的是, 本次只选取了TOC含量大于0.6%代表烃源岩的岩石样品进行统计建模。统计建模结果表明, 松辽盆地青山口组油页岩的TOC含量分别与含油率和产烃潜量呈现明显的正相关性, 相关系数(R2)均超过0.93, 说明TOC含量与含油率和产烃潜量均具有良好的相关性(图2)。由此可知, 松辽盆地油页岩含油率为3.5%分别对应TOC含量为5.0%、产烃潜量为35.0 mg HC/g rock。因此, 在松辽盆地把油页岩含油率大于3.5%、TOC含量大于5.0%、产烃潜量大于35.0 mg HC/g rock作为识别油页岩的标准, 从而可以在松辽盆地较好地运用含油率、TOC含量与产烃潜量3个参数实现对油页岩的地球化学定量识别与评价。
油页岩的地球物理识别技术因其快速有效性而备受关注(谭廷栋, 1988; Passey et al., 1990; 贺君玲等, 2006; El Sharawy和Gaafar, 2012; 张佳佳等, 2012; Jia et al., 2012), 其评价技术也因业内对页岩油气的广泛关注而得以发展。油页岩作为极富有机质的细粒沉积, 其地球物理识别更具有代表性。油页岩的地球物理识别主要从单井测井评价与空间地震评价2方面入手, 首要解决的是地球物理响应机理问题, 这是确定地球物理参数是否具有明确地质意义的前提。
1.2.1 油页岩岩石物理建模及地球物理响应
油页岩地球物理综合解释的任务主要把特定地球物理参数反演成识别油页岩的地质信息, 为此需要建立地球物理解释模型, 从而确立地球物理响应与地质参数之间的关系(孙建孟和王永刚, 2001; 贾建亮等, 2015)。常规岩石物理模型主要针对纯砂岩或泥质砂岩模型, 而非常规富含有机质的泥页岩模型研究相对较少(Meyer and Nederlof, 1984; 赵彦超, 1990)。
在油页岩岩石物理建模时, 建议把油页岩固体部分划分为非黏土矿物、黏土矿物与有机质3种组分, 并把其作为3个宏观物理量研究岩石每部分对地球物理响应值的贡献。地球物理响应的贡献不单纯看各固体部分的物理性质, 更要注重不同岩石组分所构成的孔隙结构及流体性质所起作用。大量模拟实验表明, 不同孔径大小的孔隙并非独立存在, 而是以树枝状等不同的方式相互连接, 构成了沉积岩石内部微观的流体运移通道(Nelson, 2009), 从而影响岩石的地球物理响应特征。黏土矿物和有机质赋存与作用方式的多样性使富有机质泥页岩孔隙度的定义变得极其复杂(Kuila et al., 2014)。Passey 等(2010)和Kuila等(2014)把富有机质泥页岩总孔隙度定义为泥页岩饱含流体的含量, 主要包括自由水、微毛细管水、黏土束缚水和烃类流体的孔隙, 而有效孔隙度被进一步定义为自由水、微毛细管水与烃类流体的孔隙度之和。前人研究表明, 富有机质泥页岩孔隙中占优势的是中孔(孔径2~100 nm)和微孔(孔径小于2 nm)(黄振凯等, 2013; 杨峰等, 2013), 即地球物理意义上的微毛细管孔隙与黏土水孔隙, 但不能忽略自由水孔隙的存在。Li 和Mo(1999)、莫修文等(2001)等认为完全饱水岩石的导电是自由水、微毛细管水与黏土束缚水3部分起作用, 在渗流特性上微毛细管水与黏土束缚水一样, 在通常压力下无法流动和产出, 而在导电性上微毛细管水与自由水不同于黏土束缚水的附加导电性。因此, 合理的油页岩宏观物理组分划分, 对于建立具有地质意义的油页岩岩石物理模型至关重要。
油页岩的地球物理综合响应是各岩石组分及孔隙流体某物理性质的贡献总和, 即某物理量的加权平均值, 而异常地球物理响应是油页岩地球物理识别的关键(贾建亮等, 2015)。常规砂岩储集层地球物理识别主要利用地层孔隙度及其孔隙流体类型变化引起特定地球物理参数(如电阻率、声波、中子、密度和相对波阻抗等)相对异常的响应来判断(李勇根和徐胜峰, 2008), 而油页岩地球物理识别同时要考虑有机质高度富集对各个地球物理参数响应所起作用。因此, 在地层条件一定的前提下, 油页岩的岩石组合特征、微观孔隙发育及其孔隙流体性质将成为重点影响岩石地球物理响应的因素(贾建亮等, 2015)。
油页岩岩石组分的测井响应与弹性参数分析表明, 有机质与黏土矿物具有低地层密度、高声波时差、高中子孔隙度、高自然伽马和高电阻率的测井响应特征以及低体积密度和低纵波速度的弹性参数特征, 而非黏土矿物(如碎屑矿物、碳酸盐岩矿物和重矿物等)表现出与有机质相异的地球物理响应特征(谭廷栋, 1988; 洪有密, 1993)。因此, 油页岩的地球物理异常响应主要受有机质丰度、黏土矿物作用与微观孔隙发育程度影响, 从而把油页岩岩石物理模型简化成4种有效组成部分: 非黏土矿物、黏土矿物、有机质和总孔隙度, 据此建立油页岩等效岩石物理模型极其地球物理响应方程(Passey et al., 2010; 贾建亮, 2012; 贾建亮等, 2015)。实践表明, 松辽盆地青山口组与嫩江组油页岩较高的有机碳含量、黏土矿物含量以及束缚水孔隙等特征(贾建亮, 2012), 使其具有显著的测井参数与地震波阻抗响应, 从而为其测井与地震识别奠定基础。
1.2.2 油页岩的测井识别技术
综合地质与地球化学方法识别油页岩有一定的局限性, 而高垂向分辨率的测井信息能为油页岩识别提供更加全面高效的资料支持, 油页岩的测井识别技术对预测结果的可靠性起到至关重要的作用。油页岩的测井识别主要利用测井资料的高纵向分辨率特征以及油页岩富含有机质的异常测井响应特征, 建立特征测井参数与油页岩有机碳含量之间的定量关系模型, 并进一步计算出油页岩段有机碳含量的连续分布。
油页岩或富有机质泥页岩的测井识别从提出利用单参数识别到现今已经发展了70多年, 特别近年来随着对页岩油气关注度的不断提高, 如何快速有效地识别并评价优质泥页岩是广大科研工作者致力研究的方向之一。油页岩或富有机质泥页岩的测井识别技术发展主要可以划分成以下几个阶段: (1)单测井参数识别, 主要利用在特定地区具有特征指示意义的测井参数定性判别或定量计算有机碳含量高低。例如: 利用有机质和黏土矿物对放射性元素铀的高吸附性来建立与有机碳含量的关系(Beers, 1945; Swanson, 1960); 通过有机质与自然伽马参数之间的密切关系识别富有机质泥页岩(Schmoker, 1981; 谭廷栋, 1988; Hertzog et al., 1989); 还有利用密度测井方法确定富有机质泥页岩中有机碳含量(Schmoker, 1979; Schmoker and Hester, 1983); 也有通过声波测井参数、电阻率测井参数在特定地区进行富有机质泥页岩的定性识别评价(Autric, 1985; Hussain, 1987)。(2)多测井参数识别, 利用多种与有机质测井响应密切相关指标通过不同算法实现判定富有机质泥页岩, 避免因单参数受多地质因素叠加影响引起预测结果的不稳定性。例如: 利用双测井参数交汇图建立二元一次方程来判别烃源岩和非烃源岩(Meyer and Nederlof, 1984; 赵彦超, 1990); 通过多种测井参数建立多元线性回归方程定量评价富有机质泥页岩(Fertl and Chilingar, 1988; Kamel and Mabrouk, 2003; 徐思煌和朱义清, 2010); 还有利用人工神经网络法模拟预测烃源岩以解决测井曲线和TOC间的隐式问题(Huang and Williamson, 1996; Kamali and Mirshady, 2004; Kadkhodaie-Ilkhci et al., 2009; 张晗等, 2017)。(3)△ logR识别技术, 由Passey等(1990)提出, 无疑是多测井参数评价烃源岩中最具代表性及广泛应用的一种方法, 主要是通过电阻率和三孔隙度测井指标的重叠来定量评价富有机质泥页岩, 该模型建立具有明确的物理意义而不只是简单的线性回归, 尽管Δ logR模型在叠合系数求取、基线值选取、有机碳含量背景值、成熟度参数方面存在不足(胡慧婷等, 2011; 卢双舫等, 2016), 但不影响后人对Δ logR模型的积极改进与推广应用(许晓宏等, 1998; 张志伟和张龙海, 2000; 朱振宇等, 2003; Kamali and Mirshady, 2004; 贺君玲等, 2006; 霍秋立等, 2011; El Sharawy and Gaafar, 2012), 特别是近年卢双舫团队等在优化Δ logR模型与推广应用方面做了大量积极有效工作(胡慧婷等, 2011; 黄文彪等; 2014; 刘超等, 2015; 卢双舫等, 2016; 陈海峰等, 2017), 不断优化的Δ logR模型使有机碳含量预测精度进一步提高(卢双舫等, 2016)。需要特别指出的是, 上述富有机质泥页岩的测井识别技术应用进程中, 绝大多数典型实例中都有油页岩层段发育(例如Passey et al., 1990)。油页岩因其极富有机质的岩石学特征使其具备更加显著的测井响应, 从而为富有机质泥页岩测井识别技术的应用与推广提供更加有力的证据。因此, 贺君玲等(2006)、Jia等(2012)和朱建伟等(2012)和贾建亮(2012)等利用多元统计法、Δ logR模型开展油页岩的测井识别时获得了较好的应用效果。
尽管如此, 测井曲线易受井眼环境等非地层因素的影响而不能准确地反映原状地层信息, 故不受地层外因素影响的地震属性组合(特别是缺乏测井资料地区)结合有效测井信息来综合评价富有机质泥页岩将受到关注, 从而试图利用测井— 地震多属性逐步回归分析定量估算有机碳含量(贾建亮, 2012)。测井— 地震多属性方法把Δ logR模型(或改进模型)预测的TOC含量作为目标曲线, 采用步— 进递归法从众多测井— 地震属性中按属性对TOC含量作用的显著程度优选影响显著的测井— 地震属性组合, 然后利用多元逐步回归分析方法与褶积因子算法对优选的测井— 地震多属性进行信息融合获得TOC含量(Jia et al., 2012; 贾建亮, 2012)。在松辽盆地临江南三维地震研究区的实践应用表明, 多属性预测TOC相比Δ logR预测TOC含量与测试TOC含量具有更好对应关系(图3), 预测结果含有更多反映TOC含量变化的属性信息, TOC含量背景值在垂向偏离零线, 基线变化合理而更符合原状地层有机质的垂向变化规律(Jia et al., 2012; 贾建亮, 2012)。因此, 测井— 地震多属性方法在基础资料综合运用与优选融合属性算法两方面均具备提高单井TOC预测精度的优势, 从而使其在具有丰富地球物理资料特别是地震资料的地区具有广泛的应用前景。
1.2.3 油页岩的地震识别技术
尽管测井识别油页岩能获得精度相对较高、垂向分辨率较高的有机质变化曲线, 但却无法解决有机质连续横向变化及分布预测问题。因此, 利用空间密集的地震资料及技术(地震反演和地震属性分析等)进行三维有机质空间定量刻画, 是实现油页岩或富有机质泥页岩三维空间定量识别评价的一个必然发展趋势。
受地震分辨精度及其可靠性的影响, 油页岩或富有机质泥页岩的地震识别的研究相对较薄弱。近年随着页岩油气勘探的需要, 富有机质泥页岩的地震评价不管是基于叠后地震资料还是叠前地震资料都得以发展。油页岩或富有机质泥页岩的地震识别技术发展主要可以划分为以下几个阶段: (1)叠后地震资料定性识别, 主要在低勘探程度盆地应用地震相法、岩性速度法及地震属性法对优质烃源岩进行定性早期预测(于建国等, 2005; 刘震等, 2007; 张寒和朱光有, 2007; 曹强等, 2009; Lø seth et al., 2009)。(2)叠后地震资料定量评价, 主要在勘探程度相对较高盆地, 结合测井与地球化学等资料, 基于岩石物理建模利用地震多属性分析或地震反演方法来定量评价烃源岩(Lø seth et al., 2011; 金强等, 2012)。(3)叠前地震资料定量评价, 主要通过叠前弹性参数反演与叠前扩展弹性参数反演直接获得弹性参数、密度和自然伽马数据体, 利用实测TOC含量建立与各类地震敏感参数关系估算TOC含量(陈祖庆, 2014; 侯华星等, 2016, 2017; 何沂等, 2018)。然而, 受地震资料分辨率以及地震响应强度的影响, 一般烃源岩的地震识别在界面反射型数据体中往往难以获得理想的效果, 而有机质高度富集的油页岩更易形成显著的界面反射特征(贾建亮等, 2015), 从而使油页岩的地震识别比烃源岩的更具有典型代表性。张佳佳等(2012)、Jia 等(2012)、贾建亮等(2015)等利用叠后地震正演模拟、地震多属性分析以及地震反演方法等识别与评价油页岩获得了较显著的应用效果。
针对松辽盆地临江南青山口组和嫩江组的未成熟油页岩, 结合单井估算TOC含量数据, 利用叠后地震资料正演— 反演方法识别与评价油页岩的研究结果表明, 有机质丰度和孔隙发育程度影响油页岩的地震响应, 地震资料可分辨最小油页岩厚度为6.4 m(贾建亮等, 2015), 波阻抗反演结果只能定性评价油页岩发育的空间展布特征, 而有机质与含油率反演结果不仅能预测油页岩发育的空间展布特征, 还可以评价油页岩发育的品质(图4)(Jia et al., 2012; 贾建亮, 2012), 且地震预测与钻孔测试结果吻合度较高(相关系数R2=0.796)(贾建亮等, 2015), 从而实现油页岩的地震定量识别与评价。因此, 油页岩的地震识别方法能够精细刻画地质体中油页岩的空间展布规律, 从而使其在地震资料品质较高且油页岩发育的地区具有更好的适应性。
油页岩成矿作用研究是进行油页岩地球物理识别、资源评价与成矿富集规律研究的基础, 油页岩的形成必须具备丰富的有机质来源、良好的有机质保存条件以及后期的矿藏保存条件, 归结起来即是古构造、古地理、古湖泊环境和古气候条件等地质因素综合作用的结果(刘招君等, 2009a)。
根据板块的研究成果, 中国含油页岩盆地可简单地归纳为4类, 即裂谷盆地、前陆盆地、走滑盆地和克拉通盆地。不同盆地类型的不同演化阶段其沉积环境有所差异, 不同的沉积环境造成了有机质来源、保存条件和有机质埋藏的差异, 由此导致不同类型油页岩的沉积。因此, 油页岩的沉积环境是影响有机质富集的直接关键因素。借助露头、钻井岩心、测井和地震等资料, 通过含油页岩岩系的沉积环境研究, 可以确定有利于油页岩发育的沉积相, 分析其沉积环境特征, 为进一步的资源预测及研究有机质富集时的古湖泊环境提供依据。古湖泊的沉积物填充样式和沉积相空间展布规律主要受控于其层序地层的基准面变化, 基准面又是古气候、古构造、沉积物供给速率和湖平面变化的综合响应结果, 最终影响可容纳空间变化(Vail, 1987; 刘招君等, 2002), 故含油页岩岩系层序地层研究是确定油页岩形成环境及展布规律的关键所在。
通过研究典型盆地含油页岩岩系层序地层与沉积特征(刘招君等, 2011; 孙平昌等, 2011; 贾建亮等; 2014; 徐进军, 2015)发现, 油页岩层底界面往往可以作为准层序界面(刘招君等, 2011)(图5-q), 油页岩主要发育在三级层序水进体系域和高水位体系域的深湖、半深湖及湖沼(与沼泽有成因联系的浅湖)环境。深湖腐泥型油页岩主要赋存于大中型拗陷盆地中, 往往具有分布范围广、空间展布稳定、品质和厚度中等的特点。在一些小型断陷盆地中也存在此种类型的矿床, 具有分布范围小、横向差异大、品质良好以及厚度巨大的特点; 半深湖腐殖腐泥型油页岩在断陷与拗陷盆地中均有发育, 其特点是分布范围小、空间上品质和厚度差异明显; 湖沼腐泥腐殖型油页岩主要沉积于断陷盆地中, 且与煤共生。往往在1个准层序中, 下部发育油页岩, 上部发育煤, 油页岩具有分布局限、品质好和厚度薄的特点(刘招君等, 2016)。
目前, 中国已发现的绝大多数油页岩均为湖相沉积, 油页岩作为一种典型的高富含有机质的细粒湖相沉积物, 其古湖泊环境的研究是揭示油页岩有机质富集机制及成矿环境的一种重要手段。古湖泊环境研究主要包括古湖泊生产力、湖泊水体性质和湖水分层等。
2.2.1 古湖泊生产力
古湖泊生产力对油页岩有机质的丰度具有至关重要的影响, 在已发现的湖泊中富有机质沉积物的形成常伴随高古湖泊初始生产力(尹秀珍, 2008; 刘招君等, 2010; Meyers et al., 2012), 并认为表层水的生产力比底部水体的保存条件更为重要(汪品先, 1991)。古湖泊初始生产力研究内容一般包括有机质的来源、母质类型及初始湖泊生产力等。
有机质来源包括内源和外源2种, 其中内源有机质主要来源于湖泊自生生物, 如藻类、菌类、浮游生物和水生植物等, 而外源有机质主要来源于高等植物, 如松柏科、杉科和蕨类植物等。判别有机质来源通常应用地球化学方法, 包括有机显微组分、生物标志化合物、有机碳同位素和某些元素含量变化等(图5-a至5-e)。
生物标志化合物中, 不同碳数的正烷烃可以较准确地判断有机质母质类型(Brassell et al., 1978)。一般小于C20的短链正烷烃源自algaes(藻类)和microorganisms(微生物), C21-25的中链正烷烃为Sphagnum(水藓)或aquatic macrophytes(泥炭藓), 而长链的C27-31正烷烃则来自于陆源高等植物的输入(Bechtel et al., 2007)(图6), 值得注意的是, 一些未成熟油页岩中较高含量的长链的C27-31正烷烃也可能来自于微藻类或细菌(Bechtel et al., 2012)。甾烷化合物中的C27-Sterane主要为菌藻类来源, C29-Sterane为陆源高等植物输入为主, 而C28-Sterane则是水生生物与陆源高等植物的混合来源产物(许世红等, 2002; Peters et al., 2005)。
此外, 伽马蜡烷(Gammacerane)、芳基类异戊二烯(Arylisoprenoids)、17α (β )藿烷等可指示细菌来源的有机质(Peters et al., 1993)(图5-c)。
研究发现, 不同亚相油页岩有机质来源有所差别。半深湖— 深湖亚相油页岩中有机质主要以层状藻和结构藻为主(图7-a, 7-b), 含少量的孢子体(图7-c)、镜质体和惰质体等陆源有机质; 与煤伴生的浅湖亚相油页岩有机质以孢子体和角质体等陆源高等植物为主(图7-d, 7-e), 含少量的藻类和菌类(图7-f)。
古湖泊生产力的定量恢复可采用Mü ller和Suess(1979)提出的有机碳法计算古海洋生产力的公式, 该公式的提出主要针对缺少陆地有机质供应的海相环境, 但对于中国中、新生代陆相湖盆, 需剔除陆源有机质供给的影响。此外, 刘惠民等(2003)建立了陆相湖泊生产力与有机碳之间的关系, 推导了陆相湖泊生产力的计算公式。应用上述2种方法分别计算了松辽盆地上白垩统青山口组一段和嫩江组二段以及桦甸盆地始新统桦甸组沉积时期古湖泊生产力(孟庆涛, 2010; 徐进军, 2015), 结果表明含油页岩段古湖泊生产力均明显高于上、下泥岩段, 表明高的湖泊生产力有利于优质油页岩的发育(图5-e)。
2.2.2 古湖泊水体性质
古湖泊水体性质主要包括湖水的盐度、氧化还原性和矿化度等。通常, TOC/S值(Berner and Raiswell, 1983)、生标参数如伽马蜡烷指数GI(伽马蜡烷/C30霍烷)(Peters et al., 2005)、MTTC值[三— 甲基化烷基色瞒(Tri-MTTC)与总甲基化烷基色瞒(MTTCs)的比值](Sinninghe Damsté et al., 1987)和某些元素含量及比值如B、B/Ga、 Sr/Ba、Ca/Mg、Sr/Ca等可以反映水体的盐度变化(邓宏文和钱凯, 1993)(图5-f至5-h)。水体的含氧量往往与湖水分层有密切的关系, 底部水体的缺氧环境往往是分层湖泊的标志之一(刘招君等, 2010)。目前, 对湖泊水体分层的研究大多处于定性阶段, 湖水的密度分层, 可以是热分层, 也可以是盐度分层。一般来说, 盐度分层比较稳定, 在一些含油页岩湖泊中常见, 有利于有机质的保存; 介形虫类等底栖生物化石组合及其壳体中氧同位素、微量元素比值、沉积物色泽、自生矿物、有机碳含量、变价铁元素指标、遗迹化石和生物扰动构造等均可指示底水含氧量与湖水分层, 此外, 芳基类戊异二烯、β -胡萝卜烷等生标参数也常指示强还原条件导致的盐度分层(Peters et al., 2005)(图5-j至5-m)。
古水体盐度和氧化还原性与油页岩的有机质丰度有一定的对应关系。通常, 缺氧还原环境有利于有机质的保存, 并且水体强的还原性质和盐度更有利于优质油页岩形成(图8)。
古气候通过降雨量改变古湖盆的可容纳空间, 影响古湖泊中的生物发育、物源供给和水体性质, 是控制有机质富集变化和油页岩成矿的重要因素(刘招君等, 2012; Meng et al., 2018)。前人研究表明, 温暖湿润的古气候是油页岩中有机质富集的非常有利的条件之一(Meng et al., 2012, 2018)。
古气候的研究方法很多, 黏土矿物及其组合特征(Deconinck et al., 2000)、化学蚀变指数CIA(Nesbitt and Wilson, 1992)、某些元素比值如Sr/Cu、Fe/Mn、Rb/Sr、Th/U(王随继等, 1997; Johann et al., 2006)等可以定性恢复古气候演化(图6-n至6-p)。笔者所在团队分别利用上述方法研究了准噶尔盆地中二叠统、柴达木盆地中侏罗统、老黑山盆地下白垩统、松辽盆地上白垩统及桦甸盆地始新统油页岩沉积时期古气候变化及其对油页岩成矿的控制作用。
与此同时, 通过植物化石作为气候代用指标定量重建古气候模型受到国内外学者的重视, 近年来, 在一些新生代含油页岩盆地古气候研究中得到广泛应用(史冀忠等, 2008; 孟庆涛等, 2016; Meng et al., 2012, 2018; Quan et al., 2012)。利用该方法, 可以定量计算油页岩沉积时期的古温度、降雨量等参数, 对油页岩的分布及成矿富集规律具有较好的指示作用(孟庆涛, 2010; Meng et al., 2012, 2018)。
1)油页岩识别的关键技术从基于岩石学的地质识别到精确的地球化学指标识别, 综合地球物理资料的二次开发利用从基于岩石物理建模响应的垂向高分辨率测井技术快速识别到地震技术三维空间定量刻画, 历经定性到定量、二维到三维识别的发展过程, 不断发展的测试技术、创新的理论方法以及高品质的地球物理资料不断提升油页岩预测识别的精度。
2)油页岩成矿作用研究的实质是讨论古地理、古湖泊环境和古气候条件等地质因素综合作用下有机质来源与保存条件的变化, 其研究方法从传统的层序地层学、沉积学、古湖泊学逐渐发展到分子地球化学, 由定性向半定量— 定量发展, 精细刻画有机质聚集过程, 有助于较好地预测油页岩展布及富集规律。
文中所介绍的油页岩成矿作用研究中的一些关键方法和技术, 是笔者所在团队近10余年来对中国陆相典型含油页岩盆地研究的成果总结, 但中国油页岩成因类型多样, 除陆相成因油页岩外, 还存在海陆交互相和海相成因油页岩, 其识别特征及成因机制与陆相成因油页岩存在明显差异, 仍有待进一步深入研究与总结。因此, 中国油页岩成矿作用研究尚处于起步阶段, 油页岩地质学家仍需不断努力探索。
作者声明没有竞争性利益冲突.
1 |
|
2 |
|
3 |
|
4 |
|
5 |
|
6 |
|
7 |
|
8 |
|
9 |
|
10 |
|
11 |
|
12 |
|
13 |
|
14 |
|
15 |
|
16 |
|
17 |
|
18 |
|
19 |
|
20 |
|
21 |
|
22 |
|
23 |
|
24 |
|
25 |
|
26 |
|
27 |
|
28 |
|
29 |
|
30 |
|
31 |
|
32 |
|
33 |
|
34 |
|
35 |
|
36 |
|
37 |
|
38 |
|
39 |
|
40 |
|
41 |
|
42 |
|
43 |
|
44 |
|
45 |
|
46 |
|
47 |
|
48 |
|
49 |
|
50 |
|
51 |
|
52 |
|
53 |
|
54 |
|
55 |
|
56 |
|
57 |
|
58 |
|
59 |
|
60 |
|
61 |
|
62 |
|
63 |
|
64 |
|
65 |
|
66 |
|
67 |
|
68 |
|
69 |
|
70 |
|
71 |
|
72 |
|
73 |
|
74 |
|
75 |
|
76 |
|
77 |
|
78 |
|
79 |
|
80 |
|
81 |
|
82 |
|
83 |
|
84 |
|
85 |
|
86 |
|
87 |
|
88 |
|
89 |
|
90 |
|
91 |
|
92 |
|
93 |
|
94 |
|
95 |
|
96 |
|
97 |
|
98 |
|
99 |
|
100 |
|