In conventional fi nite diff erence numerical simulation of seismic waves,regular grids in Cartesian coordinates are used to divide the calculated region.When simulating seismic wave fi elds under an irregular surface...In conventional fi nite diff erence numerical simulation of seismic waves,regular grids in Cartesian coordinates are used to divide the calculated region.When simulating seismic wave fi elds under an irregular surface,such grids are unsuitable to realize the free boundary condition.They also easily generate false scattered waves at the corners of the grids owing to the approximation of the stepped grids.These issues affect the simulation accuracy.This study introduces an orthogonal body-fitted grid generation technique in computational fl uid dynamics for generating grids in transversely isotropic(TI)media under an irregular surface.The fi rst-order velocity-stress equation in curvilinear coordinates is calculated using the optimized nonstaggered grids finite difference method.The point oscillation generated by the nonstaggered grids difference is eliminated by selective filtering.The orthogonal body-fitted grids can accurately describe the irregular surface.Further,the orthogonality of the grids allows the implementation of free boundary conditions without complicated coordinate transformation and interpolation operations.Numerical examples show that the numerical solutions obtained by this method agree well with the analytical solutions.By comparing the simulation results of the proposed method with those of the regular grid difference method,the proposed method can eff ectively eliminate the false scattered waves caused by the stepped grids under the condition of the same grid spacing.Thus,the accuracy of the numerical simulation is improved.In addition,the simulation results of the three-layer TI media model on an irregular surface show that the proposed method is also suitable for complex models.展开更多
基金supported by the National Key Research and Development Program of China (Grant No.2023YFC3206501 and 2022YFFO802600)the National Natural Science Foundation of China (Grant No.52369003,42262010 and 42374166)+6 种基金the Natural Science Foundation of Inner Mongolia Autonomous Region of China (Grant No.2023LHMS04011 and2022MS04009)the Application Technology Research and Development Project of Jungar Banner (Grant No.2023YY-18 and 2023YY-19)the First-class Academic Subjects Special Research Project of the Education Department of Inner Mongolia Autonomous Region (Grant No.YLXKZX-NND-010)the Inner Mongolia Autonomous Region Science and Technology Leading Talent Team (Grant No.2022LJRC0007)the Inner Mongolia Agricultural University Basic Research Project(BR22-12-04)the Program for Innovative Research Team in Universities of Inner Mongolia Autonomous Region (Grant No.NMGIRT2313)the Basic Scientific Research Project of Institutions of Higher(Grant No.JY20230090)。
文摘In conventional fi nite diff erence numerical simulation of seismic waves,regular grids in Cartesian coordinates are used to divide the calculated region.When simulating seismic wave fi elds under an irregular surface,such grids are unsuitable to realize the free boundary condition.They also easily generate false scattered waves at the corners of the grids owing to the approximation of the stepped grids.These issues affect the simulation accuracy.This study introduces an orthogonal body-fitted grid generation technique in computational fl uid dynamics for generating grids in transversely isotropic(TI)media under an irregular surface.The fi rst-order velocity-stress equation in curvilinear coordinates is calculated using the optimized nonstaggered grids finite difference method.The point oscillation generated by the nonstaggered grids difference is eliminated by selective filtering.The orthogonal body-fitted grids can accurately describe the irregular surface.Further,the orthogonality of the grids allows the implementation of free boundary conditions without complicated coordinate transformation and interpolation operations.Numerical examples show that the numerical solutions obtained by this method agree well with the analytical solutions.By comparing the simulation results of the proposed method with those of the regular grid difference method,the proposed method can eff ectively eliminate the false scattered waves caused by the stepped grids under the condition of the same grid spacing.Thus,the accuracy of the numerical simulation is improved.In addition,the simulation results of the three-layer TI media model on an irregular surface show that the proposed method is also suitable for complex models.
