Search

Geologic CO$_2$ storage is an important strategy for reducing greenhouse gas emissions to the atmosphere and mitigating climate change. In this process, coupling between mechanical deformation and fluid flow in fault zones is a key determinant of fault instability, induced seismicity, and CO$_2$ leakage. Using a recently developed methodology, PREDICT, we obtain probability distributions of the permeability tensor in faults from the stochastic placement of clay smears that accounts for geologic uncertainty. We build a comprehensive set of fault permeability scenarios from PREDICT and investigate the effects of uncertainties from the fault zone internal structure and composition on forecasts of CO$_2$ permanence and fault stability. To tackle the prohibitively expensive computational cost of the large number of simulations required to quantify uncertainty, we develop a deep-learning-based surrogate model capable of predicting flow migration, pressure buildup, and geomechanical responses in CO$_2$ storage operations. We also compare our probabilistic estimation of CO$_2$ leakage and fault instability with previous studies based on deterministic estimates of fault permeability. The results highlight the importance of including uncertainty and anisotropy in modeling of complex fault structures and improved management of geologic CO$_2$ storage projects.

Inspired by presentations by the authors during a workshop organized at the Princeton Center for Theoretical Science (PCTS) in January 2022, we present a perspective on some of the outstanding questions related to the "physics of the ground beneath our feet." These identified challenges are intrinsically shared with the field of Soft Matter but also have unique aspects when the natural environment is studied., Comment: Perspective Paper, 30 pages, 15 figures

A Darcy-Cahn-Hilliard model coupled with damage is developed to describe multiphase-flow and fluid-driven fracturing in porous media. The model is motivated by recent experimental observations in Hele-Shaw cells of the fluid-driven fracturing of a synthetic porous medium with tunable fracture resistance. The model is derived from continuum thermodynamics and employs several simplifying assumptions, such as linear poroelasticity and viscous-dominated flow. Two distinct phase fields are used to regularize the interface between an invading and a defending fluid, as well as the ensuing damage. The damage model is a cohesive version of a phase-field model for fracture, in which model parameters allow for control over both nucleation and crack growth. Model-based simulations with finite elements are then performed to calibrate the model against recent experimental results. In particular, an experimentally-inferred phase diagram differentiating two flow regimes of porous invasion and fracturing is recovered. Finally, the model is employed to explore the parameter space beyond experimental capabilities, giving rise to the construction of an expanded phase diagram that suggests a new flow regime.

We study fluid-induced deformation and fracture of cohesive granular media, and apply photoporomechanics to uncover the underpinning grain-scale mechanics. We fabricate photoelastic spherical particles of diameter d=2mm, and make a monolayer granular pack with tunable intergranular cohesion in a circular Hele-Shaw cell that is initially filled with viscous silicone oil. We inject water into the oil-filled photoelastic granular pack, varying the injection flow rate, defending-fluid viscosity, and intergranular cohesion. We find two different modes of fluid invasion: viscous fingering, and fracturing with leak-off of the injection fluid. We directly visualize the evolving effective stress field through the particles' photoelastic response, and discover a hoop effective stress region behind the water invasion front, where we observe tensile force chains in the circumferential direction. Outside the invasion front, we observe compressive force chains aligning in the radial direction. We conceptualize the system's behavior by means of a two-phase poroelastic continuum model. The model captures granular pack dilation and compaction with the boundary delineated by the invasion front, which explains the observed distinct alignments of the force chains. Finally, we rationalize the crossover from viscous fingering to fracturing by comparing the competing forces behind the process: viscous force from fluid injection that drives fractures, and intergranular cohesion and friction that resist fractures.

Successful deployment of geological carbon storage (GCS) requires an extensive use of reservoir simulators for screening, ranking and optimization of storage sites. However, the time scales of GCS are such that no sufficient long-term data is available yet to validate the simulators against. As a consequence, there is currently no solid basis for assessing the quality with which the dynamics of large-scale GCS operations can be forecasted. To meet this knowledge gap, we have conducted a major GCS validation benchmark study. To achieve reasonable time scales, a laboratory-size geological storage formation was constructed (the "FluidFlower"), forming the basis for both the experimental and computational work. A validation experiment consisting of repeated GCS operations was conducted in the FluidFlower, providing what we define as the true physical dynamics for this system. Nine different research groups from around the world provided forecasts, both individually and collaboratively, based on a detailed physical and petrophysical characterization of the FluidFlower sands. The major contribution of this paper is a report and discussion of the results of the validation benchmark study, complemented by a description of the benchmarking process and the participating computational models. The forecasts from the participating groups are compared to each other and to the experimental data by means of various indicative qualitative and quantitative measures. By this, we provide a detailed assessment of the capabilities of reservoir simulators and their users to capture both the injection and post-injection dynamics of the GCS operations.

We study the characteristics of fluid-fluid displacement in simple mixed-wet porous micromodels numerically using a dynamic pore network model. The porous micromodel consists of distinct water-wet and oil-wet regions, whose fractions are systematically varied to yield a variety of displacement patterns over a wide range of capillary numbers. We find that the impact of mixed-wettability is most prominent at low capillary number, and it depends on the complex interplay between wettability fraction and the intrinsic contact angle of the water-wet regions. For example, the fractal dimension of the displacement pattern is a monotonically increasing function of wettability fraction in flow cells with strongly water-wet clusters, but it becomes non-monotonic with respect to wettability fraction in flow cells with weakly water-wet clusters. Additionally, mixed-wettability also manifests itself in the injection-pressure signature, which exhibits fluctuations especially at low wettability fraction. Specifically, preferential filling of water-wet regions leads to reduced effective permeability and higher injection pressure, even at vanishingly small capillary numbers. Finally, we demonstrate that scaling analyses based on a weighted average description of the overall wetting state of the mixed-wet system can effectively capture the variations in observed displacement pattern morphology.

The accuracy and robustness of numerical models of geologic CO$_2$ sequestration are almost never quantified with respect to direct observations that provide a ground truth. Here, we conduct CO$_2$ injection experiments in meter-scale, quasi-2D tanks with porous media representing stratigraphic sections of the subsurface, and combine them with numerical simulations of those experiments. We evaluate (1) the value of prior knowledge of the system, expressed in terms of ex-situ measurements of the tank sands' multiphase flow properties (local data), to obtain an accurate simulation; and (2) the predictive capability of the matched numerical models, when applied to different settings. We match three different simulation models-each with access to an increasing level of local data-to a CO$_2$ injection experiment in tank 1 (89.7$\times$47$\times$1.05 cm). Matching is based on a quantitative comparison of CO$_2$ migration at different times from timelapse image analysis. Next, we simulate a different injection scenario in tank 1, and, finally, a different injection scenario in tank 2 (2.86$\times$1.3$\times$0.019 m), which represents an altogether different stratigraphic section. Our models can qualitatively match the CO$_2$ plume migration and convective mixing of the experimental truth. Quantitatively, simulations are accurate during the injection phase but their performance decreases with time. Using local data reduces the time required to history match. The predictive capability of matched models, however, is found to be similar. The sand-water-CO$_{2(\mathrm{g})}$ system is very sensitive to effective permeability and capillary pressure changes; where heterogeneous structures are present, accurate deterministic estimates of CO$_2$ migration are unlikely.

Carbon, capture, and storage (CCS) is an important bridging technology to combat climate change in the transition towards net-zero. The FluidFlower concept has been developed to visualize and study CO$_2$ flow and storage mechanisms in sedimentary systems in a laboratory setting. Meter-scale multiphase flow in two geological geometries, including normal faults with and without smearing, is studied. The experimental protocols developed to provide key input parameters for numerical simulations are detailed, including an evaluation of operational parameters for the FluidFlower benchmark study. Variability in CO$_2$ migration patterns for two different geometries is quantified, both between 16 repeated laboratory runs and between history-matched models and a CO$_2$ injection experiment. The predicative capability of a history-matched model is then evaluated in a different geological setting., Comment: manuscript + supplementary material

We propose a solution strategy for parameter identification in multiphase thermo-hydro-mechanical (THM) processes in porous media using physics-informed neural networks (PINNs). We employ a dimensionless form of the THM governing equations that is particularly well suited for the inverse problem, and we leverage the sequential multiphysics PINN solver we developed in previous work. We validate the proposed inverse-modeling approach on multiple benchmark problems, including Terzaghi's isothermal consolidation problem, Barry-Mercer's isothermal injection-production problem, and nonisothermal consolidation of an unsaturated soil layer. We report the excellent performance of the proposed sequential PINN-THM inverse solver, thus paving the way for the application of PINNs to inverse modeling of complex nonlinear multiphysics problems.

It is well-known that wettability exerts fundamental control over multiphase flow in porous media, which has been extensively studied in uniform-wet porous media. In contrast, multiphase flow in porous media with heterogeneous wettability (i.e., mixed-wet) is less well-understood, despite its common occurrence. Here, we study the displacement of silicone oil by water in a mostly oil-wet porous media patterned with discrete water-wet clusters that have precisely controlled wettability. Surprisingly, the macroscopic displacement pattern varies dramatically depending on the details of wettability alteration -- the invading water preferentially fills strongly water-wet clusters but encircles weakly water-wet clusters instead, resulting in significant trapping of the defending oil. We explain this counter-intuitive observation with pore-scale simulations, which reveal that the fluid-fluid interfaces at mixed-wet pores resemble an S-shaped saddle with mean curvatures close to zero. We show that incorporation of the capillary entry pressures at mixed-wet pores into a dynamic pore-network model reproduces the experiments. Our work demonstrates the complex nature of wettability control in mixed-wet porous media, and it presents experimental and numerical platforms upon which further insights can be drawn.

Physics-Informed Neural Networks (PINNs) have received increased interest for forward, inverse, and surrogate modeling of problems described by partial differential equations (PDE). However, their application to multiphysics problem, governed by several coupled PDEs, present unique challenges that have hindered the robustness and widespread applicability of this approach. Here we investigate the application of PINNs to the forward solution of problems involving thermo-hydro-mechanical (THM) processes in porous media, which exhibit disparate spatial and temporal scales in thermal conductivity, hydraulic permeability, and elasticity. In addition, PINNs are faced with the challenges of the multi-objective and non-convex nature of the optimization problem. To address these fundamental issues, we: (1)~rewrite the THM governing equations in dimensionless form that is best suited for deep-learning algorithms; (2)~propose a sequential training strategy that circumvents the need for a simultaneous solution of the multiphysics problem and facilitates the task of optimizers in the solution search; and (3)~leverage adaptive weight strategies to overcome the stiffness in the gradient flow of the multi-objective optimization problem. Finally, we apply this framework to the solution of several synthetic problems in 1D and~2D., Comment: 21 pages, 7 figures

Physics-informed neural networks (PINNs) have received significant attention as a unified framework for forward, inverse, and surrogate modeling of problems governed by partial differential equations (PDEs). Training PINNs for forward problems, however, pose significant challenges, mainly because of the complex non-convex and multi-objective loss function. In this work, we present a PINN approach to solving the equations of coupled flow and deformation in porous media for both single-phase and multiphase flow. To this end, we construct the solution space using multi-layer neural networks. Due to the dynamics of the problem, we find that incorporating multiple differential relations into the loss function results in an unstable optimization problem, meaning that sometimes it converges to the trivial null solution, other times it moves very far from the expected solution. We report a dimensionless form of the coupled governing equations that we find most favourable to the optimizer. Additionally, we propose a sequential training approach based on the stress-split algorithms of poromechanics. Notably, we find that sequential training based on stress-split performs well for different problems, while the classical strain-split algorithm shows an unstable behaviour similar to what is reported in the context of finite element solvers. We use the approach to solve benchmark problems of poroelasticity, including Mandel's consolidation problem, Barry-Mercer's injection-production problem, and a reference two-phase drainage problem. The Python-SciANN codes reproducing the results reported in this manuscript will be made publicly available at https://github.com/sciann/sciann-applications.

Deep learning has been the most popular machine learning method in the last few years. In this chapter, we present the application of deep learning and physics-informed neural networks concerning structural mechanics and vibration problems. Demonstration problems involve de-noising data, solution to time-dependent ordinary and partial differential equations, and characterizing the system's response for a given data.

The Physics-Informed Neural Network (PINN) framework introduced recently incorporates physics into deep learning, and offers a promising avenue for the solution of partial differential equations (PDEs) as well as identification of the equation parameters. The performance of existing PINN approaches, however, may degrade in the presence of sharp gradients, as a result of the inability of the network to capture the solution behavior globally. We posit that this shortcoming may be remedied by introducing long-range (nonlocal) interactions into the network's input, in addition to the short-range (local) space and time variables. Following this ansatz, here we develop a nonlocal PINN approach using the Peridynamic Differential Operator (PDDO)---a numerical method which incorporates long-range interactions and removes spatial derivatives in the governing equations. Because the PDDO functions can be readily incorporated in the neural network architecture, the nonlocality does not degrade the performance of modern deep-learning algorithms. We apply nonlocal PDDO-PINN to the solution and identification of material parameters in solid mechanics and, specifically, to elastoplastic deformation in a domain subjected to indentation by a rigid punch, for which the mixed displacement--traction boundary condition leads to localized deformation and sharp gradients in the solution. We document the superior behavior of nonlocal PINN with respect to local PINN in both solution accuracy and parameter inference, illustrating its potential for simulation and discovery of partial differential equations whose solution develops sharp gradients.

In this paper, we introduce SciANN, a Python package for scientific computing and physics-informed deep learning using artificial neural networks. SciANN uses the widely used deep-learning packages Tensorflow and Keras to build deep neural networks and optimization models, thus inheriting many of Keras's functionalities, such as batch optimization and model reuse for transfer learning. SciANN is designed to abstract neural network construction for scientific computations and solution and discovery of partial differential equations (PDE) using the physics-informed neural networks (PINN) architecture, therefore providing the flexibility to set up complex functional forms. We illustrate, in a series of examples, how the framework can be used for curve fitting on discrete data, and for solution and discovery of PDEs in strong and weak forms. We summarize the features currently available in SciANN, and also outline ongoing and future developments.

We present the application of a class of deep learning, known as Physics Informed Neural Networks (PINN), to learning and discovery in solid mechanics. We explain how to incorporate the momentum balance and constitutive relations into PINN, and explore in detail the application to linear elasticity, and illustrate its extension to nonlinear problems through an example that showcases von~Mises elastoplasticity. While common PINN algorithms are based on training one deep neural network (DNN), we propose a multi-network model that results in more accurate representation of the field variables. To validate the model, we test the framework on synthetic data generated from analytical and numerical reference solutions. We study convergence of the PINN model, and show that Isogeometric Analysis (IGA) results in superior accuracy and convergence characteristics compared with classic low-order Finite Element Method (FEM). We also show the applicability of the framework for transfer learning, and find vastly accelerated convergence during network re-training. Finally, we find that honoring the physics leads to improved robustness: when trained only on a few parameters, we find that the PINN model can accurately predict the solution for a wide range of parameters new to the network---thus pointing to an important application of this framework to sensitivity analysis and surrogate modeling., Comment: 19 pages, 14 figures

We develop a novel `moving capacitor' dynamic network model to simulate immiscible fluid-fluid displacement in porous media. Traditional network models approximate the pore geometry as a network of fixed resistors, directly analogous to an electrical circuit. Our model additionally captures the motion of individual fluid-fluid interfaces through the pore geometry by completing this analogy, representing interfaces as a set of moving capacitors. By incorporating pore-scale invasion events, the model reproduces, for the first time, both the displacement pattern and the injection pressure signal under a wide range of capillary numbers and substrate wettabilities. We show that at high capillary numbers the invading patterns advance symmetrically through viscous fingers. In contrast, at low capillary numbers the flow is governed by the wettability-dependent fluid-fluid interactions with the pore structure. The signature of the transition between the two regimes manifests itself in the fluctuations of the injection pressure signal.

We investigate tracer transport on random discrete fracture networks that are characterized by the statistics of the fracture geometry and hydraulic conductivity. While it is well known that tracer transport through fractured media can be anomalous and particle injection modes can have major impact on dispersion, the incorporation of injection modes into effective transport modelling has remained an open issue. The fundamental reason behind this challenge is that---even if the Eulerian fluid velocity is steady---the Lagrangian velocity distribution experienced by tracer particles evolves with time from its initial distribution, which is dictated by the injection mode, to a stationary velocity distribution. We quantify this evolution by a Markov model for particle velocities that are equidistantly sampled along trajectories. This stochastic approach allows for the systematic incorporation of the initial velocity distribution and quantifies the interplay between velocity distribution and spatial and temporal correlation. The proposed spatial Markov model is characterized by the initial velocity distribution, which is determined by the particle injection mode, the stationary Lagrangian velocity distribution, which is derived from the Eulerian velocity distribution, and the spatial velocity correlation length, which is related to the characteristic fracture length. This effective model leads to a time-domain random walk for the evolution of particle positions and velocities, whose joint distribution follows a Boltzmann equation. Finally, we demonstrate that the proposed model can successfully predict anomalous transport through discrete fracture networks with different levels of heterogeneity and arbitrary tracer injection modes.

Self-organization and pattern formation in network-organized systems emerges from the collective activation and interaction of many interconnected units. A striking feature of these non-equilibrium structures is that they are often localized and robust: only a small subset of the nodes, or cell assembly, is activated. Understanding the role of cell assemblies as basic functional units in neural networks and socio-technical systems emerges as a fundamental challenge in network theory. A key open question is how these elementary building blocks emerge, and how they operate, linking structure and function in complex networks. Here we show that a network analogue of the Swift-Hohenberg continuum model---a minimal-ingredients model of nodal activation and interaction within a complex network---is able to produce a complex suite of localized patterns. Hence, the spontaneous formation of robust operational cell assemblies in complex networks can be explained as the result of self-organization, even in the absence of synaptic reinforcements. Our results show that these self-organized, local structures can provide robust functional units to understand natural and socio-technical network-organized processes., Comment: 11 pages, 4 figures