水平层状介质中的快速两点间射线追踪方法

提出了一种快速准确的水平层状介质中两点间的射线追踪方法. 该方法对任意复杂的水平层状速度结构、 任意大小的入射角都能稳定、 快速、 准确地收敛. 将本方法分别与伪弯曲法和Kim方法进行了比较. 结果表明: 本方法比伪弯曲法的计算效率和计算精度更高, 但是只适用于一维速度结构; 在入射角不大的情况下, Kim方法与本方法效果相当. 在大入射角的情况下, Kim方法失效, 但本方法能快速稳定地收敛到真解. 本文提出的射线追踪方法为很多实际地震问题提供了强有力的工具, 但适用范围仅限于水平层状介质.      

复杂速度模型的地震交切定位方法

交切法是最基本的地震定位方法之一,该定位方法具有稳定性强、速度快的优点。由于速度模型假设为均匀或水平均匀介质,与实际介质有较大偏差,因此定位精度较低,特别是震源深度。为克服传统交切法的缺点,我们对其进行了改进。在改进的方法中,震源轨迹不再假设为圆形或双曲线形,而是利用最小走时树射线追踪方法进行精确计算。数值模型计算表明,本文提出的改进方法,可用于复杂速度模型的地震定位。     

基于射线追踪的微地震多波场正演模拟

对于微地震正演模拟,本文以射线追踪的原理为基础,对两点间的射线追踪问题进行了研究,应用二分算法、改进二分算法和微变网格算法对水平层状匀速模型、弯曲层状匀速模型和复杂地质模型进行射线追踪,使得计算效率和适用范围都得到了很大的改善.文中对每种算法误差范围和计算效率进行了对比验证,对于不同的地质模型,选用合适的算法才能在计算速度和精度上得到双重保证,最后正演模拟了多波三分量记录.在模型建立上引入了超薄层概念,并在前人模拟的直达波、透射波、反射波基础上拟了折射波,使正演模拟的多波场信息更丰富.文中的应用实例及模型结果表明:与二分法相比,改进二分法能够对弯曲界面进行射线追踪,并能保证结果的精度.弯曲层状模型中,改进二分法与微变网格法相比计算速度有显著提高,能够应用到资料的反演中.     

运用改进的地震定位交切法重定位2001—2010年云南地震

传统地震定位方法利用震源轨迹确定震源位置,但基于均匀或横向均匀介质模型必然导致定位误差。为此对传统方法进行改进,发展适用于三维复杂地壳速度模型的地震定位交切法。利用最小走时树射线追踪技术,以离散方式准确计算三维复杂地壳速度模型中的震源轨迹,将震源定位于震源轨迹交汇的密集点。将该方法应用于云南地区地震重定位,得到较高定位精度。     

测井约束地震波形反演的同伦摄动法

测井数据和地震数据是地震勘探中两种最重要的资料. 测井约束地震波形反演是在非线性波形反演的基础上,利用已知测井资料详细的垂直分辨能力和地震资料均匀密集的水平采样特点, 通过迭代反演来求取一个具有较高分辨率的速度参数.本文建立了测井约束反演模型,研究了测井约束下地震波形反演的同伦摄动求解方法.同伦摄动法作为一种新的、求解数学物理中各种非线性问题的有效方法,具有计算速度快、计算精度高的优点.这对于提高反演的精度和效率是十分有益的. 为了表征该方法的有效性和稳定性,分别对水平层状介质模型和逆冲断层带模型进行了数值模拟,并与Landweber迭代法相对比,结果表明该算法具有更好的收敛性,能够取得更为满意的反演效果.     

水平层状VTI介质高效射线追踪计算方法

准确、快速计算水平层状VTI(vertical transverse isotropic)介质的走时和射线路径是实现VTI介质地震叠前时间偏移和微地震监测速度建模的基础.水平层状VTI介质射线追踪的常用算法包括打靶法和最短路径法.打靶法原理简单、实现容易,但在长偏移距时存在微小初射角度变化导致射线路径巨大扰动的问题;最短路径法原理直观,但需要沿射线群速度方向(群角)计算走时,计算特定群角对应的群速度值需要先搜索对应的相速度角度(相角),显著降低了计算效率.本文综合打靶法和最短路径法的优点,从震源同时出射一系列射线,以这些射线到达界面的交点构成稀疏的网格节点,自适应加密实现检波点周围网格节点间距高于精度要求,采用插值算法获得检波点的走时和动态网格方法获得不位于网格节点上的检波点的射线路径.本文方法严格计算已知相角对应的群角和群速度,未使用弱各向异性群速度近似公式,适用于任意强度各向异性VTI介质qP、qSV和qSH直达波和反射波射线追踪;以相角确定震源出射射线,不用遍历群角和群速度对应关系.同传统最短路径方法的数值实验对比表明本文方法具有高精度和高效性,非常适合于需要多次正演计算的地震叠...     

层状TI介质中微地震定位和各向异性速度结构同时反演

微地震监测被广泛应用于非常规油气资源的水力压裂作业、油藏描绘和水驱前缘监测工程中.微地震定位采用的初始速度模型一般是基于地震测井记录和射孔数据建立,该速度模型的不准确性易引起定位误差.为降低这种定位误差,本文发展了一种微地震定位和各向异性速度结构同时反演的方法.研究对象为1-D的层状TI介质,其中对称轴方向任意.利用改进的分区多步最短路径算法计算qP、qSV和qSH波的到达时间和射线路径,结合共轭梯度法求解带约束的阻尼最小二乘问题.数值模拟结果表明,该算法能同时进行各向异性速度结构模型(每层的Thomsen参数和界面深度)和微震震源参数(空间坐标和发震时刻)的反演,并且对随机噪声不敏感,有利于实际工程应用.     

吴忠,灵武地区地壳中上部速度结构

应用射线追踪方法进行震源参数和速度结构联合反演,对宁夏吴忠、灵武地区地壳中上部Ｐ波速度结构进行了研究,茯结果与目前该区域地震定位所用的速度模型差异较大,与穿研究区的人工地震剖面相比,随深度的增加两者的速度值趋于一致。     

一种基于插值技术高精度计算稀疏网格地震定位中震源轨迹的方法

区域和全球地震定位越来越多地基于更接近实际的横向非均匀速度模型.速度模型主要来自于地震体波层析成像结果,分辨率不是很高.这样,模型宜于以稀疏网格剖分以减少计算时间和计算机内存需求.当剖分的模型单元较大时,基于射线追踪技术计算复杂介质地震定位中震源轨迹的方法——选取震源轨迹所经过单元节点(位于单元中心,称为震源轨迹节点)为轨迹参考点,利用最小走时树射线追踪方法计算绝对残差场中连接轨迹参考点的射线路径作为震源轨迹——计算结果误差较大,难以满足精定位需要.针对该问题,本文对其进行了改进:不将震源轨迹节点作为轨迹参考点,而是基于插值技术计算每个轨迹节点其法线点对(即该节点与其周围残差正负极性不同的相邻节点组成的点对中梯度绝对值最大的那对)间残差为零的点作为震源轨迹参考点.算例表明:和原方法相比,改进方法计算的震源轨迹更为精细,计算精度提高数十(线性插值)至数百倍(非线性插值),而计算效率基本保持在同一数量级,使利用震源轨迹进行直观、快速和高精度的区域或全球地震事件定位成为可能;壳幔界面反射纵波(PmP)对震源的约束和直达纵波(Pg)相似;同一台站PmP-Pg波到时差约束的震源轨迹对震源深度有很好的约束.     

TI介质射线追踪及其高斯束成像应用

射线追踪是一种高频近似前提下快速有效的波场近似计算方法,传统的基于弹性参数的各向异性介质运动学和动力学射线追踪方程,求解过程中需要处理Christoffel方程的特征值问题,因而计算效率比较低.为了解决这一问题,本文通过引入相速度和群速度,对运动学和动力学追踪方程进行修改和简化,有效地提高了各向异性介质射线追踪算法的计算效率.另外,我们将该算法应用到各向异性偏移中,实现了共炮域TI介质高斯束叠前深度偏移方法.VTI介质Hess模型和TTI介质洼陷模型的试算结果说明了该方法的正确性和有效性.