Wave Propagation In Viscoelastic Media Research Articles

Analysis of P-wave propagation in filled jointed rock mass with viscoelastic properties

Stress wave propagation in filled jointed rock mass is an important basis for the analysis of rock dynamic stability. The presence of small cracks in rock mass often causes the rock to exhibit viscoelastic behavior, resulting in the amplitude attenuation and time delay of the stress wave. In view of the incomplete understanding of stress wave propagation in viscoelastic media, this study presents an analysis on stress wave propagation in filled jointed rock mass with viscoelastic properties based on the time-domain recursive analysis method. By introducing the quality factors of stress wave in viscoelastic medium, the propagation equation of plane P-wave in filled jointed rock mass with viscoelasticity is established. Subsequently, the transmission and reflection characteristics of P-wave are accordingly analyzed. The results indicate that the transmission and reflection coefficients of P-wave through the filled joint decrease obviously due to the attenuation effect of rock viscoelasticity on wave amplitude, but the coefficients gradually stabilize as the quality factor increases. As the propagation distance of the P-wave increases in viscoelastic rock masses, the transmission and reflection coefficients decrease gradually. The normal stiffness of joint contact interface has a decisive influence on the wave propagation. Specifically, a higher normal stiffness leads to stronger transmission and poorer reflection. The larger the thickness of filled joint, the more obvious the attenuation of P-wave becomes. Additionally, the viscoelasticity of filling layer further aggravates the influence of joint thickness on the attenuation of P-wave. The study comprehensively takes into account the effect of viscoelasticity of both the rock mass and the filling layer on transmission and reflection of waves. The research results provide certain theoretical support for the dynamic response and stability analysis of rock mass containing structural planes.

Research on the forward-modeling method of viscoelastic Chapman extension in partially saturated rocks

The propagation of seismic waves in fractured rocks is greatly influenced by the fracture system and fluid content. Seismic anisotropy varies with frequency and fluid saturation. Previous theories on frequency-dependent anisotropy (FDA) are mostly limited to the assumption of single-phase fluid, whereas almost all reservoirs are usually partially or fully saturated with one or more fluids. The inelasticity, heterogeneity, and anisotropy exhibited by real earth media seriously challenge the traditional theory of uniform fully elastic media. Therefore, based on the Chapman elastic theory, we consider the viscoelastic model, the relative mobility of saturated fluids, and the coupled effects of commonly independently considered squirt flow and patch effects on the frequency-dependent anisotropic seismic wave propagation in viscoelastic media. By calculating the frequency-dependent anisotropic elastic coefficients of two-phase immiscible fluid-saturated fractured rocks, a new viscoelastic Chapman extension model for FDA in partially saturated rocks, which includes the unified seismic wave propagation effects of squirt flow and patch effects, is constructed. The effect of relative permeability is significant, and compared to the case of complete saturation, the effective fluid mobility in partially saturated rocks may be lower, which may lead to nonmonotonic changes in elastic modulus and water saturation. Numerical experiments are conducted on our viscoelastic Chapman extension model to discuss the coupling effects of squirt flow and patch effects on the FDA of viscoelastic media in cracked rocks with partially saturated reservoirs, under the conditions of no fractures and the existence of fractures. The simulation results of the P-wave modulus, P-wave velocity, and anisotropy parameters confirm the validity of our method and model. We link the anisotropy of viscoelastic media with fluid flow parameters in fractures, which is conducive to improving the understanding of frequency-dependent seismic anisotropy in partially saturated rocks with fractures and the combination of seismology and reservoir engineering.

Distributed-order time-fractional wave equations

Distributed-order time-fractional wave equations appear in the modeling of wave propagation in viscoelastic media. The material characteristics of the medium are modeled through constitutive functions or distributions in the distributed-order constitutive law. In this work we propose to take positive Radon measures for the constitutive "functions". First, we derive a thermodynamical restriction on the constitutive measures which is easy to check, and therefore suitable for applications. Then we prove that the setting with measures in combination with the derived thermodynamical restriction guarantee existence and uniqueness of solutions for the distributed-order fractional wave equation. We further discuss the support and regularity of the fundamental solution, and conclude with a discussion on wave velocities.

Spectral solutions for a class of nonlinear wave equations with Riesz fractional based on Legendre collocation technique

A numerical investigation is presented in this work for a class of Riesz space-fractional nonlinear wave equations (MD-RSFN-WEs). The presence of a spatial Laplacian of fractional order, stated by fractional Riesz derivatives, is taken into consideration by the model. The fractional wave equation governs mechanical diffusive wave propagation in viscoelastic medium with power-law creep and, as a result, gives a physical understanding of this equation within the context of dynamic viscoelasticity. To deal with the independent variables, a totally spectral collocation approach is used. Our approach has shown to be more precise, efficient, and practical for the present model. The findings demonstrated that the spectral scheme is exponentially convergent.

Nonlinear least-squares estimation of shear elasticity and shear viscosity in viscoelastic media

Rapid simulations of three-dimensional (3D) shear wave propagation in viscoelastic media with a Kelvin-Voigt model are enabled by Green’s functions accelerated with graphics processing unit parallelization and high-performance computing resources. Cross-correlation analysis of the wave motion yields an estimate of shear elasticity, where the errors in the estimated values increase with rising shear viscosity. To reduce errors in the shear elasticity estimates, a nonlinear least-squares routine that accounts for the effects of propagation, attenuation, and dispersion is applied to 3D simulated shear wave data in viscoelastic media. The nonlinear least-squares approach also yields an estimate of the shear viscosity when applied to simulated particle velocity waveforms. Post-processing, cross-correlation analysis, and the nonlinear least-squares routine are evaluated on a desktop computer using MATLAB. Estimation of viscoelastic parameters is performed at the focal depth for several combinations of shear elasticity and shear viscosity. The errors in the estimated shear elasticities obtained from the nonlinear least-squares routine and cross-correlation are determined as a function of the input shear elasticity and shear viscosity values. The estimates provided by the nonlinear least-squares routine for shear elasticity and shear viscosity are shown to approach the simulation input values for each combination of viscoelastic parameters.

Convergence Study of Variational Space-Time Coupled Least-Squares Frameworks in Simulation of Wave Propagation in Viscoelastic Medium

There are many novel applications of space-time decoupled least squares and Galerkin formulations that simulate wave propagation through different types of media. Numerical simulation of stress wave propagation through viscoelastic medium is commonly carried out using the space-time decoupled Galerkin weak form in site response problem, etc. In a recent investigation into accuracy of this formulation in simulating elastic wave propagation, it was shown that the diffusive and dispersive errors are greatly reduced when space-time coupled least squares formulation is used instead in variational form. This paper investigates convergence characteristics of both formulations. To this end, two test cases, which are site response and impact models for viscoelastic medium, are employed, one dominated by dispersive and the other by diffusive numerical error. Convergence studies reveal that, compared to the commonly used space-time decoupled Galerkin and the coupled least squares formulation has much lower numerical errors, higher rates of convergence, and ability to take far larger time increments in the evolution process. In solving such models, coefficient matrices would require updating after each time step, a process that can be very costly. However large time steps allowed by cLs are expected to be a significant feature in reducing the overall computational cost.

General nonlocal Kelvin–Voigt viscoelasticity: application to wave propagation in viscoelastic media

General nonlocal Kelvin–Voigt viscoelasticity: application to wave propagation in viscoelastic media

Mathematical modeling of wave propagation in viscoelastic media with the fractional Zener model

The question of interest for the presented study is the mathematical modeling of wave propagation in dissipative media. The generalized fractional Zener model in the case of dimension d (d=1,2,3) is considered. This work is devoted to the mathematical analysis of such model: existence and uniqueness of the strong and weak solution and energy decay result which guarantees the wave dissipation. The existence of the weak solution is shown using a priori estimates for solutions which are also presented.

$Q$ Estimation by Combining ISD With LSR Method Based on Shaping-Regularized Inversion

The quality factor $Q$ is an indispensable parameter for studying wave propagation in viscoelastic media. $Q$ can not only be implemented for improving the quality of wave records but can also be used for directly indicating frequency-dependent anomalies induced by fluids, so it is widely used in seismic exploration and clinical medicine. $Q$ estimation here refers to the extraction of $Q$ information from seismic data; it has aroused lots of attention but is still somewhat controversial due to the limitations of existing methods. In this letter, combined with a logarithmic spectral ratio (LSR) algorithm, we have introduced a sparse-constrained inversion spectral decomposition (ISD) method for average- $Q$ estimation (LSR-ISD), and have used shaping regularization to solve for the spectrum ratio. Then, through regularized linear inversion, average- $Q$ was converted to an interval- $Q$ value. Finally, we have applied this method to synthetic data and field data. Numerical examples and field data application demonstrate that the proposed method produces a series of results with high resolution and good stability.

Q full-waveform inversion based on the viscoacoustic equation

Presently, most full-waveform inversion methods are developed for elastic media and ignore the effect of attenuation. The calculation of the quality factor Q is based on velocity parameter inversion under the assumption of a given Q-model that is obtained by tomographic inversion. However, the resolution of the latter is low and cannot reflect the amplitude attenuation and phase distortion during wave propagation in viscoelastic media. Thus, a Q waveform inversion method is proposed. First, we use standard linear body theory to describe attenuation and then we derive the simplified viscoacoustic equation that characterizes amplitude attenuation and phase distortion. In comparison with conventional equations, the simplified equation involves no memory variables and therefore requires less memory during computation. Moreover, the implementations of the attenuation compensation are easier. The adjoint equation and the corresponding gradient equation with respect to either L2-norm or the zero-lag cross-correlation objective function are then derived and the regularization strategy for overcoming the instability during numerical solution of the adjoint equation is proposed. The Q waveform inversion is developed using the limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) iteration method for known velocity. To alleviate the dependence of the waveform inversion on the initial model and overcome cycle skipping to some extent, we adopt multiscale analysis. Furthermore, anti-noise property and double-parameter inversion are assessed based on the results of numerical modeling.

Shear wave propagation in viscoelastic media: validation of an approximate forward model

Many approaches to elastography incorporate shear waves; in some systems these are produced by acoustic radiation force (ARF) push pulses. Understanding the shape and decay of propagating shear waves in lossy tissues is key to obtaining accurate estimates of tissue properties, and so analytical models have been proposed. In this paper, we reconsider a previous analytical model with the goal of obtaining a computationally straightforward and efficient equation for the propagation of shear waves from a focal push pulse. Next, this model is compared with an experimental optical coherence tomography (OCT) system and with finite element models, in two viscoelastic materials that mimic tissue. We find that the three different cases—analytical model, finite element model, and experimental results—demonstrate reasonable agreement within the subtle differences present in their respective conditions. These results support the use of an efficient form of the Hankel transform for both lossless (elastic) and lossy (viscoelastic) media, and for both short (impulsive) and longer (extended) push pulses that can model a range of experimental conditions.

A high-order accurate five-equations compressible multiphase approach for viscoelastic fluids and solids with relaxation and elasticity

A novel Eulerian approach is proposed for numerical simulations of wave propagation in viscoelastic media, for application to shocks interacting with interfaces between fluids and solids. We extend the five-equations multiphase interface-capturing model, based on the idea that all the materials (gases, liquids, solids) obey the same equation of state with spatially varying properties, to incorporate the desired constitutive relation; in this context, interfaces are represented by discontinuities in material properties. We consider problems in which the deformations are small, such that the substances can be described by linear constitutive relations, specifically, Maxwell, Kelvin–Voigt or generalized Zener models. The main challenge lies in representing the combination of viscoelastic, multiphase and compressible flow. One particular difficulty is the calculation of strains in an Eulerian framework, which we address by using a conventional hypoelastic model in which an objective time derivative (Lie derivative) of the constitutive relation is taken to evolve strain rates instead. The resulting eigensystem is analyzed to identify wave speeds and characteristic variables. The spatial scheme is based on a solution-adaptive formulation, in which a discontinuity sensor discriminates between smooth and discontinuous regions. To compute the convective fluxes, explicit high-order central differences are applied in smooth regions, while a high-order finite-difference Weighted Essentially Non-Oscillatory (WENO) scheme is used at discontinuities (shocks, material interfaces and contacts). The numerical method is verified in a comprehensive fashion using a series of smooth and discontinuous (shocks and interfaces), one- and two-dimensional test problems.

Subordination in a Class of Generalized Time-Fractional Diffusion-Wave Equations

Abstract Motivated by recently proposed generalizations of the diffusion-wave equation with the Caputo time fractional derivative of order α ∈ (1, 2), in the present survey paper a class of generalized time-fractional diffusion-wave equations is introduced. Its definition is based on the subordination principle for Volterra integral equations and involves the notion of complete Bernstein function. Various members of this class are surveyed, including the distributed-order time-fractional diffusion-wave equation and equations governing wave propagation in viscoelastic media with completely monotone relaxation moduli.

Application of a forward model of axisymmetric shear wave propagation in viscoelastic media to shear wave elastography

A simple but general solution of Navier's equation for axisymmetric shear wave propagation in a homogeneous isotropic viscoelastic medium is presented. It is well-suited for use as a forward model for some acoustic radiation force impulse based shear wave elastography applications because it does not require precise knowledge of the strength of the source, nor its spatial or temporal distribution. Instead, it depends on two assumptions: (1) the source distribution is axisymmetric and confined to a small region near the axis of symmetry, and (2) the propagation medium is isotropic and homogeneous. The model accounts for the vector polarization of shear waves and exactly represents geometric spreading of the shear wavefield, whether spherical, cylindrical, or neither. It makes no assumption about the frequency dependence of material parameters, i.e., it is material-model independent. Validation using measured shear wavefields excited by acoustic radiation force in a homogeneous gelatin sample show that the model accounts for well over 90% of the measured wavefield "energy." An optimal fit of the model to simulated shear wavefields with noise in a homogeneous viscoelastic medium enables estimation of both the shear storage modulus and shear wave attenuation to within 1%.

Simulating Seismic Wave Propagation in Viscoelastic Media with an Irregular Free Surface

In seismic numerical simulations of wave propagation, it is very important for us to consider surface topography and attenuation, which both have large effects (e.g., wave diffractions, conversion, amplitude/phase change) on seismic imaging and inversion. An irregular free surface provides significant information for interpreting the characteristics of seismic wave propagation in areas with rugged or rapidly varying topography, and viscoelastic media are a better representation of the earth’s properties than acoustic/elastic media. In this study, we develop an approach for seismic wavefield simulation in 2D viscoelastic isotropic media with an irregular free surface. Based on the boundary-conforming grid method, the 2D time-domain second-order viscoelastic isotropic equations and irregular free surface boundary conditions are transferred from a Cartesian coordinate system to a curvilinear coordinate system. Finite difference operators with second-order accuracy are applied to discretize the viscoelastic wave equations and the irregular free surface in the curvilinear coordinate system. In addition, we select the convolutional perfectly matched layer boundary condition in order to effectively suppress artificial reflections from the edges of the model. The snapshot and seismogram results from numerical tests show that our algorithm successfully simulates seismic wavefields (e.g., P-wave, Rayleigh wave and converted waves) in viscoelastic isotropic media with an irregular free surface.

A constant fractional-order viscoelastic wave equation and its numerical simulation scheme

Efficient modeling schemes currently exist to handle the spatially variable-order fractional Laplacians in the fractional Laplacian viscoacoustic wave equation. The simplest approach is to change the spatially variable-order fractional Laplacians into a linear combination of several constant fractional-order Laplacians. We generalize the constant fractional-order scheme to a spatially variable fractional-order viscoelastic wave equation and develop an almost-equivalent constant fractional-order viscoelastic wave equation. Our constant fractional-order scheme avoids the simulation error introduced by directly averaging the spatially varying fractional order; thus, our scheme simulates seismic wave propagation in viscoelastic media with sharp [Formula: see text] contrasts well. The fast Fourier transform is used in the approximation of the fractional Laplacians, which improves the spectral accuracy in space. Several simulation examples are performed to verify that the numerical solution of a homogeneous [Formula: see text] model obtained by solving our constant fractional-order viscoelastic wave equation agrees well with that obtained by solving the original viscoelastic wave equation. The numerical simulations for spatially varying [Formula: see text] models obtained by the new wave equation are more straightforward than those currently in use and match the reference solutions obtained by accurate, but inefficient, methods. This match of simulation results verifies the accuracy of our viscoelastic wave equation.

Stability for the damped wave equation with neutral delay

Of concern is a wave equation with a distributed neutral delay. We prove that, despite the destructive nature of delays in general, solutions may go back to the equilibrium state in an exponential manner as time goes to infinity. Reasonable conditions on the distributed neutral delay are established. This type of problems appear in the study of wave propagation in viscoelastic media and in acoustic wave propagation. It is not well studied so far.

Coherent wave propagation in viscoelastic media with mode conversions and pair-correlated scatterers

The influence of correlation between scatterers on coherent waves propagation is studied in the case of a viscoelastic medium hosting a random configuration of either spherical or cylindrical scatterers. A distinction is made between the hole correction and the additional disturbances to the pair correlation function beyond the excluded volume via a radial and concentration dependent Ursell function. The effect of the Ursell function on the effective wavenumber is shown to be of order 3 in concentration and order 2 in scattering, and the corresponding formulas generalize those of Caleap et al. (2012) for an ideal fluid host medium. The whole order 3 in concentration is calculated; its other part is of order 3 in scattering. Both parts of the order 3 in concentration are the sum of two terms, one related to mode conversions, the other not. The numerical study is performed mostly for aluminum spheres in epoxy, which is a rather illustrative situation of the different phenomena that participate to the coherent propagation. The Ursell function effect is enhanced at low frequency, while counteracted partly at higher frequency, by the other term of order 3 in concentration. The most visible effects of both terms are on the attenuation. The Ursell term related to mode conversions is larger than the one with no mode conversions included in the low frequency regime.

Blow-up of a hyperbolic equation of viscoelasticity with supercritical nonlinearities

We investigate a hyperbolic PDE, modeling wave propagation in viscoelastic media, under the influence of a linear memory term of Boltzmann type, and a nonlinear damping modeling friction, as well as an energy-amplifying supercritical nonlinear source:{utt−k(0)Δu−∫0∞k′(s)Δu(t−s)ds+|ut|m−1ut=|u|p−1u, in Ω×(0,T),u(x,t)=u0(x,t), in Ω×(−∞,0], where Ω is a bounded domain in R3 with a Dirichlét boundary condition. The relaxation kernel k is monotone decreasing and k(∞)=1. We study blow-up of solutions when the source is stronger than dissipations, i.e., p>max⁡{m,k(0)}, under two different scenarios: first, the total energy is negative, and the second, the total energy is positive with sufficiently large quadratic energy. This manuscript is a follow-up work of the paper [30] in which Hadamard well-posedness of this equation has been established in the finite energy space. The model under consideration features a supercritical source and a linear memory that accounts for the full past history as time goes to −∞, which is distinct from other relevant models studied in the literature which usually involve subcritical sources and a finite-time memory.

Statistical correlation of fractional oscillator response by complex spectral moments and state variable expansion

The statistical characterization of the oscillator response with non-integer order damping under Gaussian noise represents an important challenge in the modern stochastic mechanics. In fact, this kind of problem appears in several issues of different type (wave propagation in viscoelastic media, Brownian motion, fluid dynamics, RLC circuit, etc.). The aim of this paper is to provide a stochastic characterization of the stationary response of linear fractional oscillator forced by normal white noise. In particular, this paper shows a new method to obtain the correlation function by exact complex spectral moments. These complex quantities contain all the information to describe the random processes but in the considered case their analytical evaluation needs some mathematical manipulations. For this reason such complex spectral moment characterization is used in conjunction with a fractional-order state variable analysis. This kind of analysis permits to find the exact expression of complex spectral moments, and the correlation function by using the Mellin transform. Moreover, the proposed approach provides an analytical expression of the response variance of the fractional oscillator. Capability and efficiency of the present method are shown in the numerical examples in which correlation and variance of fractional oscillator response are found and compared with those obtained by Monte Carlo simulations.