一种用于颗粒多孔介质的三维分形渗透率模型，该模型考虑了孔隙的粗糙度因素

《Journal of Contaminant Hydrology》：A three-dimensional fractal permeability model for granular porous media incorporating pore roughness

三维孔隙介质渗透率模型构建及表面粗糙度效应研究，提出基于球形颗粒簇的三维等效毛细管束模型，改进Hagen-Poiseuille方程引入表面粗糙度参数，实验验证显示模型在低渗透介质中预测精度达85%以上，突破传统二维分形模型的局限性。

　　
本文针对低渗透率多孔介质渗流特性建模中的关键科学问题，提出了一种融合三维孔隙结构特征与表面粗糙度效应的解析分形渗透率模型。研究团队通过构建球状颗粒簇的三维孔隙结构模型，创新性地将复杂多孔介质等效为毛细管束网络体系，在继承传统分形模型优势的基础上，首次系统性地纳入孔隙表面粗糙度的影响因素。该成果为解决低渗透介质渗流预测难题提供了新的理论框架和技术路径。

在研究背景方面，低渗透率地质介质（如煤层、页岩、致密砂岩等）的渗流行为具有显著的非线性特征。传统Kozeny-Carman模型等经典理论基于理想化平滑孔隙表面和二维结构假设，难以准确表征纳米尺度孔隙结构的三维分形特征及其表面粗糙度对渗流阻力的影响。实验研究表明，当孔隙直径进入微纳米尺度（通常小于1微米）时，表面粗糙度导致的流动路径收缩效应可使渗透率下降达2-3个数量级，这对传统模型形成了严峻挑战。

为突破现有模型的局限性，研究团队建立了多维度创新体系：首先在微观结构表征方面，采用球状颗粒簇模型模拟真实多孔介质的三维空间排布。通过控制颗粒尺寸分布和空间堆积方式，实现了孔隙连通性、曲折度与表面粗糙度的协同描述。其次在等效模型构建方面，创造性提出将非均质孔隙空间等效为三维毛细管束网络。这种等效转换既保留了原结构的多尺度分形特征，又使得渗流力学方程的推导成为可能，有效解决了传统二维模型无法表征三维空间流动特性的问题。

在表面粗糙度量化方面，研究引入了分形几何中的参数化方法。通过建立粗糙度要素的尺寸分布函数，将表面微观形貌特征转化为可计算的工程参数。特别值得注意的是，粗糙度对渗透率的影响被分解为局部流动截面缩减效应和整体路径曲折度增强效应两个独立维度，这种双参数耦合机制突破了以往单一粗糙度修正模式的理论局限。

模型验证阶段采用分层验证策略：首先通过参数敏感性分析确定关键变量，包括孔隙率φ、分形维度Df、曲折度τav、相对粗糙度ε等核心参数。然后选取不同岩心样本（砂岩、页岩、煤层）进行交叉验证，结果显示预测值与实验数据在10^-15至10^-8 mD量级范围内均具有良好吻合性（R2>0.92）。特别在致密砂岩样本中，当孔隙直径小于500nm时，模型对表面粗糙度参数的敏感性分析准确预测了渗透率的非线性衰减规律。

该研究在理论层面实现了三突破：其一，构建了首个融合三维孔隙结构、表面粗糙度及多尺度分形特征的统一理论框架，突破了传统模型在空间维度和表面效应方面的双重局限；其二，建立了基于分形粗糙度参数的等效流动截面积修正公式，将表面粗糙度从定性描述提升为定量分析工具；其三，创新性地将Hagen-Poiseuille定律拓展到三维非均质介质，推导出考虑粗糙度影响的修正渗透率表达式。

工程应用方面，研究提出的模型展现出显著优势。在非常规油气储层预测中，通过输入实际岩心的三维CT扫描数据，可精确计算分形维度（Df=2.63-2.89）和平均曲折度（τav=5.2-8.7），较传统模型提高预测精度达37%。在碳酸盐岩缝洞型储层评价中，该模型成功揭示了纳米级孔隙表面粗糙度对超临界流体渗透率的影响规律，为开发低伤害压裂技术提供了理论支撑。特别是在页岩气开发领域，模型对微纳米孔隙中"分子级"渗流行为的预测精度达到实验误差范围以内（<8%）。

未来研究可沿三个方向深化：首先在微观结构表征方面，建议引入分子动力学模拟技术，精确解析表面粗糙度与流体-岩石相互作用机制；其次在模型拓展方面，可考虑将非球形颗粒的形状因子纳入等效模型，进一步提升复杂孔隙体系的描述精度；最后在工程应用层面，需要建立不同地质条件下模型参数的快速标定方法，开发配套的工程软件实现模型参数的自动化提取和渗透率预测。

该研究的重要启示在于：低渗透率多孔介质的渗流行为本质上是孔隙结构几何特征、表面微观形貌与流体动力学特性共同作用的结果。当孔隙尺寸进入纳米尺度范围时，表面粗糙度的影响权重从传统模型的30%提升至65%以上，这一发现彻底改变了人们对于多孔介质渗流机制的理解。研究建立的参数化分析框架，为后续开发智能压裂控制算法、纳米级孔隙材料设计等提供了重要理论工具。