Parameter-robust discretization and preconditioning of Biot's consolidation model

Biot's consolidation model in poroelasticity has a number of applications in science, medicine, and engineering. The model depends on various parameters, and in practical applications these parameters ranges over several orders of magnitude. A current challenge is to design discretization techniques and solution algorithms that are well behaved with respect to these variations. The purpose of this paper is to study finite element discretizations of this model and construct block diagonal preconditioners for the discrete Biot systems. The approach taken here is to consider the stability of the problem in non-standard or weighted Hilbert spaces and employ the operator preconditioning approach. We derive preconditioners that are robust with respect to both the variations of the parameters and the mesh refinement. The parameters of interest are small time-step sizes, large bulk and shear moduli, and small hydraulic conductivity.Comment: 24 page

Robust preconditioners for a new stabilized discretization of the poroelastic equations

In this paper, we present block preconditioners for a stabilized discretization of the poroelastic equations developed in [C. Rodrigo, X. Hu, P. Ohm, J. Adler, F. Gaspar, and L. Zikatanov, Comput. Methods Appl. Mech. Engrg., 341 (2018), pp. 467-484]. The discretization is proved to be well-posed with respect to the physical and discretization parameters and thus provides a framework to develop preconditioners that are robust with respect to such parameters as well. We construct both norm-equivalent (diagonal) and field-of-value-equivalent (triangular) preconditioners for both the stabilized discretization and a perturbation of the stabilized discretization, which leads to a smaller overall problem after static condensation. Numerical tests for both two-and three-dimensional problems confirm the robustness of the block preconditioners with respect to the physical and discretization parameters

A Study of Interwell Interference and Well Performance in Unconventional Reservoirs Based on Coupled Flow and Geomechanics Modeling with Improved Computational Efficiency

Completion quality of tightly spaced horizontal wells in unconventional reservoirs is important for hydrocarbon recovery efficiency. Parent well production usually leads to heterogeneous stress evolution around parent wells and at infill well locations, which affects hydraulic fracture growth along infill wells. Recent field observations indicate that infill well completions lead to frac hits and production interference between parent and infill wells. Therefore, it is important to characterize the heterogeneous interwell stress/pressure evolutions and hydraulic fracture networks. This work presents a reservoir-geomechanics-fracturing modeling workflow and its implementation in unconventional reservoirs for the characterization of interwell stress and pressure evolutions and for the modeling of interwell hydraulic fracture geometry. An in-house finite element model coupling fluid flow and geomechanics is first introduced and used to characterize production-induced stress and pressure changes in the reservoir. Then, an in-house complex fracture propagation model coupling fracture mechanics and wellbore/fracture fluid flow is used for the simulation of hydraulic fractures along infill wells. A parallel solver is also implemented in a reservoir geomechanics simulator in a separate study to investigate the potential of improving computational efficiency. Results show that differential stress (DS), parent well fracture geometry, legacy production time, bottomhole pressure (BHP) for legacy production, and perforation cluster location are key parameters affecting interwell fracture geometry and the occurrence of frac hits. In general, transverse infill well fractures are obtained in scenarios with large DS and small legacy producing time/BHP. Non-uniform parent well fracture geometry leads to frac hits in certain cases, while the assumption of uniform parent well fracture half-lengths in the numerical model could not capture the phenomenon of frac hits. Perforation cluster locations along infill wells do not play an important role in determining whether an infill well hydraulic fracture is transverse, while they are important for the occurrence of frac hits. In addition, the implementation of a parallel solver, PETSc, in a fortran-based simulator indicates that an overall speedup of 14 can be achieved for simulations with one million grid blocks. This result provides a reference for improving computational efficiency for geomechanical simulation involving large matrices using finite element methods (FEM)

Iterative Solutions of Large-Scale Soil-Structure Problems

Ph.DDOCTOR OF PHILOSOPH

Preconditioners for Soil-Structure Interaction Problems with Significant Material Stiffness Contrast

Ph.DDOCTOR OF PHILOSOPH

A Study of Interwell Interference and Well Performance in Unconventional Reservoirs Based on Coupled Flow and Geomechanics Modeling with Improved Computational Efficiency

Completion quality of tightly spaced horizontal wells in unconventional reservoirs is important for hydrocarbon recovery efficiency. Parent well production usually leads to heterogeneous stress evolution around parent wells and at infill well locations, which affects hydraulic fracture growth along infill wells. Recent field observations indicate that infill well completions lead to frac hits and production interference between parent and infill wells. Therefore, it is important to characterize the heterogeneous interwell stress/pressure evolutions and hydraulic fracture networks. This work presents a reservoir-geomechanics-fracturing modeling workflow and its implementation in unconventional reservoirs for the characterization of interwell stress and pressure evolutions and for the modeling of interwell hydraulic fracture geometry. An in-house finite element model coupling fluid flow and geomechanics is first introduced and used to characterize production-induced stress and pressure changes in the reservoir. Then, an in-house complex fracture propagation model coupling fracture mechanics and wellbore/fracture fluid flow is used for the simulation of hydraulic fractures along infill wells. A parallel solver is also implemented in a reservoir geomechanics simulator in a separate study to investigate the potential of improving computational efficiency. Results show that differential stress (DS), parent well fracture geometry, legacy production time, bottomhole pressure (BHP) for legacy production, and perforation cluster location are key parameters affecting interwell fracture geometry and the occurrence of frac hits. In general, transverse infill well fractures are obtained in scenarios with large DS and small legacy producing time/BHP. Non-uniform parent well fracture geometry leads to frac hits in certain cases, while the assumption of uniform parent well fracture half-lengths in the numerical model could not capture the phenomenon of frac hits. Perforation cluster locations along infill wells do not play an important role in determining whether an infill well hydraulic fracture is transverse, while they are important for the occurrence of frac hits. In addition, the implementation of a parallel solver, PETSc, in a fortran-based simulator indicates that an overall speedup of 14 can be achieved for simulations with one million grid blocks. This result provides a reference for improving computational efficiency for geomechanical simulation involving large matrices using finite element methods (FEM)

Domain Decomposition And Time-Splitting Methods For The Biot System Of Poroelasticity

In this thesis, we develop efficient mixed finite element methods to solve the Biot system of poroelasticity, which models the flow of a viscous fluid through a porous medium along with the deformation of the medium. We study non-overlapping domain decomposition techniques and sequential splitting methods to reduce the computational complexity of the problem. The solid deformation is modeled with a mixed three-field formulation with weak stress symmetry. The fluid flow is modeled with a mixed Darcy formulation. We introduce displacement and pressure Lagrange multipliers on the subdomain interfaces to impose weakly the continuity of normal stress and normal velocity, respectively. The global problem is reduced to an interface problem for the Lagrange multipliers, which is solved by a Krylov space iterative method. We study both monolithic and split methods. For the monolithic method, the cases of matching and non-matching subdomain grid interfaces are analyzed separately. For both cases, a coupled displacement-pressure interface problem is solved, with each iteration requiring the solution of local Biot problems. For the case of matching subdomain grids, we show that the resulting interface operator is positive definite and analyze the convergence of the iteration. For the non-matching subdomain grid case, we use a multiscale mortar mixed finite element (MMMFE) approach. We further study drained split and fixed stress Biot splittings, in which case we solve separate interface problems requiring elasticity and Darcy solves. We analyze the stability of the split formulations. We also use numerical experiments to illustrate the convergence of the domain decomposition methods and compare their accuracy and efficiency in the monolithic and time-splitting settings. Finally, we present a novel space-time domain decomposition technique for the mixed finite element formulation of a parabolic equation. This method is motivated by the MMMFE method, where we split the space-time domain into multiple subdomains with space-time grids of different sizes. Scalar Lagrange multiplier (mortar) functions are introduced to enforce weakly the continuity of the normal component of the mixed finite element flux variable over the space-time interfaces. We analyze the new method and numerical experiments are developed to illustrate and confirm the theoretical results

MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications

Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described

Generalized averaged Gaussian quadrature and applications

A simple numerical method for constructing the optimal generalized averaged Gaussian quadrature formulas will be presented. These formulas exist in many cases in which real positive GaussKronrod formulas do not exist, and can be used as an adequate alternative in order to estimate the error of a Gaussian rule. We also investigate the conditions under which the optimal averaged Gaussian quadrature formulas and their truncated variants are internal