格子波尔兹曼方法及其在研究岩石流通特性中的应用.pdf

格子波尔兹曼方法及其在研究岩石流通特性中的应用 朱益华陶果 (中国石油大学(北京)) 摘要：在岩石物理的研究领域，计算岩石物理已经成为一种必然的发展趋势。 与传统的计算流体动力学方法(如FDM、FEM等)相比，格子波尔兹曼方法 (LBM)具有算法简单、精度高、压力可以直接求解、能够模拟具有复杂边界条 件的问题以及适合于并行计算等优点，能够模拟形状极其复杂的真实孔隙岩石中 的两相流动和成岩过程，由岩心的SEM、CT照片来相当准确地计算绝对和相对 渗透率、毛细管压力曲线等。这在非均质复杂油气藏与各物理场关系的研究中， 是一项有广阔应用前景的开创性研究工作，对于建立符合我国地质特点的测井新 方法和解释模型具有重要的意义。 关键词：格子气自动机(LGA)格子波尔兹曼(LB)计算岩石物理数字 岩心数值模拟 在岩石物理的研究领域，数值模拟方法在过去的10年里取得了令人瞩目的 进展。岩石物理的研究对象是由形状极其复杂的孔隙、岩石骨架及多相填充物 组成的体系。其主要任务是建立电磁场、声场、核物理场及核磁共振响应与岩 石孔隙、骨架体积和单相及多相流体饱和度、渗透率等参数的关系。正如量子 场的模拟对于材料科学的进展和分子尺度的模拟对于生物学的研究一样，透彻 地理解和精确地模拟孔隙尺度的物理和化学变化机理对于岩石物理的研究也是 至关重要的，特别是对于那些非线性、多物理场耦合的问题。传统的模拟油气 藏条件下的岩心实验常常是非常困难和昂贵的。在这种宏观实验过程中岩心的 微观孔隙结构、组分及流体分布状态不可被直接观察和控制。另外，目前的实 验设备还难于达到需要的测量范围(例如，准确测量孔隙度小于6％、饱和度 小于20％的岩石物理参数)。而数值模拟技术一开始就显现出这些方面的优势。 型。他使用规则的二维或三维网格结构来代表孑L隙结构。这种简单的网络结构 有一些特定的孔隙结构信息，比如，连通性及孔隙大小分布。同时，提供了一 一62— 种简单的易实现的流体计算模型。网络由一些节点构成，这些节点描述了孔隙 网络模型的最大限制就是需要将孔隙空间分成离散的孔隙和孔喉。在真实的孔 隙空间中，孔隙和孑L喉并无截然的区别。因此，在区别它们时很随意。所以， 网络模型不是真实孑L隙结构的反映，这些模型与实际岩石中孑L隙结构和流体分 布相差太大，理论计算难以适应极其复杂的实际岩石孑L隙结构，计算结果有很 大的局限性，常常得不到令人满意的结果。 一些其他领域内常用的、针对宏观物理性质的经典方法，如有限差分 很精确，但需要非常繁杂的网格剖分方法和特殊的库函数，需要大量地存储中 间结果，因此一般用于似稳场问题。FDM虽然容易计算，能用于时变场问题， 常常将复杂的孑L隙结构简化。所以，用这些常规的方法模拟孔隙尺度下的多相 流问题是非常困难的(Rothman，1988)[s,93。 岩石物理的研究一直需要一种高效、可靠并且可以处理复杂的微孔隙几何 就随之应运丽生。 当代统计物理研究认为，众多宏观复杂系统的整体复杂性，可表现为大量 运动规律十分简单的基本单体之间相互作用的结果。从微观机制出发研究宏观 复杂系统的运动规律，也是一种值得探索的科学研究方法。按照这个思路，法 et 化，建立直接模拟物理系统演化的格子气模型。与从连续系统的偏微分方程出 发的传统数值方法正好相反，该方法的根本思想是对数学物理问题重新建立离 一63— 散模型。物理系统的宏观特性直接由模型的演化方程来计算，从而导出某些宏 观的物理参数等。另外，从计算的角度来看，格子气自动机模型按纯布尔运动 演化，其算子是局域化的，具有完全并行性。因此它又是一类求解数值问题的 非数值大规模并行算法。这样，格子气自动机方法不仅克服了传统数值方法中 的主要矛盾——方法的稳定性和数值舍人误差，而且对各种常见边界问题的处 理也十分灵活方便。 该方法引起了物理学家、计算机科学家与数学家们的浓厚兴趣，致使格子 气方法无论是理论分析还是计算机模拟都取得了许多进展，关于格子气的应用 领域日益扩大，相关的文献也层出不穷。