The coupling between subsurface flow and geomechanical deformation is critical in the assessment of the environmental impacts of groundwater use, underground liquid waste disposal, geologic storage of carbon dioxide, and exploitation of shale gas reserves. In particular, seismicity induced by fluid injection and withdrawal has emerged as a central element of the scientific discussion around subsurface technologies that tap into water and energy resources. Here we present a new computational approach to model coupled multiphase flow and geomechanics of faulted reservoirs. We represent faults as surfaces embedded in a three-dimensional medium by using zero-thickness interface elements to accurately model fault slip under dynamically evolving fluid pressure and fault strength. We incorporate the effect of fluid pressures from multiphase flow in the mechanical stability of faults and employ a rigorous formulation of nonlinear multiphase geomechanics that is capable of handling strong capillary effects. We develop a numerical simulation tool by coupling a multiphase flow simulator with a mechanics simulator, using the unconditionally stable fixed-stress scheme for the sequential solution of two-way coupling between flow and geomechanics. We validate our modeling approach using several synthetic, but realistic, test cases that illustrate the onset and evolution of earthquakes from fluid injection and withdrawal.

Viscous fingering is a well-known hydrodynamic instability that sets in when a less viscous fluid displaces a more viscous fluid [1]. When the two fluids are miscible, viscous fingering introduces disorder in the velocity field and exerts a fundamental control on the rate at which the fluids mix. We present a fluid dynamics video of the mixing process in a viscously unstable flow, generated from a high-resolution numerical simulation using a computational strategy that is stable for arbitrary viscosity ratios. We develop a two-equation dynamic model of concentration variance and mean dissipation rate to quantify the degree of mixing in such a displacement process [2]. The model reproduces accurately the evolution of these two quantities as observed in highresolution numerical simulations and captures the nontrivial interplay between channeling and creation of interfacial area as a result of viscous fingering.

In this paper, we present a computational framework for the simulation of coupled flow and reservoir geomechanics. The physical model is restricted to Biot's theory of single-phase flow and linear poroelasticity, but is sufficiently general to be extended to multiphase flow problems and inelastic behavior. The distinctive technical aspects of our approach are: (1) the space discretization of the equations. The unknown variables are the pressure, the fluid velocity, and the rock displacements. We recognize that these variables are of very different nature, and need to be discretized differently. We propose a mixed finite element space discretization, which is stable, convergent, locally mass conservative, and employs a single computational grid. To ensure stability and robustness, we perform an implicit time integration of the fluid flow equations. (2) The strategies for the solution of the coupled system. We compare different solution strategies, including the fully coupled approach, the usual (conditionally stable) iteratively coupled approach, and a less common unconditionally stable sequential scheme. We show that the latter scheme corresponds to a modified block Jacobi method, which also enjoys improved convergence properties. This computational model has been implemented in an objectoriented reservoir simulator, whose modular design allows for further extensions and enhancements. We show several representative numerical simulations that illustrate the effectiveness of the approach.

Viscous fingering is a well-known hydrodynamic instability that sets in when a less viscous fluid displaces a more viscous fluid. When the two fluids are miscible, viscous fingering introduces disorder in the velocity field and exerts a fundamental control on the rate at which the fluids mix. Here we analyze the characteristic signature of the mixing process in viscously unstable flows, by means of high-resolution numerical simulations using a computational strategy that is stable for arbitrary viscosity ratios. We propose a reduced-order model of mixing, which, in the spirit of turbulence modeling and in contrast with previous approaches, recognizes the fundamental role played by the mechanical dissipation rate. The proposed model captures the nontrivial interplay between channeling and creation of interfacial area as a result of viscous fingering.

Fluid mixing plays a fundamental role in many natural and engineered processes, including groundwater flows in porous media, enhanced oil recovery, and microfluidic lab-on-a-chip systems. Recent developments have explored the effect of viscosity contrast on mixing, suggesting that the unstable displacement of fluids with different viscosities, or viscous fingering, provides a powerful mechanism to increase fluid-fluid interfacial area and enhance mixing. However, existing studies have not incorporated the effect of medium heterogeneity on the mixing rate. Here, we characterize the evolution of mixing between two fluids of different viscosity in heterogeneous porous media. We focus on a practical scenario of divergent-convergent flow in a quarter five spot geometry prototypical of well-driven groundwater flows. We study by means of numerical simulations the impact of permeability heterogeneity and viscosity contrast on the breakthrough curves and mixing efficiency, and we rationalize the nontrivial mixing behavior that emerges from the competition between the creation of fluid-fluid interfacial area and channeling.

Seismicity induced by fluid injection and withdrawal has emerged as a central element of the scientific discussion around subsurface technologies that tap into water and energy resources. Here we present the application of coupled flow‐geomechanics simulation technology to the post mortem analysis of a sequence of damaging earthquakes (Mw = 6.0 and 5.8) in May 2012 near the Cavone oil field, in northern Italy. This sequence raised the question of whether these earthquakes might have been triggered by activities due to oil and gas production. Our analysis strongly suggests that the combined effects of fluid production and injection from the Cavone field were not a driver for the observed seismicity. More generally, our study illustrates that computational modeling of coupled flow and geomechanics permits the integration of geologic, seismotectonic, well log, fluid pressure and flow rate, and geodetic data and provides a promising approach for assessing and managing hazards associated with induced seismicity.

We study mixing of two fluids of different viscosity in a microfluidic channel or porous medium. We show that the synergetic action of alternating injection and viscous fingering leads to a dramatic increase in mixing efficiency at high Péclet numbers. Based on observations from high-resolution simulations, we develop a theoretical model of mixing efficiency that combines a hyperbolic mixing model of the channelized region ahead and a mixing-dissipation model of the pseudosteady region behind. Our macroscopic model quantitatively reproduces the evolution of the average degree of mixing along the flow direction and can be used as a design tool to optimize mixing from viscous fingering in a microfluidic channel.

The interaction between seasonally-induced non-tectonic and tectonic deformation along the Himalayan plate boundary remains debated. Here, we propose that tectonic deformation along this plate boundary can be significantly influenced by the deformation induced by the non-tectonic hydrological loading cycles. We explore seasonal mass oscillations by continental water storage in Southeast Asia and Himalayan arc region using continuous Global Positioning System measurements and satellite data from the Gravity Recovery and Climate Experiment. We suggest that the substantially higher transient displacements above the base of the seismogenic zone indicate a role of changes in aseismic slip rate on the deep megathrust that may be controlled by seasonal hydrological loading. We invoke modulation of aseismic slip on the megathrust down-dip of the seismogenic zone due to a fault resonance process induced by the seasonal stress changes. This process modulates mid-crustal ramp associated micro-seismicity and influences the timing of Central Himalayan earthquakes.