3D elastic waveform modeling with an optimized equivalent staggered-grid finite-difference method

Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.

A staggered-grid high-order finite-difference modeling for elastic wave field in arbitrary tilt anisotropic media

The paper presents a staggered-grid any even-order accurate finite-difference scheme for two-dimensional (2D), three-component (3C), first-order stress-velocity elastic wave equation and its stability condition in the arbitrary tilt anisotropic media; and derives a perfectly matched absorbing layer (PML) boundary condition and its stag- gered-grid any even-order accurate difference scheme in the 2D arbitrary tilt anisotropic media. The results of nu- merical modeling indicate that the modeling precision is high, the calculation efficiency is satisfactory and the absorbing boundary condition is better. The wave-front shapes of elastic waves are complex in the anisotropic media, and the velocity of qP wave is not always faster than that of qS wave. The wave-front triplication of qS wave and its events in both reflected domain and propagated domain, which are not commonly hyperbola, is a common phenomenon. When the symmetry axis is tilted in the TI media, the phenomenon of S-wave splitting is clearly observed in the snaps of three components and synthetic seismograms, and the events of all kinds of waves are asymmetric.

Stability of High-Order Staggered-Grid Schemes for 3D Elastic Wave Equation in Heterogeneous Media

In this paper, we firstly derive the stability conditions of high-order staggered-grid schemes for the three-dimensional (3D) elastic wave equation in heterogeneous media based on the energy method. Moreover, the plane wave analysis yields a sufficient and necessary stability condition by the von Neumann criterion in homogeneous case. Numerical computations for 3D wave simulation with point source excitation are given.

Two exact first-order k-space formulations for low-rank viscoacoustic wave propagation on staggered grids

Wave propagation in the viscoacoustic media is physically dispersive and dissipated.Completely excluding the numerical dispersion error from the physical dispersion in the viscoacoustic wave simu-lation is indispensable to understanding the intrinsic property of the wave propagation in attenuated media for the petroleum exploration geophysics.In recent years,a viscoacoustic wave equation char-acterized by fractional Laplacian gains wide attention in geophysical community.However,the first-order form of the viscoacoustic wave equation,often solved by the conventional staggered-grid pseu-dospectral method,suffers from the numerical dispersion error in time due to the low-order finite-difference approximation.It is challenging to completely eliminate the error because the viscoacoustic wave equation contains two temporal derivatives,which stem from the time stepping and the amplitude attenuation terms,respectively.To tackle the issue,we derive two exact first-order k-space viscoacoustic formulations that can fully exclude the numerical error from the physical dispersion.For the homoge-neous case,two formulations agree with the viscoacoustic analytical solution very well and have the same efficiency.For the heterogeneous case,our second k-space formulation is more efficient than the first one because the second formulation significantly reduces the number of the wavenumber-space mixed-domain operators,which are the expensive part of the viscoacoustic k-space simulation.Nu-merical cases validate that the two first-order k-space formulations are effective and efficient alternatives to the current staggered-grid pseudospectral formulation for the viscoacoustic wave simulation.

Track-before-detect for bistatic radar based on pseudo-spectrum

Track-Before-Detect(TBD) is an efficient method to detect dim targets for radars. Conventional TBD usually follows an approximate motion model of the target, which may cause an inaccurate integration of the target energy. A TBD technique on basis of pseudo-spectrum in mixed coordinates adopting an accurate motion model for bistatic radar system is developed in this paper.The predicted position in bistatic polar plane is derived according to a nonlinear function that exactly describes the constant Cartesian velocity motion. Then around the predicted position, a pseudo-spectrum is formulated and its samples are accumulated to the integration frame for energy integration. The evolution of the state and the procedure of accumulation of the target energy are derived elaborately. The superior performance of the proposed method is demonstrated by some simulations.

A stable staggered-grid finite-difference scheme for acoustic modeling beyond conventional stability limit

Staggered-grid finite-difference(SGFD)schemes have been widely used in acoustic wave modeling for geophysical problems.Many improved methods are proposed to enhance the accuracy of numerical modeling.However,these methods are inevitably limited by the maximum Courant-Friedrichs-Lewy(CFL)numbers,making them unstable when modeling with large time sampling intervals or small grid spacings.To solve this problem,we extend a stable SGFD scheme by controlling SGFD dispersion relations and maximizing the maximum CFL numbers.First,to improve modeling stability,we minimize the error between the FD dispersion relation and the exact relation in the given wave-number region,and make the FD dispersion approach a given function outside the given wave-number area,thus breaking the conventional limits of the maximum CFL number.Second,to obtain high modeling accuracy,we use the SGFD scheme based on the Remez algorithm to compute the FD coefficients.In addition,the hybrid absorbing boundary condition is adopted to suppress boundary reflections and we find a suitable weighting coefficient for the proposed scheme.Theoretical derivation and numerical modeling demonstrate that the proposed scheme can maintain high accuracy in the modeling process and the value of the maximum CFL number of the proposed scheme can exceed that of the conventional SGFD scheme when adopting a small maximum effective wavenumber,indicating that the proposed scheme improves stability during the modeling.

Theoretical analysis and numerical simulation of acoustic waves in gas hydrate-bearing sediments

Based on Carcione-Leclaire model,the time-splitting high-order staggered-grid finite-difference algorithm is proposed and constructed for understanding wave propagation mechanisms in gas hydrate-bearing sediments.Three compressional waves and two shear waves,as well as their energy distributions are investigated in detail.In particular,the influences of the friction coefficient between solid grains and gas hydrate and the viscosity of pore fluid on wave propagation are analyzed.The results show that our proposed numerical simulation algorithm proposed in this paper can effectively solve the problem of stiffness in the velocity-stress equations and suppress the grid dispersion,resulting in higher accuracy compared with the result of the Fourier pseudospectral method used by Carcione.The excitation mechanisms of the five wave modes are clearly revealed by the results of simulations.Besides,it is pointed that,the wave diffusion of the second kind of compressional and shear waves is influenced by the friction coefficient between solid grains and gas hydrate,while the diffusion of the third compressional wave is controlled by the fluid viscosity.Finally,two fluid-solid(gas-hydrate formation)models are constructed to study the mode conversion of various waves.The results show that the reflection,transmission,and transformation of various waves occur on the interface,forming a very complicated wave field,and the energy distribution of various converted waves in different phases is different.It is demonstrated from our studies that,the unconventional waves,such as the second and third kinds of compressional waves may be converted into conventional waves on an interface.These propagation mechanisms provide a concrete wave attenuation explanation in inhomogeneous media.

Multiframe weak target track-before-detect based on pseudo-spectrum in mixed coordinates

Traditional multiframe Track-Before-Detect(TBD)may incur adverse integration loss resulting from model mismatch in sensor coordinates.Its suboptimal integration strategy may cause target envelope degradation.To address these issues,a pseudo-spectrum-based multiframe TBD in mixed coordinates is proposed firstly.The data search for energy integration is conducted based on an accurate model in the x-y plane while target energy is integrated based on pseudo-spectrum in sensor coordinates.The algorithm performance is improved since the model mismatch is eliminated,and the pseudo-spectrum based integration facilitates well maintained target envelope.The detailed multiframe integration procedure and theoretical target integrated envelope are derived.Secondly,to cope with the unknown target velocity,a velocity filter bank based on pseudo-spectrum in mixed coordinates is proposed.The effect of velocity mismatch on algorithm performance is analyzed and an efficient method for filter bank design is presented.Thirdly,a parameter estimation method using characteristics of integrated envelope is presented for improved target polar position and Cartesian velocity estimation.Finally,numerical results are provided to demonstrate the effectiveness of the proposed method.

On the Volume Conservation of the Immersed Boundary Method

The immersed boundary(IB)method is an approach to problems of fluid-structure interaction in which an elastic structure is immersed in a viscous incompressible fluid.The IB formulation of such problems uses a Lagrangian description of the structure and an Eulerian description of the fluid.It is well known that some versions of the IB method can suffer from poor volume conservation.Methods have been introduced to improve the volume-conservation properties of the IB method,but they either have been fairly specialized,or have used complex,nonstandard Eulerian finite-difference discretizations.In this paper,we use quasi-static and dynamic benchmark problems to investigate the effect of the choice of Eulerian discretization on the volume-conservation properties of a formally second-order accurate IB method.We consider both collocated and staggered-grid discretization methods.For the tests considered herein,the staggered-grid IB scheme generally yields at least a modest improvement in volume conservation when compared to cell-centered methods,and in many cases considered in this work,the spurious volume changes exhibited by the staggered-grid IB method are more than an order of magnitude smaller than those of the collocated schemes.We also compare the performance of cell-centered schemes that use either exact or approximate projection methods.We find that the volumeconservation properties of approximate projection IB methods depend strongly on the formulation of the projection method.When used with the IB method,we find that pressure-free approximate projection methods can yield extremely poor volume conservation,whereas pressure-increment approximate projection methods yield volume conservation that is nearly identical to that of a cell-centered exact projection method.