In this article, we present a time-dependent dual-porosity–Navier–Stokes model with four interface conditions, including Beavers–Joseph interface condition, to describe a coupling system of complex porous media and conduit networks. This system has many applications, such as the flow simulation problems for a multistage fractured horizontal wellbore with suitable boundary/interface conditions and complex geometries. For this coupling system, a decoupled finite element method is proposed based on Robin type transmission conditions. In order to avoid the iteration for the traditional domain decomposition, the interface information, which is needed for the Robin type transmission conditions at the current time step, is predicted directly from the numerical solution of the previous time steps. The stability and convergence of this Robin-type decoupled finite element method are rigorously analyzed. A series of analysis techniques are utilized to analyze various components of this complicated multi-physics problem term by term, especially for the Beavers–Joseph interface condition, the interaction terms of dual-porosity model, and the nonlinear advection of Navier–Stokes equation. Moreover, a pair of Robin parameters are considered in the scheme and analysis, and their effects on the robustness for the case with low permeability and low storativity are numerically investigated. Numerical experiments are provided to validate the convergence of the proposed algorithm and illustrate the features of the application to flow problems around multistage fractured horizontal wellbore.

在本文中，我们提出了一个具有四种界面条件（包括 Beavers-Joseph 界面条件）的瞬态双孔隙度 - Navier-Stokes 模型，以描述复杂多孔介质和管道网络的耦合系统。该系统有许多应用，例如具有合适边界/界面条件和复杂几何形状的多级压裂水平井筒的流动模拟问题。针对该耦合系统，提出了基于Robin型传输条件的解耦有限元方法。为了避免传统域分解的迭代，直接从前一时间步的数值解中预测当前时间步的Robin型传输条件所需的界面信息。严格分析了这种Robin型解耦有限元方法的稳定性和收敛性。一系列的分析技术被用来逐项分析这个复杂的多物理场问题的各个组成部分，特别是对于 Beavers-Joseph 界面条件、双孔隙模型的相互作用项和 Navier-Stokes 方程的非线性平流。此外，在方案和分析中考虑了一对Robin参数，并数值研究了它们对低渗透性和低储存性情况下稳健性的影响。数值实验验证了所提出算法的收敛性，并说明了该算法在多级压裂水平井眼周围流动问题中的应用特点。一系列分析技术被用来逐项分析这个复杂的多物理场问题的各个组成部分，特别是对于 Beavers-Joseph 界面条件、双孔隙模型的相互作用项和 Navier-Stokes 方程的非线性平流。此外，在方案和分析中考虑了一对Robin参数，并数值研究了它们对低渗透性和低储存性情况下稳健性的影响。数值实验验证了所提出算法的收敛性，并说明了该算法在多级压裂水平井眼周围流动问题中的应用特点。一系列分析技术被用来逐项分析这个复杂的多物理场问题的各个组成部分，特别是对于 Beavers-Joseph 界面条件、双孔隙模型的相互作用项和 Navier-Stokes 方程的非线性平流。此外，在方案和分析中考虑了一对Robin参数，并数值研究了它们对低渗透性和低储存性情况下稳健性的影响。数值实验验证了所提出算法的收敛性，并说明了该算法在多级压裂水平井眼周围流动问题中的应用特点。双孔隙度模型的相互作用项和 Navier-Stokes 方程的非线性平流。此外，在方案和分析中考虑了一对Robin参数，并数值研究了它们对低渗透性和低储存性情况下稳健性的影响。数值实验验证了所提出算法的收敛性，并说明了该算法在多级压裂水平井眼周围流动问题中的应用特点。双孔隙度模型的相互作用项和 Navier-Stokes 方程的非线性平流。此外，在方案和分析中考虑了一对Robin参数，并数值研究了它们对低渗透性和低储存性情况下稳健性的影响。数值实验验证了所提出算法的收敛性，并说明了该算法在多级压裂水平井眼周围流动问题中的应用特点。