孔隙网络中砂岩和流体模拟的三维孔隙尺度建模讲述.docVIP

孔隙网络中砂岩和流体模拟的三维孔隙尺度建模 Stig Bakke and Pal-Ericφren 王少依 翻译；杨晓宁 校对 摘要： 本文提出了一种新方法，用来建立真实均质和非均质三维孔隙尺度的砂岩模型。此方法的本质是通过模拟砂岩形成的主要地质学过程——砂粒的沉积、压实及成岩作用——建立实际砂岩体的类比模型。模型所需的数据来源于实际砂岩薄片的图像分析。 砂岩模型在X, Y和Z方向上的空间连续性是通过基于独立比例侵入渗透算法确定的。最终得到的空间连续性椭圆函数，可以用来对砂岩模型进行非均质性描述。非均质性分析表明，与垂向相比，压实作用在水平方向上能够更迅速的减少空间连续性。 孔隙空间所呈现的网络结构和几何形态，是直接在经过充分表征的砂岩模型中应用各种3-D图像分析算法得到的。从一个强水润湿的Bentheimer 砂岩薄片数据中得到的3-D孔隙网络被用作一个双相网络流动模拟系统的输入数据。模拟得到的砂岩模型的传导特性和实验得到的结果具有良好的一致性。 前言 高渗透率空隙或孔隙的存在是储集岩区别于其他岩石的特征。储集岩的孔隙高度不规则，其包含的孔隙空间网络由较大孔隙（孔隙本身）和连通这些孔隙的较窄孔隙（孔隙喉道）组成。孔隙空间及其相应颗粒基质的结构和几何形态决定了岩石的一些宏观特性，如绝对渗透率、相对渗透率、毛管压力、地层因素以及电阻率系数等。 由微观结构和孔隙尺度的物理现象预测岩石的宏观特性问题，一直是被广泛研究的主题。此项研究中得到普遍应用的一个工具是网络流体模拟系统（网络模型）。该方法要求对发生在孔隙尺度的物理过程的细致理解，以及对孔隙空间形态的完整描述。通过在微模型中识别出的物理现象，此步骤已经在简单或理想多孔介质的两相流体中得到成功应用。最近，许多三相流体孔隙尺度的物理现象的公开，导致了三相流体网路模型的发展。由于对真实多孔介质三维孔隙结构的精确描述仍然比较困难，所以针对真实多孔介质的技术扩展也十分复杂。关于储集岩孔隙结构的信息经常借助于压汞数据和薄片图像分析。Yanuka 和 Wardlaw 对表征孔隙结构的方法进行了很好的概述。压汞数据提供了孔隙喉道大小分布的静态信息，或者精确点儿说，是被标定的孔隙喉道大小中可能被侵入的体积的分布。对薄片图像进行的传统的孔隙几何形态分析，一般基于Delesse 定律，其认为薄片中的2-D观察可以代表3-D样本。但是，这只适用于标量，例如面积和体积方面的考虑。因此，这个定律适用于孔隙度和矿物成分的测定，而不适用于孔隙和喉道的几何形态测定，因为这些是向量。 为了在薄片图像中表征孔隙或喉道的几何形态，有必要定义2-D孔隙元素中哪些是孔隙，哪些是喉道。进一步，我们需要知道每一个孔隙元中流体的流向，才能测量垂直于流向的平面上孔隙喉道的几何形态。当前，没有方法可以直接的从薄片中提取这些信息。因此，传统孔隙结构分析得到的数据，不足以直接提供对真实岩石复杂孔隙网络的完全3-D描述。 高级技术例如微CT, 连续切片结合BSE切片成像，为岩石的3-D孔隙结构提供了详细的描述。然而这些技术是昂贵而又低效的。利用Finney 构建的完全定义的等球粒密封体，Bryant 经过模拟砂岩形成过程的结果，建造了3-D孔隙网络。重建的多孔介质的方法也已经被用来表征3-D孔隙结构。多孔介质建模也成为了许多化学工程领域的研究课题。本文提出了一种新的建立真实而又充分表征3-D孔隙网络模型的方法。此方法基于砂岩形成的主要的地质学过程结果的数字模型。模型砂岩的孔隙空间网络的表示，通过图像分析技术，从相应的矿物基质网络中分离出来。流体流模拟直接在不规则的孔隙空间进行，然后计算传导特性和现有的实验数据进行对比。 砂岩的3-D过程/随机建模 砂岩是一系列复杂的地质学和流体动力学过程的最终产物，其过程可以简述如下：开始是含石英岩石被侵蚀，然后是砂粒通过空气、水和冰川的搬运，最后在沉积盆地沉积。沉积的颗粒可能会经过再次或多次侵蚀，搬运和沉积的循环而被重新改造。经过压实作用和多种成岩作用的过程，最终形成了砂岩。 砂岩样本及其岩相学参数是其经历的所有地质学和流体动力学过程导致的结果，这些过程影响了沉积盆地及其环境。我们并不想建立这些极度复杂过程的详细动力学模型，而是试图通过使用薄片分析得到的岩相学数据，模拟砂岩形成的主要过程的结果，从而重建3-D砂岩。 砂岩重建过程包含三个主要的建模阶段：（ⅰ）沉积阶段（砂粒沉积），（ⅱ）压实阶段，（ⅲ）成岩阶段。 （ⅰ）沉积阶段 沉积过程的建模开始于粒度分布曲线的产生，而粒度分布曲线来源于真实砂岩薄片的BSE图像。我们利用图像分析技术测定图像中所有砂粒的直径。目前，所有的碎屑砂粒被当作球形的石英颗粒。复杂的地质学和流体动力学过程影响了砂床上每一个单独砂粒的沉积，为了模拟其随机性，