横向各向同性介质中地震波走时模拟

横向各向同性介质是地球内部广泛分布的一种各向异性介质.针对这种介质,我们对各向同性介质的最小走时树走时模拟方法进行了推广,推广后的方法可适用于非均匀、对称轴任意倾斜的横向各向同性介质模型.为保证计算效率,最小走时树的构建采用了一种子波传播区域随地震波传播动态变化的改进算法.对于弱各向异性介质,我们使用了一种新的地震波群速度近似表示方法,该方法基于用射线角近似表示相角的思想,对3种地震波（qP, qSV和qSH）均有较好的精度.应用本文地震波走时模拟方法对均匀介质、横向非均匀介质模型进行了计算,并将后者结果与弹性波方程有限元方法的模拟结果进行了对比,结果表明两者符合得很好.本文方法可用于横向各向同性介质的深度偏移及地震层析成像的深入研究.     

起伏地表条件下各向异性地震波最短路径射线追踪

在地震波正反演研究中,考虑起伏地表和地震各向异性具有非常重要的理论意义和实际应用价值.本文在前人研究的基础上,将最短路径追踪算法引入到起伏地表各向异性介质模型的地震波走时计算中.模型剖分时,整体模型划分成正方形单元,起伏边界附近以不规则网格逼近,进而采用非规则节点布置实现非规则网格处的最短路径计算.追踪计算中采用Sena群速度近似公式,得到各向异性地震波的走时,实现了复杂地表情况下各向异性介质模型中地震波的射线追踪.理论模型计算结果显示,本文方法能够可靠地应用于复杂各向异性介质模型,具有较高的计算精度.     

层状各向异性介质转换波克希霍夫叠前时间偏移

在克希霍夫叠前时间偏移处理中,地震波走时的计算方法是决定大偏移距地震资料成像品质的重要因素.在常规的三维转换波各向异性叠前时间偏移公式中,走时的计算是基于等效单层各向异性介质的非双曲线方法.用这种方法处理的成像道集,在偏移/深度比超过一定阈值后,成像道集中的反射同相轴将出现过偏现象,这种偏移不平的同相轴将影响偏移叠加的最佳响应,使得偏移成像波组呈低频化特征,最终降低三维转换波偏移成像质量.我们采用层状介质的走时计算方法代替常规算法,并且利用了常规方法的转换波各向异性偏移速度模型.基于层状介质的算法能够提高大偏移距转换波走时计算精度,克服中浅地层大偏移距远道成像道集中反射同相轴逐渐上翘的问题.两个地区的三维转换波资料处理结果证实,基于层状各向异性介质的转换波克希霍夫叠前时间偏移方法,明显改善了反射成像剖面的连续性和分辨率,提高成像剖面构造的可解释性.     

求解二次方程的二维VTI介质qSV波和qSH波走时快速扫描算法

地震波走时计算在观测系统设计、偏移成像、速度模型走时反演和地震定位等方面起到重要作用.各向异性广泛存在于地球介质中,影响地震波传播的振幅和走时,忽略各向异性的影响将对成像、反演以及地震定位等造成一定的误差.因此对于高分辨率成像和反演,走时计算中考虑各向异性十分重要.快速扫描法不需要存储和追踪波前面信息,在各向异性初至波...     

二维复杂黏弹性各向异性介质中的地震射线追踪

实际地层中传播的地震波普遍存在速度各向异性和能量衰减现象,因此基于弹性介质假设条件下的地震射线正反演算法具有一定的局限性,而研究黏弹性各向异性介质中的地震波传播规律可为揭示地下结构提供更加可靠的理论依据.射线追踪技术是揭示高频地震波传播特性的有效手段之一,然而绝大多数研究仅限于弹性介质.针对黏弹性各向异性介质,本文首先给出了一种射线速度和振幅衰减的计算方法,然后结合改进型最短路径算法,在实空间内实现了计算复杂介质模型中地震波(qP,qSV,qSH)的射线路径和传播走时以及能量衰减(虚走时).该算法适用于复杂黏弹性各向异性介质.误差分析结果显示,实走时的最大相对误差小于0.13％,虚走时的最大相对误差小于0.55％,表明该算法具有较高的计算精度.     

各向异性弹性波动方程多分量联合叠后逆时偏移

为了更好地描述地下介质的各向异性特性,实现准确的波场成像,开展了各向异性介质逆时偏移研究。采用高阶交错网格有限差分法推导出各向异性弹性波动方程叠后逆时深度偏移方程,证明各向异性介质向各向同性介质转化的条件,并研究Thom son参数对弹性波场的影响,合成了各向同性介质和各向异性介质多分量零偏移距剖面,实现了多分量模拟记录联合叠后逆时偏移。计算结果表明,能够准确地实现多分量波场叠后逆时成像,使地质层位中的断层、断点等目标成像清晰准确,且偏移成像精度较高,因此各向异性逆时偏移法是一种高精度的深度域偏移成像方法,可以指导实际天然地震资料和人工地震资料的数据处理。     

各向同性介质长偏移距地震同相轴动校正

传统二阶动校正方法基于较小最大偏移距与目标层深度比和地震波沿直线传播假设,进行长偏移距地震资料处理时,这些假设不再成立.高阶项动校正公式能提高长偏移动校正精度,文中对几种典型的高阶项动校正方法进行了比较,并提出了优化四阶、优化六阶动校正方法.模型计算表明,高阶项动校正方法能取得较常规动校正方法好的动校正结果,但并非阶数越高动校正精度就越高;在纵向速度变化剧烈时,高阶动校正或优化高阶动校正方法一般不能适用于最大偏移距与目标层深度大于3.5的地震反射同相轴,优化四阶和优化六阶动校正公式由于考虑了无穷大偏移距的影响,具有更稳定、更加精确动校正效果,适合于实际的各向同性长偏移距地震资料处理.     

曲线坐标系因式分解程函方程及其走时计算

地震波走时广泛应用于静校正、层析成像、Kirchhoff偏移成像、地震定位等研究.复杂地表条件是影响走时计算精度的重要因素.近年来,发展的曲线坐标系程函方程为精细刻画起伏地表条件下的地震波走时场特征提供了新的思路.然而,基于有限差分程函方程的求解方法不可避免地受到震源奇异性的影响,即震源附近波前的曲率较大,此时使用平面波近似假设的差分格式会导致较大误差.而震源误差会随着波前的传播到达整个计算区域,从而影响整个区域的求解精度.针对该问题,本文借鉴因式分解的思想,推导建立了曲线坐标系因式分解程函方程,并针对性地发展了其数值求解方法,从根源上解决了复杂模型走时计算中的震源奇异性问题.数值实例表明因式分解法能够有效降低震源误差,显著提高起伏地表走时计算的精度和效率,为起伏地表地震波走时计算提供更佳的选择,在复杂模型的地震资料处理中展现出广泛的应用前景.     

各向异性介质弹性多波高斯束正演模拟

射线类正演方法以其高效性和灵活性的特点,被广泛应用于地震勘探中.然而,普通射线类正演方法存在焦散区、阴影区和多值走时等问题,计算精度不够理想.为此,本文在传统的各向同性介质声波高斯束正演的基础上,推导了各向异性介质运动学和动力学射线追踪方程,发展了各向异性介质弹性多波射线追踪算法;并将该算法应用到高斯束正演模拟中,实现了一种各向异性介质弹性多波高斯束正演模拟方法.各向异性VTI介质断块和VTI介质复杂构造模型试算的结果表明:本文研究的方法能够对各向异性介质构造进行正演模拟,在保证计算精度的前提下,具有较高的计算效率,模型试算的结果说明了方法的有效性和正确性.     

VTI介质中P-SV波转换点与各向异性参数关系

对于转换波地震勘探中的转换点位置这个重要问题,提出转换点位置不仅与纵横波速度比,偏移距深度比以及源检距有关,还与地下介质的各向异性的性质有关,计算了忽略地下介质的各向异性影响对转换点的确定带来的严重误差,从而影响地下地质体的精确成像.通过对层状VTI介质中的转换点近似方程的推导过程,引入该方程不同于传统方程的导出是对层状各向同性介质而言,该方程通过引入各向异性参数,使我们对转换波可以有进一步的认识,拓展了转换波处理中各向异性的应用.该方程对于偏移距深度比小于3.0的情况是比较准确的,这对于大偏移距转换波勘探具有实际意义.