Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.     

初至折射波静校正技术在复杂探区地震数据处理中的应用

山区及黄土塬地区近地表结构复杂，地形起伏大、低降速带横向变化剧烈，地震资料处理中的静校正难度加大。长波长静校正容易引起构造形态畸变，短波长静校正影响地震剖面信噪比。针对此类问题，从长波长静校正问题及短波长地表非一致性问题分析入手，提出部分炮检距叠加识别长波长静校正技术及分段统计逼近地表非一致性校正技术。该技术在复杂地表区与其它相关技术综合应用，可有效解决长、短波长静校正问题。     

斯通利波和渗透率的关系分析

理论和实际资料证明,在渗透性地层,低频斯通利波出现的频散和能量衰减与地层渗透率有关,因此可利用斯通利波的传播速度和能量衰减等特性来估算地层渗透率。斯通利波的能量和渗透率负相关。因为斯通利波能量代表了斯通利波的有关信息,因此可以认定斯通利波和渗透率负相关。     

地震波混合阶褶积算法模拟

基于正反傅立叶变换,提出了地震波模拟的混合阶褶积算法。该方法原理简单、易于实现,结合了有限差分与伪谱法的优点,具有较高的精度和计算效率,适用于地震波场的正演计算,同时还给出了二维地震波场的理论计算实例。计算结果表明,此算法模拟结果正确、精度高、速度快、能适应较为复杂地质模型,并且易于推广到各向异性介质中去。     

地面地震层析技术在金属矿勘查中的试验研究

由于金属矿区地质构造的复杂性和复杂的地表条件, 适合于探测层状介质的地震反射方法技术在金属矿勘查中的应用受到了限制。为提高地震记录的横向分辨率, 探测复杂的地质构造及隐伏岩体等, 需研究相应的地震方法技术, 地面地震层析技术在金属矿勘查中的试验研究就是为解决上述问题而开展的。本文就非线性地面地震层析技术的基本原理进行了阐述, 并就该方法在金属矿勘查中的试验效果进行了讨论。