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.展开更多
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 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.
基金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.
