scholarly journals Carleman estimates for the wave equation in heterogeneous media with non-convex interface

2022 ◽  
Vol 311 ◽  
pp. 1-28
Author(s):  
Lucie Baudouin ◽  
Pamela Godoy ◽  
Alberto Mercado
Keyword(s):  
Wave Equation ◽  
Heterogeneous Media ◽  
Carleman Estimates

An analytic signal based accurate time-domain visco-acoustic wave equation from the Constant-Q theory

Geophysics ◽  
2021 ◽  
pp. 1-58
Author(s):  
Hongwei Liu ◽  
Yi Luo
Keyword(s):  
Wave Equation ◽  
Acoustic Wave ◽  
Time Domain ◽  
Seismic Waves ◽  
Heterogeneous Media ◽  
Relaxation Property ◽  
Analytic Signal ◽  
Acoustic Wave Equation ◽  
Real Component ◽  
Frequency Components

We present a concise time-domain wave equation to accurately simulate wave propagation in visco-acoustic media. The central idea behind this work is to dismiss the negative frequency components from a time-domain signal by converting the signal to its analytic format. The negative frequency components of any analytic signal are always zero, meaning we can construct the visco-acoustic wave equation to honor the relaxation property of the media for positive frequencies only. The newly proposed complex-valued wave equation (CWE) represents the wavefield with its analytic signal, whose real part is the desired physical wavefield, while the imaginary part is the Hilbert transform of the real component. Specifically, this CWE is accurate for both weak and strong attenuating media in terms of both dissipation and dispersion and the attenuation is precisely linear with respect to the frequencies. Besides, the CWE is easy and flexible to model dispersion-only, dissipation-only or dispersion-plus-dissipation seismic waves. We have verified these CWEs by comparing the results with analytical solutions, and achieved nearly perfect matching. Except for the homogeneous Q media, we have also extended the CWEs to heterogeneous media. The results of the CWEs for heterogeneous Q media are consistent with those computed from the nonstationary operator based Fourier Integral method and from the Standard Linear Solid (SLS) equations.

scholarly journals Modelling damped acoustic waves by a dissipation-preserving conformal symplectic method

2017 ◽  
Vol 473 (2199) ◽  
pp. 20160798 ◽  
Author(s):  
Wenjun Cai ◽  
Huai Zhang ◽  
Yushun Wang
Keyword(s):  
Wave Equation ◽  
Acoustic Wave ◽  
Acoustic Waves ◽  
Conservation Law ◽  
Heterogeneous Media ◽  
Numerical Dispersion ◽  
Discrete Version ◽  
Symplectic Method ◽  
Ode System ◽  
Strang Splitting

We propose a novel stable and efficient dissipation-preserving method for acoustic wave propagations in attenuating media with both correct phase and amplitude. Through introducing the conformal multi-symplectic structure, the intrinsic dissipation law and the conformal symplectic conservation law are revealed for the damped acoustic wave equation. The proposed algorithm is exactly designed to preserve a discrete version of the conformal symplectic conservation law. More specifically, two subsystems in conjunction with the original damped wave equation are derived. One is actually the conservative Hamiltonian wave equation and the other is a dissipative linear ordinary differential equation (ODE) system. Standard symplectic method is devoted to the conservative system, whereas the analytical solution is obtained for the ODE system. An explicit conformal symplectic scheme is constructed by concatenating these two parts of solutions by the Strang splitting technique. Stability analysis and convergence tests are given thereafter. A benchmark model in homogeneous media is presented to demonstrate the effectiveness and advantage of our method in suppressing numerical dispersion and preserving the energy dissipation. Further numerical tests show that our proposed method can efficiently capture the dissipation in heterogeneous media.

scholarly journals Simulation of elastic wave propagation across fractures using a nodal discontinuous Galerkin method—theory, implementation and validation

2019 ◽  
Vol 219 (3) ◽  
pp. 1900-1914 ◽  
Author(s):  
T Möller ◽  
W Friederich
Keyword(s):  
Wave Propagation ◽  
Wave Equation ◽  
Elastic Wave ◽  
Discontinuous Galerkin ◽  
Heterogeneous Media ◽  
Rheological Behaviour ◽  
Displacement Discontinuity ◽  
Signal Frequency ◽  
Elastic Wave Equation ◽  
Numerical Flux

SUMMARY An existing nodal discontinuous Galerkin (NDG) method for the simulation of seismic waves in heterogeneous media is extended to media containing fractures with various rheological behaviour. Fractures are treated as 2-D surfaces where Schoenberg’s linear slip or displacement discontinuity condition is applied as an additional boundary condition to the elastic wave equation which is in turn implemented as an additional numerical flux within the NDG formulation. Explicit expressions for the new numerical flux are derived by considering the Riemann problem for the elastic wave equation at fractures with varying rheologies. In all cases, we obtain further first order differential equations that fully describe the temporal evolution of the particle velocity jump at the fracture. Our flux formulation allows to separate the effect of a fracture from flux contributions due to simple welded interfaces enabling us to easily declare element faces as parts of a fracture. We make use of this fact by first generating the numerical mesh and then building up fractures by selecting appropriate element faces instead of adjusting the mesh to pre-defined fracture surfaces. The implementation of the new numerical fluxes into NDG is verified in 1-D by comparison to an analytical solution and in 2-D by comparing the results of a simulation valid in 1-D and 2-D. Further numerical examples address the effect of fracture systems on seismic wave propagation in 1-D and 2-D featuring effective anisotropy and coda generation. Finally, a study of the reflective and transmissive behaviour of fractures indicates that reflection and transmission coefficients are controlled by the ratio of signal frequency and relaxation frequency of the fracture.

scholarly journals Effective models for the multidimensional wave equation in heterogeneous media over long time and numerical homogenization

2016 ◽  
Vol 26 (14) ◽  
pp. 2651-2684 ◽  
Author(s):  
Assyr Abdulle ◽  
Timothée Pouchon
Keyword(s):  
Wave Equation ◽  
Wave Equations ◽  
Heterogeneous Media ◽  
Multiscale Model ◽  
Bloch Wave ◽  
Time Interval ◽  
Numerical Homogenization ◽  
Effective Equation ◽  
Long Time ◽  
Effective Models

A family of effective equations that capture the long time dispersive effects of wave propagation in heterogeneous media in an arbitrary large periodic spatial domain [Formula: see text] is proposed and analyzed. For a wave equation with highly oscillatory periodic spatial tensors of characteristic length [Formula: see text], we prove that the solution of any member of our family of effective equations is [Formula: see text]-close to the true oscillatory wave over a time interval of length [Formula: see text] in a norm equivalent to the [Formula: see text] norm. We show that the previously derived effective equation in [T. Dohnal, A. Lamacz and B. Schweizer, Bloch-wave homogenization on large time scales and dispersive effective wave equations, Multiscale Model. Simulat. 12 (2014) 488–513] belongs to our family of effective equations. Moreover, while Bloch wave techniques were previously used, we show that asymptotic expansion techniques give an alternative way to derive such effective equations. An algorithm to compute the tensors involved in the dispersive equation and allowing for efficient numerical homogenization methods over long time is proposed.

Multiscale modeling of acoustic wave propagation in 2D media

Geophysics ◽  
2014 ◽  
Vol 79 (2) ◽  
pp. T61-T75 ◽  
Author(s):  
Richard L. Gibson ◽  
Kai Gao ◽  
Eric Chung ◽  
Yalchin Efendiev
Keyword(s):  
Wave Propagation ◽  
Wave Equation ◽  
Acoustic Wave ◽  
Heterogeneous Media ◽  
Coarse Grid ◽  
Basis Functions ◽  
Seismic Survey ◽  
Fine Scale ◽  
Acoustic Wave Propagation ◽  
Fine Grid

Conventional finite-difference methods produce accurate solutions to the acoustic and elastic wave equation for many applications, but they face significant challenges when material properties vary significantly over distances less than the grid size. This challenge is likely to occur in reservoir characterization studies, because important reservoir heterogeneity can be present on scales of several meters to ten meters. Here, we describe a new multiscale finite-element method for simulating acoustic wave propagation in heterogeneous media that addresses this problem by coupling fine- and coarse-scale grids. The wave equation is solved on a coarse grid, but it uses basis functions that are generated from the fine grid and allow the representation of the fine-scale variation of the wavefield on the coarser grid. Time stepping also takes place on the coarse grid, providing further speed gains. Another important property of the method is that the basis functions are only computed once, and time savings are even greater when simulations are repeated for many source locations. We first present validation results for simple test models to demonstrate and quantify potential sources of error. These tests show that the fine-scale solution can be accurately approximated when the coarse grid applies a discretization up to four times larger than the original fine model. We then apply the multiscale algorithm to simulate a complete 2D seismic survey for a model with strong, fine-scale scatterers and apply standard migration algorithms to the resulting synthetic seismograms. The results again show small errors. Comparisons to a model that is upscaled by averaging densities on the fine grid show that the multiscale results are more accurate.

Sign in / Sign up

Export Citation Format

Share Document