共查询到20条相似文献，搜索用时 31 毫秒
1.
Reverse-time migration (RTM) directly solves the two-way wave equation for wavefield propagation; therefore, how to solve the wave equation accurately and quickly is very important for RTM. The conventional staggered-grid finite-difference (SFD) operators are usually based on the Taylor-series expansion theory. If they are used to solve wave equation on a larger frequency content, a strong dispersion will occur, which directly affects the seismic image quality. In this paper, we propose an optimal SFD operator based on least squares to solve acoustic wave equation for prestack RTM, and obtain a new antidispersion RTM algorithm that can use short spatial difference operators. The synthetic and real data tests demonstrate that the least squares SFD (LSSFD) operator can mitigate the numerical dispersion, and the acoustic RTM using the LSSFD operator can effectively improve image quality comparing with that using the Taylor-series expansion SFD (TESFD) operator. Moreover, the LSSFD method can adopt a shorter spatial difference operator to reduce the computing cost. 相似文献
2.
3.
Zhang Zhifu Meng Xiaohong Liu Hong Chen Jingbo Li YouMing 《应用地球物理》2006,3(4):218-225
An accurate and wide-angle one-way propagator for wavefield extrapolation is an important topic for research on wave-equation prestack depth migration in the presence of large and rapid velocity variations. Based on the optimal separable approximation presented in this paper, the mixed domain algorithm with forward and inverse Fourier transforms is used to construct the 3D one-way wavefield extrapolation operator. This operator separates variables in the wavenumber and spatial domains. The phase shift operation is implemented in the wavenumber domain while the time delay for lateral velocity variation is corrected in the spatial domain. The impulse responses of the one-way wave operator show that the numeric computation is consistent with the theoretical value for each velocity, revealing that the operator constructed with the optimal separable approximation can be applied to lateral velocity variations for the case of small steps. Imaging results of the SEG/EAGE model and field data indicate that the new method can be used to image complex structure. 相似文献
4.
DAN KOSLOFF ANIBAL QUEIROZ FILHO EKKEHART TESSMER ALFRED BEHLE 《Geophysical Prospecting》1989,37(4):383-394
We present a new rapid expansion method (REM) for the time integration of the acoustic wave equation and the equations of dynamic elasticity in two spatial dimensions. The method is applicable to spatial grid methods such as finite differences, finite elements or the Fourier method. It is based on a Chebyshev expansion of the formal solution to the appropriate wave equation written in operator form. The method yields machine accuracy yet it is faster than methods based on temporal differencing. Its disadvantages are that it does not apply to all types of material rheology, and it can also require much storage when many snapshots and time sections are desired. Comparisons between numerical and analytical solutions for simple acoustic and elastic problems demonstrate the high accuracy of the REM. 相似文献
5.
Staggered-grid finite-difference (SGFD) schemes have been used widely in seismic modeling. The spatial difference coefficients of the SGFD scheme are generally determined by a Taylor-series expansion (TE) method or optimization methods. However, high accuracy is hardly guaranteed both at small and large wavenumbers by using these conventional methods. We propose a new optimal SGFD scheme based on combining TE and minimax approximation (MA) for high accuracy modeling. The optimal spatial SGFD coefficients are calculated by applying a combination of TE and MA to the dispersion relation, where the implementation of the MA method is based on a Remez algorithm. We adopt the optimal SGFD coefficients to solve first-order spatial derivatives of the elastic wave equations and then perform numerical modeling. Dispersion analyses and seismic modeling show the advantage of the proposed optimal method. The optimal SGFD scheme has greater accuracy than the TE-based SGFD scheme for the same spatial difference operator length. In addition, the optimal SGFD scheme can also adopt a shorter operator length to achieve the high accuracy reducing the computational cost. 相似文献
6.
M. D. Ruiz-Medina J. M. Angulo R. Fernández-Pascual 《Stochastic Environmental Research and Risk Assessment (SERRA)》2007,21(3):273-281
A wavelet-based orthogonal decomposition of the solution to stochastic differential/pseudodifferential equations of parabolic
type is derived in the cases of random initial conditions and random forcing. The family of spatiotemporal models considered
can represent anomalous diffusion processes when the spatial operator involved is a fractional or multifractional pseudodifferential
operator. The results obtained are applied to the generation of the sample paths of Gaussian spatiotemporal random fields
in the family studied. 相似文献
7.
O. HOLBERG 《Geophysical Prospecting》1987,35(6):629-655
Conventional finite-difference operators for numerical differentiation become progressively inaccurate at higher frequencies and therefore require very fine computational grids. This problem is avoided when the derivatives are computed by multiplication in the Fourier domain. However, because matrix transpositions are involved, efficient application of this method is restricted to computational environments where the complete data volume required by each computational step can be kept in random access memory. To circumvent these problems a generalized numerical dispersion analysis for wave equation computations is developed. Operators for spatial differentiation can then be designed by minimizing the corresponding peak relative error in group velocity within a spatial frequency band. For specified levels of maximum relative error in group velocity ranging from 0.03% to 3%, differentiators have been designed that have the largest possible bandwidth for a given operator length. The relation between operator length and the required number of grid points per shortest wavelength, for a required accuracy, provides a useful starting point for the design of cost-effective numerical schemes. To illustrate this, different alternatives for numerical simulation of the time evolution of acoustic waves in three-dimensional inhomogeneous media are investigated. It is demonstrated that algorithms can be implemented that require fewer arithmetic and I/O operations by orders of magnitude compared to conventional second-order finite-difference schemes to yield results with a specified minimum accuracy. 相似文献
8.
Vladislav Zheligovsky 《地球物理与天体物理流体动力学》2013,107(5):489-540
We consider stability of regimes of hydromagnetic thermal convection in a rotating horizontal layer with free electrically-conducting boundaries, to perturbations involving large spatial and temporal scales. Equations governing the evolution of weakly nonlinear mean perturbations are derived under the assumption that the α-effect is insignificant in the leading-order (e.g. due to a symmetry of the system). The mean-field equations generalise the standard equations of hydromagnetic convection: New terms emerge – a second-order linear operator representing the combined eddy diffusivity and quadratic terms associated with the eddy advection. If the perturbed CHM regime is nonsteady and insignificance of the α-effect in the system does not rely on the presence of a spatial symmetry, the combined eddy diffusivity operator also involves a nonlocal pseudodifferential operator. If the perturbed CHM state is almost symmetric, α-effect terms appear in the mean-field equations as well. Near a point of a symmetry-breaking bifurcation, cubic nonlinearity emerges in the equations. All the new terms are in general anisotropic. A method for evaluation of their coefficients is presented; it requires solution of a significantly smaller number of auxiliary problems than in a straightforward approach. 相似文献
9.
A new time integration technique for use in forward modelling programmes is introduced. The technique presents an alternative to second-order temporal differencing. It is based on a Chebyshev expansion of the formal evolution operator to the spatially discretized wave equation. The computational effort in forward modelling based on the new technique is about the same as in methods based on temporal differencing. However, machine accuracy can be obtained. The implementation of the technique to solve the acoustic wave equation in two spatial dimensions is described. Finally, an example of using the technique to solve a problem of wave propagation in a single layer is presented. 相似文献
10.
伪谱法是一种高效、高精度计算非均匀介质地震波传播的数值方法,但是由于它的微分算子的全局性,使得该方法不适用于分散内存的并行计算.本文将有限差分算子的局部性和伪谱法算子的高效、高精度相结合,发展基于两种方法的伪谱/有限差分混合方法.该方法在一个空间坐标方向上利用交错网格高阶有限差分算子,在另外的空间坐标方向上利用交错网格伪谱法算子,既保留了后者的高效、高精度优势,又便于在PC集群上实现并行计算.对二维模型的计算显示,混合方法能有效处理介质不连续面,在保证伪谱法计算精度的情况下,提供了一种并行计算的可能途径. 相似文献
11.
A short convolutional differentiator (CD) for computing second spatial derivatives in the acoustic wave equation is presented. This differentiator is obtained by tapering the inverse Fourier transform of the band-limited Fourier spectrum of the second-derivative operator. This new filter has been applied to seismogram computations for inhomogeneous media and results are compared with the conventional high-order finite-difference (FD) and Fourier schemes. The operator can be progressively shortened at the model edges to reduce boundary artefacts. The CD method is superior to the conventional FD operator and comparable with the Fourier method in accuracy but faster to run. A strategy to reduce computation time by 20%, which exploits the localized nature of the operator, is given. The method is illustrated using simple 2D models. 相似文献
12.
本文在原有位场分层分离技术的基础上,提出了动态改进型插值切割算子.通过理想模型实验和谱分析,证明了改进算子对异常源产生的重力场有更好的分离效果.使用新方法对安徽省五河地区布格重力异常进行了分离和分析,分离后的浅部地球物理特征与地表地质调查结果的吻合度高,从而验证了改进算子的有效性和准确性.在结合基础地质资料和分离后的地球物理特征的基础上,进一步分析了五河地区构造格架、区域岩浆岩及红层盆地的空间分布和形态特征,并获取了对研究区内郯庐断裂系各条断裂的延伸、形态和相互关系的新认识. 相似文献
13.
《Advances in water resources》1988,11(3):123-126
Operator-splitting techniques are applied to convective-diffusive transport problems in porous media. The convection is treated by applying a modified method of characteristics to time-step along the characteristics of the convective part of the flow. The nonsymmetry in the spatial operator is addressed via a Petrov-Galerkin method which uses a test function to achieve stability through a balancing of the remaining convection, the diffusion, and any possible reaction terms. The use of time-stepping along characteristics allows the use of large time-steps in a stable but accurate fashion. If local phenomena are important, self-adaptive local grid refinement techniques can be coupled with the operator splitting. 相似文献
14.
M. D. Ruiz-Medina J. M. Angulo 《Stochastic Environmental Research and Risk Assessment (SERRA)》2002,16(4):241-266
In this paper, a class of spatio-temporal processes with first-order autoregressive temporal structure and functional spatio-temporal
interaction is introduced. The spatial second-order regularity is allowed to change over time and is characterized in terms
of fractional Sobolev spaces. The associated filtering problem is considered, assuming that observations are defined by spatial
linear functionals of the process of interest, being affected by additive noise. Conditions under which a stable solution
to this problem is obtained are studied. A functional least-squares linear estimate fusion method is derived to calculate
this solution A multiscale finite-dimensional approximation to the problem is obtained from the wavelet-based orthogonal expansions
of the time cross-section spatial processes, which allows the numerical inversion of the linear operator involved. 相似文献
15.
SAMUEL H. GRAY 《Geophysical Prospecting》1992,40(5):565-571
Integral migration techniques perform a sum over an aperture of input traces to obtain output at a single point. The length of the aperture is limited by a spatial Nyquist criterion, which typically prohibits imaging very steep dips at very high frequencies without generating severe migration artifacts (migration operator aliasing). For time-domain Kirchhoff migration, this can be a fatal shortcoming. The standard way to address this problem is to interpolate traces spatially before migration. This reduces the trace spacing, thereby increasing the frequency content which can be migrated without aliasing at steep dips. An alternative remedy to the operator aliasing problem is to modify the phase response of the Kirchhoff migration operator. This operator is frequency-selective across the migration aperture: it passes all temporal frequencies of the input traces in the innermost portion of the aperture (referring to the shallow dips), and gradually cuts out the higher frequencies as it approaches the outer portion of the aperture. Thus, while all frequencies of the input data contribute to the shallow-dip portion of the migrated image, only the permissible low frequencies of the input data contribute to imaging the steepest dips. Using a simple realization of a frequency-selective Kirchhoff migration operator, this technique is illustrated on a synthetic data set involving greater than vertical dips. 相似文献
16.
V. A. Zheligovsky 《Izvestiya Physics of the Solid Earth》2006,42(12):1051-1067
The problem of weakly nonlinear stability with respect to large-scale perturbations in 3-D convective magnetohydrodynamic (MHD) states in which the α-effect is absent or insignificant (e.g., because the system has symmetry relative to a center or a vertical axis) is examined. It is assumed that the MHD state whose stability is studied is free from large spatiotemporal scales and is insensitive to perturbations with the same small spatial scale as in the state under study. The equations for mean perturbation fields derived by asymptotic methods generalize the standard equations of magnetohydrodynamics (the Navier-Stokes and magnetic induction equations). A combined eddy diffusion operator, generally anisotropic and not necessarily negative definite, and additional quadratic terms similar to advective terms arise in the inferred generalized equations. 相似文献
17.
在数值模拟中,隐式有限差分具有较高的精度和稳定性.然而,传统隐式有限差分算法大多由于需要求解大型矩阵方程而存在计算效率偏低的局限性.本文针对一阶速度-应力弹性波方程,构建了一种优化隐式交错网格有限差分格式,然后将改进格式由时间-空间域转换为时间-波数域,利用二范数原理建立目标函数,再利用模拟退火法求取优化系数.通过对均匀模型以及复杂介质模型进行一阶速度-应力弹性波方程数值模拟所得单炮记录、波场快照分析表明:这种优化隐式交错网格差分算法与传统的几种显式和隐式交错网格有限差分算法相比不但降低了计算量,而且能有效的压制网格频散,使弹性波数值模拟的精度得到有效的提高. 相似文献
18.
19.
In this paper, we develop a new nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique for numerically solving two-dimensional acoustic wave equations. We first split two-dimensional acoustic wave equation into the local one-dimensional equations and transform each of the split equations into a Hamiltonian system. Then, we use both a nearly analytic discrete operator and a central difference operator to approximate the high-order spatial differential operators, which implies the symmetry of the discretized spatial differential operators, and we employ the partitioned second-order symplectic Runge–Kutta method to numerically solve the resulted semi-discrete Hamiltonian ordinary differential equations, which results in fully discretized scheme is symplectic unlike conventional nearly analytic symplectic partitioned Runge–Kutta methods. Theoretical analyses show that the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique exhibits great higher stability limits and less numerical dispersion than the nearly analytic symplectic partitioned Runge–Kutta method. Numerical experiments are conducted to verify advantages of the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique, such as their computational efficiency, stability, numerical dispersion and long-term calculation capability. 相似文献
20.
本文提出了一种全新的基于等效Q值的时域反Q滤波算法,其允许等效Q值在垂向上随时间连续变化,在空间上存在弱变化;将加权最小平方方法优化设计思想引入到时域反Q补偿短算子设计当中,给出最优时域短算子设计,将大量的频率域乘法工作转化为少量的时域褶积运算;采取表驱动方案,将短算子的构建与反Q补偿运算相剥离,极大地提升了计算效率;提出了一种新的稳定性控制方法,其既保证算法具有良好的稳定性,又满足短算子设计精度的要求.数值计算表明：时域反Q滤波算法可以取得与频域算法相同的补偿效果,并保证算法具备良好的稳定性和较高的计算效率. 相似文献