基于改进交错网格差分结构的地震波场正演模拟OA北大核心CSTPCD

Finite difference(FD)discretization schemes are widely used in seismic wave numerical simulation.Compared with the conventional temporal second-order FD,temporal high-order FDs are more accurate and sta-ble.Implicit FD is more numerical than the explicit FD.Combining the advantages of temporal high-order meth-od and implicit FD,this paper develops an improved staggered grid FD(SFD)stencil to solve the three-dimen-sional first-order acoustic wave equation.In this approach,temporal high-order discretization is applied for com-pensating the accuracy of the conventional temporal second-order FD,and the spatial derivatives are approxima-ted by the implicit SFD.We further provide Taylor-expansion(TE)-based coefficients for temporal discretization scheme,and estimate implicit spatial coefficients using TE and least-square.Compared with the existing SFD methods,the improved method can adopt different grid intervals and diffrent operator lengths in different spatial directions.While ensuring simulation accuracy,it is more flexible.Numerical analyses and modeling examples have proven that the developed method can generate high-accuracy simulation results with shorter operator lengths.Therefore,it can be regeared as an effective tool for seismic wavefield simulation.