Abstract:We study single-phase flow in a fractured porous medium at a macroscopic scale that allows to model fractures individually. The flow is governed by Darcy's law in both fractures and porous matrix. We derive a new mixed-dimensional model, where fractures are represented by $(n-1)$-dimensional interfaces between $n$-dimensional subdomains for $n\ge 2$. In particular, we suggest a generalization of the model in [22] by accounting for asymmetric fractures with spatially varying aperture. Thus, the new model is particularly convenient for the description of surface roughness or for modeling curvilinear or winding fractures. The wellposedness of the new model is proven under appropriate conditions. Further, we formulate a discontinuous Galerkin discretization of the new model and validate the model by performing two- and three-dimensional numerical experiments.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the single - phase flow problem in fractures with different apertures in porous media. Specifically, the researchers aim to develop a new hybrid - dimensional model to describe the geometric characteristics of fractures and their influence on fluid flow more accurately.
### Problem Background
In many practical applications, such as geothermal energy, enhanced oil recovery, groundwater flow, nuclear waste disposal, and carbon sequestration, the existence of fractures significantly affects or even dominates the fluid flow in porous media. Fractures usually have extremely thin apertures but extensive extensions, forming narrow heterogeneous regions, which pose challenges to classical continuum modeling. Traditional models usually assume that fractures are planar and have a constant aperture, which cannot accurately describe the situation of surface roughness or curved fractures.
### Research Objectives
To solve the above problems, the researchers proposed a new hybrid - dimensional model, in which fractures are represented as (n - 1)-dimensional interfaces, and the porous matrix is an n - dimensional sub - domain (\(n\geq2\)). This new model is particularly suitable for describing asymmetric fractures with spatially varying apertures, and thus can better describe surface roughness or curved fractures.
### Main Contributions
1. **Derivation of the New Model**: Starting from a full - dimensional decomposition system, the researchers derived a new model considering spatially varying apertures through the weak - form averaging method and approximate curvilinear integration.
2. **Model Validation**: The effectiveness and accuracy of the new model were verified through two - dimensional and three - dimensional numerical experiments.
3. **Mathematical Proof**: It was proved that the new model is well - posed under appropriate conditions, and the discontinuous Galerkin discretization method was introduced.
### Formulas and Methods
The key to the new model lies in introducing the effective quantities averaged over fractures and simplifying the problem by decomposing the integral within the fracture into a surface integral along Γ and a line integral in the normal direction. In addition, jump and average operators were introduced to handle the pressure and flux on both sides of the fracture:
\[
\llbracket f \rrbracket(t):=f|_{\Gamma_2}(t)-f|_{\Gamma_1}(t)
\]
\[
\{\{ f \}\}(t):=\frac{1}{2}(f|_{\Gamma_1}(t)+f|_{\Gamma_2}(t))
\]
For the effective pressure \(p_\Gamma\) within the fracture and the total source term \(q_\Gamma\), they are defined as follows:
\[
p_\Gamma(t):=\frac{1}{d(t)}\int_{-d_1(t)}^{d_2(t)}p_f(\eta,t)d\eta
\]
\[
q_\Gamma(t):=\int_{-d_1(t)}^{d_2(t)}q_f(\eta,t)d\eta
\]
Through these formulas and methods, the researchers successfully established a hybrid - dimensional model capable of handling complex fracture geometries, thus providing a powerful tool for research in related fields.