Elastic reverse time migration based on vector wavefield decomposition

Prestack elastic reverse time migration ( RTM) requires multicomponent seismic data .But for multi-component elastic Kirchhoff migration , there is a limitation that ray theory no longer applies if thegeology be-comes complicated .In this paper, the authors have created a new 2D migration context for isotropic , elastic RTM, which included decomposition of the elastic source and receiver wavefields into P and S wave vectors by decoupled elastodynamic extrapolation , which retained the same stress and particle velocity components as the input data .Then we appliedsource-normalized crosscorrelation imaging condition in elastic reverse time migra-tion to compensate the energy of deep strata .We found that the resulting images were nearly identical to the ve-locity model , and the resolution has been improved .Our method is a wavefielddecomposition based on vector , and we can alsoavoid the problem of polarity reversal of converted shear wave imaging .It proved the applicabili-ty of the method proposed in our paper .     

基于f-x EMD波场分离的绕射多次波压制方法

自由表面多次波压制（SRME）方法是海洋地震数据处理中的关键步骤。当存在崎岖海底条件时,绕射波十分发育,常规SRME算法很难将复杂绕射多次波有效消除。对压制效果较差的多次波进行成因分析,首先基于f-x EMD波场分离算法,将全数据分解为反射部分和绕射部分,再对其进行多次波预测,获得绕射预测多次波子集：反射-绕射预测多次波、绕射-反射预测多次波和绕射-绕射预测多次波。然后,改进常规SRME算法的迭代过程,在SRME的迭代过程中使用与原数据中多次波更加匹配的绕射多次波子集进行多次波压制。与常规SRME算法的对比结果表明,改进的SRME算法在理论数据和实际数据处理中对多次波的压制都取得更好的效果。     

逆时偏移计算中的边界处理分析及应用

在地震资料的处理中,逆时偏移方法可以实现复杂构造高精度成像,但计算量和存储量两大问题影响了该方法的实际应用.逆时偏移算法中的边界处理方式的优化可以大大减少存储量,本文讨论了记录波场边界信息的吸收边界和随机边界这两种方法的计算效率和成像效果,提出了吸收边界中只记录边界范围内未衰减单层波场的方法并通过数值实验验证了其可行性,选择了不同的边界策略应用于模型数据和实际资料的处理.结果表明:记录单层波场边界信息的吸收边界的成像效果同传统存储波场历史的方式几乎无差别,但要额外存储每个时刻的波场边界信息;随机边界不需要额外存储波场信息,但会带来边界漫反射影响和计算区域的增加;记录单层波场边界可以明显减少存储量,并且不影响成像效果.     

2013年芦山MS7.0地震的震前及临震应变异常

距离芦山地震震中70 km的姑咱地震台YRY四分量钻孔应变仪2013年4月16—19日记录到8次幅度达到十倍固体潮幅的张性应变阶跃, 这在该仪器6年多运行期间尚属首次. 这些应变异常与芦山地震在时间、 空间上相关; 应变异常的主压应力方向与震前宝兴水压致裂地应力测量主压应力方向相近, 并与震前5个月中长达3个月的应变异常主压应力方向一致; 应变阶跃起始与远震震波同时到达. 综上认为, 4月16—19日的应变异常有可能是芦山地震发震断裂在破裂错动前地层蠕滑的反映.      

一种加权抛物Radon变换地震波场重建方法

基于带限正反抛物Radon变换的基本原理和计算方法,笔者给出了一种加权抛物Radon变换（PRT）算法,该方法可应用于叠前地震波场恢复和重建,如缺失偏移距地震数据的插值外推和反假频重采样,其权系数完全由数据域来确定,也就是说在迭代计算过程中,该权系数也在变化。对于均匀采样的地震数据,笔者给出了近偏移距外推和波场的反假频重采样的计算结果,试算表明笔者给出的加权PRT方法可有效地应用于地震波场重建中,并具有算法稳健、计算精度高和效率高的特点。     

井间地震物理模拟研究

井间地震资料具有极高的分辨率,但是,其波场十分复杂,特别是有效一次反射波场能量较弱,且往往被管波等强相干干扰淹没,波场的识别与分离难度较大.地震物理模拟技术是认识复杂地震波场的有效手段.我们通过单一地质体模型、复杂地质体模型和真实井间模拟模型等多种逼真地质模型的物理模拟,揭示了井间地震观测下直达波、透射波、反射波、折射波,多次波和导波等各类波型的特征.物理模拟试验表明,在炮检域平面上分析直达波的信噪比与能量分布,可以使处理人员对旅行时的误差有更直观的认识.所以识别初至与拾取旅行时间时,不仅要在共震源点道集(CSG)上分析拾取,而且还要在共炮检距道集(COG)上分析拾取,这样就提高了直达波识别的可靠性与拾取时间的精度.     

时间域全波场各向异性弹性参数反演

从各向异性弹性波的有限元正演方程出发 ,导出了反问题中时间域雅可比矩阵求解的计算公式 .它具有与时间域有限元正演方程相同的表达形式 ,故可通过有限元正演计算来获得雅可比矩阵 .研究了有限元正演算法的效率和精度、吸收边界条件等方面的问题 ,以提高反演系统的效率和精度 .在此基础上 ,实现了叠前全波场各向异性弹性参数反演 .计算表明 ,在初始模型偏离真实模型较大的情况下 ,层状模型和横向不均匀模型的反演结果均能准确地收敛到真实模型上 .