过去几十年来，人们对裂隙多孔介质流动的数值模拟越来越感兴趣。迄今为止，大多数研究都集中在二维或伪三维计算模型上，其中三维复杂结构对渗流的影响尚未完全解决。本工作提出了一种三维非均质多孔介质渗流建模方法。能够模拟随机离散裂隙和夹杂物这种复杂的结构。提出一种网格策略对复杂域进行离散化。特别地，为了保证网格的一致性，对相交单元进行了处理。然后进行数值离散，其中包括裂隙、包裹体和围岩基质的通量相互作用。对三维裂隙介质水力特性进行数值试验分析。首先，将数值解与嵌入离散断裂模型的结果进行对比，验证了所开发的框架。接下来分析裂缝的方位、开度和半径对流体流动和等效渗透率张量的影响。研究了非均质和均质介质中压力分布的变化。最后，研究了含复杂结构介质的水力特性，揭示了裂缝与夹杂物水力特性的差异。

图1

可视化

(上行)和网格划分(下行)几种典型模式

图2 基于DARSim的EDFM网格和仿真结果

过去几十年来，人们对裂隙多孔介质流动的数值模拟越来越感兴趣。迄今为止，大多数研究都集中在二维或伪三维计算模型上，其中三维复杂结构对渗流的影响尚未完全解决。本工作提出了一种三维非均质多孔介质渗流建模方法。能够模拟随机离散裂隙和夹杂物这种复杂的结构。提出一种网格策略对复杂域进行离散化。特别地，为了保证网格的一致性，对相交单元进行了处理。然后进行数值离散，其中包括裂隙、包裹体和围岩基质的通量相互作用。对三维裂隙介质水力特性进行数值试验分析。首先，将数值解与嵌入离散断裂模型的结果进行对比，验证了所开发的框架。接下来分析裂缝的方位、开度和半径对流体流动和等效渗透率张量的影响。研究了非均质和均质介质中压力分布的变化。最后，研究了含复杂结构介质的水力特性，揭示了裂缝与夹杂物水力特性的差异。

图1

可视化

(上行)和网格划分(下行)几种典型模式

图2 基于DARSim的EDFM网格和仿真结果