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.

An adaptive finite-difference method for seismic traveltime modeling based on 3D eikonal equation

3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.

Viscoacoustic prestack reverse time migration based onthe optimal time-space domain high-order finite-difference method

Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.

Perfect plane-wave source for a high-order symplectic finite-difference time-domain scheme

The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order symplectic finite- difference time-domain （SFDTD） scheme for the first time. By splitting the fields on one-dimensional grid and using the nature of numerical plane-wave in finite-difference time-domain （FDTD）, the identical dispersion relation can be obtained and proved between the one-dimensional and three-dimensional grids. An efficient plane-wave source is simulated on one-dimensional grid and a perfect match can be achieved for a plane-wave propagating at any angle forming an integer grid cell ratio. Numerical simulations show that the method is valid for SFDTD and the residual field in SF region is shrinked down to -300 dB.

Numerical simulation of seismic wavefields in TTI media using the rotated staggered-grid compact finite-difference scheme

Numerical simulation in transverse isotropic media with tilted symmetry axis（TTI） using the standard staggered-grid finite-difference scheme（SSG）results in errors caused by averaging or interpolation. In order to eliminate the errors, a method of rotated staggered-grid finite-difference scheme（RSG） is proposed. However, the RSG brings serious numerical dispersion. The compact staggered-grid finite-difference scheme（CSG） is an implicit difference scheme, which use fewer grid points to suppress dispersion more effectively than the SSG. This paper combines the CSG with the RSG to derive a rotated staggered-grid compact finite-difference scheme（RSGC）. The numerical experiments indicate that the RSGC has weaker numerical dispersion and better accuracy than the RSG.

Finite-difference calculation of traveltimes based on rectangular grid

To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.

Optimization of the seismic processing phase-shift plus finite-difference migration operator based on a hybrid genetic and simulated annealing algorithm

Although the phase-shift seismic processing method has characteristics of high accuracy, good stability, high efficiency, and high-dip imaging, it is not able to adapt to strong lateral velocity variation. To overcome this defect, a finite-difference method in the frequency-space domain is introduced in the migration process, because it can adapt to strong lateral velocity variation and the coefficient is optimized by a hybrid genetic and simulated annealing algorithm. The two measures improve the precision of the approximation dispersion equation. Thus, the imaging effect is improved for areas of high-dip structure and strong lateral velocity variation. The migration imaging of a 2-D SEG/EAGE salt dome model proves that a better imaging effect in these areas is achieved by optimized phase-shift migration operator plus a finite-difference method based on a hybrid genetic and simulated annealing algorithm. The method proposed in this paper is better than conventional methods in imaging of areas of high-dip angle and strong lateral velocity variation.

A Finite-Difference Approach to the Time-Dependent Mild-Slope Equation

A finite-difference approach is used to develop a time-dependent mild-slope equation incorporating the effects of bottom dissipation and nonlinearity. The Enler predietor-corrector method and the three-point finite-difference method with varying spatial steps are adopted to discretize the time derivatives and the two-dimensional horizontal ones, respectively, thus leading both the time and spatial derivatives to the second-order accuracy. The boundary conditions for the present model are treated on the basis of the general conditions for open and fixed boundaries with an arbitrary reflection coefficient and phase shift. Both the linear and nonlinear versions of the numerical model are applied to the wave propagation and transformation over an elliptic shoal on a sloping beach, respectively, and the linear version is applied to the simulation of wave propagation in a fully open rectangular harbor. From comparison of numerical results with theoretical or experimental ones, it is found that they are in reasonable agreement.

Uniform stable conformal convolutional perfectly matched layer for enlarged cell technique conformal finite-difference time-domain method

Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer（CPML） is given for solving the truncation problem of the open port when the enlarged cell technique conformal finite-difference time-domain（ECT-CFDTD） method is used to simulate the wave propagation inside a perfect electric conductor（PEC） waveguide.The algorithm has the same numerical stability as the ECT-CFDTD method.For the long-time propagation problems of an evanescent wave in a waveguide,several numerical simulations are performed to analyze the reflection error by sweeping the constitutive parameters of the integral-based conformal CPML.Our numerical results show that the integral-based conformal CPML can be used to efficiently truncate the open port of the waveguide.

A wavelet finite-difference method for numerical simulation of wave propagation in fluid-saturated porous media

In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.

Poroelastic finite-difference modeling for ultrasonic waves in digital porous cores

Scattering attenuation in short wavelengths has long been interesting to geophysicists. Ultrasonic coda waves, observed as the tail portion of ultrasonic wavetrains in laboratory ultrasonic measurements, are important for such studies where ultrasonic waves interact with smallscale random heterogeneities on a scale of micrometers, but often ignored as noises because of the contamination of boundary reflections from the side ends of a sample core. Numerical simulations with accurate absorbing boundary can provide insight into the effect of boundary reflections on coda waves in laboratory experiments. The simulation of wave propagation in digital and heterogeneous porous cores really challenges numerical techniques by digital image of poroelastic properties, numerical dispersion at high frequency and strong heterogeneity, and accurate absorbing boundary schemes at grazing incidence. To overcome these difficulties, we present a staggered-grid high-order finite-difference （FD） method of Biot＇s poroelastic equations, with an arbitrary even-order （2L） accuracy to simulate ultrasonic wave propagation in digital porous cores with strong heterogeneity. An unsplit convolutional perfectly matched layer （CPML） absorbing boundary, which improves conventional PML methods at grazing incidence with less memory and better computational efficiency, is employed in the simulation to investigate the influence of boundary reflections on ultra- sonic coda waves. Numerical experiments with saturated poroelastic media demonstrate that the 2L FD scheme with the CPML for ultrasonic wave propagation significantly improves stability conditions at strong heterogeneity and absorbing performance at grazing incidence. The boundary reflections from the artificial boundary surrounding the digital core decay fast with the increase of CPML thick- nesses, almost disappearing at the CPML thickness of 15 grids. Comparisons of the resulting ultrasonic coda Qsc values between the numerical and experimental ultrasonic S waveforms for a cylindrical rock sample demonstrate that the boundary reflection may contribute around one-third of the ultrasonic coda attenuation observed in laboratory experiments.

Investigation of three-pulse photon echo in thick crystal using finite-difference time-domain method

This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the three-pulse photon echo＇s amplitude and efficiency is analyzed with the Maxwell-Bloch equations solved by finite-difference timedomain method.We demonstrate that the amplitude of three-pulse echo will increase with the increasing of thickness and the optimum thickness to generate three-pulse photon echo is 0.3 cm for Tm^（3＋）：YAG when the attenuation of the input pulse is taken into account.Meanwhile,we find the expression 0.09 exp（α＇L）,which is previously employed to describe the relationship between echo＇s efficiency and thickness,should be modified as 1.3 · 0.09 exp（2.4 ·α＇L） with the propagation of echo considered.

Acoustic finite-difference modeling beyond conventional Courant-Friedrichs-Lewy stability limit:Approach based on variable-length temporal and spatial operators

Conventional finite-difference(FD)methods cannot model acoustic wave propagation beyond Courant-Friedrichs-Lewy(CFL)numbers 0.707 and 0.577 for two-dimensional(2D)and three-dimensional(3D)equal spacing cases,respectively,thereby limiting time step selection.Based on the definition of temporal and spatial FD operators,we propose a variable-length temporal and spatial operator strategy to model wave propagation beyond those CFL numbers while preserving accuracy.First,to simulate wave propagation beyond the conventional CFL stability limit,the lengths of the temporal operators are modified to exceed the lengths of the spatial operators for high-velocity zones.Second,to preserve the modeling accuracy,the velocity-dependent lengths of the temporal and spatial operators are adaptively varied.The maximum CFL numbers for the proposed method can reach 1.25 and 1.0 in high velocity contrast 2D and 3D simulation examples,respectively.We demonstrate the effectiveness of our method by modeling wave propagation in simple and complex media.

Finite-difference model of land subsidence caused by cluster loads in Zhengzhou,China

Groundwater exploitation has been regarded as the main reason for land subsidence in China and thus receives considerable attention from the government and the academic community.Recently,building loads have been identified as another important factor of land subsidence,but researches in this sector have lagged.The effect of a single building load on land subsidence was neglected in many cases owing to the narrow scope and the limited depth of the additional stress in stratum.However,due to the superposition of stresses between buildings,the additional stress of cluster loads is greater than that of a single building load under the same condition,so that the land subsidence caused by cluster loads cannot be neglected.Taking Shamen village in the north of Zhengzhou,China,as an example,a finite-difference model based on the Biot consolidation theory to calculate the land subsidence caused by cluster loads was established in this paper.Cluster loads present the characteristics of large-area loads,and the land subsidence caused by cluster loads can have multiple primary consolidation processes due to the stress superposition of different buildings was shown by the simulation results.Pore water migration distances are longer when the cluster loads with high plot ratio are imposed,so that consolidation takes longer time.The higher the plot ratio is,the deeper the effective deformation is,and thus the greater the land subsidence is.A higher plot ratio also increases the contribution that the deeper stratigraphic layers make to land subsidence.Contrary to the calculated results of land subsidence caused by cluster loads and groundwater recession,the percentage of settlement caused by cluster loads in the total settlement was 49.43%and 55.06%at two simulated monitoring points,respectively.These data suggest that the cluster loads can be one of the main causes of land subsidence.

The Best Finite-Difference Scheme for the Helmholtz Equation

The best finite-difference scheme for the Helmholtz equation is suggested. A method of solving obtained finite-difference scheme is developed. The efficiency and accuracy of method were tested on several examples.

Optical simulation of in-plane-switching blue phase liquid crystal display using the finite-difference time-domain method

The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell＇s equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.

Optimization of a global seventh-order dissipative compact finite-difference scheme by a genetic algorithm

A global seventh-order dissipative compact finite-difference scheme is optimized in terms of time stability. The dissipative parameters appearing in the boundary closures are assumed to be different, resulting in an optimization problem with several parameters determined by applying a generic algorithm. The optimized schemes are analyzed carefully from the aspects of the eigenvalue distribution, the ε-pseudospectra, the short time behavior, and the Fourier analysis. Numerical experiments for the Euler equations are used to show the effectiveness of the final recommended scheme.

A spherical higher-order finite-difference time-domain algorithm with perfectly matched layer

A higher-order finite-difference time-domain（HO-FDTD） in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spherical resonator is modeled in order to demonstrate the advantage of this scheme over the finite-difference time-domain（FDTD） and the multiresolution time-domain（MRTD） schemes with respect to memory requirements and CPU time. Moreover, the Berenger＇s perfectly matched layer（PML） is derived for the spherical HO-FDTD grids, and the numerical results validate the efficiency of the PML.

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.

Numerical Solution of a Problem of Thermal Stresses of a Magnetothermoelastic Cylinder with Rotation by Finite-Difference Method

The present article deals with the investigation thermal stress of a magnetothermoelastic cylinder subjected to rotation,open or closed circuit,thermal and mechanical boundary conditions.The outer and inner surfaces of the cylinder are subjected to both mechanical and thermal boundary conditions.A The transient coupled thermoelasticity in an infinite cylinder with its base abruptly exposed to a heat flux of a decaying exponential function of time is devised solve by the finite-difference method.The fundamental equations’system is solved by utilizing an implicit finite-difference method.This current method is a second-order accurate in time and space;it is also unconditionally stable.To illustrate the present model’s efficiency,we consider a suitable material and acquire the numerical solution of temperature,displacement components,and the components of stresses with time t and through the radial of an infinite cylinder.The results indicate that the effect of coupled thermoelasticity,magnetic field,and rotation on the temperature,stresses,and displacement is quite pronounced.In order to illustrate and verify the analytical developments,the numerical solution of partial differential equations,stress components,displacement components and temperature is carried out and computer simulated results are presented graphically.This study is helpful in the development of piezoelectric devices.