气枪激发信号传播的谱元法数值模拟研究

为了研究气枪激发信号的波场,本文利用谱元法对双相介质中波的传播做了数值模拟,分析了波的传播特征.本文主要做了以下工作:(1)研究了使用谱无法(SEM)模拟孔隙弹性介质中波的传播,模拟结果表明,采用谱元法能有效解决双相介质的波场传播模拟问题.(2)验证了Biot理论中慢纵波的存在.双相介质中存在明显的慢纵波,流相波场的慢...     

分段光滑曲线边界波动方程数值模拟研究

矩形网格有限差分法在地震波传播数值模拟方面具有计算速度快的显著优势,但该方法在处理复杂边界问题上存在着效率低的严重缺陷.本文针对分段光滑曲线边界定义了尖点处的一种正则导数,给出了矩形网格情形分段光滑曲线网格边界点法向导数的一种插值计算方法.采用矩形网格有限差分法对复杂边界地球介质模型进行地震波场数值模拟,并采用波场系列快照技术揭示地震波在起伏地表和复杂介质中的传播规律.模拟结果表明：法向导数插值计算方法为矩形网格有限差分法处理复杂边界提供了有效途径,采用波场系列快照技术可以清晰地展现地震波在反射界面的反射和透射规律、在尖点的绕射规律以及在自由表面的直达波和多次反射规律.     

A discontinuous Galerkin method for seismic wave propagation in coupled elastic and poroelastic media

Numerical simulation in coupled elastic and poroelastic media is important in oil and gas exploration. However, the interface between elastic and poroelastic media is a challenge to handle. In order to deal with the coupled model, the first-order velocity–stress wave equations are used to unify the elastic and poroelastic wave equations. In addition, an arbitrary high-order discontinuous Galerkin method is used to simulate the wave propagation in coupled elastic–poroelastic media, which achieves same order accuracy in time and space domain simultaneously. The interfaces between the two media are explicitly tackled by the Godunov numerical flux. The proposed forms of numerical flux can be used efficiently and conveniently to simulate the wave propagation at the interfaces of the coupled model and handle the absorbing boundary conditions properly. Numerical results on coupled elastic–poroelastic media with straight and curved interfaces are compared with those from a software that is based on finite element method and the interfaces are handled by boundary conditions, demonstrating the feasibility of the proposed scheme in dealing with coupled elastic–poroelastic media. In addition, the proposed method is used to simulate a more complex coupled model. The numerical results show that the proposed method is feasible to simulate the wave propagation in such a media and is easy to implement.     

TTI介质的交错网格伪P波正演方法

研究了三维弱各向异性近似下,利用伪P波(伪纵波)模拟弹性波场P分量在倾斜对称轴的横向各向同性(TTI)介质中的传播过程,并对比了分别基于弹性Hooke定律、弹性波投影和运动学色散方程所建立的三种二阶差分伪P波方程的正演特点.目前这些伪P波方程数值计算主要采用规则网格差分,但是规则网格在TTI模拟中有低效率、低精度以及不稳定的缺点.为了提高计算的精度,本文构建出相应方程的交错网格有限差分格式.通过对比伪P波方程在三维TTI介质中不同的数值模拟的表达形式,本文认为基于色散方程所建立的伪P波方程在模拟弹性波中P波传播的过程中具有最小的噪声.本文分析不同的各向同性对称轴空间角度的频散特征,并引入适当的横波速度维持计算的稳定.二维模型算例表明,本文提出的交错网格正演算法可以得到稳定光滑的伪P波正演波场.使用本文交错网格算法对二维BP TTI模型的逆时偏移也具有较稳定的偏移结果.     

用于地震波场模拟的变网格边长声格固体模型

提出一种新的用于地震波场模拟的变网格边长声格固体模型（Phononiclatticesolidwithvariousgridlength,简称PLSVL模型）．该模型通过改变不同介质中网格的边长来体现介质速度的变化,粒子在所有网格中的运动均保持每一步运行一个网格边长．详细推导了该模型的Boltzmann方程,证明该方程及其散射项的表达式与PLS模型相同,因而从该方程出发,可寻出宏观变量所满足的波动方程．新模型在同种网格内部的传输过程不存在误差,在界面产生的误差不会随时间而积累．本文同时给出了一个理论模型的波场模拟结果．     

深度均匀采样梯形网格有限差分地震波场模拟方法

由于重力引起的岩石压实效应,一般来说,地震波传播速度由浅入深整体逐渐增大.梯形坐标系设计可耦合速度由浅入深逐渐增大的变化,该坐标系中均匀网格采样所对应的物理直角坐标系网格由浅入深逐渐增大,也即浅部低速区对应细网格,深部高速区对应粗网格.在梯形坐标系表征波动方程后利用有限差分求解,本文实现一种深度均匀采样、横向采样间隔随深度增加逐渐线性增大的有限差分地震波模拟方法.梯形坐标系波动方程离散后,仍采用常规均匀网格有限差分算法对其求解.由于横向网格大小由浅入深线性增加,本方法可避免不同大小网格区域过渡所产生的虚假反射.梯形坐标系波场模拟浅层精度高,深层横向响应范围广,可有效减少有限差分网格数量.本文提出的方法是在更广义的坐标系下利用有限差分求解波动方程,正交坐标系仅为该梯形坐标系之特例.本文旨在为大速度动态范围深地高效高精度地震波场模拟提供一种思路.     

A Hamiltonian Particle Method with a Staggered Particle Technique for Simulating Seismic Wave Propagation

We present a Hamiltonian particle method (HPM) with a staggered particle technique for simulating seismic wave propagation. In the conventional HPM, physical variables, such as particle displacement and stress, are defined at the center, i.e., at the same position, of each particle. As most seismic simulations using finite difference methods (FDM) are practiced with staggered grid techniques, we know the staggered alignment of space variables could improve the numerical accuracy. In the present study, we hypothesized that staggered technique could improve the numerical accuracy also in the HPM and tested the hypothesis. First, we conducted a plane wave analysis for the HPM with the staggered particles in order to verify the validity of our strategy. The comparison of grid dispersion in our strategy with that in the conventional one suggests that the accuracy would be improved dramatically by use of the staggered technique. It is also observed that the dispersion of waves is dependent on the propagation direction due to the difference in the average spacing of the neighboring two particles for the same parameters, as is usually observed in FDM with a rotated staggered grid. Next, we compared the results from the conventional Lamb’s problem using our HPM with those from an analytical approach in order to demonstrate the effectiveness of the staggered particle technique. Our results showed better agreement with the analytical solutions than those from HPM without the staggered particles. We conclude that the staggered particle technique would be a method to improve the calculation accuracy in the simulation of seismic wave propagation.     

Pseudo‐spectral method using rotated staggered grid for elastic wave propagation in 3D arbitrary anisotropic media

Staggering grid is a very effective way to reduce the Nyquist errors and to suppress the non‐causal ringing artefacts in the pseudo‐spectral solution of first‐order elastic wave equations. However, the straightforward use of a staggered‐grid pseudo‐spectral method is problematic for simulating wave propagation when the anisotropy level is greater than orthorhombic or when the anisotropic symmetries are not aligned with the computational grids. Inspired by the idea of rotated staggered‐grid finite‐difference method, we propose a modified pseudo‐spectral method for wave propagation in arbitrary anisotropic media. Compared with an existing remedy of staggered‐grid pseudo‐spectral method based on stiffness matrix decomposition and a possible alternative using the Lebedev grids, the rotated staggered‐grid‐based pseudo‐spectral method possesses the best balance between the mitigation of artefacts and efficiency. A 2D example on a transversely isotropic model with tilted symmetry axis verifies its effectiveness to suppress the ringing artefacts. Two 3D examples of increasing anisotropy levels demonstrate that the rotated staggered‐grid‐based pseudo‐spectral method can successfully simulate complex wavefields in such anisotropic formations.     

Wave dispersion studies in granular media by analytical and analytical–numerical methods

The phenomenon of wave dispersion in dry sand is studied both by purely analytical studies and by analytical–numerical experiments on the basis of gradient elastic and viscoelastic material models. These material models are employed in order to simulate the microstructural characteristics of dry sand. The analytical studies treat the material body as a one-dimensional (for the viscoelastic case) and three-dimensional (for the gradient elastic case) and for both material models provide explicit expressions for the velocity of propagation of harmonic compressional (P) and shear (S) waves. These velocities are found to be functions of frequency, i.e., dispersive. The analytical–numerical studies treat the material as a one-dimensional one and try to simulate P and S wave propagation along the axial direction of cylindrical dry sand specimens. Thus, a sinusoidal pulse with a specific frequency is applied at one end of the specimen and the response is determined at some other point by solving a transient dynamic boundary value problem with the aid of a numerical Laplace transform. This analytical–numerical experiment is repeated for various frequencies. Thus, one determines the velocities of P and S waves as functions of frequency, thereby proving again that wave propagation in dry sand is dispersive.     

地震波数值模拟的褶积微分算子法与伪谱法的对比分析

将基于Forsyte广义正交多项式的褶积微分算子法运用于复杂非均匀介质地震波场模拟中,并将计算结果与伪谱法计算结果进行分析比较。通过二者的计算时间对比发现:在同样的计算条件下,褶积微分算子法的采样时间始终小于伪谱法,这是其进行地震波数值模拟的一个明显优势。通过波场快照的对比,褶积微分算子法的模拟结果与伪谱法数值模拟结果的频散效应相当,可为地震波场的值计算提供一种新的选择。     

错格实数傅里叶变换微分算子及其在非均匀介质波动传播研究中的应用

提出一种新的数值微分运算方法，即错格实数傅里叶变换微 分法. 该方法的运算速度 比错格复数傅里叶变换数值微分解法快0.33倍；因为该微分算法在整个微分运算过程中保留 了奈奎斯特分量，使得它比普通分格的实数傅里叶变换数值微分算法的精度高，稳定性好. 将该方法和Cagniard De Hoop解析法在求解半无限空间地震波动的问题中进行比较，结果 表明，新微分法的精度和解析方法的精度相同. 在非均匀介质中的地震波传播数值模拟的结 果表明，该方法是一种研究非均匀介质中地震波传播问题的有效的数值微分方法.     

低速地层裸眼井中弹性波的传播

本文在文献[4]的基础上分析了低速地层（地层横波速度低于井内流体速度）裸眼井中由点源激发的弹性波的传播。指出在这种情况下井内传播的波仅有共振纵波、斯通利波和衰减的反射波,共振横波和伪瑞利波将不再出现。由数值计算得出了相速度和群速度频散曲线,给出了势函数的振幅表达式和反射波随井轴的衰减规律。     

黏滞声波地震波场模拟的SG-ONADM混合方法

交错网格方法（SG）和最优近似解析离散化方法（ONADM）是两类典型的地震波场数值模拟方法.两类方法各有其优势,相对于ONADM方法,SG方法在单个时间层内的计算更为简单;相对于SG方法,ONADM方法可以在较大空间步长条件下有效压制数值频散.结合两种方法的优势,本文提出了一种新的地震波场模拟方法（SG-ONADM）.该方法对控制方程中的一阶偏导数采用SG方法给出的一阶偏导数近似公式,高阶偏导数采用ONADM方法给出的高阶偏导数逼近公式.理论分析及数值算例表明,SG-ONADM方法保留了两种方法的优势,不仅能在较大空间步长条件下有效压制数值频散,同时具有较低的内存需求量;在对同一计算区域进行波场模拟时,SG-ONADM方法的计算效率要高于SG方法和ONADM方法.最后,我们使用SG-ONADM方法进行黏滞声波波场模拟,研究了黏滞声波在复杂介质中的传播.     

基于牛顿迭代法求解群速度的地震波射线追踪模拟

地下介质中普遍存在着各向异性,当前基于各向异性的地震波射线追踪多是在弱各向异性介质中进行且采用群速度近似表示方法,这些近似方法在强各项异性介质中会导致很大误差而无法真正模拟地震波的传播规律。根据地下普遍存在各向异性的事实和地震波基本传播规律,提出利用牛顿迭代法高效求解群速度,基于Paraview平台自动化构建三维地质模型,采用最短路径法进行地震波射线追踪模拟及可视化,实现对复杂三维地质的速度不均匀性和各向异性的表达,为三维地质模型的构建和地震波射线追踪模拟及可视化提供一种新思路,并以华北克拉通山西断陷带北部局部区域为例进行研究。结果表明,该方法能够减少由各向异性对地震波传播模拟造成的影响,清晰表达了研究区地质结构和各向异性特点,在对复杂三维地质结构的解读中能够较好应用。     

复杂地表边界元-体积元波动方程数值模拟

复杂近地表引起来自深部构造的地震反射信号振幅和相位的异常变化,是影响复杂近地表地区地震资料品质的主要原因.本文采用边界元-体积元方法,通过求解含复杂地表的波动积分方程,来模拟地震波在复杂近地表构造中的传播.其中,边界元法模拟地形起伏和表层地质结构对地震波传播的影响;体积元法模拟起伏地表下非均质低降速层的影响.与其他数值...     

Mesh grading approach for wave propagation in high velocity-contrast media

A mesh grading approach based on investigated lump method has been presented for simulating wave propagation in high velocity-contrast media. Unstructured fine grids are used to discretize the low wave-velocity medium in order to ensure the accuracy of numerical computation, and unstructured coarse grids are used for the high wave-velocity medium in order to substantially reduce the computational cost. On the interface, one coarse grid can match the fine grids of arbitrary odd number. The key feature of the proposed method is the constructions of investigated lumps on the interfaces of media. The transition zone, which is commonly used in the discontinuous grid scheme based on the staggered-grid finite-difference method, will not be used any more. Moreover, the computational instability that the discontinuous grid schemes frequently encountered does not arise in the proposed method. The comparisons with the analytical solutions and the application in studying the effects of sedimentary basin demonstrated that the mesh grading approach is a valid, accurate, convenient and flexible algorithm in simulating wave propagations in high velocity-contrast media with irregular interfaces.     

六分量弹性波场数值模拟与分析

地震波传播过程中,质点的振动不仅包括三个独立的平移部分,还包括三个独立的旋转部分.本文基于一阶速度-应力弹性波方程,采用分裂完全匹配层(SPML)的吸收边界条件,推导了时间导数二阶精度和空间导数高阶精度的交错网格有限差分格式的弹性波速度与应力各分量计算公式,模拟了各向同性介质中均匀模型和层状模型下的六分量波场,并对二维各向同性层状模型下的三个分量地震记录做高分辨率线性拉东变换得到各自的频散能谱.数值模拟分析结果表明:(1)旋转分量的能量要比平动分量弱的多;(2)在平动分量上,面波能量强,频率低,反射P波能量较强,反射S波能量稍弱;在旋转分量上,反射P波能量很弱,S波能量强;(3)与平动分量相比,旋转分量的频散能谱效果更好,能看到基阶和完整的高阶面波,即旋转分量能反映更多的地下介质信息.     

地震应力波的传播与岩石的动态响应

本文用 Split-Hopkinson 压力杆装置研究了应力波通过破裂岩石试件的衰减现象。给出了灰砂岩、石灰岩两种岩石动态的破裂强度与应力应变关系.结合受应力作用的岩样切片的微观观察,对应力裂纹、裂纹的动态分叉现象进行了分析讨论.      

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

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

用于声波方程数值模拟的时间-空间域有限差分系数确定新方法

声波方程数值模拟已广泛应用于理论地震计算,同时构成了地震逆时偏移成像技术的基础.对于有限差分法而言,在满足一定的稳定性条件时,普遍存在着因网格化而形成的数值频散效应.如何有效地缓解或压制数值频散是有限差分方法研究的关键所在.为精确求解空间偏导数,相继发展了高阶差分格式优化方法和伪谱方法.近期,为更好地缓解数值频散,提出了时间-空间域有限差分方法,该方法采用了泰勒展开近似方法来确定有限差分格式系数,因而只能保证在一定的小范围内很好的拟合波场传播规律.为进一步压制数值频散效应,本文引入了时间-空间域特定波数点满足频散关系的方法,根据震源、波速和网格间距确定波数范围,同时考虑了多个传播角度,然后建立方程确定了相应的有限差分格式系数,使得差分系数能在更大范围符合波场传播规律.通过频散分析和正演模拟,验证了本文方法的有效性.