Viscoelastic Wave Research Articles

On cost-efficient parallel iterative solvers for 3D frequency-domain seismic multisource viscoelastic anisotropic wave modeling

Solving large sparse linear systems in 3D frequency-domain seismic wave modeling, especially in viscoelastic anisotropic media, poses significant challenges due to the increasing number of discrete moduli and nonzero elements in the linear system matrix. The computational load surpasses that of acoustic or viscoacoustic media, making it even more challenging when dealing with multisource problems. Popular scientific tools for solving a linear system, such as multifrontal massively parallel sparse direct solver (MUMPS), structured matrix package (STRUMPACK), and portable extensible toolkit for scientific computation (PETSc), can be used, but their applicability to our specific problem has not been comprehensively evaluated. Our study aims to tackle the challenges in solving large, sparse, complex-valued symmetric linear systems with multiple right-hand side vectors for 3D frequency-domain seismic wave modeling. We have adopted the preconditioned conjugate gradient iterative algorithms as the foundation for our research, introducing two highly cost-effective parallel iterative solvers: the parallel symmetric successive overrelaxation conjugate gradient (P-SSORCG) and the parallel incomplete Cholesky conjugate gradient (P-ICCG). These novel solvers are subjected to a comprehensive comparative analysis against well-established scientific tools, such as MUMPS, STRUMPACK, and PETSc, in the context of 3D frequency-domain seismic wave modeling. We show their promising performances in a practical 3D SEG/EAGE overthrust model and demonstrate that the grouped P-SSORCG offers an efficient alternative to parallel direct solvers, particularly in situations wherein computational resources are limited.

General energy decay estimates for a variable coefficient viscoelastic semilinear wave equation with logarithmic source term and acoustic boundary conditions

In this paper, we consider a variable coefficient viscoelastic semilinear wave equation with logarithmic source term and acoustic boundary conditions in domains with non‐locally reacting boundary. By potential well method, we get the global existence of solution and through several important lemmas, we prove the general energy decay estimate of the system.

Fractional-order velocity-dilatation-rotation viscoelastic wave equation and numerical solution based on a constant-Q model

Viscoelastic wave equations based on the constant- Q (CQ) model can accurately describe the amplitude dissipation and phase distortion of waves in anelastic media. However, only three velocity or displacement components can be obtained directly by solving such equations. Starting from the time-domain second-order displacement viscoelastic wave equation, we derived the decoupled P- and S-wave displacement vector viscoelastic wave equation by using the polarization difference of P- and S-wave propagation in isotropic media. The equation can be transformed into the velocity-dilatation-rotation viscoelastic wave equation containing the first-order temporal derivative and fractional Laplacian operators, which can be solved directly by using the staggered-grid finite-difference and pseudospectral methods. We use the low-rank decomposition method to approximate the derived mixed space-wavenumber domain fractional Laplacian operators for modeling wave propagation in heterogeneous attenuating media. We also demonstrated the precision of our equation by comparing the numerical solutions with the analytical solutions. Furthermore, compared with the conventional velocity-stress viscoelastic wave equation, experimental results demonstrate that our equation can separate the pure P and S waves from the mixed wavefield during wavefield continuation. In addition, it can be separated into an equation containing predominantly an amplitude attenuation or a phase distortion term.

On a coupled system of viscoelastic wave equation of infinite memory with acoustic boundary conditions

This work deals with a coupled system of viscoelastic wave equation of infinite memory with mixed Dirichlet-Neumann boundary conditions. The coupling is via the acoustic boundary conditions on a portion of the boundary. The semigroup theory is used to show the well-posedness and regularity of the initial and boundary value problem. Moreover, we investigate exponential stability of the system taking into account Gearhart-Prüss’s theorem. Mathematics Subject Classification (2010): 35A01, 74B05, 93D15. Received 11 June 2021; Accepted 13 October 2021

Energy behavior for Sobolev solutions to viscoelastic damped wave models with time‐dependent oscillating coefficient

AbstractIn this work, we study the asymptotic behavior of the structurally damped wave equations arising from the viscoelastic mechanics. We are particularly interested in the complicated interaction of the time‐dependent oscillating coefficients on the Dirichlet Laplacian operator and the structurally damped terms. On the one hand, by the application of WKB analysis, we explore the asymptotic energy estimates of the wave equations influenced by four types of oscillating mechanisms. On the other hand, in order to prove the optimality of the energy estimates for the critical cases, typical coefficients and initial Cauchy data will be constructed to show the lower bound of the energy growth rate by the application of instability arguments.

Thermodynamically extended symplectic numerical simulation of viscoelastic, thermal expansion and heat conduction phenomena in solids

AbstractSymplectic numerical schemes for reversible dynamical systems predict the solution reliably over large times as well, and are a good starting point for extension to schemes for simulating irreversible situations like viscoelastic wave propagation and heat conduction coupled via thermal expansion occuring in rocks, plastics, biological samples etc. Dissipation error (artificial nonpreservation of energies and amplitudes) of the numerical solution should be as small as possible since it should not be confused with the real dissipation occurring in the irreversible system. In addition, the other well-known numerical artefact, dispersion error (artificial oscillations emerging at sharp changes), should also be minimal to avoid confusion with the true wavy behavior. The continuum thermodynamical aspects (respect for balances with fluxes, systematic constitutive relationships between intensive quantities and fluxes, the second law of thermodynamics with positive definite entropy production, and the spacetime-based kinematic viewpoint) prove valuable for obtaining such extended schemes and for monitoring the solutions. Generalizing earlier works in this direction, here, we establish and investigate such a numerical scheme for one-dimensional viscoelastic wave propagation in the presence of heat conduction coupled via thermal expansion, demonstrating long-term reliability and the applicability of thermodynamics-based quantities in supervising the quality of the solution.

Amplitude-preserving P-/S-wavefield separation with the discontinuous Galerkin method on unstructured meshes

The decoupling of P and S waves is an essential prerequisite for elastic reverse time migration, which effectively reduces crosstalk artifacts, but most wavefield separation algorithms are implemented on uniform rectangular grids. We develop an amplitude- and phase-preserving P- and S-wavefield separation approach on unstructured meshes, which can effectively decompose the original elastic wavefield into P and S wavefields. The isotropic case is considered. With the aid of viscoelastic theory, we choose to attenuate P or S waves and preserve the other wave mode, so as to achieve wavefield decomposition. Viscoelastic wave equations are first reformulated as decoupling wave equations with a selective strong attenuation. We then use the discontinuous Galerkin (DG) method to simulate decoupling P or S wavefield propagation on triangular and tetrahedral meshes. We adopt a quadrature-free DG approach and the arbitrary mesh is mapped into the reference mesh for numerical calculation, where no additional volume and surface integrations are involved. The amplitude and phase information of this vector decomposition agrees with that of the original elastic data. Four numerical examples are used to demonstrate the superior performance of this vector decomposition algorithm. The isotropic example indicates the applicability and correctness of our scheme and the second example displays superiority in handling strong velocity contrasts. The third example exhibits mesh flexibility in dealing with complex structures, such as caves, faults, and undulating surfaces. The last example indicates the effectiveness of our developed algorithm extended to a 3D case.

On the numerical resolution of the second-order viscoacoustic and viscoelastic anisotropic wave equations using the recursive convolution method

Accurate and efficient seismic wave modeling is fundamental for the high-resolution seismic full-waveform inversion and correct interpretation of seismic data. Based on the viscoelastic mechanism of the Generalized Standard Linear Solid and the constitutive relationship between stress and strain, we demonstrate a generalized recursive formula using the Taylor series expansion to directly compute the temporal convolution in the constitutive equation, rather than solving auxiliary differential equations of the memory variables in the conventional method. New sets of the second order viscoacoustic and viscoelastic anisotropic wave equations are established with the new recursive temporal convolutions, which completely replaces the traditional auxiliary differential equations of the memory variables. The error analysis shows the proposed method can offer better accuracy than the existing recursive methods, and is very competitive with the widely used memory variable method. The numerical results comparing the different methods show the accuracy and feasibility of the generalized recursive convolution method to simulate the viscoacoustic and viscoelastic anisotropic waves.

Stabilization of the Viscoelastic Wave Equation with Variable Coefficients and a Delay Term in Nonlocal Boundary Feedback

Stabilization of the Viscoelastic Wave Equation with Variable Coefficients and a Delay Term in Nonlocal Boundary Feedback

An Accurate and Efficient Recursive Convolution Method to Simulate Viscoacoustic and Viscoelastic Waves

An Accurate and Efficient Recursive Convolution Method to Simulate Viscoacoustic and Viscoelastic Waves

General decay of a variable coefficient viscoelastic wave equation with the logarithmic nonlinearity, not necessarily decreasing kernel and acoustic boundary conditions

General decay of a variable coefficient viscoelastic wave equation with the logarithmic nonlinearity, not necessarily decreasing kernel and acoustic boundary conditions

Qualitative analysis of the asymptotic behavior for a viscoelastic wave equation in the presence of acoustic and fractional conditions combined with nonlinear distributed delay in boundary feedback

Qualitative analysis of the asymptotic behavior for a viscoelastic wave equation in the presence of acoustic and fractional conditions combined with nonlinear distributed delay in boundary feedback

Decay rate of solutions to the Cauchy problem for a coupled system of viscoelastic wave equations with a strong delay in Rⁿ

"Using weighted spaces, we establish a general decay rate properties of solutions as T → ∞ for a coupled system of viscoelastic wave equations in Rⁿ under some conditions on g1, g2, φ. We exploit a density function to introduce weighted spaces for solutions and using an appropriate Lyapunov function. Keywords: Lyapunov function, relaxation function, density, decay rate, weighted spaces."

High-Accuracy Simulation of Rayleigh Waves Using Fractional Viscoelastic Wave Equation

The propagation of Rayleigh waves is usually accompanied by dispersion, which becomes more complex with inherent attenuation. The accurate simulation of Rayleigh waves in attenuation media is crucial for understanding wave mechanisms, layer thickness identification, and parameter inversion. Although the vacuum formalism or stress image method (SIM) combined with the generalized standard linear solid (GSLS) is widely used to implement the numerical simulation of Rayleigh waves in attenuation media, this type of method still has its limitations. First, the GSLS model cannot split the velocity dispersion and amplitude attenuation term, thus limiting its application in the Q-compensated reverse time migration/full waveform inversion. In addition, GSLS-model-based wave equation is usually numerically solved using staggered-grid finite-difference (SGFD) method, which may result in the numerical dispersion due to the harsh stability condition and poses complexity and computational burden. To overcome these issues, we propose a high-accuracy Rayleigh-waves simulation scheme that involves the integration of the fractional viscoelastic wave equation and vacuum formalism. The proposed scheme not only decouples the amplitude attenuation and velocity dispersion but also significantly suppresses the numerical dispersion of Rayleigh waves under the same grid sizes. We first use a homogeneous elastic model to demonstrate the accuracy in comparison with the analytical solutions, and the correctness for a viscoelastic half-space model is verified by comparing the phase velocities with the dispersive images generated by the phase shift transformation. We then simulate several two-dimensional synthetic models to analyze the effectiveness and applicability of the proposed method. The results show that the proposed method uses twice as many spatial step sizes and takes 0.6 times that of the GSLS method (solved by the SGFD method) when achieved at 95% accuracy.

A Class of Viscoelastic Wave Equations with Exponential Source and the Nonlinear Strong Damping

A Class of Viscoelastic Wave Equations with Exponential Source and the Nonlinear Strong Damping

General Stability for the Viscoelastic Wave Equation with Nonlinear Time-Varying Delay, Nonlinear Damping and Acoustic Boundary Conditions

This paper is focused on energy decay rates for the viscoelastic wave equation that includes nonlinear time-varying delay, nonlinear damping at the boundary, and acoustic boundary conditions. We derive general decay rate results without requiring the condition a2>0 and without imposing any restrictive growth assumption on the damping term f1, using the multiplier method and some properties of the convex functions. Here we investigate the relaxation function ψ, namely ψ′(t)≤−μ(t)G(ψ(t)), where G is a convex and increasing function near the origin, and μ is a positive nonincreasing function. Moreover, the energy decay rates depend on the functions μ and G, as well as the function F defined by f0, which characterizes the growth behavior of f1 at the origin.

Existence and General Decay of Solutions for a Weakly Coupled System of Viscoelastic Kirchhoff Plate and Wave Equations

In this paper, a weakly coupled system (by the displacement of symmetric type) consisting of a viscoelastic Kirchhoff plate equation involving free boundary conditions and the viscoelastic wave equation with Dirichlet boundary conditions in a bounded domain is considered. Under the assumptions on a more general type of relaxation functions, an explicit and general decay rate result is established by using the multiplier method and some properties of the convex functions.

Initial boundary value problem for a viscoelastic wave equation with Balakrishnan–Taylor damping and a delay term: decay estimates and blow-up result

In this paper, we study the initial boundary value problem for the following viscoelastic wave equation with Balakrishnan–Taylor damping and a delay term where the relaxation function satisfies g'(t)leq -xi (t)g^{r}(t), tgeq 0, 1leq r< frac{3}{2}. The main goal of this work is to study the global existence, general decay, and blow-up result. The global existence has been obtained by potential-well theory, the decay of solutions of energy has been established by introducing suitable energy and Lyapunov functionals, and a blow-up result has been obtained with negative initial energy.

Attenuation characteristics of shock waves in drilling and blasting based on viscoelastic wave theory

This study investigates the propagation and attenuation characteristics of shock waves produced by drilling and blasting. The focus is on a novel method for calculating the attenuation of the shock wave pressure. Based on the propagation of shock waves in a rod with a variable section and a Maxwell model for rock, theoretical solutions for the attenuation of the blast-induced shock wave pressure were derived and the main factors influencing the attenuation speed were analyzed. The attenuation law of shock waves for siltstone, limestone, and granite, with common borehole diameters of 76–150 mm and initial shock pressures of 1–6 GPa, was analyzed through numerical simulations. The results show that the shock wave falls off as the square root of the relative distance. The attenuation speed of the shock wave is mainly related to the borehole size, rock properties, and shock wave intensity. The attenuation speed increases with increasing initial shock pressure on the borehole wall, but decreases with increasing borehole radius. Compared with soft rock, the shock wave pressure attenuates faster in hard rock. Combining theoretical analysis and numerical simulations, an exponential attenuation coefficient was introduced to the attenuation formula. Statistical analysis shows that the exponential attenuation coefficient is linearly proportional to the initial shock pressure. Furthermore, a new method for calculating the attenuation of shock wave pressure was proposed, reflecting the comprehensive influence of rock characteristics and shock wave intensity. The new method is an improvement over existing methods, which only consider the effect of rock properties on shock wave attenuation. The reliability of the new method was verified against the results of field tests. The findings provide a theoretical basis for the optimization of blasting parameters and assessments of blast-resistance for structures.

General decay and blow-up of solution for a viscoelastic wave equation with nonlinear degenerate damping and source terms

We are concerned with the initial-boundary value problem for a nonlinear viscoelastic problem with the degenerate damping and the nonlinear source term in this paper. Under appropriate assumptions on the initial data, memory kernel and parameters, we establish several results regarding the global solution, blow-up properties and general decay in the case of positive initial energy. We overcome the difficulties brought by positive initial energy by using the potential well theory.