在地球深处，巨大的油气储层往往并不像充满孔隙的砂岩那样如湿海绵般存在，而是被困在复杂的裂缝和缝隙网络之中。这些裂缝岩石为试图预测流体如何在其中流动的工程师们带来了独特的挑战。在储层工程领域，两种无形的力量支配着油和水的流动：毛细管压力和相对渗透率。毛细管压力是使流体保持在微小空间内的力量，而相对渗透率则描述了当另一种流体已经存在时，一种流体流动的难易程度。几十年来，模拟这些裂缝系统的标准方法一直依赖于一个便捷但有缺陷的假设：即这些裂缝非常宽且均匀，以至于毛细管压力实际上为零，并且流体以完美线性、可预测的方式相互滑动。这种简化虽然让计算变得容易，但近期的证据表明这在物理上是错误的，会导致对可采油量及其流动速度的预测出现显著误差。
伊朗石油技术大学的一个研究小组开发了一种全新的、统一的方法来描述这些裂缝岩石中的流体流动，该方法摒弃了这些旧有的捷径。他们没有通过任意数字进行猜测或拟合实验数据曲线，而是构建了一个从裂缝的实际物理形状出发的模型。他们认识到没有任何裂缝是完美光滑或均匀的；岩壁之间的间隙（被称为开度）在不同位置的变化极大。通过将这种变化视为一种特定的统计模式，他们创建了一个将裂缝几何形状与内部流体行为直接联系起来的框架。他们的工作提供了一种清晰的、基于物理的方法，仅通过一个关键测量值——裂缝开度的标准差——即可计算出裂缝的毛细管压力和相对渗透率。这个描述裂缝宽度变化程度的单一数字，取代了以往模型中困扰已久的多个混乱的拟合参数。
研究人员首先将裂缝构想为一个网络，而非单一的平坦间隙，在这个网络中，开口大小沿其长度和宽度不断变化。他们确定这些变化遵循特定的统计分布，类似于人群中身高的差异。利用这种分布，他们推导出了数学关系，展示了水和油如何在裂缝内分布。较小的开口倾向于留住润湿流体（如水），而较大的开口则允许非润湿流体（如油）通过。通过整合这些物理现实，他们生成了一组“类型曲线”。这些曲线充当了一张通用地图，使工程师只需了解裂缝宽度的变化程度，就能预测其行为。如果宽度非常一致，流动方式便呈现一种状态；如果宽度变化显著，流动方式则表现得不同。这种方法消除了旧有的、过度简化的“X型”曲线，因为那些曲线假设无论裂缝多么复杂，流体都会沿直线运动。
为了证明其模型的有效性，研究小组将其与广泛的现有数据进行了对比测试，包括理论研究和真实的实验室实验。他们将预测结果与数十项研究的数据进行了比较，涵盖了从光滑的实验室裂缝到粗糙的天然岩石表面。在裂缝物理属性已知的情况下，他们的模型以极高的精度匹配了观测数据。在物理属性未知的情况下，他们使用了一种匹配技术来寻找单一的最佳拟合变化参数，随后成功预测了不同条件下的流体行为。结果是一致的：新模型能够高精度地复制复杂的实验数据，其表现通常优于依赖多个经验常数的旧方法。至关重要的是，该模型捕捉到了裂缝在压力下的变化。随着上方岩层的重量挤压岩石，裂缝会闭合，其宽度变化也会随之改变。研究人员展示了他们的模型可以追踪这些机械变化，调整流动预测以反映裂缝新的、更紧凑的几何结构。
这项工作的意义不仅在于对计算方法的改进。通过提供一个依赖于单一具有物理意义的参数而非一系列拟合数字的模型，研究人员提供了一个既简单又稳健的工具。这意味着工程师现在可以直接将这些更真实的流动描述集成到用于管理油田的大规模计算机模拟中。他们不再依赖于“裂缝是没有任何内部阻力的开放通道”这一过时的假设，而是可以模拟流体在实际且不均匀的岩石几何结构中进行的复杂“舞蹈”。该研究证实，裂缝毛细管压力并非为零，且相对渗透率并非简单的直线关系；它们是复杂的、取决于几何形状的现象，必须予以考虑才能理解储层的真实潜力。虽然该模型是针对单个裂缝开发的，但作者指出，该框架可以扩展到处理整个裂缝网络以及更复杂的三相流，从而为理解地球破碎地壳中的隐藏动力学提供一条更清晰的路径。
技术摘要：具有物理一致性的裂缝毛细管压力型曲线与相对渗透率
问题陈述
准确模拟裂缝储层中的流体流动仍然是一个重大挑战，因为这涉及到复杂的裂缝网络几何形状以及复杂的基质-裂缝相互作用。虽然裂缝毛细管压力 (Pc) 和相对渗透率 (kr) 是控制相干涉和驱替效率的基本组成关系，但传统的油藏工程实践通常依赖于过于简化的假设。具体而言，数值模拟器经常假设裂缝中毛细管压力为零，且相对渗透率表现为线性（X型曲线）行为。近期的实验和建模证据表明，这些假设缺乏坚实的物理基础，可能导致性能预测出现显著误差。此外，现有模型往往依赖于多个缺乏明确物理意义或具有特定案例性的经验“拟合参数”，限制了其在不同裂缝几何形状下的预测能力。在裂缝微观结构特征（特别是孔隙度分布）与这些流动参数之间存在直接联系的环节上，目前仍存在关键性的研究空白。
研究方法
本研究通过将裂缝孔隙度分布与基本流动力学相结合，开发了一个统一的、基于物理机制的裂缝 Pc 和 kr 模型框架。核心研究方法包括以下步骤：
- 统计特征表征： 将裂缝系统概念化为一个包含完全连通网络的单位体积，其中单个裂缝表现出空间变化的孔隙度。研究采用对数正态分布来表征裂缝的孔隙度、宽度和长度，并使用正态分布来表征裂缝的方向性。孔隙度标准差 (σb) 被确定为主要的控制参数。
- 相饱和度推导： 通过计算在由毛细管压力阈值定义的特定孔隙度范围内的相体积，推导出润湿相和非润湿相的可动饱和度和不可动饱和度。这些范围对应于不可逆润湿相、残留非润湿相以及总相饱和度。
- 型曲线开发：
- 毛细管压力： 将无量纲润湿相饱和度公式化为有效饱和度和无量纲裂缝毛细管压力 (PcDf) 的函数。通过引入对数正态孔隙度分布，利用误差函数 (erf) 推导出排驱和吸吮过程的裂缝毛细管压力型曲线。这些曲线取决于单个物理参数 σb 以及无量纲边界压力 (PcDminf, PcDmaxf)。
- 相对渗透率： 本模型通过整合用于导管流体流动的 Hagen-Poiseuille 方程与用于孔隙介质的达西定律来推导裂缝相对渗透率。该集成过程考虑了孔隙度分布以及排驱和吸吮过程中的相饱和度。模型中还引入了曲折度系数，其指数通过与实验数据的验证确定。
- 验证与匹配： 将模型与广泛的理论框架（如 Pruess 和 Tsang, Zhou 和 Wheater, Liu 等）以及实验数据集（如 Reitsma 和 Kueper, Huo 等, Chen, Diomampo）进行对比验证。
- 当 σb 已知时，模型直接生成曲线以供比较。
- 当 σb 未知时，采用型曲线匹配程序。将实验测得的 Pc 与饱和度关系在对数-对数坐标系下绘制，并与开发的无量次元型曲线进行匹配，以提取最优的 σb 及其他参数。
- 定量指标： 使用均方误差 (MSE) 和平均绝对相对误差 (AARE) 来评估模型的准确性。
主要贡献
- 统一的基于孔隙度的模型： 本研究提出了一个单一框架，直接从裂缝孔隙度分布推导出ทั้ง 裂缝毛细管压力和相对渗透率，消除了对任意经验拟合系数的需求。
- 单一物理参数： 模型仅依赖于裂缝孔隙度标准差 (σb) 作为唯一输入参数。该参数具有明确的物理可解释性，代表了裂缝的几何异质性，这与需要多个特定案例参数的传统模型形成了对比。
- 型曲线方法论： 引入的裂缝毛细管压力型曲线实现了双重应用：(1) 从已知的孔隙度分布预测 Pc 和 kr；(2) 通过型曲线匹配，从实验 Pc 数据中表征裂缝孔隙度分布。
- 数学简洁性与鲁棒性： 推导出的方程利用标准的统计函数（误差函数），旨在实现与商业油藏数值模拟器的无缝集成。
结果
开发的模型在多种裂缝类型和流动条件下均表现出卓越的鲁棒性和有效性：
- 毛细管压力： 模型成功复制了理论模型（如 Zhou 和 Wheater），实现了近乎完美的吻合（MSE 在 10−9 至 10−10 范围内），并与实验数据（如 Huo 等, Reitsma 和 Kueper）高度匹配（AARE 通常 < 20%）。值得注意的是，模型捕捉到了物理趋势，即随着有效应力的增加，优化的 σb 会随之增加，这反映了裂缝闭合和异质性的增加。
- 相对渗透率： 模型准确预测了排驱和吸吮过程的相对渗透率曲线。它与理论极限（如 λ→∞ 时的 Brooks 和 Corey 模型）以及实验数据集（如 Chen, Diomampo）表现出强一致性。模型成功地将相干涉从次要实验变量中解耦出来，提供了比传统基于饱和度图表更严谨的验证。
- 参数效率： 在所有验证案例中，该模型在仅使用单一参数 σb 的情况下，达到了与复杂现有模型相当或更优的准确度。
意义与主张
本文声称，所提出的框架为确定裂缝毛细管压力和相对渗透率提供了一种可靠且具有成本效益的方法，可以替代高成本、耗时的实验室实验。通过将模型建立在裂缝孔隙度分布的物理现实而非经验曲线拟合之上，本研究为油藏工程师提供了一个具有物理一致性的工具。
作者强调，模型的数学结构便于直接应用于商业油藏模拟器，从而实现对裂缝系统更真实的规模化模拟。研究结论指出，裂缝毛细管压力既不是零也不是常数；裂缝相对渗透率也显著偏离线性“X型曲线”；必须使用基于孔隙度的模型来解释这些现象，以准确捕捉相干涉。这项工作弥合了复杂解析解与实际油藏工程应用之间的鸿沟，提供了一种无需面对任意拟合参数歧义的标准方法来表征裂缝介质。
每周获取最佳 electrical engineering 论文。
受到斯坦福、剑桥和法国科学院研究人员的信赖。
请查收邮箱确认订阅。
出了点问题，再试一次？
无垃圾邮件，随时退订。