skip to main content

Title: 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):
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Advances in water resources
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. 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
  4. 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
  5. null (Ed.)
    Abstract. The continuum of behavior that emerges during fracturenetwork development in crystalline rock may be categorized into threeend-member modes: fracture nucleation, isolated fracture propagation, andfracture coalescence. These different modes of fracture growth producefracture networks with distinctive geometric attributes, such as clusteringand connectivity, that exert important controls on permeability and theextent of fluid–rock interactions. To track how these modes of fracturedevelopment vary in dominance throughout loading toward failure and thushow the geometric attributes of fracture networks may vary under theseconditions, we perform in situ X-ray tomography triaxial compressionexperiments on low-porosity crystalline rock (monzonite) under upper-crustalstress conditions. To examine the influence of pore fluid on the varyingdominance of the three modes of growth, we perform two experiments undernominally dry conditions and one under water-saturated conditions with 5 MPa ofpore fluid pressure. We impose a confining pressure of 20–35 MPa and thenincrease the differential stress in steps until the rock failsmacroscopically. After each stress step of 1–5 MPa we acquire athree-dimensional (3D) X-ray adsorption coefficient field from which weextract the 3D fracture network. We develop a novel method of trackingindividual fractures between subsequent tomographic scans that identifieswhether fractures grow from the coalescence and linkage of several fracturesor from the propagation of a single fracture. Throughout loading in all ofthe experiments, the volume of preexisting fractures is larger than that ofnucleating fractures, indicating that the growth of preexisting fracturesdominates the nucleation of new fractures. Throughout loading until close tomacroscopic failure in all of the experiments, the volume of coalescingfractures is smaller than the volume of propagating fractures, indicatingthat fracture propagation dominates coalescence. Immediately precedingfailure, however, the volume of coalescing fractures is at least double thevolume of propagating fractures in the experiments performed at nominallydry conditions. In the water-saturated sample, in contrast, although thevolume of coalescing fractures increases during the stage preceding failure,the volume of propagating fractures remains dominant. The influence ofstress corrosion cracking associated with hydration reactions at fracturetips and/or dilatant hardening may explain the observed difference infracture development under dry and water-saturated conditions. 
    more » « less