Stress and permeability variations around a wellbore and in the reservoir are of much interest in petroleum and geothermal reservoir development. Water injection causes significant changes in pore pressure, temperature, and stress in hot reservoirs, changing rock permeability. In this work, two- and three-dimensional finite element methods were developed to simulate coupled reservoirs with damage mechanics and stress-dependent permeability. The model considers the influence of fluid flow, temperature, and solute transport in rock deformation and models nonlinear behavior with continuum damage mechanics and stress-dependent permeability. Numerical modeling was applied to analyze wellbore stability in swelling shale with two- and three-dimensional damage/fracture propagation around a wellbore and injection-induced microseismic events. The finite element method (FEM) was used to solve the displacement, pore pressure, temperature, and solute concentration problems. Solute mass transport between drilling fluid and shale formation was considered to study salinity effects. Results show that shear and tensile failure can occur around a wellbore in certain drilling conditions where the mud pressure lies between the reservoir pore pressure and fracture gradient. The fully coupled thermo-poro-mechanical FEM simulation was used to model damage/fracture propagation and microseismic events caused by fluid injection. These studies considered wellbore geometry in small-scale modeling and point-source injection, assuming singularity fluid flux for large-scale simulation. Damage mechanics was applied to capture the effects of crack initiation, microvoid growth, and fracture propagation. The induced microseismic events were modeled in heterogeneous geological media, assuming the Weibull distribution functions for modulus and permeability. The results of this study indicate that fluid injection causes the effective stress to relax in the damage phase and to concentrate at the interface between the damage phase and the intact rock. Furthermore, induced-stress and far-field stress influence damage propagation. Cold water injection causes the tensile stress and affects the initial fracture and fracture propagation, but fracture initiation pressure and far-field stress are critical to create a damage/fracture plane, which is normal to the minimum far-field stress direction following well stimulation. Microseismic events propagate at both well scale and reservoir-scale simulation; the cloud shape of a microseismic event is affected by permeability anisotropy and far-field stress, and deviatoric horizontal far-field stress especially contributes to the localization of the microseismic cloud.

In this research we consider developing a reservoir simulator capable of simulating complex coupled poromechanical processes on massively parallel computers. A variety of problems arising from petroleum and environmental engineering inherently necessitate the understanding of interactions between fluid flow and solid mechanics. Examples in petroleum engineering include reservoir compaction, wellbore collapse, sand production, and hydraulic fracturing. In environmental engineering, surface subsidence, carbon sequestration, and waste disposal are also coupled poromechanical processes. These economically and environmentally important problems motivate the active pursuit of robust, efficient, and accurate simulation tools for coupled poromechanical problems. Three coupling approaches are currently employed in the reservoir simulation community to solve the poromechanics system, namely, the fully implicit coupling (FIM), the explicit coupling, and the iterative coupling. The choice of the coupling scheme significantly affects the efficiency of the simulator and the accuracy of the solution. We adopt the fixed-stress iterative coupling scheme to solve the coupled system due to its advantages over the other two. Unlike the explicit coupling, the fixed-stress split has been theoretically proven to converge to the FIM for linear poroelasticity model. In addition, it is more efficient and easier to implement than the FIM. Our computational results indicate that this approach is also valid for multiphase flow. We discretize the quasi-static linear elasticity model for geomechanics in space using the continuous Galerkin (CG) finite element method (FEM) on general hexahedral grids. Fluid flow models are discretized by locally mass conservative schemes, specifically, the mixed finite element method (MFE) for the equation of state compositional flow on Cartesian grids and the multipoint flux mixed finite element method (MFMFE) for the single phase and two-phase flows on general hexahedral grids. While both the MFE and the MFMFE generate cell-centered stencils for pressure, the MFMFE has advantages in handling full tensor permeabilities and general geometry and boundary conditions. The MFMFE also obtains accurate fluxes at cell interfaces. These characteristics enable the simulation of more practical problems. For many reservoir simulation applications, for instance, the carbon sequestration simulation, we need to account for thermal effects on the compositional flow phase behavior and the solid structure stress evolution. We explicitly couple the poromechanics equations to a simplified energy conservation equation. A time-split scheme is used to solve heat convection and conduction successively. For the convection equation, a higher order Godunov method is employed to capture the sharp temperature front; for the conduction equation, the MFE is utilized. Simulations of coupled poromechanical or thermoporomechanical processes in field scales with high resolution usually require parallel computing capabilities. The flow models, the geomechanics model, and the thermodynamics model are modularized in the Integrated Parallel Accurate Reservoir Simulator (IPARS) which has been developed at the Center for Subsurface Modeling at the University of Texas at Austin. The IPARS framework handles structured (logically rectangular) grids and was originally designed for element-based data communication, such as the pressure data in the flow models. To parallelize the node-based geomechanics model, we enhance the capabilities of the IPARS framework for node-based data communication. Because the geomechanics linear system is more costly to solve than those of flow and thermodynamics models, the performance of linear solvers for the geomechanics model largely dictates the speed and scalability of the coupled simulator. We use the generalized minimal residual (GMRES) solver with the BoomerAMG preconditioner from the hypre library and the geometric multigrid (GMG) solver from the UG4 software toolbox to solve the geomechanics linear system. Additionally, the multilevel k-way mesh partitioning algorithm from METIS is used to generate high quality mesh partitioning to improve solver performance. Numerical examples of coupled poromechanics and thermoporomechanics simulations are presented to show the capabilities of the coupled simulator in solving practical problems accurately and efficiently. These examples include a real carbon sequestration field case with stress-dependent permeability, a synthetic thermoporoelastic reservoir simulation, poroelasticity simulations on highly distorted hexahedral grids, and parallel scalability tests on a massively parallel computer. / text