摘要
本发明公开一种粗糙裂隙入侵‑逾渗两相流的模拟方法，采用平面网格单元对粗糙裂隙开度场进行划分，建立入侵‑逾渗两相流数值模型，将驱替过程模拟为准静态的界面推进过程，在单位时间步内具有最小毛细压力阈值的单元被入侵驱替；采用Young‑Laplace方程，通过计算平面内曲率半径和开度方向曲率半径，确定每一待驱替单元的毛细压力阈值；采用二叉树数据结构，记录每个单元的驱替过程信息，并在驱替的单个时间步内，更新变化的二叉树数据结构；结合深度搜索算法和宽度搜索算法，判别和找寻毛细捕获的单元位置信息。本发明通过考虑平面内曲率对入侵驱替过程的影响，并结合高效的数据结构与搜索算法，实现了粗糙裂隙中入侵‑逾渗驱替过程的精确、高效模拟。
