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.

An improved convolution perfectly matched layer for elastic second-order wave equation

A convolution perfectly matched layer(CPML)can efficiently absorb boundary reflection in numerical simulation.However,the CPML is suitable for the first-order elastic wave equation and is difficult to apply directly to the second-order elastic wave equation.In view of this,based on the first-order CPML absorbing boundary condition,we propose a new CPML(NCPML)boundary which can be directly applied to the second-order wave equation.We first systematically extend the first-order CPML technique into second-order wave equations,neglecting the space-varying characteristics of the partial damping coefficient in the complex-frequency domain,avoiding the generation of convolution in the time domain.We then transform the technique back to the time domain through the inverse Fourier transform.Numerical simulation indicates that the space-varying characteristics of the attenuation factor have little influence on the absorption effect and increase the memory at the same time.A number of numerical examples show that the NCPML proposed in this study is effective in simulating elastic wave propagation,and this algorithm is more efficient and requires less memory allocation than the conventional PML absorbing boundary.

P-and S-wavefield simulations using both the firstand second-order separated wave equations through a high-order staggered grid finite-difference method

In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.

THE PARAMETER PERTURBATION METHOD ON ELASTIC WAVE EQUATION IN INHOMOGENEOUS MEDIUM

In this paper, the medium parameters of the elastic wave equation in inhomogeneous medium are rewritten by introducing the referential variables and the perturbational variables, and the wave equation whose sources are the medium parameter perturbational term in homogeneous medium is obtained. By using the Green function theory, the integral equation of the perturbational parameters is obtained. Then the displacement field in homogeneous medium is considered the result of the first iteration, and the displacement field is solved by this integral equation. When the perturbations of medium parameters are about 50 percent, this method can solve the displacement field effectively. from the analysis of the numerical results, the characteristics of wave field in inhomogeneous medium are obtained. The results conform with the local principles of wave function in inhomogeneous medium.

Finite element equations and numerical simulation of elastic wave propagation in two-phase anisotropic media

Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical solution of finite element equations is given.Finally,Properties of elastic wave propagation are observed and analyzed through FEM modeling.

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.

SCATTERING OF HARMONIC ANTI-PLANE SHEAR WAVES BY AN INTERFACE CRACK IN MAGNETO-ELECTRO-ELASTIC COMPOSITES

The dynamic behavior of an interface crack in magneto-electro-elastic composites under harmonic elastic anti-plane shear waves is investigated for the permeable electric boundary conditions. By using the Fourier transform, the problem can be solved with a pair of dual integral equations in which the unknown variable was the jump of the displacements across the crack surfaces. To solve the dual integral equations, the jump of the displacements across the crack surface was expanded in a series of Jacobi polynomials. Numerical examples were provided to show the effect of the length of the crack, the wave velocity and the circular frequency of the incident wave on the stress, the electric displacement and the magnetic flux intensity factors of the crack. From the results, it can be obtained that the singular stresses in piezoelectric/piezomagnetic materials carry the same forms as those in a general elastic material for anti-plane shear problem.

DYNAMIC BEHAVIOR OF TWO PARALLEL SYMMETRY CRACKS IN MAGNETO-ELECTRO-ELASTIC COMPOSITES UNDER HARMONIC ANTI-PLANE WAVES

The dynamic behavior of two parallel symmetry cracks in magneto-electro-elastic composites under harmonic anti-plane shear waves is studied by Schmidt method. By using the Fourier transform, the problem can be solved with a pair of dual integral equations in which the unknown variable is the jumps of the displacements across the crack surfaces. To solve the dual integral equations, the jumps of the displacements across the crack surface were expanded in a series of Jacobi polynomials. The relations among the electric filed, the magnetic flux and the stress field were obtained. From the results, it can be obtained that the singular stresses in piezoelectric/piezomagnetic materials carry the same forms as those in a general elastic material for the dynamic anti-plane shear fracture problem. The shielding effect of two parallel cracks was also discussed.

RAYLEIGH LAMB WAVES IN MICROPOLAR ISOTROPIC ELASTIC PLATE

The propagation of waves in a homogeneous isotropic micropolar elastic cylindrical plate subjected to stress free conditions is investigated. The secular equations for symmetric and skew symmetric wave mode propagation are derived. At short wave limit, the secular equations for symmetric and skew symmetric waves in a stress free circular plate reduces to Rayleigh surface wave frequency equation. Thin plate results are also obtained. The amplitudes of displacements and microrotation components are obtained and depicted graphically. Some special cases are also deduced from the present investigations. The secular equations for symmetric and skew symmetric modes are also presented graphically.

Theconvolutionaldiferentiatormethodfornu┐mericalmodelingofacousticandelasticwave┐fields

Based on the techniques of forward and inverse Fourier transformation, the authors discussed the design scheme of ordinary differentiator used and applied in the simulation of acoustic and elastic wavefields in isotropic media respectively. To compress Gibbs effects by truncation effectively, Hanning window is introduced in. The model computation shows that, the convolutional differentiator method has the advantages of rapidity, low requirements of computer′s inner storage and high precision, which is a potential method of numerical simulation.

THE SCATTERING OF HARMONIC ELASTIC ANTI-PLANE SHEAR WAVES BY A FINITE CRACK IN INFINITELY LONG STRIP USING THE NON-LOCAL THEORY

In this paper, the scattering of harmonic anti-plane shear wavesby a finite crack in infinitely long strip is studied using thenon-local theory. The Fourier transform is applied and a mixedboundary value problem is formulated. Then a set of dual integralequations is solved using the Schmidt method instead of the first orthe second integral equation method. A one-dimensional non-localkernel is used instead of a two-di- mensional one for the anti-planedynamic problem to obtain the stress occurring at the crack tips.Contrary to the classical elasticity solution, it is found that nostress singularity is present at the crack tip. The non-local dynamicelastic solutions yield a finite hoop stress at the crack tip, thusallowing for a fracture criterion based on the maximum dynamic stresshypothesis. The finite hoop stress at the crack tip depends on thecrack length, the width of the strip and the lattice parameters.

A NUMERICAL METHOD OF MODELING ELASTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA

The velocity-stress finite-difference method is adopted to simulate the elastic wave propa- gation in azimuthal anisotropic media.The difference grids are completely staggered in the numerical im- plementation.To reduce the computational work,the absorbin8 boundary conditions for anisotropic media are introduced first and the corner points are specially treated.Examples show that more accurate results can be obtained from the modeling algorithm,which cost much less computational time than the conven- tional methods.Therefore,the algorithm has broad application prospects in engineering.

Frequency-domain Elastic wave Simulation Based on the Nonoverlapping Domain Decomposition Method

A new wave simulation technique for the elastic wave equation in the frequency domain based on a no overlapping domain decomposition algorithm is investigated. The boundary conditions and the finite difference discrimination of the elastic wave equation are derived. The algorithm of no overlapping domain decomposition method is given. The method solves the elastic wave equation by iteratively solving sub problems defined on smaller sub domains. Numerical computations both for homogeneous and inhomogeneous media show the effectiveness of the proposed method. This method can be used in the full-waveform inversion.

Using Refined Theory to Studied Elastic Wave Scattering and Dynamic Stress Concentrations in Plates with Two Cutouts

In this paper, based on complex variables and conformal mapping methods, using the refined dynamic equation of plates, elastic wave scattering and dynamic stress concentrations in plates with two cutouts were studied. Applying the orthogonal function expansion method, the problem to be solved can be reduced into the solution of a set of infinite algebraic equations. According to free boundary conditions, numerical results of dynamic moment concentration factors in thick plates with two circular cutouts analyze that: there will be more complex interaction changes between two-cutout situation than single cutout situation. In the case of low frequency or high frequency and thin plate, the hole-spacing in the absence of coupling interactions was larger or smaller. The numerical results and method can be used to analyze the dynamics and strength of plate-like structures.

SCATTERING OF ANTI-PLANE SHEAR WAVES BY A SINGLE CRACK IN AN UNBOUNDED TRANSVERSELY ISOTROPIC ELECTRO-MAGNETO-ELASTIC MEDIUM

A theoretical treatment of the scattering of anti-plane shear(SH) waves is provided by a single crack in an unbounded transversely isotropic electro-magneto-elastic medium. Based on the differential equations of equilibrium, electric displacement and magnetic induction intensity differential equations, the governing equations for SH waves were obtained. By means of a linear transform, the governing equations were reduced to one Helmholtz and two Laplace equations. The Cauchy singular integral equations were gained by making use of Fourier transform and adopting electro-magneto impermeable boundary conditions. The closed form expression for the resulting stress intensity factor at the crack was achieved by solving the appropriate singular integral equations using Chebyshev polynomial. Typical examples are provided to show the loading frequency upon the local stress fields around the crack tips. The study reveals the importance of the electro-magneto-mechanical coupling terms upon the resulting dynamic stress intensity factor.

The Quasi-Static Approximation of Heat Waves in Anisotropic Thermo-Elastic Media

The equilibrium equations of anisotropic media, coupled to the heat conduction equations, are studied here based on the standard spaces of the physical presentation, in which an new thermo-elastic model based on the second law of thermodynamics is induced. The uncoupled heat wave equation for anisotropic media is deduced. The results show that the equation of heat wave is of the properties of dissipative waves. In final part of this paper, we discuss the propagation behaviour of heat waves for transversely isotropic media.

INTERACTION OF ELASTIC WAVES WITH A PERIODIC ARRAY OF COLLINEAR INPLANE CRACKS

A boundary integral equation method is applied to the study of the interaction of plane elastic waves with a periodic array of collinear inplane cracks.Numerical results are presented for the dynamic stress in- tensity factors.The effects of the wave type,wave frequency,wave incidence angle,and crack spacing on the dynamic stress intensity factors are analyzed in detail.

Elastic Full Waveform Inversion Based on the Trust Region Strategy

In this paper, we investigate the elastic wave full-waveform inversion (FWI) based on the trust region method. The FWI is an optimization problem of minimizing the misfit between the observed data and simulated data. Usually, the line search method is used to update the model parameters iteratively. The line search method generates a search direction first and then finds a suitable step length along the direction. In the trust region method, it defines a trial step length within a certain neighborhood of the current iterate point and then solves a trust region subproblem. The theoretical methods for the trust region FWI with the Newton type method are described. The algorithms for the truncated Newton method with the line search strategy and for the Gauss-Newton method with the trust region strategy are presented. Numerical computations of FWI for the Marmousi model by the L-BFGS method, the Gauss-Newton method and the truncated Newton method are completed. The comparisons between the line search strategy and the trust region strategy are given and show that the trust region method is more efficient than the line search method and both the Gauss-Newton and truncated Newton methods are more accurate than the L-BFGS method.

New Periodic Wave Solutions and Their Interaction for （2＋1）-dimensional KdV Equation

的新二倍地周期的波浪解决方案的一个类(2+1 ) 维的 KdV 方程被介绍适当 Jacobi 获得椭圆形的函数和 Weierstrass 在一般解决方案的椭圆形的函数(包含二个任意的函数) 借助于 multilinear 变量分离途径得到了为(2+1 ) 维的 KdV 方程。限制盒子被考虑，一些局部性的刺激被导出，例如 dromion， mul-tidromions， dromion-antidromion， multidromions-antidromions 等等。dromion-antidromion 和 multidromions-antidromions 的一些答案在一个方向是周期的但是在另外的方向局部性。这些答案的相互作用性质，数字地被学习，表明他们中的一些是无弹性的，一些是完全有弹性的。而且，这些结果被设想。