scholarly journals A discontinuous Galerkin method for wave propagation in orthotropic poroelastic media with memory terms

2019 ◽  
Vol 397 ◽  
pp. 108865 ◽  
Author(s):  
Jiangming Xie ◽  
M. Yvonne Ou ◽  
Liwei Xu
Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. T163-T174 ◽  
Author(s):  
Jonás D. De Basabe ◽  
Mrinal K. Sen ◽  
Mary F. Wheeler

We have formulated and implemented a discontinuous Galerkin method (DGM) for elastic wave propagation that allows for discontinuities in the displacement field to simulate fractures or faults. The approach is based on the interior-penalty formulation of DGM, and the fractures are simulated using the linear-slip model, which is incorporated into the weak formulation by including an additional term that is similar to the penalty term but uses the fracture compliance instead of an arbitrary penalty parameter. We have calibrated our results against an analytic solution of fracture-induced anisotropy for a set of elongated horizontal fractures, and we have evaluated numerical examples that simulate the reflection and transmission of waves at a fracture and at fracture interface waves. This method can further be used with models containing intersecting fractures and multiple fracture sets in 2D or 3D domains.


Geophysics ◽  
2011 ◽  
Vol 76 (4) ◽  
pp. N13-N28 ◽  
Author(s):  
Bastien Dupuy ◽  
Louis De Barros ◽  
Stephane Garambois ◽  
Jean Virieux

Biphasic media with a dynamic interaction between fluid and solid phases must be taken into account to accurately describe seismic wave amplitudes in subsurface and reservoir geophysical applications. Consequently, the modeling of the wave propagation in heteregeneous porous media, which includes the frequency-dependent phenomena of the fluid-solid interaction, is considered for 2D geometries. From the Biot-Gassmann theory, we have deduced the discrete linear system in the frequency domain for a discontinuous finite-element method, known as the nodal discontinuous Galerkin method. Solving this system in the frequency domain allows accurate modeling of the Biot wave in the diffusive/propagative regimes, enhancing the importance of frequency effects. Because we had to consider finite numerical models, we implemented perfectly matched layer techniques. We found that waves are efficiently absorbed at the model boundaries, and that the discretization of the medium should follow the same rules as in the elastodynamic case, that is, 10 grids per minimum wavelength for a P0 interpolation order. The grid spreading of the sources, which could be stresses or forces applied on either the solid phase or the fluid phase, did not show any additional difficulties compared to the elastic problem. For a flat interface separating two media, we compared the numerical solution and a semianalytic solution obtained by a reflectivity method in the three regimes where the Biot wave is propagative, diffusive/propagative, and diffusive. In all cases, fluid-solid interactions were reconstructed accurately, proving that attenuation and dispersion of the waves were correctly accounted for. In addition to this validation in layered media, we have explored the capacities of modeling complex wave propagation in a laterally heterogeneous porous medium related to steam injection in a sand reservoir and the seismic response associated to a fluid substitution.


Sign in / Sign up

Export Citation Format

Share Document