基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移

时间域常Q黏声波方程,由于含分数阶时间导数项,数值求解需要大量内存,计算效率低,不利于地震偏移的实施.通过一系列近似,可将该方程简化为介质频散效应和衰减效应解耦的分数阶拉普拉斯算子黏声波方程,数值求解内存需求少,计算效率高.本文采用交错网格有限差分逼近时间导数,改进的伪谱法计算空间导数,PML吸收边界去除边界反射,对该方程进行数值离散和地震正演模拟,开展地震数据的黏声介质逆时偏移,实现波场逆时延拓过程中同时完成频散校正和衰减补偿.改善深层构造的成像精度,数值结果表明,基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移可大幅度提高地震模拟计算效率,偏移剖面明显优于常规声波偏移剖面,极大改善深层构造的成像品质.     

最小二乘傅立叶有限差分偏移

一般偏移算法是用反演算子通过解析方法求解.最小二乘偏移方法采用另一种思路,即采用数值方法,通过解一个线性离散反问题来索求解.这样我们试着寻找一个模型匹配地震数据并能表现出其某些特点来得到偏移图像.最小二乘法能减少偏移赝像,得到更精确的偏移效果.Kirchhoff算子在最小二乘偏移方法中应用较广,但需要较多的迭代次数,而且具有Kirchhoff偏移的缺点.本文把最小二乘方法运用到基于波长延拓的波动方程偏移方法中,为提高最小二乘偏移的效率,可采用效率较高的正传播算子和反传播算子.我们利用效率较高,能适应剧烈横向变速的傅立叶有限差分正传播和反传播算子来做叠后最小二乘偏移.数值实例表明,通过少数的共轭梯度法迭代,就能得到与真实模型差别很微小的偏移效果.对于傅式变换我们采用了数值软件FFTW,其变换速度比常规FFT算法一般要快六倍以上,进一步提高了效率.本文算法很容易在并行机上实现,这些特点在处理大型数据时大有裨益.     

三维复杂构造中地震波模拟的单程波方法

复杂构造中单程波与双程波方法模拟结果的比较表明,就地震勘探中主要关心的一次反射波而言,单程波算法已具有足够的精度. 使用单程波方程将极大地减少数值计算的计算量,同时对介质的几何和物理参数建模也降低了要求. 单程波算法可视为深度偏移的“逆运算”,这样可以很好地借用已知的深度偏移方法及其程序系统. 基于计算效率和计算精度的双重考虑,本文在介质速度结构较复杂时采用显式短算子波场延拓方法,而在介质速度结构相对简单时采用分裂步相移法. 反射系数的计算中考虑了其随入射角的变化.     

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

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

The One-Way Wave Equation: A Full-Waveform Tool for Modeling Seismic Body Wave Phenomena

The study of seismic body waves is an integral aspect in global, exploration and engineering scale seismology, where the forward modeling of waves is an essential component in seismic interpretation. Forward modeling represents the kernel of both migration and inversion algorithms as the Green’s function for wavefield propagation and is also an important diagnostic tool that provides insight into the physics of wave propagation and a means of testing hypotheses inferred from observational data. This paper introduces the one-way wave equation method for modeling seismic wave phenomena and specifically focuses on the so-called operator-root one-way wave equations. To provide some motivation for this approach, this review first summarizes the various approaches in deriving one-way approximations and subsequently discusses several alternative matrix narrow-angle and wide-angle formulations. To demonstrate the key strengths of the one-way approach, results from waveform simulation for global scale shear-wave splitting modeling, reservoir-scale frequency-dependent shear-wave splitting modeling and acoustic waveform modeling in random heterogeneous media are shown. These results highlight the main feature of the one-way wave equation approach in terms of its ability to model gradual vector (for the elastic case) and scalar (for the acoustic case) waveform evolution along the underlying wavefront. Although not strictly an exact solution, the one-way wave equation shows significant advantages (e.g., computational efficiency) for a range of transmitted wave three-dimensional global, exploration and engineering scale applications.     

裂步法最小二乘偏移

最小二乘偏移的核心思想是用正向传播算子的逆算子代替其共轭转置算子进行偏移.它克服了常规偏移方法用共轭转置算子偏移的缺点.本文详细地描述了裂步法最小二乘偏移的原理,实现算法,并通过两个模型展示了裂步法最小二乘偏移的可行性和有效性.并且通过研究发现最小二乘较容易被扩展到基于其它波场延拓算子的最小二乘偏移.     

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

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

The Numerical Modeling of 3-D Elastic Wave Equation Using a High-Order,Staggered-Grid, Finite Difference Scheme

This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, temporal-and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial resuits for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.     

图法及其在Toeplitz矩阵分解中的应用

实现了基于图法的稀疏正定系统的求解，并在此基础上实现了具有Toeplitz结构的大型稀疏矩阵的快速LU分解，在基于波边方程的地震数据处理如地震波场模拟和叠前深度偏移等隐式方法中，拉普拉氏算子或亥姆霍兹算子的快速分解是这些方法能否实现的关键，在螺旋边界条件下，这些算子的表示矩阵是具有Toeplitz结构的正定厄密矩阵，可以通过本文方法实现快速分解。     

基于照明补偿的单程波最小二乘偏移

最小二乘偏移是一种基于反射地震数据与地下反射率间线性关系而建立起来的地震数据线性反演方法,相比常规偏移成像具有更好的保幅性能.本文提出了一种基于照明补偿的单程波最小二乘偏移方法,首先利用单程波方程的稳定Born近似广义屏波场传播算子构建反射地震数据与地下反射率间的线性算子,然后再应用线性最优化方法求解最小二乘偏移所对应的线性反问题.在迭代求解最优化问题的过程中,以地震波场的地下照明强度作为迭代反演的预条件算子加快迭代的收敛速度.单程波传播过程中考虑了速度分界面产生的透射效应,并用单极震源代替常规偏移中的偶极震源.把本文提出的方法应用于层状理论模型和Marmosi模型地震数据的数值试验中均取得了理想的结果.     

波场模拟中的数值频散分析与校正策略

波动方程有限差分法正演模拟,对认识地震波传播规律、进行地震属性研究、地震资料地质解释、储层评价等,均具有重要的理论和实际意义.但有限差分法本身固有存在着数值频散问题,数值频散在正演模拟中是一种严重的干扰,会降低波场模拟的精度与分辨率.针对TI介质波场模拟的交错网格有限差分方法,本文从空间网格离散、时间网格离散和算子近似等三个方面对其产生的数值频散进行了分析,并结合其他学者的研究成果给出了TI介质波场模拟中压制数值频散的方法与策略：在已知介质频散关系时,对差分算子可实施算子校正;通过提高差分方程的阶数来提高波场模拟精度;采用流体力学中守恒式方程的通量校正传输方法来压制波场模拟中的数值频散;在实际正演模拟时,采用交错网格高阶有限差分方程,不仅在空间上采用高阶差分,而且在时间上也要采用高阶差分,否则只在单一方向上（空间或时间）提高方程的阶数对压制数值频散也不会取得理想的效果.     

Acoustic Reverse-time Migration using Optimal Staggered-grid Finite-difference Operator Based on Least Squares

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.     

非零偏VSP弹性波叠前逆时深度偏移技术探讨

非零偏VSP地震资料是一种多分量资料,处理非零偏VSP资料,弹性波叠前逆时深度偏移技术无疑是最适合的处理技术.本文从二维各向同性介质的弹性波波动方程出发,研究了对非零偏VSP资料进行叠前逆时深度偏移的偏移算法,讨论了逆时传播过程中的边值问题和数值频散问题及其相应的解决方案;采用求解程函方程计算得到地下各点的地震波初至时间作为成像时间,实现了非零偏VSP资料的叠前逆时深度偏移.最后进行了模型试算和非零偏VSP地震资料的试处理,结果表明该方法不受地层倾角限制,较适用于高陡构造地区或介质横向速度变化较大地区的非零偏VSP地震资料处理.     

Poynting and polarization vectors based wavefield decomposition and their application on elastic reverse time migration in 2D transversely isotropic media

With the progress in computational power and seismic acquisition, elastic reverse time migration is becoming increasingly feasible and helpful in characterizing the physical properties of subsurface structures. To achieve high-resolution seismic imaging using elastic reverse time migration, it is necessary to separate the compressional (P-wave) and shear (S-wave) waves for both isotropic and anisotropic media. In elastic isotropic media, the conventional method for wave-mode separation is to use the divergence and curl operators. However, in anisotropic media, the polarization direction of P waves is not exactly parallel to the direction of wave propagation. Also, the polarization direction of S-waves is not totally perpendicular to the direction of wave propagation. For this reason, the conventional divergence and curl operators show poor performance in anisotropic media. Moreover, conventional methods only perform well in the space domain of regular grids, and they are not suitable for elastic numerical simulation algorithms based on non-regular grids. Besides, these methods distort the original wavefield by taking spatial derivatives. In this case, a new anisotropic wave-mode separation scheme is developed using Poynting vectors. This scheme can be performed in the angle domain by constructing the relationship between group and polarization angles of different wave modes. Also, it is performed pointwise, independent of adjacent space points, suitable for parallel computing. Moreover, there is no need to correct the changes in phase and amplitude caused by the derivative operators. By using this scheme, the anisotropic elastic reverse time migration is more efficiently performed on the unstructured mesh. The effectiveness of our scheme is verified by several numerical examples.     

Acoustic finite-difference modeling beyond conventi-onal Courant-Friedrichs-Lewy stability limit: App-roach based on variable-length temporal and spatial operators

Conventional finite-difference (FD) methods cannot model acoustic wave propagation beyond Courant-Friedrichs-Lewy (CFL) numbers 0.707 and 0.577 for two-dimensional (2D) and three-dimensional (3D) equal spacing cases, respectively, thereby limiting time step selection. Based on the definition of temporal and spatial FD operators, we propose a variable-length temporal and spatial operator strategy to model wave propagation beyond those CFL numbers while preserving accuracy. First, to simulate wave propagation beyond the conventional CFL stability limit, the lengths of the temporal operators are modified to exceed the lengths of the spatial operators for high-velocity zones. Second, to preserve the modeling accuracy, the velocity-dependent lengths of the temporal and spatial operators are adaptively varied. The maximum CFL numbers for the proposed method can reach 1.25 and 1.0 in high velocity contrast 2D and 3D simulation examples, respec-tively. We demonstrate the effectiveness of our method by modeling wave propagation in simple and complex media.     

Principle and application of high density spatial sampling in seismic migration

To avoid spatial aliasing problems in broad band high resolution seismic sections, I present a high density migration processing solution. I first analyze the spatial aliasing definition for stack and migration seismic sections and point out the differences between the two. We recognize that migration sections more often show spatial aliasing than stacked sections. Second, from wave propagation theory, I know that migration output is a new spatial sampling process and seismic prestack time migration can provide the high density sampling to prevent spatial aliasing on high resolution migration sections. Using a 2D seismic forward modeling analysis, I have found that seismic spatial aliasing noise can be eliminated by high density spatial sampling in prestack migration. In a 3D seismic data study for Daqing Oilfield in the Songliao Basin, I have also found that seismic sections obtained by high-density spatial sampling (10 × 10 m) in prestack migration have less spatial aliasing noise than those obtained by conventional low density spatial sampling (20 × 40 m) in prestack migration.     

用于波场成像的谱法LU分解

地震波场模拟和偏移成像等有限差分隐格式算法中的重要环节，是实现亥姆霍兹算子表示矩阵H的快速求逆运算. 在螺旋边界条件下，H具有Toeplitz结构的正定厄密矩阵，其快速求逆可由谱法LU分解实现. 本文分析了谱法LU分解对提高计算速度的原理及特点，并着重讨论了在不同类型的介质模型中，采用谱法分解矩阵H时带来的数值误差、误差的分布及其对波场计算的影响. 研究结果表明，对均匀介质而言，矩阵H各列具有相同的非零元素分布，谱法LU分解的误差在吸收边界条件下，不影响波场模拟和成像计算；但对于非均匀介质模型，矩阵H各列具有不同的非零元素分布，谱法LU分解的误差随介质不均匀性程度的增大而增大，势必影响非均匀介质中波场计算. 在波场模拟和成像等有限差分隐格式算法中，采用谱法LU分解完成矩阵求逆时，必须考虑到并尽量减少该方法的误差对波场计算的影响.     

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

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

Seismic wavefield modeling based on time-domain symplectic and Fourier finite-difference method

Seismic wavefield modeling is important for improving seismic data processing and interpretation. Calculations of wavefield propagation are sometimes not stable when forward modeling of seismic wave uses large time steps for long times. Based on the Hamiltonian expression of the acoustic wave equation, we propose a structure-preserving method for seismic wavefield modeling by applying the symplectic finite-difference method on time grids and the Fourier finite-difference method on space grids to solve the acoustic wave equation. The proposed method is called the symplectic Fourier finite-difference (symplectic FFD) method, and offers high computational accuracy and improves the computational stability. Using acoustic approximation, we extend the method to anisotropic media. We discuss the calculations in the symplectic FFD method for seismic wavefield modeling of isotropic and anisotropic media, and use the BP salt model and BP TTI model to test the proposed method. The numerical examples suggest that the proposed method can be used in seismic modeling of strongly variable velocities, offering high computational accuracy and low numerical dispersion. The symplectic FFD method overcomes the residual qSV wave of seismic modeling in anisotropic media and maintains the stability of the wavefield propagation for large time steps.     

A probabilistic model for ghost delay time estimation based on recording geometry

In order to deconvolve the ghost response from marine seismic data, an estimate of the ghost operator is required. Typically, this estimate is made using a model of in‐plane propagation, i.e., the ray path at the receiver falls in the vertical plane defined by the source and receiver locations. Unfortunately, this model breaks down when the source is in a crossline position relative to the receiver spread. In this situation, in‐plane signals can only exist in a small region of the signal cone. In this paper, we use Bayes' theory to model the posterior probability distribution functions for the vertical component of the ray vector given the known source–receiver azimuth and the measured inline component of the ray vector. This provides a model for the ghost delay time based on the acquisition geometry and the dip of the wave in the plane of the streamer. The model is fairly robust with regard to the prior assumptions and controlled by a single parameter that is related to the likelihood of in‐plane propagation. The expected values of the resulting distributions are consistent with the deterministic in‐plane model when in‐plane likelihood is high but valid everywhere in the signal cone. Relaxing the in‐plane likelihood to a reasonable degree radically simplifies the shape of the expected‐value surface, lending itself for use in deghosting algorithms. The model can also be extended to other plane‐wave processing problems such as interpolation.