玛湖凹陷二叠系风城组页岩油储层微观孔隙结构精细表征_邵广辉.pdf

收稿日期 基金项目中国石油天然气股份有限公司重大科技专项“陆相中高成熟度页岩油勘探开发关键技术研究与应用”（）。第一作者邵广辉（），男，工程师，现主要从事测井储层评价方面的研究工作，。邵广辉，高衍武，蔺敬旗，等 玛湖凹陷二叠系风城组页岩油储层微观孔隙结构精细表征 长江大学学报（自然科学版），（）：，（），（）：玛湖凹陷二叠系风城组页岩油储层微观孔隙结构精细表征邵广辉，高衍武，蔺敬旗，吴伟，杨帆，肖华，王先虎，彭宇，曾蕊 中国石油集团测井有限公司地质研究院，陕西 西安 中国石油集团测井有限公司新疆分公司，新疆 克拉玛依 中国石油集团测井有限公司制造公司，陕西 西安 中国石油集团测井有限公司长庆分公司，陕西 西安 摘要为了研究玛湖凹陷二叠系风城组页岩油储层微观孔隙结构，利用扫描电镜、核磁共振、高压压汞及低温氮气吸附种手段对页岩样品开展实验分析，获取了不同实验方式下的微观孔隙结构特征和定量评价参数，分析了各实验方法在孔隙特征定性和定量评价上的利弊。在此基础上，利用二次拟合方法构建了孔隙结构的全孔径曲线，实现了微观孔隙结构全孔径联合精细表征。结果显示，风城组页岩油储层整体孔隙度低、孔隙类型复杂多样；样品的全孔径分布曲线呈多峰分布特征，孔径分布跨度大，的孔体积主要在孔径 的孔隙中，说明页岩样品的孔隙主要以小孔及微孔为主；在孔径小于 的区间内，孔隙大小与总有机碳（）含量及黏土含量关系密切。上述研究结果为研究区页岩油储层微观孔隙结构评价提供了更为精确的方案。关键词玛湖凹陷；页岩油；微观孔隙结构；核磁共振；低温氮气吸附 中图分类号 文献标志码 文章编号 （），：，；，；，：；长江大学学报（自然科学版）年 第 卷 第期 （），DOI:10.16772/ki.1673-1409.20220921.002近年来，玛湖凹陷二叠系风城组页岩油作为准噶尔盆地重要的非常规接替资源，已成为勘探开发的热点。已有研究成果表明，研究区储层微观孔隙结构复杂，微米、纳米级孔隙均有发育，孔径分布范围广，页岩油储层物性及可动油丰度受孔隙结构影响大。对于页岩油储层孔隙结构评价，前人通常运用低温氮气吸附、核磁共振及高压压汞等方法进行孔隙结构特征分析，利用单一实验或多种实验结果单独分析页岩油储层孔隙结构特征。但受实验原理的影响，每种实验所得出的孔隙结构特征在表征上各有利弊，单一实验方法往往只能体现储层一部分孔隙结构特征，无法全面表征页岩油储层孔隙结构及其分布范围，对于原油赋存状态多样的页岩油储层，各类孔隙都存在原油的分布，只评价小孔或大孔的孔隙结构不足以满足储层吸附油和游离油丰度的评价，对试油及勘探开发造成不确定性。准确表征微观孔隙结构是进行页岩油储层含油性评价的关键。为此，笔者以玛湖凹陷二叠系风城组页岩油储层为例，以实验方法原理为基础，逐一分析实验结果并判断其在页岩孔隙结构表征方面上的优缺点，最终优选实验结果开展多方法孔隙结构联合表征，以期建立页岩油全孔径精细表征方法，为正确认识储层孔隙结构特征及后期页岩油资源的开发提供依据。样品的基本特征为了增加不同测试结果的对比性，将同一深度点上的页岩样品切分成平行样，分别开展矿物成分、地化特征和气测孔隙度测试，结果如表所示。由表可知，页岩样品总有机碳（）含量分布范围为 ，氦气孔隙度分布范围为 。矿物成分以石英与白云石为主，长石含量次之，部分样品发育填充堵塞孔隙的方解石，黏土平均含量为 ，其他矿物成分含量较少。石英对储层孔隙有一定的支撑作用，长石和白云石易发生溶蚀形成次生孔隙，三者整体含量高的样本，其孔隙发育情况要优于其他样品。表页岩样品基本信息 样品编号深度氦气孔隙度 含量不同矿物及其含量石英长石方解石白云石黄铁矿黏土 平均值 实验结果 扫描电镜利用 型扫描电镜仪器对岩样进行抽提清洗，在样品表面上用细聚焦电子束由点到行进行扫描，最终利用探测器对图像进行放大处理，通过不同放大倍数下的观察，从而对微米级、亚微米级和纳米级孔隙进行观察和分析。图像分析结果显示，研究区页岩样品的孔隙分为有机孔和无机孔两大类。有机孔是页岩有机组分在生烃过程中在有机质内部或边缘所产生的孔隙，是页岩油储层重要的孔隙类型。如图所示，通过不同倍数观测到有机孔形态多样，多以近圆形和椭圆形为主，局部发育狭长型和不规则状有机质孔。研究区页岩热演化程度偏低，并不发育蜂窝状密集分布的有机质孔隙，多为孤立的形式分布在有机质中，孔径以纳米级和亚微米级孔隙为主，少数狭长型有机孔可达到微米级。第 卷 第期邵广辉 等：玛湖凹陷二叠系风城组页岩油储层微观孔隙结构精细表征图页岩样品有机孔形态特征 无机孔主要发育粒间孔、粒内孔、黏土层间孔（见图）。在长石、石英等矿物颗粒之间或与有机质的交接边缘处通常发育粒间孔，以近卵形、三角形或不规则状为主，多数孔径可达到数十至百纳米，普遍具有较好的孔隙连通性，是页岩储层主要的渗流通道。粒内孔发育在矿物颗粒的内部，多数以点状分布，尺寸较小。黏土层间孔是黏土矿物之间的孔隙，一般呈狭缝状，定向平行分布，长度可达。由于岩石在生长过程中矿物晶体堆积疏松，造成岩样局部在黄铁矿、黏土等矿物晶体之间发育少量晶间孔。其次，由于成岩过程中的压实、脱水收缩和重结晶作用，在无机矿物基质中发育微裂缝，具有一定的延展性，缝宽在 范围内，局部微裂缝连接形成缝网，这种缝网在页岩油运移的过程中起到了改善渗流环境的突出作用。最后，长石与白云石均有溶蚀作用，对储层起到增孔的作用，其中长石是溶蚀主体，白云石次之。图页岩样品无机孔形态特征 核磁共振根据核磁共振原理，在均匀的静磁场中，岩样饱和流体时的氢质子由于静磁场发生极化而产生磁矢量。同时，氢质子被频率脉冲激发，产生核磁共振现象。在去除频率脉冲后，可以得到幅度随时间衰减长江大学学报（自然科学版）年月的信号。由于储层孔隙中流体的氢质子与孔径大小呈正相关性，实验得出的横向弛豫时间（）谱图能够展现出储层孔隙的分布特征。越短，孔隙体积越小，反之孔隙体积越大，。对研究区页岩样品进行长时间的彻底洗油并烘干，将页岩样品模拟饱和原始地层矿化度的水，利用 核磁共振分析仪最终得到饱和水页岩样品的谱。研究区页岩样品谱可以分为两种典型类型：类和类（见图）。两种类型均为双峰特征，说明样品至少发育两种孔隙类型，谱图的左侧峰为微小孔特征峰，右侧峰为中大孔特征峰。类样品谱峰的位置在 之间，谱峰靠左，说明该类样品微小孔发育较好，中大孔少量发育；类样品谱峰的位置在 之间，说明微小孔隙不发育，黏土及毛细管束缚水含量少，该样品以发育中大孔为主。结合全岩矿物、扫描电镜及地化实验分析，类样品谱整体黏土及有机质含量较高，发育多种孔隙类型，孔隙结构复杂，有机孔占其样品孔隙类型的主导位置；类样品谱则以无机孔为主，主要发育剩余粒间孔及溶蚀孔，中大孔隙发育，整体物性优于类样品，但生油能力较弱。图页岩样品核磁共振典型双峰谱图 高压压汞利用 型压汞仪，设置 为实验最高压力，测量研究区页岩样品的进汞和退汞曲线。利用进汞、退汞曲线，结合 理论，可计算出样品的孔喉大小、孔径分布及连通性等参数。研究区页岩样品整体排驱压力大，平台段斜率大，进汞退汞效率低（见图）。进汞饱和度平均为 ，进汞饱和度与压力呈线性关系，平台段斜率大，分选较差，退汞曲线随着压力下降，退汞饱和度并没有匹配的变化趋势；当压力下降至极低时，样品退汞才出现下降趋势，最终退汞效率较低表明页岩样品的孔隙中仍滞留了大量残留汞。上述进汞和退汞曲线的特征共同表明了研究区页岩样品主要以较小的孔隙网络为主，整体喉道细小，大部分孔隙不连通，孤立孔占比大，这种结构不利于孔隙内部流体的运移，整体属于微孔细喉储层。由于页岩样品整体进汞饱和度过低，压汞曲线反映的并不是真实的孔隙结构，只能用来定性分析判断储层孔隙结构的特征。低温氮气吸附首先选取洗油前后的残余样品，粉碎为 目，碎样分别在 、下抽真空至压力为 ，再分别通过 、环境下烘干 ，将样品中过剩水分和杂质气体进行去除，通过低温氮吸附设备 在临界温度下测量氮气的吸附量和脱附量，最终得到洗油前后的低温氮气吸附、脱附曲线及样品孔径分布等实验数据。第 卷 第期邵广辉 等：玛湖凹陷二叠系风城组页岩油储层微观孔隙结构精细表征图页岩样品压汞曲线特征 通过低温氮气等温吸附脱附曲线形态的变化规律，可以判断出特定区间孔隙的分布特征。吸附等温线在低（为氮气分压；为液氮温度下氮气的饱和蒸气压；为相对压力）区域，等温线呈上凸形态，在拐点处完成单层吸附，随着相对压力的增大，开始形成第二层吸附；在高区域，受孔隙结构影响，吸附质发生毛细管凝聚，产生吸附滞后，形成滞后环，吸附曲线末端斜率较大，说明储层还有更大的孔隙存在（见图）。根据国际理论化学与应用化学联合会（）分类标准，研究区页岩样品等温吸附脱附曲线类型属于型，滞后环类型属于、型，而对应的孔隙形态分别为墨水瓶型和平行板状狭缝型。图为洗油前后孔体积变化率与孔径关系曲线，可以发现，样品洗油后具有更大的孔体积，且在 和黏土较为发育的样品中更为明显，表明样品发育大量微孔，同时也反映出有机质对页岩微孔的发育具有主要贡献，洗油后释放的孔隙原先可能被油沥青所占据，导致洗油前后孔体积变化明显。图低温氮气等温吸附脱附曲线图孔体积变化率与孔径的关系曲线 实验结果分析某种程度上，孔隙度的大小是孔隙结构差异的定量表现形式。为了精细描述扫描电镜、核磁共振、高压压汞、低温氮气吸附种实验手段评价页岩油孔隙结构的差异性，对其孔隙度测量值进行深入分析。氦气是惰性气体，常温下在岩石中不发生吸附现象，分子较小，能够进入大部分孔隙中，测量结果与真实值最为接近。本次研究选用氦气孔隙度近似作为岩样真实孔隙度。各实验获取孔隙度结果如表所示。）扫描电镜实验显示，研究区页岩油储层实际孔径分布范围从纳米级到微米级，但扫描电镜受视域影响，无法在同一视域下兼顾大小尺度孔隙，更偏向定性认识。）由表可知，核磁孔隙度平均值是氦气孔隙度平均值的近，约 的孔隙度没有被探测，其原因一方面是核磁实验受取样时间及地层压力释放的影响，孤立孔隙中的流体发生逸散；另一方面水分子要大于氦气分子，导致核磁样品的微小孔隙难以完全饱和水，减少了样品孔隙度的测量值。长江大学学报（自然科学版）年月表不同实验方法孔隙度测量结果 样品编号氦气孔隙度核磁孔隙度压汞孔隙度液氮孔隙度 平均值 ）对于高压压汞实验，由于页岩样品整体孔隙体积较小，孤立孔发育，孔隙连通性差，排驱压力大，并且压汞实验受仪器最大压力的影响，无法将样品的孔隙完全冲注。因此，汞只进入了一部分连通性相对较好的孔隙中，导致压汞实验测量孔隙度与氦气法测量偏差较大，与氦气孔隙度相比近 的孔隙度被严重低估，其结果不能完全反映微小孔发育样品的孔隙特征。）低温氮气吸附实验的测量原理对于检测纳米级孔隙以及孤立不连通孔隙更加具有优势，能更好地反映微小孔隙的结构特征，但其对大孔隙反映不够灵敏，液氮孔隙度平均值只占氦气孔隙度平均值的 。综上分析认为，采用单一的传统实验手段不能够充分表征页岩油储层的孔隙结构，储层物性的定量评价误差较大。孔隙特征多尺度联合表征为了从整体上展现页岩油储层的孔隙特征，首先利用核磁共振谱与低温氮气吸附实验的相关性，采用线性方式转换 ，得出核磁共振实验的孔径分布曲线。图为编号 、的样品不同实验方法得出的孔径分布曲线，可以发现，低温氮气吸附在表征孔径小于 的孔隙特征时，能更精细地显示微小孔径的变化，而核磁共振和高压压汞在表征中、大孔隙时效果明显；在 的区间内，随着孔隙变小，高压压汞孔径分布逐渐偏离核磁共振，且孔径越小分布差异越大，这因为高压压汞实验主要针对连通孔喉发育的样品，测量的孔隙度及孔径分布图主要反映的是大孔部分，而对于孤立孔发育、孔喉配比差、以窄颈孔隙为主的页岩样品，高压压汞测量结果低于真实孔隙度，所以中大孔采用核磁共振实验结果表征。图页岩样品不同实验方法测量孔径分布图 第 卷 第期邵广辉 等：玛湖凹陷二叠系风城组页岩油储层微观孔隙结构精细表征表多尺度孔隙度计算结果 样品编号全尺度孔隙度 平均值 综上分析，利用低温氮气吸附和核磁共振的孔径分布曲线拼接拟合，可以实现页岩油储层孔隙的多尺度表征，其多尺度孔隙度如表所示，实验样品的全孔径精细表征图如图所示。在拟合过程中分段求取液氮孔隙度和核磁孔隙度，再依据低温氮气吸附和核磁共振两种实验孔径分布曲线的拼接位置将孔隙度进行累加，得到实验样品整体的孔隙度大小，即全尺度孔隙度，其值与氦气孔隙度基本符图页岩样品全孔径精细表征图 合（见表、表），表明拼接后的孔径分布曲线代表了页岩油储层整体的孔隙分布。分析图发现，研究区页岩油储层孔径分布具有以下个特征：孔径分布范围广，可横跨个数量级，除个别大孔较发育的样品，孔径大于 的孔隙占总体的，的孔隙孔径在 之间，说明页岩样品的孔隙主要以小孔及微孔为主，这与页岩样品粒度细、成熟度高，在成岩过程中普遍经历强烈压实的减孔作用有关，；在孔径小于 范围内，黏土的 含量越高，孔径的主峰位置越靠左，表明有机孔主要集中在小于 孔径的孔隙中，有机孔在页岩的储集空间中占有特殊的一部分。结论）玛湖凹陷二叠系风城组页岩油储层孔隙结构复杂，孔隙度低，孔喉细小，连通性差，孔径分布跨度大，传统的单一实验方法不能满足孔隙结构精细评价要求。）多尺度孔隙结构联合表征结果显示，研究区储层孔隙系统完整，样品的全孔径分布曲线呈多峰，孔径主要集中在 之间，小孔隙与 含量、黏土含量关系密切，整体纳米级孔隙发育，少数样品发育微米级以上孔隙。参考文献金之钧，朱如凯，梁新平，等当前陆相页岩油勘探开发值得关注的几个问题 石油勘探与开发，（）：，（）：姚兰兰 页岩油储层微观孔隙结构特征评价及渗流机理研究 北京：中国科学院大学，：，许琳，常秋生，冯玲丽，等 准噶尔盆地玛湖凹陷二叠系风城组页岩油储层特征及控制因素 中国石油勘探，（）：，（）：，长江大学学报（自然科学版）年月 ，（）：王剑，周路，靳军，等准噶尔盆地吉木萨尔凹陷芦草沟组页岩油储层孔隙结构、烃类赋存及其与可动性关系 石油实验地质，（）：，（）：吴松涛，朱如凯，崔景伟，等非常规储层孔隙结构表征：思路、思考与展望 地质论评，（）：，（）：，：，：李廷微，姜振学，许辰璐，等沾化凹陷沙三下亚段陆相页岩储层微纳米孔隙结构特征 石油科学通报，（）：，（）：陶金雨，申宝剑，潘安阳桂中坳陷北部下石炭统页岩储层微观孔隙结构特征 地质论评，（）：，（）：，（）：，：王子萌，蒋裕强，付永红，等 基于核磁共振表征渝西地区五峰组龙一亚段页岩储层孔隙结构及非均质性 地球科学，（）：，（）：汤天知，李庆峰，赵小青，等基于电成像与核磁共振测井的古龙页岩油储层有效性评价 大庆石油地质与开发，（）：，（）：，（）：（），（）：刘焯，欧阳传湘，王长权，等基于截止值的幂函数构建毛管压力曲线 波谱学杂志，（）：，（）：，：王迪，陈敏，王健伟，等 基于和毛管压力分类转换的渗透率计算新方法 地球物理学进展，（）：，（）：李晋宁泥页岩储层孔隙结构表征和连通方式研究 南京：南京大学，：，张盼盼，刘小平，关铭，等沧东凹陷孔二段低熟页岩纳米孔隙特征及主控因素 特种油气藏，（）：，（）：编辑龚丹第 卷 第期邵广辉 等：玛湖凹陷二叠系风城组页岩油储层微观孔隙结构精细表征