-
南海西部海域目前已发现大量高温高压致密砂岩气藏，作为一种天然的多孔介质气藏，其储层本身具有低孔、低—特低渗透的特征，给这类气田的商业开发带来了巨大挑战。储层孔隙结构是微观岩石物理研究的核心内容[1]，对提高气藏产能和采收率有着积极的作用，而致密砂岩储层在经历构造、沉积、成岩等作用过程中形成更为复杂的微观孔隙结构。随着大量致密油气藏的发现和开发，众多学者借助不同实验方法，从微观角度对致密砂岩储层孔隙结构与渗流规律进行了研究[2-15]，但传统的压汞、铸体薄片、扫描电镜等评价手段难以实现对这类储层内部结构的精细表征，而且经过压汞实验的岩样无法二次利用，对海上高成本的取心样品损耗很大。此外，基于达西定律所建立的渗流方程，只能对流体的运动特征进行宏观表征，无法对岩心孔隙结构（微米尺度）进行精细刻画[16]，更无法实现岩心微-纳米尺度的微观渗流机理研究。近年来快速发展的数字岩心技术[17-21]，在一定程度上实现了对储层内部微观结构的定量表征，尤其在微观层面对低渗透致密砂岩储层内部不同尺度的孔喉空间分布特征有着更为细致的刻画，从而进一步揭示了这类储集体岩石物理属性的变化规律和微观渗流机理。
-
笔者以莺歌海盆地乐东 10 区高温高压气田致密砂岩储层为研究对象，开展岩石组构定量分析，借助多尺度 Micro-CT 扫描技术建立三维孔隙结构模型，实现对致密砂岩储层数字岩心的两相渗流过程模拟，明确致密砂岩储层微观渗流机理，以期为乐东10区高温高压气田的合理开发提供理论基础。
-
1 区域地质概况
-
莺歌海盆地总体形态呈北北西走向的长条形，由东南部的莺歌海凹陷和西北部的河内凹陷组成，两凹陷间为临高凸起带分隔。根据构造变形特征，进一步将莺歌海盆地划分为莺东斜坡带、临高凸起带、中央坳陷和莺西斜坡带 4 个一级构造单元（图1）。乐东10区气田位于莺歌海盆地莺东斜坡，主要含气层位黄流组为大型水道中细砂岩储层，砂岩厚度大、分选好。岩心孔隙度主要为7.7%～14.5%，中值为 10%；渗透率主要为 0.05～10.58 mD，中值为 0.51 mD，整体属于特低孔、特低渗透储层，局部发育物性相对较好的“甜点”储层。气藏埋深为3 916.0～4 322.3 m，地层压力系数高，为 2.174～2.306，原始地层温度为 181.25～192.45℃，为高温超高压弹性气驱岩性气藏。
-
图1 乐东10 区构造位置
-
Fig.1 Structure location at study area
-
2 岩石组构定量分析
-
2.1 矿物组分特征
-
选取研究区 2口井 7块致密砂岩样品进行数字岩心分析，样品的常规物性实验测试显示，孔隙度为 7.7%～11.3%，渗透率为 0.059～0.268 mD。采用纳米级扫描电镜能谱（EDS）和矿物标准库的岩石定量测试系统进行矿物组分及含量测试，并与具备更高分辨率且可表征密度分布的背散射电子成像 （BSE）对照，可获得精度远高于 EDS 分辨率的矿物分布，准确区分粒径为几十纳米的矿物。测试结果表明（表1），岩石类型为长石质石英砂岩，主要矿物成分为石英、长石和碳酸盐与黏土胶结物，构成了致密砂岩孔隙空间发育的物质基础。长石以较为稳定的钾长石、钠长石为主，基性斜长石消耗殆尽。碳酸盐胶结物含量较高，主要发育方解石、白云石和铁白云石。黏土矿物含量较低，以伊利石为主，另有少量绿泥石、高岭石和蒙皂石。碎屑颗粒之间呈现-凹凸接触。
-
2.2 孔隙分布特征
-
基于样品背散射电子成像下的岩石矿物与孔隙发育特征定量分析，乐东10区致密砂岩储层孔隙空间发育程度主要受压实作用、胶结作用和溶蚀作用等成岩作用影响，高温高压环境对孔隙发育存在一定的改善。整体上储层孔隙空间以次生孔隙为主，原生孔隙为辅，其中次生孔隙以长石粒内溶孔和铸膜孔为主，原生孔隙受强烈胶结作用的影响所占比例较低。粒间溶孔与铸膜孔孔隙尺寸最大，这部分孔隙整体分布较为分散，无法形成有效的孔隙网络（图2a）；长石颗粒内部发育较多溶解作用形成的粒内溶孔，保存较完整的长石颗粒中可见解理缝发育，部分解理缝内部有其他矿物填充（图2b，2c）；胶结物主要为方解石和白云石，其内部发育较多晶间孔（图2d）；黏土矿物及云母主要作为粒内/粒间充填矿物存在（图2e，2f），由于黏土矿物间的转换作用以及孔隙压力差异，黏土矿物及云母内晶间孔尺寸差异较大。
-
图2 背散射电子成像下样品孔隙空间发育特征
-
Fig.2 Characteristics of pore space in BSE imaging
-
基于对扫描电镜图像的二值化处理，提取出样品中的孔隙并进行定量分析。从孔隙半径分布（图3）可以看出，孔隙半径分布曲线呈现双峰或三峰的特征，说明孔隙发育呈现多尺度特征。其中粒间溶孔及铸膜孔孔隙半径主要集中在 40~300 μm 左右，这部分孔隙尺寸较大，是孔隙发育的主要构成部分。部分残余粒间孔、粒内溶孔及胶结物溶孔的孔隙半径主要集中在 3~40 μm，为次级孔隙发育的主要范围。胶结物内的晶间孔孔隙半径最小，这部分孔隙分布较广，且数量较多，其孔隙半径主要集中在 0.070~3 μm。孔隙结构整体上表现为大孔细喉的特征，孔隙半径中等、喉道细，孔隙连通性较差。
-
3 基于数字岩心的孔隙结构模型建立
-
数字岩心技术是以真实岩样为基础，通过一系列的图像处理技术和数值算法将岩心数字化[22]，从而在孔隙尺度上较真实地反映岩石微观孔隙结构。目前，业界常用的数字岩心构建方法主要有物理实验方法和数值重建方法[18，23-26]，前者是借助高精度光学仪器直接获取岩心的平面图像，结果较为准确、直观，但成本高、耗时长，常用的有CT扫描法、聚焦离子束成像法和序列成像组合法等；后者是借助岩心二维切片，利用各种数学算法间接实现基于数字岩心的孔隙结构三维构建，主要包含随机法和过程模拟法，这种方法简便易行、建模速度快，但得到的孔隙结构参数与实际岩心往往存在一定的差异。
-
图3 不同类型孔隙对应的半径分布范围
-
Fig.3 Distribution of different types of pore radius
-
本次研究采用 CT 扫描法，利用 FEI Heliscan 微米CT采集预处理后柱塞样品的真实三维投影数据，选取合适的图像重建方法将投影数据重建岩心灰度图像，最后采用图像二值分割方法分离灰度图像中的孔隙空间和岩石骨架后得到岩心内部孔隙结构的真实数字化模型。考虑到样品中孔隙发育的多尺度特征，采用 2 种不同扫描精度对样品内部不同尺寸范围的孔隙进行捕捉。分别以 4.7 μm 扫描精度对直径为 8 mm 的岩心进行图像扫描，以 0.65 μm扫描精度对直径为2 mm的岩心进行图像扫描。
-
3.1 CT扫描数据处理
-
图像扫描及滤波处理 结合孔隙半径分布特征，首先对直径为 8 mm 的样品进行 4.7 μm 尺度的扫描得到岩心整体三维图像，用于捕捉孔隙半径较大的孔隙信息。基于对图像的非均质性分析，优选次级孔隙较发育的区域钻取直径为 2 mm 的小柱塞并进行 0.65 μm 尺度的高分辨率精细扫描，用于捕捉孔隙半径中等的残余粒间孔、粒内溶孔及胶结物晶间孔。经过扫描获取岩心的三维灰度图像后，应用高斯及中值滤波对图像进行优化处理，便于减少扫描系统背景噪声的干扰，提高灰度图像的质量与对比度，方便下一步的图像分割与孔隙骨架提取 （图4）。
-
图4 滤波处理前后灰度图像对比
-
Fig.4 Comparison of grey images before and after filtering
-
图像二值化处理 二值化处理是一个极为重要的图像前期处理过程，核心是通过阈值的确定来提取孔隙空间与岩石骨架信息，从而分割出图像中的孔隙部分[27-28]。以岩心实验测量孔隙度为依据，通过选择合适的灰度阈值，将岩心CT灰度图像中的孔隙与岩石骨架进行分割，使分割的岩心图像中孔隙像素数占总像素数的比例等于岩心孔隙度，确保了分割结果的可靠性（图5）。分割阈值计算式为：
-
图5 岩心CT图像二值化处理前后对比
-
Fig.5 Comparison of core CT images before and after binarization
-
表征单元体提取 数字岩心表征性体积单元 （简称表征单元体，记为 REV）是岩石相关特性趋于基本稳定时的最小体积。当岩石体积小于REV时，其特性随体积变化而变化，表现出随机波动特性。当岩石体积大于或等于 REV 时，其特性趋于稳定，不会发生明显变化。同时由于岩心CT图像体积较大，为提高计算效率，通常从图像中截取代表性较强的一块区域进行研究。因此需要通过 REV 分析确定对不同孔隙结构具有代表性的最小样品或模型尺寸。
-
由于不同扫描方法的视域范围和扫描精度不同，每种图像所对应的孔隙尺度也有所差异。SEM 图像视域范围大，扫描精度高，可以表征所有尺度的孔隙；由于分辨率的限制，直径为 8 mm 柱塞样品主要捕捉粒间溶孔和铸膜孔等较大孔隙；2 mm柱塞样品扫描精度较高，而视域范围较小，主要捕捉长石溶孔、岩屑粒间孔等较小孔隙。通过对 3 种扫描方式所获取的图像不同体元/面元大小下的孔隙度误差计算，可以确定微观研究的表征单元体边长。本次研究设定误差值小于5%时样品孔隙度趋于稳定，误差分析结果表明，扫描电镜粗扫柱塞样品的直径大于 3 mm 时，所有尺度的孔隙度基本保持稳定，基于 8与 2 mm 柱塞样品的 REV 分析，其岩心尺寸分别为 1.5与 1 mm 左右，对应的孔隙度基本保持稳定，据此确定8 mm柱塞样品扫描建立的六面体数字岩心模型每条边长像素数大于 320，2 mm 柱塞样品扫描建立的六面体数字岩心模型每条边长像素数大于1 539。
-
3.2 三维孔隙结构模型建立
-
针对研究区致密砂岩样品孔隙呈现多尺度发育特征的情况，选取 2 级扫描精度对样品孔隙结构进行获取，并利用分水岭分割方法，构建2种尺度的三维孔隙结构模型（图6），其体素分别为 4.7和 0.65 μm。应用连通域分析方法对三维孔隙结构模型的连通性进行分析，0.65 μm 尺度模型可以捕捉到主要的渗流通道和相应的连通孔隙空间，而4.7 μm尺度模型捕捉的主要为相互非连通的较大一级别孔隙，该级别的孔隙发育空间对孔隙度贡献较大，对渗透率影响较小。
-
图6 2种尺度三维孔隙结构模型
-
Fig.6 3D graphs of pore structures for two scales
-
基于 2 种尺度的孔隙结构模型，利用双尺度孔隙耦合的方法计算2种尺度孔隙结构模型中孔隙同时参与的情况下岩心整体的孔渗参数。双尺度孔隙耦合方法首先对大柱塞图像进行三值化分割，分别为孔隙、包含次级孔隙的区域、不含孔隙的骨架区域，而后提取孔隙网络时将包含次级孔隙的区域作为特殊的等效喉道处理，其孔隙度及渗透率可由精细模型（0.65 μm 尺度）计算得到，而等效喉道的进出口端流量变化可由达西公式计算得到，最后根据侵入-逾渗理论对孔隙结构模型进行求解，得到模型的整体孔隙度及渗透率[29]。
-
基于单尺度孔隙结构模型进行计算后得到 4.7 μm 尺度孔隙结构模型孔隙度为 5.1%，渗透率为 0，即模型不连通；0.65 μm 尺度孔隙结构模型孔隙度为 4.79%，渗透率为 0.051 mD。利用双尺度孔隙耦合的方法计算得到岩心整体孔隙度为 9.65%，渗透率为 0.097 mD。该结果与实验室测试结果（孔隙度为11.33%，渗透率为0.10 mD）吻合程度较高。
-
为进一步验证粒间孔、铸膜孔等微米-亚毫米级孔隙对岩石孔渗参数的影响，通过减少4.7 μm尺度孔隙结构真实模型中粒间孔、铸膜孔等微米-亚毫米级的孔隙数量，得到孔隙度分别为 1.25%， 2.89%，3.5% 的 3 个模型（表2，图7）。在 0.65 μm 尺度孔隙结构模型不变的情况下，当双尺度融合孔隙结构模型中大孔隙数量上升时，整体孔隙度上升明显，但整体渗透率无明显变化。进一步说明粒间孔、铸膜孔等微米-亚毫米级孔隙对孔隙度影响明显，但由于自身无法形成有效渗流网络，对整体渗透率无明显改善效果。
-
4 渗流机理
-
4.1 基于孔隙结构模型的气水两相流模拟
-
基于数字岩心的孔隙结构模型是利用简单几何体所组成的空间网络来代表真实的复杂多孔介质的孔隙空间，利用这些简单几何体等效为真实多孔介质的孔隙和喉道，并通过拟静态流动模型模拟流体在这些简单几何体中的流动过程，以近似代替流体在真实多孔介质中的流动[30-33]。
-
利用已建立的孔隙结构模型进行两相流过程渗流模拟，从模拟结果可以看出，当只存在单尺度的微孔隙时（图8a），由于气相作为非润湿相，在毛细管力的限制下优先进入尺寸较大的喉道，占据了有利的渗流通道，从而导致水相渗透率迅速下降； 而当模型中同时存在大孔隙及微孔隙时（图8b），气相克服最粗喉道进入孔隙，首先开始充填附近连接的大孔隙空间。由于大孔隙体积较大，延缓了气相侵入喉道的时间，因此等渗点的含水饱和度向左偏移，水相渗透率下降速度明显缓于单尺度模型的。
-
图7 4套双尺度融合孔隙结构模型中的大尺度孔隙空间
-
Fig.7 Large scale pore space for four two-scale fusion pore structure models
-
图8 不同尺度孔隙结构模型气水两相流模拟结果
-
Fig.8 Simulation results of gas-water two-phase flow for pore structure models at different scales
-
4.2 微观孔喉发育模式与渗流机理
-
基于所构建的2种尺度三维孔隙结构模型连通性分析，大孔隙尺度模型中的粒间溶孔及铸膜孔互不连通，无法仅通过这部分孔隙形成有效渗流通道，而微孔隙尺度模型中的长石溶孔、岩屑粒间孔为主的微孔隙是样品的主要渗流通道。
-
小孔小喉渗流模式 在以长石溶孔、岩屑粒间孔为主的微孔隙尺度下，孔隙与喉道半径差异不大，渗流空间由微孔隙和喉道构成，渗流模式呈现小孔小喉特征（图9a）。该模式中，孔隙连通率对残余水饱和度影响较大。在孔隙结构相近的情况下，随着孔隙连通率升高，残余水饱和度逐渐降低。孔隙连通率低，表明死孔隙或孔隙盲端较多，而此类孔隙结构无法形成渗流通道，根据孔隙形态不同，残余水以水柱或水珠形式封闭于死孔隙或盲端中，使得残余水饱和度变大。在气水两相流残余水状态下，强表面张力和高毛细管力共同作用下以厚水膜形式赋存于微孔隙及喉道壁上，残余水饱和度较高。而极细喉道在高毛细管力作用下形成卡断，残余水以水柱形式堵塞喉道，残余水饱和度极高。
-
图9 不同孔隙结构模型下气水两相流模式
-
Fig.9 Gas-water two-phase flow modes for different pore structure models
-
大孔小喉渗流模式 在双尺度融合后的孔隙结构模型中，存在粒间溶孔、铸膜孔等大孔隙，渗流模式呈大孔小喉特征（图9b）。在大孔隙参与下，气水渗流模式较小孔小喉模式更为复杂。大孔隙作为占比较高的储集空间，在驱替阶段对含水饱和度影响较大。在驱替开始阶段，部分气相侵入孔隙空间，大孔隙中部分含水，残余水饱和度较高；大孔隙提供的储集空间减缓了气体侵入的进程，微孔隙及喉道部分仍然保持了水相的连通，极细喉道在高毛细管力作用下基本未被气体侵入。在残余水状态下，大孔隙中的残余水在表面张力作用下以薄水膜形式赋存于孔隙表面，残余水饱和度小；微孔隙及喉道中的残余水在强表面张力和高毛细管力共同作用下以厚水膜形式赋存于微孔隙及喉道壁上，残余水饱和度较高；极细喉道中的残余水在高毛细管力作用下形成卡断，残余水以水柱形式堵塞喉道，残余水饱和度极高。
-
5 结论
-
乐东 10 区致密砂岩储层主要发育 3 种级别的孔隙，粒间溶孔及铸膜孔的孔隙半径最大，主要集中在40~300 μm左右；部分残余粒间孔、粒内溶孔及胶结物溶孔孔隙半径中等，主要集中在3~40 μm；胶结物内的晶间孔孔隙半径最小，主要集中在 0.070~3 μm。
-
构建4.7和0.65 μm尺度的三维孔隙结构模型， 4.7 μm尺度模型反映相互非连通的较大一级孔隙，对孔隙度贡献较大，对渗透率影响较小；0.65 μm尺度模型反映主要的渗流通道和相应的连通孔隙空间。利用双尺度孔隙耦合方法进行三维孔隙结构模型的渗透率计算，模型计算结果与实验测试渗透率结果基本吻合。
-
以长石溶孔、岩屑粒间孔为主的微孔隙发育，岩心表现为小孔小喉渗流模式；以粒间溶孔、铸膜孔等大孔隙及次级孔隙发育，岩心表现为大孔小喉渗流模式。大孔隙尺度中的粒间溶孔及铸膜孔互不连通，无法通过这部分孔隙形成有效渗流通道，而微孔隙尺度中的长石溶孔、岩屑粒间孔为主的微孔隙是样品的主要渗流通道。
-
符号解释
-
abs——k 取不同值时对图像的孔隙度与岩心孔隙度的差取绝对值，%；
-
i——灰度值，无量纲；
-
Imax，Imin——图像的最大、最小灰度值，无量纲；
-
k——从Imin—Imax依次取得的孔隙分割阈值，无量纲；
-
k*——最优阈值，无量纲；
-
p（i）——灰度值为i的像素数，无量纲；
-
ϕ——岩心孔隙度，%。
-
参考文献
-
[1] 罗哲潭，王允诚.油气储集层的孔隙结构［M］.北京：科学出版社，1986.LUO Zhetan，WANG Yuncheng.Pore structure of oil and gas res⁃ ervoir［M］.Beijing：Science Press，1986.
-
[2] 孙卫，史成恩，赵惊蛰，等.X-CT 扫描成像技术在特低渗透储层微观孔隙结构及渗流机理研究中的应用——以西峰油田庄19井区长82储层为例［J］.地质学报，2006，80（5）：775-779.SUN Wei，SHI Cheng’en，ZHAO Jingzhe，et al.Application of XCT scanned image technique in the research of micro-pore tex⁃ ture and percolation mechanism in ultra-permeable oil field-tak⁃ ing an example from Chang82 Formation in the Xifeng oil field［J］.Acta Geologica Sinica，2006，80（5）：775-779.
-
[3] 关振良，谢从姣，董虎，等.多孔介质微观孔隙结构三维成像技术［J］.地质科技情报，2009，28（2）：115-121.GUAN Zhenliang，XIE Congjiao，DONG Hu，et al.3D imaging and visualization technology of micro pore structure in porous media ［J］.Geological Science and Technology Information，2009，28（2）：115-121.
-
[4] 杨浩.低渗透油层微观特征及其与宏观指标关系研究［D］.大庆：大庆石油学院，2010.YANG Hao.Research on microscopic characteristics and the rela⁃ tionship between microscopic characteristics and macro indica⁃ tors of low permeability reservoir［D］.Daqing：Daqing Petroleum Institute，2010.
-
[5] 朱维耀.特低渗透油藏有效开发渗流理论和方法［M］.北京：石油工业出版社，2010.ZHU Weiyao.Seepage theory and method for effective develop⁃ ment of ultra-low permeability reservoir［M］.Beijing：Petroleum Industry Press，2010.
-
[6] 苏娜，段永刚，于春生.微CT扫描重建低渗气藏微观孔隙结构——以新场气田上沙溪庙组储层为例［J］.石油与天然气地质，2011，32（5）：792-796.SU Na，DUAN Yonggang，YU Chunsheng，et al.Reconstruction of microscopic pore structure in low permeability gas reservoirs by micro-CT scanning：an example from the Upper Shaximiao For⁃ mation in Xinchang gas field［J］.Oil & Gas Geology，2011，32（5）：792-796.
-
[7] 张全培，吴文瑞，刘丽萍，等.鄂尔多斯盆地镇北地区延长组超低渗透储层孔隙结构及其分形特征［J］.油气地质与采收率，2020，27（3）：20-31.ZHANG Quanpei，WU Wenrui，LIU Liping，et al.Pore structure and fractal characteristics of ultra-low permeability reservoirs in Yanchang Formation in Zhenbei area，Ordos Basin［J］.Petroleum Geology and Recovery Efficiency，2020，27（3）：20-31.
-
[8] 郝乐伟，王琪，唐俊.储层岩石微观孔隙结构研究方法与理论综述［J］.岩性油气藏，2013，25（5）：123-128.HAO Lewei，WANG Qi，TANG Jun.Research progress of reser⁃ voir microscopic pore structure［J］.Lithologic Reservoirs，2013，25（5）：123-128.
-
[9] 程乐利，印森林，万友利，等.羌塘盆地北羌塘坳陷上三叠统巴贡组致密储层成岩作用与孔隙演化［J］.石油实验地质，2020，42（1）：60-68.CHENG Leli，YIN Senlin，WAN Youli，et al.Diagenesis and pore evolution of tight sandstone reservoirs in Upper Triassic Bagong Formation，North Qiangtang Depression，Qiangtang Basin［J］.Pe⁃ troleum Geology & Experiment，2020，42（1）：60-68.
-
[10] 蒋裕强，陈林，蒋婵，等.致密储层孔隙结构表征技术及发展趋势［J］.地质科技情报，2014，33（3）：63-70.JIANG Yuqiang，CHEN Lin，JIANG Chan，et al.Characterization techniques and trends of the pore structure of tight reservoirs［J］.Geological Science and Technology Information，2014，33（3）：63-70.
-
[11] 盛军，孙卫，赵婷，等.致密砂岩气藏微观孔隙结构参数定量评价——以苏里格气田东南区为例［J］.西北大学学报：自然科学版，2015，45（6）：913-924.SHENG Jun，SUN Wei，ZHAO Ting，et al.Quantitative evaluation of microscopic pore structure parameters of tight sandstone gas reservoir：A case the Southeast Area of Sulige Gas Field［J］.Jour⁃ nal of Northwest University：Natural Science Edition，2015，45（6）：913-924.
-
[12] 符勇，白玉彬，黄礼，等.靖边地区长6段致密砂岩储层孔隙结构及影响因素［J］.断块油气田，2020，27（1）：34-39.FU Yong，BAI Yubin，HUANG Li，et al.Pore structure and influ⁃ encing factors of Chang6 tight sandstone reservoir in Jingbian area ［J］.Fault-Block Oil and Gas Field，2020，27（1）：34-39.
-
[13] ROSENBERG E，LYNCH J，GUEROULT P，et al.High resolu⁃ tion3D reconstructions of rocks and composites［J］.Oil & Gas Sci⁃ ence and Technology，l999，54（4）：497-51l.
-
[14] 惠威，薛宇泽，白晓路，等.致密砂岩储层微观孔隙结构对可动流体赋存特征的影响［J］.特种油气藏，2020，27（2）：87-92.HUI Wei，XUE Yuze，BAI Xiaolu，et al.Influence of micro-pore structure on the movable fluid occurrence in tight sandstone reser⁃ voir［J］.Special Oil & Gas Reservoirs，2020，27（2）：87-92.
-
[15] 马铨峥，杨胜来，杨龙，等.准噶尔盆地吉木萨尔凹陷芦草沟组致密储层微观孔隙特征［J］.大庆石油地质与开发，2020，39（6）：13-20.MA Quanzheng，YANG Shenglai，YANG Long，et al.Characteris⁃ tics of the micro-pore in Lucaogou-Formation tight reservoir of Jimsar Sag，Junggar Basin［J］.Petroleum Geology & Oilfield De⁃ velopment in Daqing，2020，39（6）：13-20.
-
[16] 姚军，赵秀才，衣艳静，等.数字岩心技术现状及展望［J］.油气地质与采收率，2005，12（6）：52-54.YAO Jun，ZHAO Xiucai，YI Yanjing，et al.The current situation and prospect on digital core technology［J］.Petroleum Geology and Recovery Efficiency，2005，12（6）：52-54.
-
[17] 邓世冠，吕伟峰，刘庆杰，等.利用CT技术研究砾岩驱油机理 ［J］.石油勘探与开发，2014，41（3）：330-335.DENG Shiguan，LÜ Weifeng，LIU Qingjie，et al.Research on dis⁃ placement mechanism in conglomerate using CT scanning method ［J］.Petroleum Exploration and Development，2014，41（3）：330-335.
-
[18] 王平全，陶鹏，刘建仪，等.基于数字岩心的低渗透储层微观渗流机理研究［J］.非常规油气，2016，3（6）：1-5.WANG Pingquan，TAO Peng，LIU Jianyi，et al.Mechanism of mi⁃ cro seepage in low-permeability reservoirs based on digital core ［J］.Unconventional Oil & Gas，2016，3（6）：1-5.
-
[19] 张思勤，汪志明，洪凯，等.基于格子Boltzmann方法的3D数字岩心渗流特征分析［J］.测井技术，2016，40（1）：18-22.ZHANG Siqin，WANG Zhiming，HONG Kai，et al.Analysis of seepage characteristics of 3D digital core based on Boltzmann method［J］.Well Logging Technology，2016，40（1）：18-22.
-
[20] 李骞，张钰祥，李滔，等.基于数字岩心建立的评价碳酸盐岩完整孔喉结构的方法——以川西北栖霞组为例［J］.油气地质与采收率，2021，28（3）：53-61.LI Qian，ZHANG Yuxiang，LI Tao，et al.A method for evaluating complete pore-throat structure of carbonate rocks based on digital cores：A case study of Qixia Formation in Northwest Sichuan［J］.Petroleum Geology and Recovery Efficiency，2021，28（3）：53-61.
-
[21] 盛军，阳成，徐立，等.数字岩心技术在致密储层微观渗流特征研究中的应用［J］.西安石油大学学报：自然科学版，2018，33（5）：83-89.SHENG Jun，YANG Cheng，XU Li，et al.Application of digital core technology in the study of microscopic seepage characteris⁃ tics of tight reservoirs［J］.Journal of Xi’an Shiyou University：Nat⁃ ural Science Edition，2018，33（5）：83-89.
-
[22] 高兴军，齐亚东，宋新民，等.数字岩心分析与真实岩心实验平行对比研究［J］.特种油气藏，2015，22（6）：93-96.GAO Xingjun，QI Yadong，SONG Xinmin，et al.Parallel compari⁃ son of digital core analysis and real core test［J］.Special Oil & Gas Reservoirs，2015，22（6）：93-96.
-
[23] 赵秀才，姚军.数字岩心建模及其准确性评价［J］.西安石油大学学报：自然科学版，2007，22（2）：16-20.ZHAO Xiucai，YAO Jun.Construction of digital core and evalua⁃ tion of its quality［J］.Journal of Xi’an Shiyou University：Natural Science Edition，2007，22（2）：16-20.
-
[24] 刘学锋，张伟伟，孙建孟.三维数字岩心建模方法综述［J］.地球物理学进展，2013，28（6）：3 066-3 072.LIU Xuefeng，ZHANG Weiwei，SUN Jianmeng.Methods of con⁃ structing 3-D digital cores：A review［J］.Progress in Geophysics，2013，28（6）：3 066-3 072.
-
[25] 屈乐，孙卫，杜环虹，等.基于CT扫描的三维数字岩心孔隙结构表征方法及应用——以莫北油田116井区三工河组为例 ［J］.现代地质，2014，28（1）：190-196.QU Le，SUN Wei，DU Huanhong，et al.Characterization technique of pore structure by 3D digital core based on CT scanning and its application：An example from Sangonghe Formation of 116 Well Field in Mobei Oilfield［J］.Geoscience，2014，28（1）：190-196.
-
[26] 林承焰，吴玉其，任丽华，等.数字岩心建模方法研究现状及展望［J］.地球物理学进展，2018，33（2）：679-689.LIN Chengyan，WU Yuqi，REN Lihua，et al.Review of digital core modeling methods［J］.Progress in Geophysics，2018，33（2）：679-689.
-
[27] 朱益华，陶果，方伟.图像处理技术在数字岩心建模中的应用 ［J］.石油天然气学报，2007，29（5）：54-57.ZHU Yihua，TAO Guo，FANG Wei.Application of image process⁃ ing technique in digital core modeling［J］.Journal of Oil and Gas Technology，2007，29（5）：54-57.
-
[28] 赵秀才.数字岩心及孔隙网络模型重构方法研究［D］.东营：中国石油大学（华东），2009.ZHAO Xiucai.Numerical rock construction and pore network ex⁃ traction［D］.Dongying：China University of Petroleum（East Chi⁃ na），2009.
-
[29] 罗瑜，王寅，王容，等.碳酸盐岩孔隙-裂缝双重网络模型构建及分析［J］.新疆石油地质，2021，42（1）：107-112.LUO Yu，WANG Yin，WANG Rong，et al.Construction and analy⁃ sis of pore-fracture network model of carbonate rock［J］.Xinjiang Petroleum Geology，2021，42（1）：107-112.
-
[30] BAKKE S，OREN P E.3-D pore-scale modelling of sandstones and flow simulations in the pore networks［J］.SPE Journal，1997，2：136-149.
-
[31] BLUNT M J.Physically-based network modeling of multiphase flow in intermediate-wet porous media［J］.Journal of Petroleum Science and Engineering，1998，20：117-125.
-
[32] 王克文，孙建孟，关继腾，等.储层岩石电性特征的逾渗网络模型［J］.石油学报，2007，28（1）：101-106.WANG Kewen，SUN Jianmeng，GUAN Jiteng，et al.Percolation network modeling of electrical properties of reservoir rock［J］.Ac⁃ ta Petrolei Sinica，2007，28（1）：101-106.
-
[33] 陶军，姚军，李爱芬，等.利用孔隙级网络模型研究油水两相流 ［J］.油气地质与采收率，2007，14（2）：74-77.TAO Jun，YAO Jun，LI Aifen，et al.Research on oil and water flow using pore-scale network model［J］.Petroleum Geology and Re⁃ covery Efficiency，2007，14（2）：74-77.
-
摘要
针对致密砂岩储层微观孔隙结构复杂导致孔隙空间展布和渗流机理认识不清的问题，以莺歌海盆地乐东10 区高温高压特低渗透致密砂岩储层为例，利用三维高精度Micro-CT扫描、纳米级电镜扫描等图像分析技术对致密砂岩储层的孔隙空间展布规律进行定量表征，建立基于数字岩心的双尺度三维孔隙结构模型，并基于双尺度孔隙耦合渗流模拟方法确定致密砂岩储层渗流机理。研究结果表明：①乐东10区致密砂岩储层主要发育3种尺度范围的孔隙，孔隙半径分别为0.070~3，3~40，40~300 μm；②在以长石溶孔、岩屑粒间孔为主的微孔隙尺度下，孔隙与喉道半径差异较小，渗流空间以微孔隙和喉道构成，表现为小孔小喉渗流模式，在气水两相渗流过程中水相渗透率下降快；③在由较大粒间溶孔、铸模孔与次级尺度孔隙共同组成的双尺度孔隙结构模型中，呈大孔小喉渗流模式，水相渗透率下降较缓。
Abstract
The complex micropore structures of tight sandstone reservoirs leads to an unclear understanding of pore space distribution and flow mechanism. Accordingly，taking the high-temperature，high-pressure，and ultra-low permeability tight sandstone reservoir at Ledong 10 area in Yinggehai Basin as an example，the pore space characteristicss of tight sand- stone reservoirs with image analysis techniques such as 3D high-precision micro CT scanning and nano electron micro- scope scanning were quantitatively characterized in this paper. A two-scale 3D pore structure model based on digital cores was constructed and the flow mechanism of tight sandstone reservoirs based on the two-scale pore coupling flow simulation was determined. The results show as following：①The tight sandstone reservoir at Ledong10 area mainly develops three scale ranges of pores，which are respectively 0.070-3 μm，3-40 μm，and 40-300 μm. ②In the micropore scale dominated by feldspar dissolved pores and rock debris intergranular pores，the difference between pore and throat radius is small. The flow space is composed of micropores and throats，which are characterized by small pores and small throats. The water phase permeability decreases rapidly in the process of gas-water two-phase flow. ③In the two-scale pore structure model composed of large intergranular dissolved pores，mold pores，and secondary scale pores，the flow model presents the charac- teristics of large pores and small throat，and the water phase permeability decreases slowly.
