一种优化的频率域三维声波有限差分模拟方法

为提高频率域有限差分（FD,finite-difference）正演模拟技术的计算精度和效率,基于旋转坐标系统的优化差分格式被广泛应用,但是只应用于正方形网格的情况.基于平均导数法（ADM）的优化差分格式,应用于正方形和长方形网格模拟.这些频率域有限差分算子,各自具有不同的差分格式和对应的优化系数求解表达式.本文基于三维声波方程发展了一种新的优化方法,只要给定FD模板形式,可直接构造频散方程,求取FD模板上各节点的优化系数.此方法的优点在于频率域FD算子的优化系数对应各个节点,可扩展优化其他格式.运用此优化方法,计算得到了不同空间采样间距比情况下27点和7点格式的优化系数.数值实验表明,优化27点格式与ADM 27点格式具有相同的精度,优化7点格式比经典的7点格式具有更小的数值频散.     

一种模拟非均匀介质中弹性波传播新的k-space方法

传统的伪谱(PS)方法,采用傅里叶变换(FT)计算空间导数具有很高的精度,每个波长仅需要两个采样点,而时间导数采用有限差分(FD)近似因而精度较低.当采用大时间步长时,由于时空精度不平衡,PS法存在不稳定性问题.原始的k-space方法可以有效地克服这些问题但是却无法适用于非均匀介质.为了提高原始k-space方法模拟非均匀介质波动方程的精度,我们提出了一种新的k-space算子族.它是用非均匀介质的变速度代替原k-space算子中的常数补偿速度构造得到,引入低秩近似可以高效求解.我们将构造的新的k-space算子应用于耦合的二阶位移波动方程,而不是交错网格一阶速度应力波动方程,使模拟弹性波的计算存储量减少.我们从数学上证明了基于二阶波动方程的k-space方法与基于一阶波动方程的k-space方法是等价的.数值模拟实验表明,与传统的PS、交错网格PS和原始的k-space方法相比,我们的新方法可以在时间和空间步长较大的均匀和非均匀介质中,为弹性波的传播提供更精确的数值解.在保持稳定性和精度的同时,采用较大的时空采样间隔,可以大大降低数值模拟的计算成本.     

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

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

基于时空域自适应高阶有限差分的声波叠前逆时偏移

叠前逆时偏移在理论上是现行偏移方法中最为精确的一种成像方法,其实现过程中的核心步骤之一是波动方程的波场延拓,而波场延拓的本质是求解波动方程,所以精确、快速地求解波动方程对逆时偏移至关重要.本文采用一种基于时空域频散关系的有限差分方法来求解声波方程,分析其频散和稳定性,实现波场数值模拟,并将分析和模拟结果与传统有限差分法进行对比.分析结果和模型数值模拟结果都表明时空域有限差分法模拟精度更高、稳定性更好.将时空域高阶有限差分法应用到叠前逆时偏移波场延拓的方程求解中,然后再利用归一化互相关成像条件成像,理论模型数据偏移处理获得了精度更高的成像.同时,在逆时偏移波场延拓的实现中,采用自适应变长度的空间差分算子求解空间导数的有限差分策略,在不影响数值模拟和成像精度的前提下,有效地提高了计算效率.     

最优化辛格式广义离散奇异核褶积微分算子地震波场模拟

将波动方程变换至Hamilton体系,构造了一种新的保结构算法,即最优化辛格式广义褶积微分算子(OSGCD). 在时间离散上,首先引入了Lie算子设计二级二阶辛格式,基于最小误差原理得到了优化的辛格式. 在空间离散上,引入广义离散奇异核褶积微分算子计算空间微分,提出了一种有效方法优化GCD并得到了稳定的算子系数. 针对本文发展的新方法,给出了OSGCD稳定性条件. 在数值实验中,将OSGCD与多种方法比较,从精度和计算效率两方面分析了OSGCD的计算优势,计算结果也表明OSGCD长时程以及非均匀介质中地震波模拟亦具有较强能力.     

各向异性介质弹性波传播的三维不规则网格有限差分方法

提出一种新的三维空间不规则网格有限差分方法,模拟具有地形构造的非均匀各向异性介质中弹性波传播过程. 该方法通过具有二阶时间精度和四阶空间精度的不规则交错网格差分算子来近似一阶弹性波动方程,与多重网格不同,无需在精细网格和粗糙网格间进行插值,所有网格点上的计算在同一次空间迭代中完成. 针对具有复杂物性参数和复杂几何特征的地层结构,使用精细不规则网格处理粗糙界面、断层和空间界面等复杂几何构造, 理论分析和数值算例表明,该方法不但节省了大量计算机内存和计算时间,而且具有令人满意的稳定性和精度.     

基于广义正交多项式褶积微分算子的地震波场数值模拟方法

地震波场模拟方法研究对于与波动现象有关的地震学问题的重要性是不言而喻的.就目前现有的各种正演算法来说,精度较高的算法（如有限元法、谱元法、高阶有限差分法等）,其计算速度较慢;计算速度较快的算法（如低阶有限差分法、付氏伪谱法等）计算精度却比较低.为了兼顾地震波场模拟的精度与速度,本文推出了一种快速的、高精度地震波场模拟方法（基于Forsyte广义正交多项式的褶积微分算子法）,该方法是以计算数学中的Forsyte广义正交多项式插值函数为基础,构建一个新的褶积微分算子,并将该算子引入到地震波动方程的一阶速度-应力方程的空间微分运算中去,采用时间交错网格有限差分算子替代普通的差分算子以匹配高精度的褶积微分算子,从而构造一种全新的地震波场数值模拟方法.该方法同时具有广义正交多项式方法的高精度和短算子低阶有限差分算法的高速度.通过对算子长度的调节及算子系数的优化,可同时兼顾波场解的全局信息与局部信息.复杂非均匀介质模型中的波场数值模拟实验证实了该方法的可行性及优越性.     

基于精细积分法的三维弹性波数值模拟(英文)

波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。     

VTI介质纯P波混合法正演模拟及稳定性分析

各向异性介质纯P波方程完全不受横波的干扰,在一定程度上可以减缓由于介质各向异性引起的数值不稳定,本文推导了具有垂直对称轴的横向各向同性(VTI)介质纯P波一阶速度-应力方程.由于纯P波方程存在一个分数形式的伪微分算子,无法直接采用有限差分法求解.针对该问题,本文采用伪谱法和高阶有限差分法联合求解波动方程,重点分析了混合法求解纯P波一阶速度-应力方程的稳定性问题,并给出了混合法求解纯P波方程的稳定性条件.数值模拟结果表明纯P波方程伪谱法和高阶有限差分混合法能够进行复杂介质的正演模拟,在强变速度、变密度的地球介质中仍然具有较好的稳定性.     

Hybrid seismic modeling based on discrete-wave number and finite-difference methods

Any calculation of seismic wave propagation comprising the seismic source, the travel path, and the receiver site in a single finite-difference (FD) model requires a considerable amount of computer time and memory. Moreover, the methods currently available for including point sources in the 2D FD calculations are far-field approximations only. Therefore we have developed a new hybrid method for treating the seismic wave fields at localized 2D near-surface structures embedded in a 1D background medium, and excited by a point source. The source radiation and propagation in the background model is solved by the discrete-wave number (DW) method, while the propagation in the local 2D structure is calculated by the FD method. The coupling between the two sets of calculations is performed on a rectangular excitation box surrounding the local structure. We show the usefulness of the method in ground-motion studies where both near-field source effects and local site effects are important. Technical problems connected with the inconsistency between the 3D source radiation and the 2D FD calculation are minor for the relatively distant in-plane point explosive sources, but are more serious for the in-plane dislocation sources.     

基于自适应优化有限差分方法的全波VSP逆时偏移

与地面地震资料相比,VSP资料具有分辨率高、环境噪声小及能更好地反映井旁信息等优点.常规VSP偏移主要对上行反射波进行成像,存在照明度低、成像范围受限等问题.为了增加照明度、拓宽成像范围、提高成像精度,本文采用直达波除外的所有声波波场数据(全波),包括一次反射波、多次反射波等进行叠前逆时偏移成像.针对逆时偏移中的四个关键问题,即波场延拓、吸收边界条件、成像条件及低频噪声的压制,本文分别采用自适应变空间差分算子长度的优化有限差分方法(自适应优化有限差分方法)求解二维声波波动方程以实现高精度、高效率的波场延拓,采用混合吸收边界条件压制因计算区域有限所引起的人工边界反射,采用震源归一化零延迟互相关成像条件进行成像,采用拉普拉斯滤波方法压制逆时偏移中产生的低频噪声.本文对VSP模型数据的逆时偏移成像进行了分析,结果表明:自适应优化有限差分方法比传统有限差分方法具有更高的模拟精度与计算效率,适用于VSP逆时偏移成像;全波场VSP逆时偏移成像比上行波VSP逆时偏移的成像范围大、成像效果好;相对于反褶积成像条件,震源归一化零延迟互相关成像条件具有稳定性好、计算效率高等优点.将本文方法应用于某实际VSP资料的逆时偏移成像,进一步验证了本文方法的正确性和有效性.     

求解二阶解耦弹性波方程的低秩分解法和低秩有限差分法

时间域的波场延拓方法在本质上都可以归结为对一个空间-波数域算子的近似.本文基于一阶波数-空间混合域象征,提出一种新的方法求解解耦的二阶位移弹性波方程.该方法采用交错网格,连续使用两次一阶前向和后向拟微分算子,推导得到了解耦的二阶位移弹性波方程的波场延拓算子.由于该混合域象征在伪谱算子的基础上增加了一个依赖于速度模型的补偿项,可以补偿由于采用二阶中心差分计算时间微分项带来的误差,有效地减少模拟结果的数值频散,提高模拟精度.然而,在非均匀介质中,直接计算该二阶的波场延拓算子,每一个时间步上需要做N次快速傅里叶逆变换,其中N是总的网格点数.为了减少计算量,提出了交错网格低秩分解方法;针对常规有限差分数值频散问题,本文将交错网格低秩方法与有限差分法结合,提出了交错网格低秩有限差分法.数值结果表明,交错网格低秩方法和交错网格低秩有限差分法具有较高的精度,对于复杂介质的地震波数值模拟和偏移成像具有重要的价值.     

基于优化有限差分和混合吸收边界条件的三维VTI介质声波和弹性波数值模拟

传统有限差分系数是通过泰勒级数展开求取的,这样导致所计算的频散曲线在大波数区域会产生较强的数值误差.针对二阶空间偏导数的显式有限差分离散,本文发展了一种新的优化差分系数方法：首先将泰勒级数展开与多点采样方法结合应用于空间频散关系,基于最大范数建立直观有效的优化目标函数,采用Remez算法求解该目标函数,从而获得最优化差分系数.利用优化有限差分方法求解三维垂直对称轴横向各向同性（VTI）介质中的声波和弹性波方程.另外,本文将二维混合吸收边界条件推广到三维VTI介质中,用于吸收人工截断边界反射;基于各向异性特征,合理调整了边界区域的速度值来提高吸收效果.考虑到三维情况下计算效率的问题,本文波场外推过程中采用图形处理器（GPU）取代传统的中央处理器（CPU）.数值精度分析表明,相比较于传统的泰勒级数展开方法,优化有限差分方法在大波数区域对频散误差的压制效果更明显.在三维均匀和修改的Hess VTI模型中的数值模拟实验证明了本文方法具有更高的精度与效率,混合吸收边界条件在三维VTI介质中具有良好的边界吸收效果.     

TOWARDS OPTIMUM ONE-WAY WAVE PROPAGATION1

Numerical wavefield extrapolation represents the backbone of any algorithm for depth migration pre- or post-stack. For such depth imaging techniques to yield reliable and interpretable results, the underlying wavefield extrapolation algorithm must propagate the waves through inhomogeneous media with a minimum of numerically induced distortion, over a range of frequencies and angles of propagation. A review of finite-difference (FD) approximations to the acoustic one-way wave equation in the space-frequency domain is presented. A straightforward generalization of the conventional FD formulation leads to an algorithm where the wavefield is continued downwards with space-variant symmetric convolutional operators. The operators can be precomputed and made accessible in tables such that the ratio between the temporal frequency and the local velocity is used to determine the correct operator at each grid point during the downward continuation. Convolutional operators are designed to fit the desired dispersion relation over a range of frequencies and angles of propagation such that the resulting numerical distortion is minimized. The optimization is constrained to ensure that evanescent energy and waves propagating at angles higher than the maximum design angle are attenuated in each extrapolation step. The resulting operators may be viewed as optimally truncated and bandlimited spatial versions of the familiar phase shift operator. They are unconditionally stable and can be applied explicitly. This results in a simple wave propagation algorithm, eminently suited for implementation on pipelined computers and on large parallel computing systems.     

一种时空域优化的高精度交错网格差分算子与正演模拟

如何有效压制数值频散是有限差分正演模拟研究中的关键问题之一.近年来,许多学者对二阶声波方程的差分算子开展了大量的优化工作,在压制频散方面取得不错的效果.一阶压强-速度方程广泛用于研究地震波在地下变密度模型中传播规律,目前针对一阶方程的优化工作大多只是在空间差分算子上展开.本文在前人研究的基础上,推导出一阶声波方程中压强场与偏振速度场之间的解析关系,据此在传统交错网格基础上给出一种高精度的显式时间递推格式,该递推格式将时间差分与空间差分算子结合在一起,并采用共轭梯度法得到精确时间递推匹配系数,实现时空差分算子的同时优化.在编程实现算法的基础上,通过频散分析与三个典型模型测试表明:本文方法能够较为有效地压制时间频散与空间频散,提高数值计算精度;同时对复杂模型也有很好适用性.     

重叠区域伪谱法计算非均匀地球介质地震波传播

伪谱法是一种高效、高精度计算非均匀介质地震波传播的数值算法,由于它的微分算子的全局性,一般认为该方法不适于并行计算. 本文介绍了并行计算非均匀介质中地震波传播的重叠区域分解算法,给出了一种基于傅里叶伪谱法的并行算法. 文中给出的算法将介质划分为相互重叠的若干区域,在各个子域上单独求解,利用重叠部分的解的传递,将各个子域连接起来,实现了伪谱法在分布式并行处理机上的计算. 文中给出了一个将二维区域分解的算例,比较了并行算法和整体算法的结果,分析了并行算法的计算精度. 结果表明,并行算法会有效降低计算时间,并且保证计算精度. 该方法在大规模三维非均匀介质的地震波场模拟方面有应用价值.      

地震波传播的哈密顿表述及辛几何算法

地震波传播过程本质上是能量在传播过程中逐步损耗直至殆尽的过程，而在实际应用中，常在无能量损耗假设下，用弹性波动方程或标量波动方程描述它.在哈密顿(Hamilton)体系表述下，地震波传播过程即为一个无限维的哈密顿系统随时间的演化过程.若不计能量损耗，波场演化过程实质上为一个单参数连续的辛变换，因而对应的数值算法应为辛几何算法.本文首先从地震波标量方程出发，给出哈密顿体系下地震波传播的表述，即任意两个时刻的波场是通过辛变换联系起来的.随后，把波场在时间和相空间离散化后，给出了用于波场计算的一些辛格式，如显式辛格式、隐式辛格式和蛙跳辛格式.并进一步讨论了有限差分格式和辛格式的异同.然后，应用显式辛格式和同阶的有限差分方法给出了同一理论速度模型下的波场和Marmousi速度模型下的单炮记录.数值结果表明，辛算法是一类可行的波场模拟的数值算法.在时间步长较小时，有限差分方法是辛算法的一个很好近似.文中的理论和方法，为地震波传播理论及实际应用研究提供了新的途径.     

Visco‐acoustic prestack reverse‐time migration based on the time‐space domain adaptive high‐order finite‐difference method

It is important to include the viscous effect in seismic numerical modelling and seismic migration due to the ubiquitous viscosity in an actual subsurface medium. Prestack reverse‐time migration (RTM) is currently one of the most accurate methods for seismic imaging. One of the key steps of RTM is wavefield forward and backward extrapolation and how to solve the wave equation fast and accurately is the essence of this process. In this paper, we apply the time‐space domain dispersion‐relation‐based finite‐difference (FD) method for visco‐acoustic wave numerical modelling. Dispersion analysis and numerical modelling results demonstrate that the time‐space domain FD method has great accuracy and can effectively suppress numerical dispersion. Also, we use the time‐space domain FD method to solve the visco‐acoustic wave equation in wavefield extrapolation of RTM and apply the source‐normalized cross‐correlation imaging condition in migration. Improved imaging has been obtained in both synthetic and real data tests. The migration result of the visco‐acoustic wave RTM is clearer and more accurate than that of acoustic wave RTM. In addition, in the process of wavefield forward and backward extrapolation, we adopt adaptive variable‐length spatial operators to compute spatial derivatives to significantly decrease computing costs without reducing the accuracy of the numerical solution.     

瞬变电磁Crank-Nicolson FDTD三维正演

时域有限差分(FDTD)方法使用Yee网格剖分电磁场的空间采样,通过时间步迭代实现电磁场数值模拟,具有内存消耗低、计算简单等特点,常用于瞬变电磁三维正演.然而,常规FDTD方法的时间迭代步长Δt受Courant-Friedrich-Lewy(CFL)条件严格限制,过多的迭代次数以及过密的采样往往导致计算速度慢、累积误差不断增大.本文提出一种不受CFL条件约束的无条件稳定隐式差分算法Crank-Nicolson FDTD(CN-FDTD)用于瞬变电磁三维正演.基于Crank-Nicolson差分方法对Maxwell方程组重新离散,空间网格仍然采用Yee元胞,时间步进采用在整时间步电场、磁场同时采样的策略,建立无条件稳定FDTD格式,突破CFL条件限制.与常规FDTD交替采样相比,CN-FDTD电场、磁场同时采样的策略构成的隐式差分格式,需要求解大型稀疏矩阵方程组.通常,瞬变电磁三维正演模型中产生的矩阵阶数往往较大,需要占用大量内存和求解时间.为解决上述问题,采用Crank-Nicolson-cycle-sweep-uniform(CNCSU-FDTD)方法近似求解CN-FDTD方程,在保...     

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

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