岩石电性物理实验的数值模拟.pdfVIP

岩石电性物理实验的数值模拟 岳文正陶果 (中国石油大学(北京)) 摘要：利用计算机建立具有不同骨架颗粒形状和不同孔隙度的孔隙介质模型， 采用格子气自动机方法模拟孔隙介质在油水饱和度不同的情况下电的传输特性， 揭示了流体饱和度与介质电阻率关系非阿尔奇现象产生的本质原因。发现饱和度 指数是饱和度和孔隙度的函数，它反映了孔隙介质导电路径的复杂程度随饱和度 变化而变化的情况，而不是传统的观点；饱和度指数是流体分布状态的反映。并 得到新的饱和度计算方法。以某油田为例，将计算的饱和度结果与实际岩心实验 的结果作对比，说明新方法的有效性。 关键词：多相流体格子气自动机电阻率增大系数饱和度指数非阿尔奇现象 石油是一种重要的战略资源。随着国家经济的发展，对石油的需求越来越 大。从经济和技术的角度，储层储量评价都是一个重要的关键技术，提高储量 评价技术无论是对国民经济的发展，还是石油储备量的合理确定，都有十分重 要的意义。 在储层含油气性评价的测井资料综合定量解释中，阿尔奇公式作为最基本 的解释关系式具有重要的作用。由于阿尔奇公式[1]是建立在纯砂岩岩石物理电 性实验基础上的半理论半经验公式，因此其在实际应用上存在一定的局限性。 此外，在岩电关系的实验研究中发现了大量的非阿尔奇现象。非阿尔奇现象是 一类非常重要的现象，对于定量的评价储层的含油气性有着至关重要的影响。 对于存在非阿尔奇现象的储层，如果利用阿尔奇公式进行含油气性评价，会造 成含油气饱和度的解释结果存在较大误差，从而影响油气层的识别。因此，研 究非阿尔奇现象的成因，建立相应的岩电关系模型是提高储层含油气定量评价 准确性，提高储量计算精度的关键。 对于非阿尔奇现象的成因及其相应的解释模型，许多学者进行了大量的研 究。李宁[2]通过研究，根据岩，tl,薄片的矿物及流体分析，建立了岩心的串并联 一48— 混合导电网络结构模型，并在此模型的基础上，通过理论分析建立了相应的岩 电关系模型，该岩电关系模型可以很容易地退化成阿尔奇公式、双水模型等， 电的机理非常复杂，非阿尔奇现象较普遍，并提出了GEM模型。 鉴于宏观岩电实验无法观测孔隙微观结构和流体、电流的流动等队岩电实 验的影响，许多人尝试用数值模拟孑L隙微观结构的方法来解释阿尔奇公式参数 的物理意义，寻找定量解释的方法。斯坦福大学的Yald4]曾利用一个三维孔隙 网络模型来模拟储层岩石的流通特性。陶果唧应用这个模型来模拟和解释在饱 步发展了Yale的模型以计算多相流体饱和孔隙介质的电传输特性。前人虽然在 数值研究上取得一定成果，然而他们的这些孔隙介质模型与实际岩石中孔隙结 构和流体分布相差太大，理论计算难以适应极其复杂的实际岩石孔隙结构，计 算结果有很大的局限性，常常得不到令人满意的结果。 本文利用格子气自动机方法对岩心的f—s。关系进行模拟，揭示非阿尔奇 现象的成因，得到饱和度计算公式。并以某油田为例，将计算结果与岩心实验 结果作对比。 1格子气自动机方法 格子气自动机(LGA)模型方法是在细胞自动机基础上演变和发展起来的。 程，并导致了其他一些不合理的模拟结果。1986年，正六边形格子的FHP模 型的提出解决了HPP模型的非旋转对称问题。 格子气自动机模型除了对空间进行离散外，还要对流体和时间进行离散。 把流体看成只有质量而没有体积的微小粒子，只能分布在节点上，沿着格线方 向运动，同时粒子的分布要遵守Pauli不相容原理。时间的离散是一种抽象离 散，把时间离散成整的时间步，每一个时间步分为两部分：粒子在每个节点上 的相互碰撞和从该节点向邻居节点的运移。 以导出达西定律。比较欧姆定律和达西定律可以发现，如果用电荷代替流体粒 一49— 子，用电场强度代替压力梯度，那么两定律式子的形式几乎是相 、／ 同的，这说明在忽略重力影响下，电的流动与其他流体的流动有 ·—号h ／＼ 着相似的规律，由此可见电流也一样满足NS方程。这样，我们 有理由认为格子气方法同样可以用来模拟电的流动。 曼1．}：竺竺三 图1中黑点为格子气体，沿箭头方向运动。在本项研究中， 模型空间