Showing papers in "arXiv: Geophysics in 2021"

U-FNO - an enhanced Fourier neural operator based-deep learning model for multiphase flow.

Abstract: Numerical simulation of multiphase flow in porous media is essential for many geoscience applications. Data-driven machine learning methods provide faster alternatives to traditional simulators by training neural network models with numerical simulation data mappings. Here we present U-FNO, a novel neural network architecture for solving multiphase flow problems with superior speed, accuracy, and data efficiency. U-FNO is designed based on the newly proposed Fourier neural operator (FNO) that learns an infinite-dimensional integral kernel in the Fourier space, which has shown excellent performance for single-phase flows. Here we extend the FNO-based architecture to a CO2-water multiphase problem, and proposes the U-FNO architecture to enhance the prediction accuracy in multiphase flow systems. Through a systematic comparison among a CNN benchmark and three types of FNO variations, we show that the U-FNO architecture has the advantages of both the traditional CNN and original FNO, providing significantly more accurate and efficient performance than previous architectures. The trained U-FNO predicts gas saturation and pressure buildup with a 6*10e4 times speed-up compared to traditional numerical simulators while maintaining similar accuracy. The trained models can act as a general-purpose simulator alternative for 2D-radial CO2 injection problems with wide ranges of permeability and porosity heterogeneity, anisotropy, reservoir conditions, injection configurations, flow rates, and multiphase flow properties.

Neural Network Attribution Methods for Problems in Geoscience: A Novel Synthetic Benchmark Dataset.

Abstract: Despite the increasingly successful application of neural networks to many problems in the geosciences, their complex and nonlinear structure makes the interpretation of their predictions difficult, which limits model trust and does not allow scientists to gain physical insights about the problem at hand. Many different methods have been introduced in the emerging field of eXplainable Artificial Intelligence (XAI), which aim at attributing the network's prediction to specific features in the input domain. XAI methods are usually assessed by using benchmark datasets (like MNIST or ImageNet for image classification), or through deletion/insertion techniques. In either case, however, an objective, theoretically-derived ground truth for the attribution is lacking, making the assessment of XAI in many cases subjective. Also, benchmark datasets for problems in geosciences are rare. Here, we provide a framework, based on the use of additively separable functions, to generate attribution benchmark datasets for regression problems for which the ground truth of the attribution is known a priori. We generate a long benchmark dataset and train a fully-connected network to learn the underlying function that was used for simulation. We then compare estimated attribution heatmaps from different XAI methods to the ground truth in order to identify examples where specific XAI methods perform well or poorly. We believe that attribution benchmarks as the ones introduced herein are of great importance for further application of neural networks in the geosciences, and for accurate implementation of XAI methods, which will increase model trust and assist in discovering new science.

Transfer learning: Improving neural network based prediction of earthquake ground shaking for an area with insufficient training data

Abstract: In a recent study (Jozinovic et al, 2020) we showed that convolutional neural networks (CNNs) applied to network seismic traces can be used for rapid prediction of earthquake peak ground motion intensity measures (IMs) at distant stations using only recordings from stations near the epicenter. The predictions are made without any previous knowledge concerning the earthquake location and magnitude. This approach differs from the standard procedure adopted by earthquake early warning systems (EEWSs) that rely on location and magnitude information. In the previous study, we used 10 s, raw, multistation waveforms for the 2016 earthquake sequence in central Italy for 915 events (CI dataset). The CI dataset has a large number of spatially concentrated earthquakes and a dense station network. In this work, we applied the CNN model to an area around the VIRGO gravitational waves observatory sited near Pisa, Italy. In our initial application of the technique, we used a dataset consisting of 266 earthquakes recorded by 39 stations. We found that the CNN model trained using this smaller dataset performed worse compared to the results presented in the original study by Jozinovic et al. (2020). To counter the lack of data, we adopted transfer learning (TL) using two approaches: first, by using a pre-trained model built on the CI dataset and, next, by using a pre-trained model built on a different (seismological) problem that has a larger dataset available for training. We show that the use of TL improves the results in terms of outliers, bias, and variability of the residuals between predicted and true IMs values. We also demonstrate that adding knowledge of station positions as an additional layer in the neural network improves the results. The possible use for EEW is demonstrated by the times for the warnings that would be received at the station PII.

An End-to-End Earthquake Detection Method for Joint Phase Picking and Association using Deep Learning

Abstract: Earthquake monitoring by seismic networks typically involves a workflow consisting of phase detection/picking, association, and location tasks. In recent years, the accuracy of these individual stages has been improved through the use of machine learning techniques. In this study, we introduce a new, end-to-end approach that improves overall earthquake detection accuracy by jointly optimizing each stage of the detection pipeline. We propose a neural network architecture for the task of multi-station processing of seismic waveforms recorded over a seismic network. This end-to-end architecture consists of three sub-networks: a backbone network that extracts features from raw waveforms, a phase picking sub-network that picks P- and S-wave arrivals based on these features, and an event detection sub-network that aggregates the features from multiple stations and detects earthquakes. We use these sub-networks in conjunction with a shift-and-stack module based on back-projection that introduces kinematic constraints on arrival times, allowing the model to generalize to different velocity models and to variable station geometry in seismic networks. We evaluate our proposed method on the STanford EArthquake Dataset (STEAD) and on the 2019 Ridgecrest, CA earthquake sequence. The results demonstrate that our end-to-end approach can effectively pick P- and S-wave arrivals and achieve earthquake detection accuracy rivaling that of other state-of-the-art approaches.

Earthquake Phase Association using a Bayesian Gaussian Mixture Model

Abstract: Earthquake phase association algorithms aggregate picked seismic phases from a network of seismometers into individual earthquakes and play an important role in earthquake monitoring. Dense seismic networks and improved phase picking methods produce massive earthquake phase data sets, particularly for earthquake swarms and aftershocks occurring closely in time and space, making phase association a challenging problem. We present a new association method, the Gaussian Mixture Model Association (GaMMA), that combines the Gaussian mixture model for phase measurements (both time and amplitude), with earthquake location, origin time, and magnitude estimation. We treat earthquake phase association as an unsupervised clustering problem in a probabilistic framework, where each earthquake corresponds to a cluster of P and S phases with hyperbolic moveout of arrival times and a decay of amplitude with distance. We use a multivariate Gaussian distribution to model the collection of phase picks for an event, the mean of which is given by the predicted arrival time and amplitude from the causative event. We carry out the pick assignment for each earthquake and determine earthquake parameters (i.e., earthquake location, origin time, and magnitude) under the maximum likelihood criterion using the Expectation-Maximization (EM) algorithm. The GaMMA method does not require the typical association steps of other algorithms, such as grid-search or supervised training. The results on both synthetic test and the 2019 Ridgecrest earthquake sequence show that GaMMA effectively associates phases from a temporally and spatially dense earthquake sequence while producing useful estimates of earthquake location and magnitude.

Reservoir Computing as a Tool for Climate Predictability Studies

Abstract: Reduced-order dynamical models play a central role in developing our understanding of predictability of climate irrespective of whether we are dealing with the actual climate system or surrogate climate-models. In this context, the Linear-Inverse-Modeling (LIM) approach, by capturing a few essential interactions between dynamical components of the full system, has proven valuable in providing insights into predictability of the full system. We demonstrate that Reservoir Computing (RC), a form of learning suitable for systems with chaotic dynamics, provides an alternative nonlinear approach that improves on the predictive skill of the LIM approach. We do this in the example setting of predicting sea-surface-temperature in the North Atlantic in the pre-industrial control simulation of a popular earth system model, the Community-Earth-System-Model so that we can compare the performance of the new RC based approach with the traditional LIM approach both when learning data is plentiful and when such data is more limited. The improved predictive skill of the RC approach over a wide range of conditions -- larger number of retained EOF coefficients, extending well into the limited data regime, etc. -- suggests that this machine-learning technique may have a use in climate predictability studies. While the possibility of developing a climate emulator -- the ability to continue the evolution of the system on the attractor long after failing to be able to track the reference trajectory -- is demonstrated in the Lorenz-63 system, it is suggested that further development of the RC approach may permit such uses of the new approach in more realistic predictability studies.

Fast Quasi-Geostrophic Magneto-Coriolis Modes in the Earth's core

Abstract: Fast changes of Earth's magnetic field could be explained by inviscid and diffusion-less quasi-geostrophic (QG) Magneto-Coriolis modes. We present a hybrid QG model with columnar flows and three-dimensional magnetic fields and find modes with periods of a few years at parameters relevant to Earth's core. For the simple poloidal magnetic field that we consider here they show a localization of kinetic and magnetic energy in the equatorial region. This concentration of energy near the equator and the high frequency make them a plausible mechanism to explain similar features observed in recent geomagnetic field observations. Our model potentially opens a way to probe the otherwise inaccessible magnetic field structure in the Earth's outer core.

On the scale dependence in the dynamics of frictional rupture: constant fracture energy versus size-dependent breakdown work.

Abstract: Potential energy stored during the inter-seismic period by tectonic loading around faults is released during earthquakes as radiated energy, heat and fracture energy. The latter is of first importance since it controls the nucleation, propagation and arrest of the seismic rupture. On one side, fracture energy estimated for natural earthquakes (breakdown work) shows a clear slip-dependence. On the other side, recent experimental studies highlighted that, fracture energy is a material property limited by an upper bound value corresponding to the fracture energy of the intact material independently of the size of the event. To reconcile these contradictory observations, we performed stick-slip experiments in a bi-axial shear configuration. We analyzed the fault weakening during frictional rupture by accessing to the near-fault stress-slip curve through strain gauge array. We first estimated fracture energy by comparing the measured strain with the theoretical predictions from Linear Elastic Fracture Mechanics and a Cohesive Zone Model. By comparing these values to the breakdown work obtained from the integration of the stress-slip curve, we show that, at the scale of our experiments, fault weakening is divided into two stages; the first one consistent with the estimated fracture energy, and a long-tailed weakening corresponding to a larger energy not localized at the rupture tip, increasing with slip. Through numerical simulations, we demonstrate that only the first weakening stage controls the rupture initiation and that the breakdown work induced by the long-tailed weakening can enhance slip during rupture propagation and allow the rupture to overcome stress heterogeneity along the fault. We conclude that the origin of the seismological estimates of breakdown work could be related to the energy dissipated in the long-tailed weakening rather than to the one dissipated near the tip.

Multiscale Data-driven Seismic Full-waveform Inversion with Field Data Study

Abstract: Seismic full-waveform inversion (FWI), which uses iterative methods to estimate high-resolution subsurface models from seismograms, is a powerful imaging technique in exploration geophysics. In recent years, the computational cost of FWI has grown exponentially due to the increasing size and resolution of seismic data. Moreover, it is a non-convex problem and can encounter local minima due to the limited accuracy of the initial velocity models or the absence of low frequencies in the measurements. To overcome these computational issues, we develop a multiscale data-driven FWI method based on fully convolutional networks (FCN). In preparing the training data, we first develop a real-time style transform method to create a large set of synthetic subsurface velocity models from natural images. We then develop two convolutional neural networks with encoder-decoder structure to reconstruct the low- and high-frequency components of the subsurface velocity models, separately. To validate the performance of our data-driven inversion method and the effectiveness of the synthesized training set, we compare it with conventional physics-based waveform inversion approaches using both synthetic and field data. These numerical results demonstrate that, once our model is fully trained, it can significantly reduce the computation time, and yield more accurate subsurface velocity models in comparison with conventional FWI.

Deep learning for low-magnitude earthquake detection on a multi-level sensor network

Abstract: Automatic detection of low-magnitude earthquakes has become an increasingly important research topic in recent years due to a sharp increase in induced seismicity around the globe. The detection of low-magnitude seismic events is essential for microseismic monitoring of hydraulic fracturing, carbon capture and storage, and geothermal operations for hazard detection and mitigation. Moreover, the detection of micro-earthquakes is crucial to understand the underlying mechanisms of larger earthquakes. Various algorithms, including deep learning methods, have been proposed over the years to detect such low-magnitude events. However, there is still a need for improving the robustness of these methods in discriminating between local sources of noise and weak seismic events. In this study, we propose a convolutional neural network (CNN) to detect seismic events from shallow borehole stations in Groningen, the Netherlands. We train a CNN model to detect low-magnitude earthquakes, harnessing the multi-level sensor configuration of the G-network in Groningen. Each G-network station contains four geophones at depths of 50, 100, 150, and 200 meters. Unlike prior deep learning approaches that use 3-component seismic records only at a single sensor level, we use records from the entire borehole as one training example. This allows us to train the CNN model using moveout patterns of the energy traveling across the borehole sensors to discriminate between events originating in the subsurface and local noise arriving from the surface. We compare the prediction accuracy of our trained CNN model to that of the STA/LTA and template matching algorithms on a two-month continuous record. We demonstrate that the CNN model shows significantly better performance than STA/LTA and template matching in detecting new events missing from the catalog and minimizing false detections.

The Landslide Velocity

Abstract: Proper knowledge of velocity is required in accurately determining enormous destructive energy of a landslide. We present the first physics-based general analytical landslide velocity model that incorporates internal deformation and external forces: net driving force and viscous resistant. The model stands as a novel non-linear advective-dissipative system where classical Voellmy and inviscid Burgers' equation are specifications. Non-linear advection and external forcing fundamentally regulate the motion which substantially enhances our understanding of a coherently deforming landslide. Since analytical solutions provide fastest, cost-effective and best rigorous answer to the problem, we construct several new/general exact analytical solutions covering wider spectrum of landslide velocity. New solutions bridge existing gap between negligibly and massively deforming landslides. This provides a novel, rapid and consistent method for efficient coupling of different types of mass transports. Mechanism of landslide advection, stretching and approaching to steady-state has been explained. Shifting, up-lifting and stretching of velocity field stem from forcing and advection. Our solution describes fascinating breaking wave and emergence of landslide folding. This happens as the solution simultaneously introduces domain propagation, velocity up-lift and non-linear advection. Domain translation and stretching depends on net driving force and viscous drag controls shock wave generation, wave breaking, folding and velocity. Landslide dynamics are architectured by advection and reigned by system forcing. Analytically obtained velocities are close to observed values, constituting a new foundation of landslide velocity. This provides the practitioners with key information in instantly/accurately estimating impact force that is important in delineating hazard zones and mitigation of landslide hazards.

Deep Bayesian inference for seismic imaging with tasks

Abstract: We propose to use techniques from Bayesian inference and deep neural networks to translate uncertainty in seismic imaging to uncertainty in tasks performed on the image, such as horizon tracking. Seismic imaging is an ill-posed inverse problem because of unavoidable bandwidth and aperture limitations, which that is hampered by the presence of noise and linearization errors. Many regularization methods, such as transform-domain sparsity promotion, have been designed to deal with the adverse effects of these errors, however, these methods run the risk of biasing the solution and do not provide information on uncertainty in the image space and how this uncertainty impacts certain tasks on the image. A systematic approach is proposed to translate uncertainty due to noise in the data to confidence intervals of automatically tracked horizons in the image. The uncertainty is characterized by a convolutional neural network (CNN) and to assess these uncertainties, samples are drawn from the posterior distribution of the CNN weights, used to parameterize the image. Compared to traditional priors, in the literature it is argued that these CNNs introduce a flexible inductive bias that is a surprisingly good fit for many diverse domains in imaging. The method of stochastic gradient Langevin dynamics is employed to sample from the posterior distribution. This method is designed to handle large scale Bayesian inference problems with computationally expensive forward operators as in seismic imaging. Aside from offering a robust alternative to maximum a posteriori estimate that is prone to overfitting, access to these samples allow us to translate uncertainty in the image, due to noise in the data, to uncertainty on the tracked horizons. For instance, it admits estimates for the pointwise standard deviation on the image and for confidence intervals on its automatically tracked horizons.

Interplay of river and tidal forcings promotes loops in coastal channel networks

Abstract: Global coastlines and their dense populations have an uncertain future due to increased flooding, storms, and human modification. The distributary channel networks of deltas and marshes that plumb these coastlines present diverse architectures, including well-studied dendritic topologies. However, the quasi-stable loops that are frequent in many coastal networks have not yet been explained. We present a model for self-organizing networks inspired by vascular biophysics to show that loops emerge when the relative forcings between rivers and tides are comparable, resulting in interplay between processes at short timescales relative to network evolution. Using field data and satellite imaging, we confirm this control on 21 natural networks. Our comparison provides the first evidence that hydrodynamic fluctuations promote loop formation in geophysical systems.

Embracing Data Incompleteness for Better Earthquake Forecasting

Abstract: We propose two new methods to calibrate the parameters of the epidemic-type aftershock sequence (ETAS) model based on expectation maximization (EM) while accounting for temporal variation of catalo...

An introduction to distributed training of deep neural networks for segmentation tasks with large seismic datasets

Abstract: Deep learning applications are drastically progressing in seismic processing and interpretation tasks. However, the majority of approaches subsample data volumes and restrict model sizes to minimise computational requirements. Subsampling the data risks losing vital spatio-temporal information which could aid training whilst restricting model sizes can impact model performance, or in some extreme cases, renders more complicated tasks such as segmentation impossible. This paper illustrates how to tackle the two main issues of training of large neural networks: memory limitations and impracticably large training times. Typically, training data is preloaded into memory prior to training, a particular challenge for seismic applications where data is typically four times larger than that used for standard image processing tasks (float32 vs. uint8). Using a microseismic use case, we illustrate how over 750GB of data can be used to train a model by using a data generator approach which only stores in memory the data required for that training batch. Furthermore, efficient training over large models is illustrated through the training of a 7-layer UNet with input data dimensions of 4096X4096. Through a batch-splitting distributed training approach, training times are reduced by a factor of four. The combination of data generators and distributed training removes any necessity of data 1 subsampling or restriction of neural network sizes, offering the opportunity of utilisation of larger networks, higher-resolution input data or moving from 2D to 3D problem spaces.

Salinity and sodicity can cause hysteresis in soil hydraulic conductivity

Abstract: Soil salinization and sodification can cause detrimental effects to soil structure, with important implications to irrigated agriculture. Knowledge of the extent to which degradation in soil structure due to salinization and sodification is reversible is still lacking, however. The objective of our study is to examine the effects of the history of solute composition on the degree of reversibility in saturated hydraulic conductivity (Ksat). We systematically investigate the effects of salinity (electrolyte concentration) and sodicity (sodium adsorption ratio) on Ksat, for three soils of varying clay content. The soil column experiments yielded hysteresis graphs, in which Ksat does not go back to its original values after initial decay. We developed indices to quantify the degree of Ksat degradation and reversibility, and found that contrary to our expectations, high susceptibility to degradation does not always correlate with low capability of rehabilitation. Measurements of soil swelling helped us discern when degradation was mainly caused by swelling or clay dispersion. Our findings underscore the need of understanding hysteresis caused by salinity and sodicity, aiming at a better management of agricultural soils for sustainable long-term use.

Predicting water flow in fully and partially saturated porous media: a new fractal-based permeability model

Abstract: Predicting the permeability of porous media in saturated and partially saturated conditions is of crucial importance in many geo-engineering areas, from water resources to vadose zone hydrology or contaminant transport predictions. Many models have been proposed in the literature to estimate the permeability from properties of the porous media such as porosity, grain size or pore size. In this study, we develop a model of the permeability for porous media saturated by one or two fluid phases with all physically-based parameters using a fractal upscaling technique. The model is related to microstructural properties of porous media such as fractal dimension for pore space, fractal dimension for tortuosity, porosity, maximum radius, ratio of minimum pore radius and maximum pore radius, water saturation and irreducible water saturation. The model is favorably compared to existing and widely used models from the literature. Then, comparison with published experimental data for both unconsolidated and consolidated samples, we show that the proposed model estimate the permeability from the medium properties very well.

Nearly constant Q models of the generalized standard linear solid type and the corresponding wave equations

Abstract: Time-domain seismic forward and inverse modeling for a dissipative medium is a vital research topic to investigate the attenuation structure of the Earth. Constant Q, also called frequency independence of the quality factor, is a common assumption for seismic Q inversion. We propose the first- and second-order nearly constant Q dissipative models of the generalized standard linear solid type, using a novel Q-independent weighting function approach. The two new models, which originate from the Kolsky model (a nearly constant Q model) and the Kjartansson model (an exactly constant Q model), result in the corresponding wave equations in differential form. Even for extremely strong attenuation (e.g., Q = 5), the quality factor and phase velocity for the two new models are close to those for the Kolsky and Kjartansson models, in a frequency range of interest. The wave equations for the two new models involve explicitly a specified Q parameter and have compact and simple forms. We provide a novel perspective on how to build a nearly constant Q dissipative model which is beneficial for time-domain large scale wavefield forward and inverse modeling. This perspective could also help obtain other dissipative models with similar advantages. We also discuss the extension beyond viscoacousticity and other related issues, for example, extending the two new models to viscoelastic anisotropy.

Tsunami hydrodynamic force on a building using a SPH real scale numerical simulation

Abstract: One of the most important aspects in tsunami studies is the wave behavior when it approaches the coast. Information on physical parameters that characterize waves is often limited because of the diffilculties in achieving accurate measurements at the time of the event. The impact of a tsunami on the coast is governed by nonlinear physics such as turbulence with spatial and temporal variability. The use of the Smoothed Particle Hydrodynamic method (SPH) presents advantages over models based on two-dimensional Shallow Waters Equations (SWE), because the assumed vertical velocity simplifies hydrodynamics in two dimensions. The study presented here reports numerical SPH simulations of the tsunami event occurred in Coquimbo (Chile) on September 16 of 2015. On the basis of the reconstruction of the physical parameters that characterized this event (flow velocities, direction and water elevations), calibrated by a reference rodel, force values on buildings located on the study coast were numerically calculated, and compared with an estimate of the Chilean Structural Design Standard. Finally, discussion and conclusions of the comparison of both methodologies are presented, including an influence analysis of the topographical detail of the model in the estimation of hydrodynamic forces.

A continuous fracture front tracking algorithm with multi layer tip elements (MuLTipEl) for a plane strain hydraulic fracture.

Abstract: The problem of a plane strain hydraulic fracture propagating in a layered formation is considered. Fracture toughness, in-situ stress, and leak-off coefficient are assumed to vary by layer, while the elastic properties are kept constant throughout the domain for simplicity. The purpose of this study is to develop a numerical algorithm based on a fixed mesh approach, which is capable to solve the above problem accurately using elements which can even be larger than the layer size. In order to do this, the concept of fictitious tip stress is first introduced for determining the fracture front location. In this technique, an additional stress is applied to the tip element with the purpose to suppress opening and to mimic width corresponding to the actual fracture front location. A theoretical basis for this concept has been established and it is further calibrated for piece-wise constant elements. Once the ability to track the crack front location is developed, the effect of layers is included by vary properties as a function of front location. Several numerical examples benchmarking the numerical solution, as well as highlighting capabilities of the algorithm to tackle multiple thin layers accurately are presented.

The effects of ENSO, climate change and human activities on the water level of Lake Toba, Indonesia: a critical literature review

Abstract: This research quantitatively and qualitatively analyzes the factors responsible for the water level variations in Lake Toba, North Sumatra Province, Indonesia. According to several studies carried out from 1993 to 2020, changes in the water level were associated with climate variability, climate change, and human activities. Furthermore, these studies stated that reduced rainfall during the rainy season due to the El Nino Southern Oscillation (ENSO) and the continuous increase in the maximum and average temperatures were some of the effects of climate change in the Lake Toba catchment area. Additionally, human interventions such as industrial activities, population growth, and damage to the surrounding environment of the Lake Toba watershed had significant impacts in terms of decreasing the water level. However, these studies were unable to determine the factor that had the most significant effect, although studies on other lakes worldwide have shown these factors are the main causes of fluctuations or decreases in water levels. A simulation study of Lake Toba's water balance showed the possibility of having a water surplus until the mid-twenty-first century. The input discharge was predicted to be greater than the output; therefore, Lake Toba could be optimized without affecting the future water level. However, the climate projections depicted a different situation, with scenarios predicting the possibility of extreme climate anomalies, demonstrating drier climatic conditions in the future. This review concludes that it is necessary to conduct an in-depth, comprehensive, and systematic study to identify the most dominant factor among the three that is causing the decrease in the Lake Toba water level and to describe the future projected water level.

Earthquake Breakdown Energy Scaling Despite Constant Fracture Energy

Abstract: In the quest to determine fault weakening processes that govern earthquake mechanics, it is common to infer the earthquake breakdown energy from seismological measurements. Breakdown energy is observed to scale with slip, which is often attributed to enhanced fault weakening with continued slip or at high slip rates, possibly caused by flash heating and thermal pressurization. However, breakdown energy varies by more than six orders of magnitude, which is physically irreconcilable with prevailing material properties. We present a dynamic model that demonstrates that breakdown energy scaling can occur despite constant fracture energy and does not require thermal pressurization or other enhanced weakening. Instead, earthquake breakdown energy scaling occurs simply due to scale-invariant stress drop overshoot, which is affected more directly by the overall rupture mode -- crack-like or pulse-like -- rather than from a specific slip-weakening relationship. Our findings suggest that breakdown energy may be used to discern crack-like earthquakes from self-healing pulses with negative breakdown energy.

Quantitative Approximation of Geothermal Potential of Bakreswar Geothermal Area in Eastern India

Abstract: Proper utilization of geothermal energy for power generation is still ignored in India even after having it's enough potential as much as the equivalent to the other nonconventional energy resources of the country. A major thrust is required in this field of technology. The source of geothermal energy is the decay of the radio-nuclei, such as Uranium, Thorium, and Potassium inside the Earth's crust apart from the primordial heat source. The noble gas 4He is also produced during the radioactive disintegration process. Therefore, measuring the amount of 4He gas generated in the terrestrial radioactive process along with some other geochemical parameters in an Indian geothermal reservoir, the potential of the reservoir can be evaluated without performing conventional and detailed geochemical & geophysical techniques. Mathematical calculations relating to the radioactive disintegration to estimate the geothermal potential of the Bakreswar geothermal area at eastern India utilizing the concept of the 4He exploration technique has been described here. The study showed that the heat energy generated by the radioactive decay of 238Th, 238U and 235U inside Bakreswar geothermal reservoir was evaluated as 38 MW only considering the He emanated from the Agni Kunda hot spring. In addition, the depth of the geothermal reservoir was also evaluated to be about 1,100 m using some geophysical characteristics of the reservoir at the study area. Furthermore, the suitable locations of deep drilling for the installation of the probable geothermal power plant are also identified by investigating the resistivity survey profile of the study area.

Growth of barchan dunes of bidispersed granular mixtures.

Abstract: Barchans are dunes of crescentic shape found on Earth, Mars and other celestial bodies, growing usually on polydisperse granular beds. In this Letter, we investigate experimentally the growth of subaqueous barchans consisting of bidisperse grains. We found that the grain distribution within the dune changes with the employed pair, and that a transient stripe appears on the dune surface. We propose that observed patterns result from the competition between fluid entrainment and easiness of rolling for each grain type, and that grains segregate with a diffusion-like mechanism. Our results provide new insights into barchan structures found in other environments.

Surrogate and inverse modeling for two-phase flow in porous media via theory-guided convolutional neural network

Abstract: The theory-guided convolutional neural network (TgCNN) framework, which can incorporate discretized governing equation residuals into the training of convolutional neural networks (CNNs), is extended to two-phase porous media flow problems in this work. The two principal variables of the considered problem, pressure and saturation, are approximated simultaneously with two CNNs, respectively. Pressure and saturation are coupled with each other in the governing equations, and thus the two networks are also mutually conditioned in the training process by the discretized governing equations, which also increases the difficulty of model training. The coupled and discretized equations can provide valuable information in the training process. With the assistance of theory-guidance, the TgCNN surrogates can achieve better accuracy than ordinary CNN surrogates in two-phase flow problems. Moreover, a piecewise training strategy is proposed for the scenario with varying well controls, in which the TgCNN surrogates are constructed for different segments on the time dimension and stacked together to predict solutions for the whole time-span. For scenarios with larger variance of the formation property field, the TgCNN surrogates can also achieve satisfactory performance. The constructed TgCNN surrogates are further used for inversion of permeability fields by combining them with the iterative ensemble smoother (IES) algorithm, and sufficient inversion accuracy is obtained with improved efficiency.

Learning by example: fast reliability-aware seismic imaging with normalizing flows.

Abstract: Uncertainty quantification provides quantitative measures on the reliability of candidate solutions of ill-posed inverse problems. Due to their sequential nature, Monte Carlo sampling methods require large numbers of sampling steps for accurate Bayesian inference and are often computationally infeasible for large-scale inverse problems, such as seismic imaging. Our main contribution is a data-driven variational inference approach where we train a normalizing flow (NF), a type of invertible neural net, capable of cheaply sampling the posterior distribution given previously unseen seismic data from neighboring surveys. To arrive at this result, we train the NF on pairs of low- and high-fidelity migrated images. In our numerical example, we obtain high-fidelity images from the Parihaka dataset and low-fidelity images are derived from these images through the process of demigration, followed by adding noise and migration. During inference, given shot records from a new neighboring seismic survey, we first compute the reverse-time migration image. Next, by feeding this low-fidelity migrated image to the NF we gain access to samples from the posterior distribution virtually for free. We use these samples to compute a high-fidelity image including a first assessment of the image's reliability. To our knowledge, this is the first attempt to train a conditional network on what we know from neighboring images to improve the current image and assess its reliability.

Fast and robust earthquake source spectra and moment magnitudes from envelope inversion

Abstract: With the present study we introduce a fast and robust method to calculate the source displacement spectra of small earthquakes on a local to regional scale. The work is based on the publicly available Qopen method of full envelope inversion which is further tuned for the given purpose. Important source parameters -- seismic moment, moment magnitude, corner frequency and high-frequency fall-off -- are determined from the source spectra by fitting a simple earthquake source model. The method is demonstrated by means of a data set comprising the 2018 West Bohemia earthquake swarm. We report moment magnitudes, corner frequencies, and centroid moment tensors inverted from short period body waves with the Grond package for all earthquakes with a local magnitude larger than 1.8. Moment magnitudes calculated by envelope inversion show a very good agreement to moment magnitudes resulting from the probabilisitc moment tensor inversion. Furthermore, source displacement spectra from envelope inversion show a good agreement with spectra obtained by multiple taper analysis of the direct onsets of body waves, but are not affected by the large scatter of the second. The seismic moments obtained with the envelope inversion scale with corner frequencies according to $M_0 \propto f_{\mathrm{c}}^{-4.7}$. Earthquakes of the present data set result in a smaller stress drop for smaller magnitudes. Self-similarity of earthquake rupture is not observed. Additionally, we report frequency-dependent site amplification at the used stations.

Simulation of flow in deformable fractures using a quasi-Newton based partitioned coupling approach

Abstract: We introduce a partitioned coupling approach for iterative coupling of flow processes in deformable fractures embedded in a poro-elastic medium that is enhanced by interface quasi-Newton (IQN) methods. In this scope, a unique computational decomposition into a fracture flow and a poro-elastic domain is developed, where communication and numerical coupling of the individual solvers are realized by consulting the open-source library preCICE. The underlying physical problem is introduced by a brief derivation of the governing equations and interface conditions of fracture flow and poro-elastic domain followed by a detailed discussion of the partitioned coupling scheme. We evaluate the proposed implementation and undertake a convergence study to compare a classical interface quasi-Newton inverse least-squares (IQN-ILS) with the more advanced interface quasi-Newton inverse multi-vector Jacobian (IQN-IMVJ) method. These coupling approaches are verified for an academic test case before the generality of the proposed strategy is demonstrated by simulations of two complex fracture networks. In contrast to the development of specific solvers, we promote the simplicity and computational efficiency of the proposed partitioned coupling approach using preCICE and FEniCS for parallel computations of hydro-mechanical processes in complex, three-dimensional fracture networks.

Compressive time-lapse seismic monitoring of carbon storage and sequestration with the joint recovery model

Abstract: Time-lapse seismic monitoring of carbon storage and sequestration is often challenging because the time-lapse signature of the growth of CO2 plumes is weak in amplitude and therefore difficult to detect seismically. This situation is compounded by the fact that the surveys are often coarsely sampled and not replicated to reduce costs. As a result, images obtained for different vintages (baseline and monitor surveys) often contain artifacts that may be attributed wrongly to time-lapse changes. To address these issues, we propose to invert the baseline and monitor surveys jointly. By using the joint recovery model, we exploit information shared between multiple time-lapse surveys. Contrary to other time-lapse methods, our approach does not rely on replicating the surveys to detect time-lapse changes. To illustrate this advantage, we present a numerical sensitivity study where CO2 is injected in a realistic synthetic model. This model is representative of the geology in the southeast of the North Sea, an area currently considered for carbon sequestration. Our example demonstrates that the joint recovery model improves the quality of time-lapse images allowing us to monitor the CO2 plume seismically.

Benchmarking Exercises for Granular Flows

Abstract: For the 2007 International Forum on Landslide Disaster Management framework, our team performed several numerical simulations on both theoretical and natural cases of granular flows. The objective was to figure out the ability and the limits of our numerical model in terms of reproduction and prediction. Our benchmarking exercises show that for almost all the cases, the model we use is able to reproduce observations at the field scale. Calibrated friction angles are almost similar to that used in other models and the shape of the final deposits is in good agreement with observation. However, as it is tricky to compare the dynamics of natural cases, these exercises do not allow us to highlight the good ability to reproduce the behavior of natural landslides. Nevertheless, by comparing with analytical solution, we show that our model presents very low numerical dissipation due to the discretization and to the numerical scheme used. Finally, in terms of mitigation and prediction, the different friction angles used for each cases figure out the limits of using such model as long as constitutive equations for granular media are not known.