An official website of the United States government
This paper is concerned with the numerical solution of compressible fluid flow in a fractured porous medium. The fracture represents a fast pathway (i.e., with high permeability) and is modeled as a hypersurface embedded in the porous medium. We aim to develop fast-convergent and accurate global-in-time domain decomposition (DD) methods for such a reduced fracture model, in which smaller time step sizes in the fracture can be coupled with larger time step sizes in the subdomains. Using the pressure continuity equation and the tangential PDEs in the fracture-interface as transmission conditions, three different DD formulations are derived; each method leads to a space-time interface problem which is solved iteratively and globally in time. Efficient preconditioners are designed to accelerate the convergence of the iterative methods while preserving the accuracy in time with nonconforming grids. Numerical results for two-dimensional problems with non-immersed and partially immersed fractures are presented to show the improved performance of the proposed methods.
more »
« less
Flow in porous media with low dimensional fractures by employing enriched Galerkin method
This paper presents the enriched Galerkin discretization for modeling fluid flow in fractured porous media using the mixed-dimensional approach. The proposed method has been tested against published benchmarks. Since fracture and porous media discontinuities can significantly influence single- and multi-phase fluid flow, the heterogeneous and anisotropic matrix permeability setting is utilized to assess the enriched Galerkin performance in handling the discontinuity within the matrix domain and between the matrix and fracture domains. Our results illustrate that the enriched Galerkin method has the same advantages as the discontinuous Galerkin method; for example, it conserves local and global fluid mass, captures the pressure discontinuity, and provides the optimal error convergence rate. However, the enriched Galerkin method requires much fewer degrees of freedom than the discontinuous Galerkin method in its classical form. The pressure solutions produced by both methods are similar regardless of the conductive or non-conductive fractures or heterogeneity in matrix permeability. This analysis shows that the enriched Galerkin scheme reduces the computational costs while offering the same accuracy as the discontinuous Galerkin so that it can be applied for large-scale flow problems. Furthermore, the results of a time-dependent problem for a three-dimensional geometry reveal the value of correctly capturing the discontinuities as barriers or highly-conductive fractures.
more »
« less
- Award ID(s):
- 1913016
- PAR ID:
- 10166992
- Date Published:
- Journal Name:
- Advances in water resources
- Volume:
- 142
- Issue:
- 103620
- ISSN:
- 0309-1708
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Unconventional shale or tight oil/gas reservoirs that have micro-/nano-sized dual-scale matrix pore throats with micro-fractures may result in different fluid flow mechanisms compared with conventional oil/gas reservoirs. Microfluidic models, as a potential powerful tool, have been used for decades for investigating fluid flow at the pore-scale in the energy field. However, almost all microfluidic models were fabricated by using etching methods and very few had dual-scale micro-/nanofluidic channels. Herein, we developed a lab-based, quick-processing and cost-effective fabrication method using a lift-off process combined with the anodic bonding method, which avoids the use of any etching methods. A dual-porosity matrix/micro-fracture pattern, which can mimic the topology of shale with random irregular grain shapes, was designed with the Voronoi algorithm. The pore channel width range is 3 μm to 10 μm for matrices and 100–200 μm for micro-fractures. Silicon is used as the material evaporated and deposited onto a glass wafer and then bonded with another glass wafer. The channel depth is the same (250 nm) as the deposited silicon thickness. By using an advanced confocal laser scanning microscopy (CLSM) system, we directly visualized the pore level flow within micro/nano dual-scale channels with fluorescent-dyed water and oil phases. We found a serious fingering phenomenon when water displaced oil in the conduits even if water has higher viscosity and the residual oil was distributed as different forms in the matrices, micro-fractures and conduits. We demonstrated that different matrix/micro-fracture/macro-fracture geometries would cause different flow patterns that affect the oil recovery consequently. Taking advantage of such a micro/nano dual-scale ‘shale-like’ microfluidic model fabricated by a much simpler and lower-cost method, studies on complex fluid flow behavior within shale or other tight heterogeneous porous media would be significantly beneficial.more » « less
-
Abstract Many geo‐engineering applications, for example, enhanced geothermal systems, rely on hydraulic fracturing to enhance the permeability of natural formations and allow for sufficient fluid circulation. Over the past few decades, the phase‐field method has grown in popularity as a valid approach to modeling hydraulic fracturing because of the ease of handling complex fracture propagation geometries. However, existing phase‐field methods cannot appropriately capture nucleation of hydraulic fractures because their formulations are solely energy‐based and do not explicitly take into account the strength of the material. Thus, in this work, we propose a novel phase‐field formulation for hydraulic fracturing with the main goal of modeling fracture nucleation in porous media, for example, rocks. Built on the variational formulation of previous phase‐field methods, the proposed model incorporates the material strength envelope for hydraulic fracture nucleation through two important steps: (i) an external driving force term, included in the damage evolution equation, that accounts for the material strength; (ii) a properly designed damage function that defines the fluid pressure contribution on the crack driving force. The comparison of numerical results for two‐dimensional test cases with existing analytical solutions demonstrates that the proposed phase‐field model can accurately model both nucleation and propagation of hydraulic fractures. Additionally, we present the simulation of hydraulic fracturing in a three‐dimensional domain with various stress conditions to demonstrate the applicability of the method to realistic scenarios.more » « less
-
In this work, we review and describe our computational framework for solving multiphysics phase-field fracture problems in porous media. Therein, the following five coupled nonlinear physical models are addressed: displacements (geo-mechanics), a phase-field variable to indicate the fracture position, a pressure equation (to describe flow), a proppant concentration equation, and/or a saturation equation for two-phase fracture flow, and finally a finite element crack width problem. The overall coupled problem is solved with a staggered solution approach, known in subsurface modeling as the fixed-stress iteration. A main focus is on physics-based discretizations. Galerkin finite elements are employed for the displacement-phase-field system and the crack width problem. Enriched Galerkin formulations are used for the pressure equation. Further enrichments using entropy-vanishing viscosity are employed for the proppant and/or saturation equations. A robust and efficient quasi-monolithic semi-smooth Newton solver, local mesh adaptivity, and parallel implementations allow for competitive timings in terms of the computational cost. Our framework can treat two- and three-dimensional realistic field and laboratory examples. The resulting program is an in-house code named IPACS (Integrated Phase-field Advanced Crack Propagation Simulator) and is based on the finite element library deal.II. Representative numerical examples are included in this document.more » « less
-
抄録 It is important to understand the long-term migration of radionuclides when considering long-lasting rock engineering projects such as the geological disposal of radioactive waste. The network of fractures and pores in a rock mass plays a major role in fluid migration as it provides pathways for fluid flow. The geometry of such a network can change due to fracture sealing by fine-grained material over extended periods of time. Groundwater commonly contains fine-grained material such as clay minerals, and it is probable that such minerals accumulate within rock fractures during groundwater flow, thereby decreasing fracture apertures and bulk permeability. It is therefore essential to conduct permeability measurements using water that includes fine-grained minerals in order to understand the evolving permeability characteristics of rock. However, this has not been studied to date in in-situ rock mass. Therefore, in the present study, we perform permeability measurements in a granite rock mass to investigate the change of permeability that occurs under the flow of water that includes clays. Our findings show that clay particles accumulate in fractures and that the permeability (hydraulic conductivity) of the granite rock mass decreases over time. The decrease was more significant in the earlier time. We conclude that the accumulation of clay minerals in the fracture decreases the permeability of a rock mass. Furthermore, we consider that the filling and closure of fractures in rock is possible under the flow of groundwater that contains clay minerals.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/https://doi.org/10.1016/j.advwatres.2020.103620
