陆相页岩气储层孔隙发育特征及其主控因素分析:以鄂尔多斯盆地长7段为例
2018-05-26冯小龙敖卫华
冯小龙,敖卫华,唐 玄
1.中国地质大学(北京)工程技术学院,北京 100083 2.中国地质大学(北京)材料科学与工程学院,北京 100083 3.中国地质大学(北京)页岩气勘查与评价国土资源部重点实验室,北京 100083
0 引言
页岩气是主要以游离态和吸附态赋存于富有机质泥页岩中的烃类气体,表现为典型的“自生自储”成藏特点。作为页岩油气资源的储层,泥页岩多具有致密、低孔低渗的特征,因此相比常规油气,页岩油气的开采难度较大[1-3]。目前,页岩气已经成为全球非常规油气勘探开发的热点。截至2016年,美国已经实现了阿巴拉契亚盆地、密歇根盆地、威利斯顿盆地等10余个盆地页岩气的商业开采[4]。我国已在四川盆地下古生界海相和中生界陆相、西北侏罗系陆相以及华北石炭系—二叠系海陆交互相等领域获得页岩气勘探的重要发现[4-8]。无论是美国还是其他等国家,目前商业化开采页岩气较为成功的地区均为海相或海陆交互相地层,而陆相页岩气的勘探开发仍处于试验和研究阶段,大规模商业开发尚未进行。
我国陆相泥页岩分布广泛。据估算,仅华北和东北的陆相页岩气地质资源量已达26.79×1012m3,其中,中—新生界陆相含气页岩是我国重要的目标层段[9]。华北地区陆相富有机质页岩以中生代延长组最为发育,且分布广泛。其中,长7段富有机质页岩是鄂尔多斯盆地中生界最主要的油源岩,同时也是含气量较高的优良储层,吸附含气量为1.97~2.24 m3/g,具有较大的页岩气资源潜力[10-11]。前人针对鄂尔多斯盆地中生界陆相页岩开展了大量研究工作,探讨了页岩的储集性、吸附特征及资源潜力[12-17],但是在气体赋存空间及主控因素方面尚未达成共识,在陆相富有机质页岩的地球化学特征及岩矿组分对储集特征的作用机理等方面有待于进一步深入认识。
本次研究选取该盆地东部Y15井长7段岩心的典型样品,通过总有机碳质量分数、显微组分分析、X射线衍射分析等测试手段,获得页岩地球化学特征和岩石学特征,结合扫描电镜(SEM)分析和低温氮等温吸附测试,获得页岩储层孔隙发育特征,从而通过相关性分析探讨页岩气储层孔隙发育的主控因素。
1 延长组页岩发育地质背景
鄂尔多斯盆地是发育在华北克拉通西部的中生代叠合型盆地。从中生代开始,华北克拉通逐步进入差异升降阶段,受滨太平洋板块向华北板块俯冲作用的影响,盆地东部在走滑-挤压作用下逐步抬升,形成了西倾的斜坡构造格局[18-19]。
上三叠统延长组在整个华北地区分布稳定。长7沉积时期,华北地区继承了晚三叠世海侵作用,随着地壳的持续稳定下沉,海侵规模逐渐增大,形成了一套大面积稳定分布的暗色泥页岩[20]。大量发育的疑源藻、蓝藻及葡萄藻微化石,代表了一组长期稳定温暖潮湿气候条件下的深水湖泊沉积[21]。
鄂尔多斯盆地长7段主要发育3种岩相:黑色页岩岩相、纹层状泥岩岩相和粉砂质泥岩岩相[22]。研究区长7段下部以纹层状泥岩和粉砂质泥岩岩相为主,上部为暗色页岩岩相,有机质类型好、丰度高。
2 样品及实验
实验样品均采自鄂尔多斯盆地X探区Y15井上三叠统延长组长7段(图1),深度为1 740~1 872 m,见表1。全岩及黏土矿物X射线衍射分析采用D8 Advance型X射线衍射仪,并分别依据SY/T5983-1994和SY/T5163-1999标准进行测试;w(TOC)采用Leco碳硫测定仪,镜质体反射率(Ro)采用蔡司Mpv-sp显微光度计测定;孔隙结构分析在TESCAN VEGA Ⅱ型场发射扫描电子显微镜(FESEM)下完成(分辨率为3.0 nm),样品实验前经过氩离子剖光和离子溅射喷镀导电膜处理;氮气吸附实验采用Quadrasorb SI型比表面分析仪在-196 ℃下完成,主要用来获取3.5~100.0nm范围内的孔径分布、比表面积及孔隙体积等参数。实验前样品在110 ℃的真空条件下经过20 h脱气处理,样品的比表面积和孔径分布分别由Brunauer-Emmett-Teller(BET)比表面积测定法和Barret-Joyner-Halenda(BJH)测定法进行计算。
图1 采样点位置示意图Fig.1 Location map of sampling sites
表1 Y15井岩心样品测试数据
BET方法基于郎格缪尔单分子层吸附理论,由等温吸附实验确定的气体分子单层吸附量对多孔介质的表面积进行计算[23-26]:
SBET=(υmNS)/cm。
(1)
式中:SBET为介质的表面积;υm为气体分子单层吸附量;N为阿伏伽德罗常数;S和c分别为氮气的横截面积和物质的量浓度;m为氮气的质量。
BJH方法基于Kelvin-Cohan公式,在压降过程中,孔隙中气体凝聚液的蒸发及吸附层的减薄过程可以表述为[27]
(2)
式中:ΔVpj(n)为第j层气体凝聚液的体积变化量;n为凝聚液的吸附层数;rkj和rpj分别为第j层气体的核部孔径和总孔径;Δtn为吸附层的减薄量;ΔV(n)为孔隙体积的增量;Scj为因凝聚液蒸发所产生的孔隙内表面积。
3 页岩储集空间发育特征
3.1 孔隙类型
页岩具有复杂的孔隙结构,孔隙的大小、形状及连通状况是油气聚集与运移的重要控制因素。目前,页岩气孔隙分类多采用国际理论与应用化学联合会(IUPAC)的分类方案[28],即:根据孔径大小划分为大孔(孔径大于50 nm)、中孔(孔径为2~50 nm)、和微孔(孔径小于2 nm)。孔径分布上,研究区包含大孔、中孔和微孔3种类型的孔隙;按成因类型,研究区页岩基质孔隙可划分为粒间孔、粒内孔、有机质孔以及微裂隙等类型[29-31]。研究采用文献[29]中页岩的孔隙分类标准,即:根据赋存状况将页岩孔隙分为矿物基质孔隙(又可细分为粒间孔、粒内孔及有机质孔)和微裂隙,并由孔隙形态做进一步细分。
有机质孔:富有机质页岩中发育大量有机质孔隙,孔隙发育特征与干酪根之间具有密切的联系。有机质孔的形成伴随干酪根的热降解和生烃过程,随着有机质成熟度的升高,有机质孔隙含量增大[32]。此外,有机质孔也包括来自于生物活动的排泄物和残留下来的生物遗体发育的微孔及纳米孔隙[33-36]。研究区页岩中有机质孔多发育于有机组分与黏土矿物的分界面(图2a),少量发育于有机质内部(图2b),孔隙形态为狭缝状或不规则的网状(图2c),孔径为0.01~0.50 μm,其形成可能与有机质生烃的体积收缩作用有关[33]。与我国南方海相页岩相比,研究区页岩样品热演化程度较低,与有机质生烃作用伴生的圆形-椭圆形残余孔隙发育较差。Jarvie[33]指出,Ro高于0.6%时才可能发育有机质孔隙。整体上,研究区页岩Ro处于0.5%左右,因此认为有机质孔隙不是研究区页岩中的主要孔隙类型。
粒内孔:页岩样品中主要发育黏土集合体矿片间孔隙和颗粒内溶蚀孔隙两种类型的粒内孔。页岩中的黏土矿物多呈絮状体的形式产出,通过单个的鳞片状结构以边-面接触或边-边接触的方式搭建成“卡片屋”的结构。Slatt等[29]指出,这种开放式的孔隙网络可能成为甲烷的渗流通道。研究区页岩中伊/蒙混层絮状体发育程度较好,混层矿物多以裙状、蜂窝状堆积,形成多角形或拉长形孔隙(图2d、e、f),孔径为0.03~1.20 μm。伊利石多发育平行板状孔隙或楔形孔隙,孔径为0.03~0.50 μm,常连通有机质孔及矿物粒间残余孔,形成连通性较高的孔隙网络。
泥页岩中的长石及碳酸盐矿物在有机质脱羧产生的酸性流体的作用下,颗粒内部及边缘常发生溶蚀形成粒内溶孔[34]。溶孔多呈次圆状-不规则状,孔径为0.6~5.0 μm,孔壁呈锯齿状。溶孔多与粒间孔或黏土矿物粒内孔相连通,形成典型的墨水瓶状孔隙。
粒间孔:原生粒间孔多被黏土矿物充填,因此粒间孔隙以残余粒间孔及黏土矿片间孔隙为主。黏土矿物的粒间孔隙类似于砂岩中的粒间孔,但在尺寸上远小于后者(黏土矿物粒间孔通常为纳米尺度上的孔隙),主要以拉长状及不规则形态产出,孔径为1.6~3.0 μm,常与黏土矿物粒内孔相连通(图2f)。
残余粒间孔多为发生充填作用的粒间孔后期被溶蚀而成,多发育于不同组分的分界面(脆性颗粒与黏土矿片之间、黏矿片与有机质之间),孔隙形态多样,多呈三角形、多边形或长条状,孔径为2.5~6.0 μm,连通性较差(图2g、h、r)。
微裂缝:微裂缝是泥页岩中重要的储集空间,同时也是页岩气的渗流通道,对烃类的生产意义重大[36-37]。微裂隙的发育受矿物组分、有机质含量、构造运动及沉积-成岩作用等多种因素控制[38]。研究区样品中微裂隙主要发育于脆性矿物内部及颗粒边缘,产状稳定(图2k、l)。裂缝面呈锯齿状,宽度为0.01~0.50 μm,长为4~15 μm。大部分裂隙均呈开启状态,大规模的溶蚀缝、脆性颗粒边缘缝和黏土矿物及有机质收缩缝常与其他孔隙相连通,构成复杂的孔隙-裂缝网络体系。
a.条带状有机质,有机质孔不发育,微裂缝;b.团块状有机质生烃裂缝;c.块状有机质收缩缝;d.伊利石粒内孔;e.黏土矿物粒内孔;f.黏土矿物粒内孔,残余粒间孔;g.残余粒间孔,溶蚀扩大;h.粒间孔;r.粒间溶孔,微裂缝;j.碳酸盐矿物粒间孔;k.溶蚀缝,粒内溶孔;l.微裂缝,粒内溶孔。除h为SEM外,其他均为FESEM。图2 鄂尔多斯盆地中生界长7段页岩孔隙类型Fig.2 Pore characteristics of Mesozoic Chang 7 Member shale in the Ordos basin
3.2 孔隙结构特征
3.2.1 氮气吸附/脱附曲线
整体上,吸附曲线形态与BET分类中的Ⅱ型吸附等温线接近,具有反“S”形的三段式的特点(图3)。吸附量在低压段(相对压力<0.3)上升较快,吸附曲线向上微凸(图3中ab段),表明微孔对氮气较强的吸附作用;中间段(相对压力为0.3~0.8)表现为多层吸附的特征,吸附量增加缓慢(图3中bc段),与低压段存在较明显的拐点B(图3右);吸附曲线后半段(相对压力>0.8)急剧上升,直至接近饱和蒸气压也未达到吸附饱和(图3中cd段),表明样品中存在一定量的中孔和大孔。
图中阴影部分为滞回环面积;G为解析量骤降对应的拐点;吸附曲线的ab段对应N2的微孔充填阶段,bc段对应微孔及中孔中的多层吸附阶段,cd段对应较大孔隙及微裂隙内的毛细凝聚阶段。图3 页岩样品氮气吸附/脱附曲线Fig.3 Nitrogen adsorption/desorption isotherms of shale samples
通常在开放的孔隙系统中,在发生毛管凝聚和蒸发时气液两相的界面不同,导致气体吸附和脱附时需要的相对压力不同,形成滞回环。滞回环形态及面积(图3中阴影部分)可以有效反映多孔介质的孔隙结构及其对吸附气体的滞留效应;通常孔隙结构越复杂、孔隙连通性越差,孔隙的滞留效应越显著,所形成的滞回环面积越大[39-40]。根据滞回环的形态特征,国际理论与应用化学联合会(IUPAC)将多孔介质分为4种类型(图4)[41]。但实际样品中往往包括不同的孔隙形态,吸附/脱附曲线受多种类型孔隙叠加的影响,实际曲线与理论曲线形态存在较大差异。根据吸附回线的形态可将样品分为两类:
Ⅰ型吸附回线(以Y15-3为代表)与IUPAC的H3型类似,兼具H2型曲线的特征,孔隙结构多为具有狭缝结构或楔形结构的平行板状孔或微裂缝,同时可能发育墨水瓶状孔。高压段出现吸附回线,吸附和脱附分支斜率均较大,表明较大孔径范围内存在开放型孔隙;随着压降曲线趋于平缓并在低压段近于重合,指示了较小孔径范围内的孔隙多为一端封闭型的不透气孔。脱附曲线存在明显的拐点G(图3左),对应相对压力在0.49左右,由Kelvin方程算得的孔隙半径r= 1.34 nm,说明孔径在1.34 nm以上的孔隙中存在开放型孔隙。
Ⅱ型吸附回线(以Y15-7为代表)与IUPAC的H1型类似,孔隙结构以两端开口的圆筒状孔为主。与第一类曲线的区别在于低压段和中等压力段曲线平缓且吸附分支和脱附分支几乎重合,当相对压力接近饱和蒸汽压时曲线斜率急剧增大并发生吸附滞后现象,属于典型的开放型圆筒状孔内的凝聚和蒸发现象。
研究区样品的滞回环面积受w(TOC)及矿物含量影响较大。具体表现为w(TOC)越高、黏土矿物质量分数越高的样品吸附滞后现象越明显,形成的滞回环面积越大;说明有机质内的墨水瓶状孔及黏土矿片间的平行板状孔对页岩孔隙结构起到重要的控制作用。结合扫描电镜对孔隙结构的分析结果来看,黏土矿物矿片间孔隙、有机质内狭缝状孔隙及基质微裂隙等提供了一端封闭或四面都开放的平行板状孔或微裂缝,而有机质孔、残余粒间孔是墨水瓶状孔隙和两端都开放的圆筒状孔隙的主要来源。
3.2.2 孔隙体积、比表面积及孔径分布
据文献[41]。图4 滞后环分类及孔隙类型Fig.4 Classification of hysteresis loops and corresponding pore structure
BJH总孔体积为0.005 4~0.017 2 cm3/g,平均0.011 8 cm3/g(表2)。其中,中孔对孔隙体积的贡献最大,可占总孔体积的76.69%;微孔体积和大孔体积分别占10.12%和13.19%(图5a)。BET比表面积为1.13~17.00 m3/g,平均8.08 m3/g(表2)。其中,中孔和微孔提供了样品中几乎全部的比表面积,分别占67.84%和30.27%,大孔对比表面积的贡献极小(图5b)。与我国南方海相高成熟页岩相比,微孔体积及比表面积所占比例均有较大程度的降低(分别由20.00%和57.60%降低至10.12%和30.27%)[42]。原因可能为研究区陆相页岩热演化程度较低,与有机质生烃作用相关的有机质纳米孔隙发育较差。
孔径分布上,研究区泥页岩样品存在2种孔隙分布类型,一种为单峰型,另一种为双峰-多峰型,两种孔径分布均以微孔及中孔为主(图6)。孔径分布结果显示Y15-1、Y15-2和Y15-4等样品孔径分布曲线在2.6~3.1 nm处存在峰值,并在孔径大于4.2 nm之后吸附量快速下降,孔隙主要集中于2.6~4.2 nm;而Y15-5和Y15-10样品的孔径峰型较为分散,峰值孔径主要分布在2.3~3.1、3.5~3.8以及4.3~5.2 nm。
结合孔径分布与吸附滞后环的形态分析发现,样品Y15-1、Y15-2及Y15-4表现为Ⅱ型吸附回线的形态,孔径分布于2.6~4.2 nm,表明2.6~4.2 nm孔径段主要为开放型的圆筒状孔;而Y15-5和Y15-10多表现为Ⅰ型吸附回线的形态,说明样品中发育2.3~3.1、3.5~3.8、4.3~5.2 nm多个孔径段的平行板状孔和墨水瓶状孔。
表2氮气吸附法孔隙结构参数
Table2Parametersofporestructurederivedfromnitrogenadsorptiondata
样品编号w(TOC)/%比表面积/(m2/g)总孔体积/(cm3/g)平均孔径/nmY15-11.239.430.009617.28Y15-23.6617.000.013913.78Y15-35.394.830.013413.07Y15-42.659.440.014215.86Y15-52.0313.000.01259.24Y15-62.937.460.005412.14Y15-75.241.130.012020.28Y15-81.242.260.007516.80Y15-96.224.830.012113.07Y15-103.6111.400.01729.01
a.总孔体积分布;b.比表面积分布。图5 鄂尔多斯盆地中生界长7段页岩孔隙结构特征Fig.5 Pore structure of Mesozoic Chang 7 Member shale in the Ordos basin
图6 鄂尔多斯盆地中生界长7段页岩孔径分布曲线Fig.6 Pore size distribution of Mesozoic Chang 7 Member shale in the Ordos basin
4 孔隙发育的主控因素
页岩是一种力学强度较高的多孔介质,其孔隙系统受有机组分和无机组分共同制约[43]。Ambrose等[44]通过对场离子显微镜(FIM)与扫描电镜图像重组得到的三维页岩图像的研究发现,页岩基质中大部分孔隙都与干酪根网络相关,因此页岩的孔隙结构与有机组分含量密切相关。通过对页岩样品的岩石热解分析及X射线衍射分析发现,页岩样品的孔隙结构主要受有机质含量及矿物组分的控制。
4.1 有机质类型及成熟度
泥页岩中有机质类型对孔隙结构具有重要的影响。前人对我国南方两套海相页岩的研究[45]表明:Ⅱ型干酪根结构复杂,热演化过程中易于生成有机质孔;相比之下,Ⅰ型干酪根有机质颗粒较小,生成的气孔较少且单个气孔的体积较小。Loucks等[46]指出:Ⅰ型和Ⅱ型干酪根以生油为主,液态烃类排出后有机质中易于形成较大的残余孔隙;而Ⅲ型干酪根以生气为主,排烃后残余有机质中形成的孔隙有限。研究区长7段沉积时期为半深湖—深湖发育区,干酪根类型以Ⅰ型和Ⅱ1型为主[47],具备大规模有机质孔隙发育的有利条件。
页岩孔隙发育特征与干酪根热演化密切相关[48-50]。模拟实验结果表明,泥页岩孔喉演化程度与成熟度正相关,随着有机质大量裂解生烃,有机质孔隙经历了内部微裂隙—蜂窝状集合体—相互连通的较大孔隙的演化过程(孔隙度由最初的0.6%增至2.0%左右)[51]。研究区泥页岩样品Ro平均为0.71%,整体上处于低成熟—成熟阶段,泥页岩孔隙演化经历了压实作用、干酪根初次裂解作用及固体沥青-油质充填作用阶段,有机质生烃过程中产生的有机质孔提供了页岩气吸附和运移的主要空间[52]。然而研究区页岩热成熟度较低,尚未发生有机质大量裂解生烃,Ro与孔隙结构参数之间相关性不明显。
4.2 有机质丰度
有机碳含量是页岩孔隙度和含气饱和度的重要影响因素,也是页岩气藏评价的重要方面。研究表明,页岩中吸附气和游离气含量均随w(TOC)的增加而增大[53]。Loucks等[54]通过对离子剖光页岩样品的扫描电镜分析首次发现并证实了有机质内部的孔隙系统。Passey等[35]和Sondergeld等[55]指出,富有机质页岩中有机质孔最多可以占据总有机质的体积50%,而此时5%的w(TOC)对岩石孔隙体积的贡献多达20%。
研究区页岩样品中除了Y15-1和Y15-6两个样品可能因显微组分的差异导致w(TOC)与孔隙结构参数的关系复杂之外,其余样品的w(TOC)与总孔体积及比表面积均呈正相关的关系(图7a),表明w(TOC)是页岩储集性能的主控因素之一。孔径分布来看,微孔和中孔的孔容均随w(TOC)的增加而增大。前人研究发现,直径小于10 nm的有机质孔隙为页岩气的主要吸附介质,多分布于与矿物伴生的有机质或分散有机质内[56],因此,w(TOC)对页岩储集特征的控制机理体现在与干酪根网络相关的有机质孔隙对微孔及中孔的控制作用上。同时,大孔体积却表现出相反的趋势(图7b),即随w(TOC)的增加而减小。原因可能为矿物基质,尤其是脆性矿物,是页岩中大孔的主要载体,w(TOC)的增加导致脆性矿物和黏土矿物的相对含量减少,因此大孔含量相对降低。
4.3 矿物组分
页岩的矿物组分受沉积环境及后期成岩作用的控制,矿物组分不仅对页岩孔隙结构起到一定的控制作用,而且影响着页岩气开采的增产措施。Hickey等[57]根据岩心矿物成分将Barnett页岩分为有机质页岩、含化石页岩、白云质菱形页岩、白云质页岩、碳酸盐结合体页岩和磷灰岩6种类型。Jarvie等[49]指出,生物成因的硅质经重结晶作用可以形成点接触的岩石框架,而成岩过程中蒙脱石向伊利石转化过程中产生的硅质可胶结陆源石英颗粒[58],这两种情况均有利于原生孔隙的保存和页岩气开采过程中增产方案的实施。
黏土矿物是页岩中重要的矿物组分,对页岩孔隙体积及比表面积的贡献较大。随着黏土矿物种类及含量的变化,矿物晶层及颗粒间空隙的大小、形态差异较大,进而对页岩的孔隙结构及吸附能力产生重要的影响。黏土矿物的孔隙类型多样,包括黏土晶层间的狭缝状和面状孔隙、黏土颗粒之间的粒间孔和溶孔等,其含量与孔隙体积及比表面积有较好的正相关性,并且不同级别孔隙的孔容随着黏土矿物含量的增加呈增大的趋势[59-61]。研究区页岩黏土矿物总质量分数为30.0%~68.5%,平均为54.5%。黏土矿物与孔隙结构参数相关性不明显,但总体上总孔体积随着黏土矿物质量分数的增加而降低,而比表面积具有相反的变化趋势(图8)。原因可能为黏土矿物对原生孔隙的充填和分割作用使孔隙结构复杂化、孔隙度降低,但以其自身较大的表面积,黏土矿物会增加泥页岩的比表面积。
通过黏土矿物质量分数与孔隙分布的研究发现,黏土矿物总质量分数与微孔孔容正相关,与中孔孔容负相关;表明研究区泥页岩中黏土矿物基质孔隙主要集中在微孔级别。吉利明等[62]对单矿物甲烷等的温吸附实验表明,蒙脱石和伊/蒙混层的孔径多集中于微孔范围内,而伊利石和绿泥石则多分布于中孔—大孔的范围。研究区泥页岩中黏土矿物以伊/蒙混层为主(质量分数为25%~65%,平均48%),因此黏土矿物总质量分数决定了微孔含量。此外,伊利石质量分数与比表面积及孔容有较好的相关性,比表面积和孔容均随着伊利石质量分数的增加而增大(图9)。前人研究表明,蒙皂石向伊利石转化过程中由体积收缩而产生微裂(孔)隙,从而增加气体的吸附及运移空间[63]。
图7 w(TOC)与孔隙结构参数之间的相关性Fig.7 Correlation between w(TOC) and pore structure parameters
图8 黏土矿物质量分数与孔隙结构参数之间的相关性Fig.8 Correlation between clay content and pore structure parameters
图9 伊利石质量分数与总孔体积及比表面积相关性Fig.9 Correlation between illite and total pore volume
图10 脆性矿物质量分数与孔隙结构参数之间的相关性Fig.10 Correlation between brittle mineral content and pore structure parameters
除了黏土矿物之外,石英、长石及碳酸盐矿物等脆性矿物也是页岩中的重要组分。脆性矿物的种类及含量不仅影响泥页岩的孔隙发育,而且对页岩气开采过程中储层的压裂改造有重要的影响[64-65]。通过对石英和长石等脆性矿物质量分数与孔容及比表面积等孔隙结构参数之间的关系研究发现,总孔体积随脆性矿物质量分数的增加而降低(图10a),而比表面积与脆性矿物质量分数之间的相关性不明显(图10b)。前人研究[63-65]表明,同沉积的石英等脆性矿物能够抑制压实过程中孔隙度的降低,有利于原生孔隙的保存。此外,石英等脆性矿物在应力作用下颗粒边缘易产生裂缝,改善页岩的孔渗条件[66-67]。但以上两种作用仅对大孔有效。比表面积及孔隙体积分布特征(图5)表明,微孔及中孔是孔隙体积及比表面积的主要贡献者。随着脆性矿物含量的增加,以有机质及黏土矿物为主要载体的微孔及中孔的含量相对降低,导致总孔体积及比表面积均减小。
5 结论
1)长7段页岩的吸附回线可分为2种类型,分别对应黏土矿物矿片间孔隙、有机质内狭缝状孔隙、基质微裂隙以及有机质孔、残余粒间孔等孔隙类型,孔隙结构为分布于2.3~3.1、3.5~3.8、4.3~5.2 nm等多个孔径段的平行板状孔或微裂缝和墨水瓶状孔以及集中分布于2.6~4.2 nm的圆筒状孔。
2)在总孔体积和比表面方面,长7段页岩较我国南方海相富有机质页岩低,分别平均为0.011 8 cm3/g和8.08 m3/g。中孔对总孔体积及比表面积的贡献最大,分别占76.69%及67.84%,其次为微孔,大孔最小。孔径分布变化较大,主要集中在1.6~4.2 nm,体现陆相页岩储层具有较强的非均质性。
3)w(TOC)是延长组页岩总孔体积及比表面积的主要控制因素,作用机制主要体现在对微孔及中孔含量的控制上;Ro和矿物组成对孔隙发育的控制作用不明显,但黏土矿物对微孔孔容贡献较大,其中伊利石质量分数与比表面积及孔容有较好的相关性。
(
):
[1] 付金华, 郭少斌, 刘新社, 等. 鄂尔多斯盆地上古生界山西组页岩气成藏条件及勘探潜力[J]. 吉林大学学报(地球科学版), 2013, 43(2):382-389.
Fu Jinhua, Guo Shaobin, Liu Xinshe, et al. Shale Gas Accumulation Condition and Exploration Potential of the Upper Paleozoic Shanxi Formation in Ordos Basin[J]. Jounal of Jilin University (Earth Science Edition), 2013, 43(2): 382-389.
[2] 卢双舫, 黄文彪, 陈方文,等. 页岩油气资源分级评价标准探讨[J]. 石油勘探与开发, 2012, 39(2): 249-256.
Lu Shuangfang, Huang Wenbiao, Chen Fangwen, et al. Classification and Evaluation Criteria of Shale Oil and Gas Resources: Discussion and Application[J]. Petroleum Exploration and Development, 2012, 39(2): 249-256.
[3] 邹才能, 杨智, 崔景伟,等. 页岩油形成机制、地质特征及发展对策[J]. 石油勘探与开发, 2013, 40(1):14-26.
Zou Caineng, Yang Zhi, Cui Jingwei, et al. Formation Mechanism, Geological Characteristics, and Development Strategy of Nonmarine Shale Oil in China[J]. Petroleum Exploration and Development, 2013, 40(1): 14-26.
[4] 王志刚. 涪陵页岩气勘探开发重大突破与启示[J]. 石油与天然气地质, 2015, 36(1): 1-6.
Wang Zhigang. Breakthrough of Fuling Shale Gas Exploration and Development and Its Inspiration[J]. Oil & Gas Geology, 2015, 36(1): 1-6.
[5] 郭彤楼. 涪陵页岩气田发现的启示与思考[J]. 地学前缘, 2016,23(1):29-43.
Guo Tonglou. Discovery and Characteristics of the Fuling Shale Gas Field and Its Enlightenment and Thinking[J]. Earth ScienceFrontiers,2016,23(1): 29-43
[6] 韩辉, 钟宁宁, 陈聪,等. 西北地区中小型盆地侏罗系陆相泥页岩的含气性[J]. 科学通报, 2014,59(9):809-815.
Han Hui, Zhong Ningning, Chen Cong, et al. The Gas Potential of Jurassic Continental Shale in the Middle-Small Basins, Northwest China[J]. Chinese Science Bulletin, 2014, 59(9): 809-815.
[7] 董大忠, 邹才能, 杨桦,等. 中国页岩气勘探开发进展与发展前景[J]. 石油学报, 2012, 33(增刊1):107-114.
Dong Dazhong, Zou Caineng, Yang Hua, et al. Progress and Prospects of Shale Gas Exploration and Development in China[J]. Acta Petrolei Sinica, 2012,33(Sup. 1): 107-114.
[8] 董大忠, 高世葵, 黄金亮,等. 论四川盆地页岩气资源勘探开发前景[J]. 天然气工业, 2014, 34(12): 1-15.
Dong Dazhong, Gao Shikui, Huang Jinliang, et al. A Discussion on the Shale Gas Exploration & Development Prospect in the Sichuan Basin[J]. Natural Gas Industry, 2014, 34(12): 1-15.
[9] 国土资源部油气资源战略研究中心. 全国页岩气资源潜力调查评价及有利区优选(2009—2011年):2014年度国土资源科学技术奖一等奖成果[J]. 科技成果管理与研究, 2014(12): 57-60.
Strategic Research Center for Oil and Gas Resources, Ministry of Land and Resources. Evaluation of Shale Gas Resource Potential and Optimization of Favorable Areas in China (2009-2011): First Prize of Science and Technology of Land and Resources in 2014[J]. Management and Research on Scientific & Technological Achievements, 2014(12): 57-60.
[10] 杨华, 邓秀芹. 构造事件对鄂尔多斯盆地延长组深水砂岩沉积的影响[J]. 石油勘探与开发, 2013, 40(5):513-520.
Yang Hua, Deng Xiuqin. Deposition of Yanchang Formation Deep-Water Sandstone Under the Control of Tectonic Events, Ordos Basin[J]. Petroleum Exploration and Development, 2013, 40(5): 513-520.
[11]Tang X, Zhang J, Jin Z, et al. Experimental Investigation of Thermal Maturation on Shale Reservoir Properties from Hydrous Pyrolysis of Chang 7 Shale, Ordos Basin[J]. Marine & Petroleum Geology, 2015, 64(Sup.1): 165-172.
[12] 杨超, 张金川, 唐玄. 鄂尔多斯盆地陆相页岩微观孔隙类型及对页岩气储渗的影响[J]. 地学前缘, 2013, 20(4): 240-250.
Yang Chao, Zhang Jinchuan, Tang Xuan. Microscopic Pore Types and Its Influence on Storage and Permeability of Continental Shale Gas, Ordos Basin[J]. Earth Science Frontiers, 2013, 20(4): 240-250.
[13] 赵靖舟, 王芮, 耳闯. 鄂尔多斯盆地延长组长7段暗色泥页岩吸附特征及其影响因素[J]. 地学前缘, 2016, 23(1): 146-153.
Zhao Jingzhou, Wang Rui, Er Chuang. Adsorption Characteristics of Chang 7 Shale from the Triassic Yanchang Formation in Ordos Basin, and Its Controlling Factor[J]. Earth Science Frontiers, 2016, 23(1):146-153.
[14] 唐玄, 张金川, 丁文龙,等. 鄂尔多斯盆地东南部上古生界海陆过渡相页岩储集性与含气性[J]. 地学前缘, 2016, 23(2): 147-157.
Tang Xuan, Zhang Jinchuan, Ding Wenlong, et al. The Reservoir Property of the Upper Paleozoic Marine-Continental Transitional Shale and Its Gas Bearing Capacity in the Southeastern Ordos Basin[J]. Earth Science Frontiers, 2016, 23(2): 147-157.
[15] 曾秋楠, 于炳松, 李昱霏. 鄂尔多斯盆地东南部延长组页岩储层特征及控制因素[J]. 特种油气藏, 2013, 20(1): 23-26.
Zeng Qiunan, Yu Bingsong, Li Yufei. Reservoir Characteristics and Their Controlling Factors of the Yanchang Formation Shale in Southeast of Ordos Basin[J]. Special Oil & Gas Reservoirs, 2013, 20(1): 23-26.
[16] 吴建国, 刘大锰, 姚艳斌. 鄂尔多斯盆地渭北地区页岩纳米孔隙发育特征及其控制因素[J]. 石油与天然气地质, 2014, 35(4): 542-550.
Wu Jianguo, Liu Dameng, Yao Yanbin. Characteristics and Controlling Factors of Nanopores in Shales in Weibei, Ordos Basin[J]. Oil and Gas Geology, 2014, 35 (4): 542-550.
[17] 曾秋楠, 周新桂, 于炳松,等. 陆相页岩气储层评价标准探讨:以延长组富有机质页岩为例[J]. 新疆地质, 2015,33(3):409-414.
Zeng Qiunan, Zhou Xingui, Yu Bingsong, et al. Evaluation Criteria of Lake Facies Shale Gas Reservoir: A Case Study of the Organic Rich Shale Developed in Yanchang Group, Ordos Basin[J]. Xinjiang Geology, 2015, 33(3): 409-414.
[18] 邓军, 王庆飞, 黄定华,等. 鄂尔多斯盆地基底演化及其对盖层控制作用[J]. 地学前缘, 2005, 12(3):91-99.
Deng Jun, Wang Qingfei, Huang Dinghua, et al. Basement Evolution of the Ordos Basin and Its Constraint on Cap Rock[J]. Earth Science Frontiers,2005, 12(3): 91-99.
[19] 赵靖舟, 王力, 孙兵华,等. 鄂尔多斯盆地东部构造演化对上古生界大气田形成的控制作用[J]. 天然气地球科学, 2010, 21(6): 875-881.
Zhao Jingzhou, Wang Li, Sun Binghua, et al. Effect of Structural Evolution on the Formation and Distribution of Upper Paleozoic Giant Gas Fields in the East Ordos Basin[J]. Natural Gas Geoscience, 2010, 21(6): 875-881.
[20] 杨华, 田景春, 王峰, 等. 鄂尔多斯盆地三叠纪延长组沉积期湖盆边界与底形及事件沉积研究[M]. 北京: 地质出版社, 2009: 90-144.
Yang Hua, Tian Jingchuan, Wang Feng, et al. Research on Lacustrine Basin Boundary and Event Deposits of Triassic Yanchang Sedimentary Period in Ordos Basin[M]. Beijing: Geology Publishing House, 2009: 90-144.
[21] 吉利明, 祝幼华. 鄂尔多斯盆地西南部甘肃西峰地区延长组孢粉组合及古气候研究[J]. 微体古生物学报, 2013, 30(4):367-378.
Ji Liming, Zhu Youhua. Sporo-Polen Assemblages and Palenoclimate of the Yanchang Formation in the Xifeng Area, Southwestern Ordos Basin, Gansu Province, NW China[J]. Acta Micropalaeontologica Sinica, 2013, 30(4): 367-378.
[22] 耳闯, 罗安湘, 赵靖舟,等. 鄂尔多斯盆地华池地区三叠系延长组长7段富有机质页岩岩相特征[J]. 地学前缘, 2016, 23(2):108-117.
Er Chuang, Luo Anxiang, Zhao Jingzhou, et al. Lithofacies Features of Organic-Rich Shale of the Triassic Yanchang Formation in Huachi Aera, Ordos Basin[J]. Earth Science frontiers, 2016, 23(2): 108-117.
[23] Brunauer S, Emmett P H, Teller E. Adsorption of Gases in Multimolecular Layers[J]. Journal of the American Chemical Society, 1938, 60(2): 309-319.
[24] Kuila U, Prasad M. Specific Surface Area and Pore-Size Distribution in Clays and Shales[J]. Geophysical Prospecting, 2013, 61(2): 341-362.
[25] ISO 9277:2010 Determination of the Specific Surface Area of Solids by Gas Adsorption: BET method[S]. Geneva: Technical Committee ISO/TC 24 and Subcommittee SC 4, 1995.
[26] Rouquerol J, Llewellyn P, Rouquerol F. Is the BET Equation Applicable to Microporous Adsorbents[J]. Studies in Surface Science & Catalysis, 2007, 160(7): 49-56.
[27] Rouquerol F, Rouquerol J, Sing K. Adsorption by Powders and Porous Solids: Principles, Methodology and Applications[M]. Marseille: Academic Press, 2014.
[28] Brunauer S, Deming L S, Deming W E, et al. On a Theory of the Van der Waals Adsorption of Gases[J]. Journal of the American Chemical Society, 1940, 62(7): 1723-1732.
[29] Slatt R M, O'Brien N R. Pore Types in the Barnett and Woodford Gas Shales: Contribution to Understanding Gas Storage and Migration Pathways in Fine-Grained Rocks[J]. Aapg Bulletin, 2011, 95(12): 2017-2030.
[30] Turner J, Turner J. Zooplankton Fecal Pellets, Ma-rine Snow and Sinking Phytoplankton Blooms: Aquat Microbial Ecol[J]. Aquatic Microbial Ecology, 2002, 27(1): 57-102.
[31] 于炳松. 页岩气储层孔隙分类与表征[J]. 地学前缘, 2013, 20(4): 211-220.
Yu Bingsong. Classification and Characterization of Gas Shale Pore System[J]. Earth Science frontiers, 2013, 20(4): 211-220.
[32] 黄磊, 申维. 页岩气储层孔隙发育特征及主控因素分析:以上扬子地区龙马溪组为例[J]. 地学前缘, 2015, 22(1):374-385.
Huang Lei, Shen Wei. Characteristics and Controlling Factors of the Formation of Pores of a Shale Gas Reservoir: A Case Study from Longmaxi Formation of the Upper Yangtze Region[J]. Earth Science frontiers, 2015, 22(1): 374-385.
[33] Jarvie D. Evaluation of Hydrocarbon Generation and Storage in Barnett Shale, for Worth Basin, Texas[R]. Texas: Humble Geochemical Services Division, 2004.
[34] 魏祥峰, 刘若冰, 张廷山,等. 页岩气储层微观孔隙结构特征及发育控制因素:以川南—黔北XX地区龙马溪组为例[J]. 天然气地球科学, 2013, 24(5):1048-1059.
Wei Xiangfeng, Liu Ruobing, Zhang Tingshan, et al. Micro-Pores Structure Characteristics and Development Control Factors of Shale Gas Reservoir: A Case of Longmaxi Formation in XX Area of Southern Sichuan and Northern Guizhou[J]. Natural Gas Geoscience, 2013, 24(5): 1048-1059.
[35] Passey Q R, Bohacs K, Esch W L, et al. From Oil-Prone Source Rock to Gas-Producing Shale Reservoir: Geologic and Petrophysical Characterization of Unconventional Shale Gas Reservoirs[C]// CPS/SPE International Oil &Gas Conference and Exhibition. Beijing: [s. n.], 2010: 131350.
[36] 黄振凯, 陈建平, 薛海涛,等. 松辽盆地白垩系青山口组泥页岩孔隙结构特征[J]. 石油勘探与开发, 2013, 40(1): 58-65.
Huang Zhenkai, Chen Jianping, Xue Haitao, et al. Microstructural Characteristics of the Cretaceous Qingshankou Formation Shale, Songliao Basin[J]. Petroleum Exploration and Development, 2013, 40(1): 58-65.
[37] Gale J F W, Reed R M, Holder J. Natural Fractures in the Barnett Shale and Their Importance for Hydraulic Fracture Treatments[J]. Aapg Bulletin, 2007, 91(4):603-622.
[38] 何龙, 郑荣才, 梁西文,等. 川东涪陵地区大安寨段裂缝控制因素及期次分析[J]. 岩性油气藏, 2014, 26(4):88-96.
He Long, Zheng Rongcai, Liang Xiwen, et al. Controlling Factors and Development Periods of Fracture of Da’anzhai Member in Fuling Area, Eastern Sichuan Basin[J]. Lithologic Reservoirs, 2014, 26(4): 88-96.
[39] 刘伟新, 俞凌杰, 张文涛,等. 川东南龙马溪组页岩微观孔隙结构特征[J]. 海洋地质与第四纪地质, 2016,36(3):127-134.
Liu Weixin, Yu Lingjie, Zhang Wentao, et al. Micro-Pore Structure of Longmaxi Shale from Southeast Sichuan Basin[J]. Marine Geology & Quaternary Geology, 2016, 36(3): 127-134.
[40] 陈萍, 唐修义. 低温氮吸附法与煤中微孔隙特征的研究[J]. 煤炭学报, 2001, 26(5):552-556.
Chen Ping, Tang Xiuyi. The Research on the Adsorption of Nitrogen in Low Temperature and Micro-Pore Properties in Coal[J]. Journal of China Coal Society, 2001, 26(5): 552-556.
[41] 严继民,张启元.吸附与聚集[M].北京:科学出版社,1979:108-120.
Yan Jimin, Zhang Qiyuan. Adsorption and Aggregation[M]. Beijing: Science Press, 1979: 108-120.
[42] 杨潇, 姜振学, 宋岩,等. 渝东南牛蹄塘组与龙马溪组高演化海相页岩全孔径孔隙结构特征对比研究[J]. 高校地质学报, 2016, 22(2):368-377.
Yang Xiao, Jiang Zhenxue, Song Yan, et al. A Comparative Study on Whole-Aperture Pore Structure Characteristics Between Niutitang and Longmaxi Formation of High-Matruity Marine Shales in Southeastern Chongqingelt[J]. Geological Journal of China Universities, 2016, 22(2): 368-377.
[43]Kang S M, Fathi E, Ambrose R J, et al. Carbon Dioxide Storage Capacity of Organic-Rich Shales[J]. Spe Journal, 2013, 16(4): 842-855.
[44] Ambrose R J, Hartman R C, Campos M D, et al. New Pore-Scale Considerations for Shale Gas in Place Calculations[C]// SPE Unconventional Gas Conference. Pittsburgh: Society of Petroleum Engineers, 2010: 131772.
[45] 聂海宽, 边瑞康, 张培先,等. 川东南地区下古生界页岩储层微观类型与特征及其对含气量的影响[J]. 地学前缘, 2014, 21(4):331-343.
Nie Haikuan, Bian Ruikang, Zhang Peixian, et al. Micro-Types and Characteristics of Shale Reservoir of the Lower Paleozoic in Southeast Sichuan Basin, and Their Effects on the Gas Content[J]. Earth Science Frontiers, 2014, 21(4): 331-343.
[46]Loucks R G, Reed R M, Ruppel S C, et al. Mo-rphology, Genesis, and Distribution of Nanometer-Scale Pores in Siliceous Mudstones of the Mississippian Barnett Shale[J]. Journal of Sedimentary Research, 2009, 79(12): 848-861.
[47] 耳闯, 赵靖舟, 王芮,等. 沉积环境对富有机质页岩分布的控制作用:以鄂尔多斯盆地三叠系延长组长7油层组为例[J]. 天然气地球科学, 2015, 26(5):823-832.
Er Chuang, Zhao Jingzhou, Wang Rui, et al. Controlling Role of Sedimentary Environment on the Distribution of Organic-Rich Shale: A Case Study of the Chang7 Member of the Triassic Yanchang Formation, Ordos Basin[J]. Natural Gas Geoscience, 2015, 26(5): 823-832.
[48] Schieber J. Common Themes in the Formation and Preservation of Porosity in Shales and Mudstones-Illustrated with Examples across the Phanerozoic[R]. Pittsburgh: SPE Unconventional Gas Conference, 2010.
[49] Jarvie D M, Jarvie B M, Weldon D, et al. Com-ponents and Processes Impacting Production Success from Unconventional Shale Resource Systems[C]// 10th Middle East Geosciences Conference and Exhibition. Manama: [s. n.], 2012: 40908.
[50] Curtis M E, Ambrose R J, Sondergeld C H, et al. Investigation of the Relationship Between Organic Porosity and Thermal Maturity in the Marcellus Shale[C]// North American Unconventional Gas Conference and Exhibition. Woodlands: Society of Petroleum Engineers, 2011: 144370.
[51] 吴松涛, 朱如凯, 崔京钢,等. 鄂尔多斯盆地长7湖相泥页岩孔隙演化特征[J]. 石油勘探与开发, 2015, 42(2):167-176.
Wu Songtao, Zhu Rukai, Hui Jinggang, et al. Characteristics of Lacustrine Shale Porosity Evolution, Triassic Chang 7 Member, Ordos Basin, NW China[J].Petroleum Exploration and Development, 2015, 42(2): 167-176.
[52] Bustin R M. Shale Gas and Shale Oil Petrology and Petrophysics[J]. International Journal of Coal Geology, 2012, 103: 1-2.
[53] Wang F P, Reed R M. Pore Networks and Fluid Flow in Gas Shales[C]// SPE Annual Technical Conference and Exhibition. New Orleans: Society of Petroleum Engineers, 2009: 124253.
[54] Loucks R G, Reed R M, Ruppel S C, et al. Mor-phology, Genesis, and Distribution of Nanometer-Scale Pores in Siliceous Mudstones of the Mississippian Barnett Shale[J]. Journal of Sedimentary Research, 2009, 79(12): 848-861.
[55] Sondergeld C H, Newsham K E, Comisky J T, et al. Petrophysical Considerations in Evaluating and Producing Shale Gas Resources[C]// SPE Unconventional Gas Conference. Pennsylvania: Society of Petroleum Engineers, 2010: 131768.
[56] 梁峰, 拜文华, 邹才能,等. 渝东北地区巫溪2井页岩气富集模式及勘探意义[J]. 石油勘探与开发, 2016, 43(3):350-358.
Liang Feng, Bai Wenhua, Zou Caineng, et al. Shale Gas Enrichment Pattern and Exploration Significance of Well Wuxi-2 in Northeast Chongqing, NE Sichuan Basin[J]. Petroleum Exploration and Development,2016, 43(3): 350-358.
[57]Hickey J J, Bo H. Lithofacies Summary of the Mississippian Barnett Shale, Mitchell 2 TP Sims Well, Wise County, Texas[J]. Aapg Bulletin, 2007, 91(4): 437-443.
[58] Thyberg B, Jahren J, Winje T, et al. From Mud to Shale: Rock Stiffening by Micro-Quartz Cementaton[J]. First Break, 2009, 22(10): 170-173.
[59] 王茂桢, 柳少波, 任拥军,等. 页岩气储层黏土矿物孔隙特征及其甲烷吸附作用[J]. 地质论评, 2015, 61(1):207-216.
Wang Maozhen, Liu Shaobo, Ren Yongjun, et al. Pore Characteristics and Methane Adsorption of Clay Minerals in Shale Gas Reservoir[J]. Geological Review, 2015, 61(1): 207-216.
[60] 杨峰, 宁正福, 张世栋,等. 基于氮气吸附实验的页岩孔隙结构表征[J]. 天然气工业, 2013, 33(4):135-140.
Yang Feng, Ning Zhengfu, Zhang Shidong, et al. Characterization of Pore Structures in Shales Through Nitrogen Adsorption Experiment[J]. Natural Gas Industry, 2013, 33(4): 135-140.
[61] 黄振凯, 陈建平, 王义军,等. 松辽盆地白垩系青山口组泥岩微观孔隙特征[J]. 石油学报, 2013, 34(1):30-36.
Huang Zhenkai, Chen Jianping, Wang Yijun, et al. Characteristics of Micropores in Mudstones of the Cretaceous Qingshankou Formation, Songliao Basin[J]. Acta Petrolei Sinica, 2013, 34(1): 30-36.
[62] 吉利明, 马向贤, 夏燕青,等. 黏土矿物甲烷吸附性能与微孔隙体积关系[J]. 天然气地球科学, 2014, 25(2): 141-152.
Ji Liming, Ma Xiangxian, Xia Yanqing, et al. Relationship Between Methane Adsorption Capacity of Clay Minerals and Micropore Volume[J]. Natural Gas Geoscience, 2014, 25(2): 141-152.
[63] 蒋裕强, 董大忠, 漆麟,等. 页岩气储层的基本特征及其评价[J]. 天然气工业, 2010, 30(10): 7-12.
Jiang Yuqiang, Dong Dazhong, Qi Lin, et al. Basic Features and Evaluation of Shale Reservoirs[J]. Natural Gas Industry, 2010, 30(10): 7-12.
[64] Ross D J K, Bustin R M. The Importance of Shale Composition and Pore Structure upon Gas Storage Potential of Shale Gas Reservoirs[J]. Marine & Petroleum Geology, 2009, 26(6): 916-927.
[65] Hill D G, Nelson C R. Reservoir Properties of the Upper Cretaceous Lewis Shale, a New Natural Gas Play in the San Juan Basin[J]. AAPG Bulletin, 2000, 84(8): 1240.
[66] 张琴, 刘洪林, 拜文华,等. 渝东南地区龙马溪组页岩含气量及其主控因素分析[J]. 天然气工业, 2013,33(5):35-39.
Zhang Qin, Liu Honglin, Bai Wenhua, et al. Shale Gas Content and Its Main Controlling Factors in Longmaxi Shales in Southeastern Chongqing[J]. Natural Gas Industry,2013, 33(5): 35-39.
[67] 梁超, 姜在兴, 杨镱婷,等. 四川盆地五峰组—龙马溪组页岩岩相及储集空间特征[J]. 石油勘探与开发, 2012, 39(6):691-698.
Liang Chao, Jiang Zaixing, Yang Yiting, et al. Characteristics of Shale Lithofacies and Reservoir Space of the Wufeng-Longmaxi Formation, Sichuan Basin[J]. Petroleum Exploration and Development, 2012, 39(6): 691-698.