A hybrid-Galerkin finite-element method for seismic wave propagation in fractured media

Author(s):  
Janaki Vamaraju ◽  
Mrinal Sen ◽  
Mary F. Wheeler ◽  
Jonas De Basabe
2020 ◽  
Vol 221 (2) ◽  
pp. 857-878 ◽  
Author(s):  
Janaki Vamaraju ◽  
Mrinal K Sen ◽  
Jonas De Basabe ◽  
Mary Wheeler

SUMMARY The discontinuous Galerkin finite element method (DGM) is a promising algorithm for modelling wave propagation in fractured media. It allows for discontinuities in the displacement field to simulate fractures or faults in a model. Our 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. On the other hand, the spectral element method (SEM) can be used to simulate elastic wave propagation in non-fractured media. SEM uses continuous basis functions which do not allow for discontinuities in the displacement field. However, the computation cost of DGM is significantly larger than SEM due primarily to increase in the number of degrees of freedom. Here we propose a hybrid Galerkin method (HGM) for elastic wave propagation in fractured media that combines the salient features of each of the algorithm resulting in significant reduction in computational cost compared to DGM. We use DGM in areas containing fractures and SEM in regions without fractures. The coupling between the domains at the interfaces is satisfied in the weak form through interface conditions. The degree of reduction in computation time depends primarily on the density of fractures in the medium. In this paper, we formulate and implement HGM for seismic wave propagation in fractured media. Using realistic 2-D/3-D numerical examples, we show that our proposed HGM outperforms DGM with reduced computation cost and memory requirement while maintaining the same level of accuracy.


2020 ◽  
Author(s):  
Yang Xu ◽  
Xiaofei Chen ◽  
Dechao Han ◽  
Wei Zhang

<p>Numerical simulation of seismic wavefield is helpful to understand the propagation law of seismic wave in complex media. In addition, accurate simulation of seismic wave propagation is of great importance for seismic inversion. The discontinuous Galerkin finite element method(DG-FEM) combines the advantages of finite element method(FEM) and finite volume method(FVM) to effectively simulate the propagation characteristics of seismic waves in complex medium.</p><p>In this study, we use the hp-adaptive DG -FEM to perform accurate simulation of seismic wave propagation in complex topography and medium, and compare the results with the analytical solution of the Generalized Reflection/Transmission(GRT) coefficient method. Furthermore, ADE CFS-PML is modified and applied to DG-FEM, which greatly reduces the impact of artificial boundaries.</p>


2011 ◽  
Vol 121-126 ◽  
pp. 4880-4884
Author(s):  
Zhao Ling Wang ◽  
Zheng Ping Liu ◽  
Chi Zhang

In the paper, two-dimensional Tunnel seismic Wave field is Simulated with finite element method, and the in the tunnel model with fault zone load Ricker wavelet source on the workface, compared the case of wave propagation according to wave field snapshot and time record, can intuitively, accurately reflect the characteristics of seismic wave propagation in tunnel seismic prediction with geological disasters such as the fault zone and so on.


Sign in / Sign up

Export Citation Format

Share Document