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...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.展开更多
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, t...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.展开更多
By the standard theory,the stable Qk+1,k−Qk,k+1/Qdck divergence-free element converges with the optimal order of approximation for the Stokes equations,but only order k for the velocity in H1-norm and the pressure in...By the standard theory,the stable Qk+1,k−Qk,k+1/Qdck divergence-free element converges with the optimal order of approximation for the Stokes equations,but only order k for the velocity in H1-norm and the pressure in L2-norm.This is due to one polynomial degree less in y direction for the first component of velocity,which is a Qk+1,k polynomial of x and y.In this manuscript,we will show by supercloseness of the divergence free element that the order of convergence is truly k+1,for both velocity and pressure.For special solutions(if the interpolation is also divergence-free),a two-order supercloseness is shown to exist.Numerical tests are provided confirming the accuracy of the theory.展开更多
This paper presents a finite-difference(FD)method with spatially non-rectangular irregular grids to simulate the elastic wave propagation.Staggered irregular grid finite difference oper- ators with a second-order time...This paper presents a finite-difference(FD)method with spatially non-rectangular irregular grids to simulate the elastic wave propagation.Staggered irregular grid finite difference oper- ators with a second-order time and spatial accuracy are used to approximate the velocity-stress elastic wave equations.This method is very simple and the cost of computing time is not much.Complicated geometries like curved thin layers,cased borehole and nonplanar interfaces may be treated with non- rectangular irregular grids in a more flexible way.Unlike the multi-grid scheme,this method requires no interpolation between the fine and coarse grids and all grids are computed at the same spatial iteration.Compared with the rectangular irregular grid FD,the spurious diffractions from'staircase' interfaces can easily be eliminated without using finer grids.Dispersion and stability conditions of the proposed method can be established in a similar form as for the rectangular irregular grid scheme.The Higdon's absorbing boundary condition is adopted to eliminate boundary reflections.Numerical simu- lations show that this method has satisfactory stability and accuracy in simulating wave propagation near rough solid-fluid interfaces.The computation costs are less than those using a regular grid and rectangular grid FD method.展开更多
基金National Natural Science Foundation of China (49894190-024) and Geophysical Prospecting Key Laboratory Foun- dation of China National Petroleum Corporation.
文摘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.
基金National Natural Science Foundation of China (49894190-024) and Geophysical Prospecting Key Laboratory Foundation of China National Petroleum Corporation.
文摘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.
文摘By the standard theory,the stable Qk+1,k−Qk,k+1/Qdck divergence-free element converges with the optimal order of approximation for the Stokes equations,but only order k for the velocity in H1-norm and the pressure in L2-norm.This is due to one polynomial degree less in y direction for the first component of velocity,which is a Qk+1,k polynomial of x and y.In this manuscript,we will show by supercloseness of the divergence free element that the order of convergence is truly k+1,for both velocity and pressure.For special solutions(if the interpolation is also divergence-free),a two-order supercloseness is shown to exist.Numerical tests are provided confirming the accuracy of the theory.
文摘This paper presents a finite-difference(FD)method with spatially non-rectangular irregular grids to simulate the elastic wave propagation.Staggered irregular grid finite difference oper- ators with a second-order time and spatial accuracy are used to approximate the velocity-stress elastic wave equations.This method is very simple and the cost of computing time is not much.Complicated geometries like curved thin layers,cased borehole and nonplanar interfaces may be treated with non- rectangular irregular grids in a more flexible way.Unlike the multi-grid scheme,this method requires no interpolation between the fine and coarse grids and all grids are computed at the same spatial iteration.Compared with the rectangular irregular grid FD,the spurious diffractions from'staircase' interfaces can easily be eliminated without using finer grids.Dispersion and stability conditions of the proposed method can be established in a similar form as for the rectangular irregular grid scheme.The Higdon's absorbing boundary condition is adopted to eliminate boundary reflections.Numerical simu- lations show that this method has satisfactory stability and accuracy in simulating wave propagation near rough solid-fluid interfaces.The computation costs are less than those using a regular grid and rectangular grid FD method.
