Search

We present a numerical method to simulate seismic wave propagation in coupled fluid and porous media. We developed a numerical finite element–based algorithm to approximate solutions to viscoacoustic and Biot's equations, considering the open pore conditions at the interfaces between both media. The algorithm architecture allows to simulate arbitrary distributions of viscoacoustic and poroelastic regions, facilitating the modelling of heterogeneous systems involving complex geometries. The algorithm includes a double parallelization scheme whose efficiency in terms of computing time and memory requirements was tested for different core distributions and mesh sizes. We validate our proposal by performing a comparison between its results and those obtained with a well‐known freely available code. We test its capabilities by studying two different scenarios with geophysical interest: a lake with an irregular bottom and a fractured porous medium. [ABSTRACT FROM AUTHOR]

The controlled source electro-magnetic (CSEM) method is increasingly used for in-land and off-shore subsurface characterization. Given its complex dependence between data and the parameters of interest, there is a crucial need for performant numerical algorithms that can simulate the CSEM response of 3-D geological structures. Here, we present two finite element (FE) algorithms for simulating the CSEM response in 3-D media with isotropic conductivity. A primary/secondary field approach is used to avoid the singularity introduced by the source. The primary field is computed semi-analytically for a horizontally layered model and different sources. The secondary field is obtained by discretizing the diffusive frequency-domain Maxwell's equations with non-conforming FE. The two numerical algorithms are specifically designed to work on distributed-memory computers: (1) an iterative procedure with domain decomposition and (2) a direct and global algorithm. We evaluate their performance by computing their speed up on parallel processors, and solving problems with realistic conductivity structures. We also compare the accuracy of the solutions with published results on canonical models. The results shown here demonstrate the functionality of the two methodologies presented for specific cases when computing 3-D CSEM solutions. [ABSTRACT FROM AUTHOR]

Numerical simulations of seismic wave propagation in fractured media are often performed in the framework of the linear slip theory (LST). Therein, fractures are represented as interfaces and their mechanical properties are characterized through a compliance matrix. This theory has been extended to account for energy dissipation due to viscous friction within fluid-filled fractures by using complex-valued frequency-dependent compliances. This is, however, not fully adequate for fractured porous rocks in which wave-induced fluid flow (WIFF) between fractures and host rock constitutes a predominant seismic attenuation mechanism. In this letter, we develop an approach to incorporate WIFF effects directly into the LST for a ID system via a complex-valued, frequency-dependent fracture compliance. The methodology is validated for a medium permeated by regularly distributed planar fractures, for which an analytical expression for the complex-valued normal compliance is determined in the framework of quasistatic poroelasticity. There is good agreement between synthetic seismograms generated using the proposed recipe and those obtained from comprehensive, but computationally demanding, poroelastic simulations. [ABSTRACT FROM AUTHOR]