非均匀TI介质P-SV波传播交错网格高阶有限差分数值模拟

给出了在非均匀横向各向同性(TI)介质情况下,四阶时间精度、高阶空间精度的一阶速度-应力P-SV波的波动方程交错网格有限差分解法.首先根据一阶速度(应力)波动方程把速度(应力)对时间的一阶和三阶导数转换为应力(速度)对空间的导数,从而在使用四阶时间精度有限差分格式计算某一时刻的波场时只需要前面两个时间步的波场值;然后在空间上采用高阶有限差分格式以提高数值模拟的精度.数值模拟结果和实测垂直地震剖面(VSP)记录符合得很好,说明该方法是可行的.     

时间-空间高阶精度矩形交错网格隐式有限差分声波正演模拟

有限差分方法因其操作简单、计算消耗低而成为地震勘探领域中最为常用的数值模拟方法之一,然而用离散的显式差分算子数值逼近地震波动方程中的连续导数容易导致数值频散,并且基于正方形网格离散形式的有限差分方法对不同地质模型的适应性较低.针对一阶变密度声波方程的数值模拟,本文发展了一种适用于矩形网格离散形式的时间高阶空间隐式有限差分格式,可以有效压制时间和空间频散,同时灵活的网格剖分增强了其应用的广泛性.基于本文矩形交错网格时间高阶空间隐式有限差分格式的时空域频散关系和变量替换的思想,首先采用泰勒级数展开方法求解不同方向的非轴上时间差分系数及轴上空间差分系数,使本文差分格式可以获得任意偶数阶时间和空间精度.为了进一步提高本文差分格式在更大波数区域的空间模拟精度,我们采用线性优化方法来求取新的轴上空间差分系数用于一阶变密度声波方程的波场迭代求解中.频散、稳定性分析及数值模拟算例表明：相比于传统十字形空间域隐式有限差分格式,本文矩形交错网格时间高阶空间隐式有限差分格式在精度、稳定性和效率方面均具有优势.     

基于新的差分结构的时-空域高阶有限差分波动方程数值模拟方法

传统的高阶有限差分波动方程数值模拟方法采用高阶差分算子近似空间偏导数,能有效抑制空间频散.然而,传统的有限差分法仅采用二阶差分算子近似时间偏导数,这使得地震波场沿时间外推的精度较低.当采用较大的时间采样间隔,传统的有限差分法模拟波场会出现明显的时间频散,甚至不稳定.本文基于新的差分结构和中心网格剖分,发展了一种空间任意偶数阶精度、时间四阶和六阶精度的时空域有限差分方法.基于对离散后的频散关系进行泰勒展开,本文推导了时空域高阶有限差分算子的差分系数.相速度分析表明时间四阶、六阶精度的差分方法能显著地减小传统时间二阶精度差分方法的时间频散.在相同的精度下与传统差分法比较,本文发展的时间四阶、六阶有限差分方法的计算效率比传统方法高.均匀和非匀均介质中的波场数值模拟实验进一步证实本文研究的时空高阶有限差分方法的优越性.     

三维波动方程时空域混合网格有限差分数值模拟方法

常规高阶和时空域高阶有限差分方法广泛应用于三维标量波动方程的数值模拟,这两种差分方法仅利用笛卡尔坐标系中的坐标轴网格点构建三维Laplace差分算子,相应的差分离散波动方程本质上仅具有2阶差分精度,模拟精度低.本文将三维笛卡尔坐标系中非坐标轴网格点分为两类:坐标平面内的非坐标轴网格点和坐标平面外的非坐标轴网格点,系统推导出了两类非坐标轴网格点构建三维Laplace差分算子的方法,进而提出了一种利用坐标轴网格点和非坐标轴网格点共同构建三维Laplace差分算子的混合网格有限差分方法,并利用时空域频散关系和泰勒展开建立差分系数方程,推导出了差分系数的通解.相比常规高阶和时空域高阶差分格式的2阶差分精度,时空域混合网格差分离散波动方程理论上能够达到任意偶数阶差分精度,模拟精度显著提高,同时稳定性更强.频散分析表明:相比常规高阶和时空域高阶差分格式,在计算效率基本相同时,时空域混合网格差分格式能更有效地减小数值频散,减弱数值各向异性,模拟精度更高;在模拟精度基本相当时,混合网格差分格式能采用更大的时间采样间隔,计算效率更高.数值模拟实例进一步验证了混合网格差分格式在提高模拟精度和计算效率方面的...     

二维地震波场的组合型紧致有限差分数值模拟

地震波场数值模拟在地球物理勘探和地震学中具有重要的支撑作用.本文将组合型紧致差分格式用于声波和弹性波方程的数值模拟中.根据泰勒级数展开和声波方程,建立了位移场时间四阶离散格式,并将组合型紧致差分格式用于位移场空间导数的求取,然后对该差分格式进行了精度分析、误差分析、频散分析和稳定性分析.理论研究结果表明：①该差分格式为时间四阶、空间六阶精度,与常规七点六阶中心差分和五点六阶紧致差分相比,具有更小的截断误差和更高的模拟精度;②每个波长仅需要5.6个采样点,且满足稳定性条件的库郎数为0.792,可以使用粗网格和较大时间步长进行计算.所以该方法具有占用内存少、计算效率高和低数值频散等优势.最后,本文进行了二维各向同性完全弹性介质的声波和弹性波方程的数值模拟,实验结果表明本文提出的方法具有更高的计算精度,能够大幅度的节约计算量和内存需求,对于三维大尺度模型问题具有更好的适应性.     

旋转网格和常规网格混合的时空域声波有限差分正演

压制数值频散,提高正演模拟精度,一直是有限差分正演模拟研究的重要内容.基于时空域频散关系的有限差分法,比基于空间域频散关系的传统有限差分法,模拟精度更高.时空域声波方程数值模拟,普遍采用常规十字交叉型高阶有限差分格式.而在频率-空间域,普遍采用旋转网格和常规网格混合的有限差分格式,有效提高了模拟精度和计算效率.本文将频率-空间域混合网格有限差分的思想引入到时空域,提出了时空域混合网格2 M+N型声波方程有限差分方法.首先推导出基于时空域频散关系的混合网格差分系数计算方法,然后进行频散分析、稳定性分析,并和传统高阶、时空域高阶有限差分法对比,结果表明:计算量相同时,新方法能有效压制数值频散,显著提高模拟精度;新方法相比传统2 M阶有限差分法,稳定性增强,与时空域2 M阶有限差分法稳定性基本相当.最后利用新方法进行均匀介质、层状介质、盐丘模型的数值模拟和盐丘模型的逆时偏移,模拟效果和成像质量进一步证实了该方法的有效性和普遍适用性.     

结合边界条件的截断高阶隐式有限差分方法(英文)

有限差分方法广泛应用于求解许多科技领域所涉及的偏微分方程,高阶显式有限差分方法通常用来提高求解精度,已经提出的高阶隐式有限差分方法和截断高阶显式有限差分方法可用来进一步提高模拟精度而不增加计算量。本文首先计算了针对常规网格上的一阶导数和二阶导数、交错网格上的一阶导数的有限差分系数,发现高阶隐式有限差分系数中存在一些小的系数。频散分析结果表明:忽略这些小的差分系数能够近似维持有限差分的精度,但是显著减小了计算量。然后,引入镜像对称边界条件来提高隐式有限差分方法的精度和稳定性,采用混合吸收边界条件来减小来自模型边界所不需要的反射。最后,给出了针对均匀和非均匀介质模型的弹性波模拟例子,表明了本文方法的优点。     

一阶弹性波方程交错网格高阶差分解法

提高计算精度和运算效率是所有波场正演方法所追求的目标,本文通过将速度 （应力）对时间的奇数阶高阶寻数转化为应力（速度）对空间的导数,运用时间和空间差分精度 均可达任意阶的高阶差分法,通过交错网格技术,对一阶速度－应力弹性波方程进行了数值求 解．波场快照以及实际模型的正演结果表明,这种求解一阶弹性波方程的高阶差分解法,和 常规的差分法相比网格频散显著减小,精度明显提高,而且可以取较大的空间步长,提高计算 效率。     

流固边界耦合介质高阶有限差分地震正演模拟方法

本文针对流固边界耦合介质提出了一种高效、稳定的正演数值模拟方法. 首先,从一阶位移-应力弹性波方程出发,基于海底流固边界的位移和应力的连续性条件,采用三次样条海底界面定量表征方法,推导出不规则海底界面下流固边界耦合介质中的地震波波动方程;其次,通过空间微分的高阶差分格式提高数值模拟的空间精度,并结合已推导的地震波波动方程,将四阶时间微分转换至高阶空间微分,进一步提高了数值模拟的时间精度;最后,在与标量波波动方程数值模拟结果对比分析的基础上,分别利用简单的水平层状模型和复杂海底模型,验证和讨论了本文提出的流固边界耦合介质高阶有限差分地震波正演模拟方法的有效性和准确性.      

基于PML边界的时间高阶伪谱法弹性波场模拟

在地震正演模拟中,传统的有限差分法采用二阶差分算子近似时间偏导数,如果时间采样间隔选取较大,模拟波场会出现时间频散,甚至导致差分算法的不稳定.本文研究了基于完美匹配层(PML)边界弹性波二阶方程时间四阶精度解法,通过对空间采用傅里叶变换提高计算精度并压制空间频散;而在时间域将位移对时间的四阶导数转化为空间的导数.与传统的有限差分法相比较,由于精度的提高,时间频散得到有效的压制,本文提出的方法可以适应较大的时间采样间隔,提高计算效率.     

Seismic modeling with an optimal staggered-grid finite-difference scheme based on combining Taylor-series expansion and minimax approximation

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.     

Optimal implicit staggered‐grid finite‐difference schemes based on the sampling approximation method for seismic modelling

We propose new implicit staggered‐grid finite‐difference schemes with optimal coefficients based on the sampling approximation method to improve the numerical solution accuracy for seismic modelling. We first derive the optimized implicit staggered‐grid finite‐difference coefficients of arbitrary even‐order accuracy for the first‐order spatial derivatives using the plane‐wave theory and the direct sampling approximation method. Then, the implicit staggered‐grid finite‐difference coefficients based on sampling approximation, which can widen the range of wavenumber with great accuracy, are used to solve the first‐order spatial derivatives. By comparing the numerical dispersion of the implicit staggered‐grid finite‐difference schemes based on sampling approximation, Taylor series expansion, and least squares, we find that the optimal implicit staggered‐grid finite‐difference scheme based on sampling approximation achieves greater precision than that based on Taylor series expansion over a wider range of wavenumbers, although it has similar accuracy to that based on least squares. Finally, we apply the implicit staggered‐grid finite difference based on sampling approximation to numerical modelling. The modelling results demonstrate that the new optimal method can efficiently suppress numerical dispersion and lead to greater accuracy compared with the implicit staggered‐grid finite difference based on Taylor series expansion. In addition, the results also indicate the computational cost of the implicit staggered‐grid finite difference based on sampling approximation is almost the same as the implicit staggered‐grid finite difference based on Taylor series expansion.     

选择性滤波同位网格有限差分法在地震波数值模拟中的应用.

引入计算空气声学领域的选择性滤波同位网格有限差分算法（SFFD法）用于二维地震波数值模拟．SFFD法使用经过优化的11点DRP同位网格差分格式,对空间一阶导数进行离散近似,同时采用选择性滤波方法来消除同位网格差分所产生的格点高频振荡,它既提高了数值模拟的精度, 又保证了求解过程的稳定性．数值实验结果表明,SFFD法能够达到O(Delta;x8, Delta;t4)阶交错网格算法同样的精度,同时该方法还具有很强的适应性,能够应用于存在着强泊松比差异的介质模型中,完整地模拟地震波传播过程中各类型的波场,并且对复杂非均匀介质的适应能力也很好．此外,由于避免了交错网格算法在曲线坐标系和一般各向异性介质的数值模拟时所需进行的复杂的插值运算, SFFD法在这些问题上也有着很好的应用前景．      

一种全局优化的隐式交错网格有限差分算法及其在弹性波数值模拟中的应用

在数值模拟中,隐式有限差分具有较高的精度和稳定性.然而,传统隐式有限差分算法大多由于需要求解大型矩阵方程而存在计算效率偏低的局限性.本文针对一阶速度-应力弹性波方程,构建了一种优化隐式交错网格有限差分格式,然后将改进格式由时间-空间域转换为时间-波数域,利用二范数原理建立目标函数,再利用模拟退火法求取优化系数.通过对均匀模型以及复杂介质模型进行一阶速度-应力弹性波方程数值模拟所得单炮记录、波场快照分析表明:这种优化隐式交错网格差分算法与传统的几种显式和隐式交错网格有限差分算法相比不但降低了计算量,而且能有效的压制网格频散,使弹性波数值模拟的精度得到有效的提高.     

地震波传播的三维伪谱和高阶有限差分混合方法并行模拟

基于交错网格伪谱法和高阶精度有限差分方法,发展了模拟非均匀介质地震波传播的三维伪谱和有限差分混合算法.该方法在两个水平方向利用交错网格伪谱算子计算空间微分,保留了该方法高效、高精度的优势,在垂直方向采用交错网格高阶精度有限差分算子实现空间微分计算．利用有限差分方法的局部性特征,将三维计算区域在垂直方向上划分为一系列子区域,并分配给不同的处理器,实现了在并行计算机集群上的三维并行计算．通过模拟算例,与离散波数法比较,检验了该算法的精度．为了检验该方法的实用性,在64个处理器上,对三维沉积盆地模型进行了67108864个网格点的并行计算,模拟的波场主频率为1.25Hz,讨论了沉积盆地深度对三维沉积盆地地面运动的影响．      

求解弹性波方程的辛RKN格式

将弹性波方程变换至Hamilton体系,构造适用于弹性波模拟的高效显式二阶辛Runge-Kutta-Nystrm(RKN)格式,运用根数理论得到此格式的阶条件方程组.通过给定系数的限定条件,得到方程的对称解.为了使时间离散误差达到极小,提出数值频率与真实频率比较,通过Taylor展开,得到关于辛系数的限定方程,求解方程组得到最小频散辛RKN格式.对比分析时间演进方程的稳定性,得到使库朗数达到极大值的限定方程,求解方程组得到最稳定辛RKN格式.发现此两种格式为同一格式.新得到的辛RKN格式不依赖于空间离散方法,为了对比的需要,选取有限差分法进行空间离散.在频散、稳定性分析中,与常见辛格式对比,从理论上分析了本文提出的格式在数值频散压制、稳定性提升等方面的优势,数值实验进一步证实了理论分析的正确性.     

ACOUSTIC MODELLING ON A GRID OF VERTICALLY VARYING SPACING1

The Fourier spectral method and high-order differencing have both been shown to be very accurate in computing spatial derivatives of the acoustic wave equation, requiring only two and three gridpoints per shortest wavelength respectively. In some cases, however, there is a lack of flexibility as both methods use a uniform grid. If these methods are applied to structures with high vertical velocity contrasts, very often most of the model is oversampled. If a complicated interface has to be covered by a fine grid for exact representation, both methods become less attractive as the homogeneous regions are sampled more finely than necessary. In order avoid this limitation we present a differencing scheme in which the grid spacings can be extended or reduced by any integer factor at a given depth. This scheme adds more flexibility and efficiency to the acoustic modelling as the grid spacings can be changed according to the material properties and the model geometry. The time integration is carried out by the rapid expansion method. The spatial derivatives are computed using either the Fourier method or a high-order finite-difference operator in the x-direction and a modified high-order finite-difference operator in the z-direction. This combination leads to a very accurate and efficient modelling scheme. The only additional computation required is the interpolation of the pressure in a strip of the computational mesh where the grid spacing changes.     

A post-processing technique and an a posteriori error estimate for the newmark method in dynamic analysis

In this paper, we present a post-processing technique and an a posteriori error estimate for the Newmark method in structural dynamic analysis. By post-processing the Newmark solutions, we derive a simple formulation for linearly varied third-order derivatives. By comparing the Newmark solutions with the exact solutions expanded in the Taylor series, we achieve the local post-processed solutions which are of fifth-order accuracy for displacements and fourth-order accuracy for velocities in one step. Based on the post-processing technique, a posteriori local error estimates for displacements, velocities and, thus, also the total energy norm error estimate are obtained. If the Newmark solutions are corrected at each step, the post-processed solutions are of third-order accuracy in the global sense, i.e. one-order improvement for the original Newmark solutions is achieved. We also discuss a method for estimating the global time integration error. We find that, when the total energy norm is used, the sum of the local error estimates will give a reasonable estimate for the global error. We present numerical studies on a SDOF and a 2-DOF example in order to demonstrate the performance of the proposed technique.     

Dispersion-relationship-preserving seismic modelling using the cross-rhombus stencil with the finite-difference coefficients solved by an over-determined linear system

Finite-difference modeling with a cross-rhombus stencil with high-order accuracy in both spatial and temporal derivatives is a potential method for efficient seismic simulation. The finite-difference coefficients determined by Taylor-series expansion usually preserve the dispersion property in a limited wavenumber range and fixed angles of propagation. To construct the dispersion-relationship-preserving scheme for satisfying high-wavenumber components and multiple angles, we expand the dispersion relation of the cross-rhombus stencil to an over-determined system and apply a regularization method to obtain the stable least-squares solution of the finite-difference coefficients. The new dispersion-relationship-preserving based scheme not only satisfies several designated wavenumbers but also has high-order accuracy in temporal discretization. The numerical analysis demonstrates that the new scheme possesses a better dispersion characteristic and more relaxed stability conditions compared with the Taylor-series expansion based methods. Seismic wave simulations for the homogeneous model and the Sigsbee model demonstrate that the new scheme yields small dispersion error and improves the accuracy of the forward modelling.