碳酸盐岩岩电特性数字岩心仿真.pdfVIP

JournalofShandcngUniversityofScienceandTechnology 岳文正等 碳酸盐岩岩电特性数字岩心仿真 Natur釉iScience 由于其在处理复杂孔隙几何结构的灵活性和方便性 ，在岩石物理领域扮演越来越重要 的角色。研究证实数 字岩心作为岩石物理实验的补充 ，是一种在 中观尺度研究岩石孑L隙传输特性 的有效手段 。 数字岩心建模可采用 图像处理技术对得到的2D铸体薄片进行处理，得到储层岩石的孔隙结构信息，用 顺序指示模拟方法进行 3D数字岩心建模。薄片的二值 图像每个节点的意义 ，可 以分为 2种 ：孔隙和骨架。 从薄片数据的图像处理中，得到了二值图像的空间统计参数——孔隙度和 自相关函数 ，然后确定随机访问每 个 网格节点路径 ，指定估计网格点的领域条件数据的个数 ，在每个节点处 ，都估算 出一个局部条件累积分布 函数。与先前沿着随机的路线模拟的点一样受薄片图像制约。估算的值作为条件保存下来后 ，沿着随机路 线处理下一个节点。当所有节点都处理后 ，用正确的空间统计 ，就能得到新 的3D数字岩心模型_g。该建模 方法得到的数字岩心能最大程度地与实际的岩石接近，重构的数字岩心不仅可以用来精确地预测物理特性， 而且还可以用来 了解不同物理响应之间的内在关系。这对岩石物理的研究有着十分重要的意义 。 (a)碳酸盐岩岩心铸体薄片 (b)数字岩心模型 图2 碳酸盐岩岩心铸体薄片和数宇岩心模型图 Fig．2 Thecorecastslicesofcarbonaterockanddigitalcoremodel 此外 ，为了检验孔隙空间结构对导 电特性的影响，还可采用大小形状不 同的骨架颗粒来构造数字岩心模 型，通常骨架的复杂程度 同样能够反映孔隙空间的复杂程度。首先可 由岩心薄片确定骨架颗粒 的粒径分布 情况 (图3(a))。根据骨架颗粒分布特征 ，采用数值方法模拟颗粒 的沉积过程 ，形成数字岩心模型，如 图 3(b)。根据需要 ，分别按照不同的孔隙度来建立模型，并按照饱和度大小 ，利用数值模拟流体 的分布特征， 对不同的孔隙结构和饱和流体 的分布情况建立相应的数字岩心模型。 matrixgrain／ixm (a)骨架颗粒粒径分布 (b)岩心模型 图 3 根据骨架颗粒粒径分布建立数字岩心模型 图 Fig．3 Thedigitalcoremodelformedwiththeradiusdistributionofmatrixgrain Imt11111t：~tt 学手 蜜然科学蹶 2第01310年卷6第月3期 VduOnl23001N1o3 数值实验过程中，模型的上下边界采用无滑移固壁边界模拟绝缘边界 。模型的左右边界设为周期性边 界 ，模拟无穷介质。具体过程可以描述为 ：如果以s出(￡)表示第 t时间步流出边界 的状态 ，s人(￡+1)表示第 ￡+1时间步流人边界的状态 ，则 S人(+1)一S出()。设置电压为0．1，初始平均粒子分布密度为0．85。 3 数值模拟结果和分析 在含油气储层饱和度评价 中，目前普遍应用的阿尔奇公式是建立在具有均匀粒 间孔隙的纯砂岩岩石物 理 电性实验基础上的半理论半经验公式 I— bS7”。 (4) 其中：n称为饱和度指数 ，对应于双对数坐标系中阿尔奇直线 的斜率；6为阿尔奇参数 ；对于同一岩心，和6 可视为常量 。 碳酸盐岩具有低孔、低渗、非均质性强等特点，岩 电物 理实验难以获得低含水饱和度情况下的电阻率资料 ，所 以 本研究通过数字岩心模拟不 同含水饱和度下 的岩 电特性 ， 建立含水饱和度计算模型。 首先研究饱和单相流体数字岩心的导电特性。根据不 同孔隙完全水饱和数字岩心的导电特性模拟结果 ，图4