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.     

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

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

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

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

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

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

利用钻孔超声成像的图像特征进行岩石风化程度分类

钻孔超声成像方法能精确地反映井壁的光滑情况, 进而判断岩石的风化发育程度, 确定风化层的厚度。阐述了钻孔超声成像的走时和振幅图像特征, 给出了1个在香港地区的应用实例。