稳定的TTI介质拟声波逆时偏移及伪横波的联合压制策略

为克服各向同性和VTI介质逆时偏移方法对复杂地质构造成像的局限,研究了TTI介质拟声波逆时偏移方法.首先从精确的TTI介质频散关系出发,引入一个各向异性控制参数σ,推导了新的二阶耦合TTI介质拟声波方程,以保证波场延拓的稳定性;然后引入波场的伪速度分量,推导了等价的一阶拟声波方程.相比于规则网格有限差分法,交错网格有限差分(SGFD)法能够有效地压制数值频散,模拟精度更高;因此利用高阶SGFD法求解TTI介质一阶拟声波方程,构建逆时偏移所需的正向和逆时波场延拓算子,并应用归一化互相关成像条件实现精确的TTI介质逆时偏移成像.最后,简单讨论了伪横波的产生机制,并给出了伪横波的联合压制策略.模型试验结果验证了方法的有效性和稳定性.     

改进的TTI介质纯P波方程正演模拟与逆时偏移

逆时偏移作为一种高精度偏移方法已成为复杂构造成像的重要技术,描述纵波独立传播的延拓方程是各向异性介质逆时偏移的一个关键问题.在对VTI介质几个经典相速度近似公式回顾的基础上,针对常用于描述纯P波的Harlan近似公式在各向异性参数ε较大情况下近似精度较低的问题,本文对Harlan公式中的非椭圆项进行了修正,在非椭圆项前添加了一个与各向异性参数ε有关的修正系数,得到了三种改进型Harlan公式,并以近似精度最高的改进式为基础,推导了TTI介质纯P波方程.针对该伪微分方程,本文利用伪谱法和有限差分法联合实现波场延拓,对于常密度二阶方程,基于中心网格实现;对于一阶应力-速度方程则基于旋转交错网格实现.通过数值试验分析了TTI介质纯P波一阶应力-速度方程的近似精度,并以一阶纯P波方程为基础进行了TTI介质逆时偏移数值模拟试验.结果表明,本文给出的方法能够较准确地描述TTI介质纯P波波场特征,可以应用至各向异性介质逆时偏移.     

带正则化校正的TTI介质qP波方程及其逆时偏移方法和应用

各向异性研究对地下介质精确成像有着重要的意义,在当前计算机硬件迅速发展及宽方位地震数据采集日益普遍的情况下,成像必须考虑介质的各向异性.逆时偏移是基于双程波动方程的较为精确的数值解的成像方法,所以相对于其他地震成像方法,它具有很大的优势,譬如不受反射界面的倾角限制、偏移速度结构合适时能够使回转波及多次波正确成像.在各向同性介质中,可使用标量波方程来模拟波场.而在各向异性介质中,P波和SV波是相互耦合的,即不存在单纯的标量波传播,通常利用能代表耦合波场中P波分量运动学特征的拟声波(qP波)进行偏移成像.本文中,我们推导出了TTI介质下qP波控制方程.该方程可采用显式有限差分格式进行求解.通过声学近似,若沿对称轴方向的剪切波速度为零,对于对称轴方向不变且ε≥δ的模型来说,可得到稳定的数值解.但对于TTI介质来说,由于沿对称轴方向各向异性参数是变化的,声学近似会引起波场传播及数值计算的不稳定.因此,我们提出了正则化有限横波的方法,很好地解决了这一问题.最后,给出了Foothill模型的测试结果及某探区实际资料试算结果,展示了采用这个方程进行复杂TTI模型正演和高质量逆时偏移成像结果,证实了该方法的正确性和实际资料应用中的有效性.     

各向异性介质Low-rank有限差分法纯qP波叠前平面波最小二乘逆时偏移

拟声波最小二乘逆时偏移是一种极具潜力的地震波成像工具,但该方法遭受各向异性拟声波近似的限制,TTI介质正演模拟不稳定、反偏移记录中遭受伪横波二次扰动及数值频散假象,另外拟声波最小二乘逆时偏移还面临计算效率低、收敛速度慢、对速度等模型参数依赖性高等问题.为了克服各向异性拟声波最小二乘逆时偏移的缺陷,在反演框架下,本文借助Low-rank有限差分算法首次提出并实现了TTI介质纯qP波线性正演模拟及纯qP波最小二乘逆时偏移;为了进一步提升反演成像效率,同时改善反演成像方法对模型参数误差的依赖性及对地震数据噪声的适应性,通过引入叠前平面波优化策略,发展了TTI介质纯qP波叠前平面波最小二乘逆时偏移成像方法.在编程实现方法的基础上,通过开展模型成像测试,展示了本方法的优势和潜力：一方面加快了反演成像效率,另一方面也提升了方法的抗噪性,同时还降低了方法对模型参数的依赖性.     

TTI介质qP波逆时偏移中伪横波噪声压制方法

在对地下复杂构造介质,特别是盐丘侧翼及岩下区域进行成像时,相对于传统的各向同性逆时偏移和VTI逆时偏移,具有倾斜对称轴的TTI逆时偏移成像效果最优.不仅反射同相轴更加的连续,而且能量得到了更好的聚焦.传统的各向异性介质全弹性波RTM的计算量大且计算效率低.由于目前仍以纵波勘探为主,因此TTI逆时偏移qP波波动方程的选取显得尤为重要.为了提高计算效率,采用将沿着对称轴方向的横波速度设为零的方法,简化得到qP波波动方程.然而,这样会引入一种严重影响成像效果的低速度、低振幅的qSV波人为干扰.本文建立了qP波方程的完全匹配层控制方程,而后借助于辅助波场采用一种高效的压制伪横波噪声传播的方法,通过模型测试验证了该方法的有效性.     

TTI介质伪解析解耦波动方程

地球介质中最常见的一种各向异性介质就是倾斜横向各向同性(tilted transversely isotropic,简称TTI)介质,目前TTI介质的正演和逆时偏移是地震勘探工业界研究的前沿热点,TTI介质中的高精度地震波正演和逆时偏移必须要考虑TTI各向异性的影响,否则会使得模拟和成像结果不够准确.本文基于Tsvankin提出的TI介质相速度精确频散方程,针对常规TTI介质解耦波动方程时间精度的不足引入伪解析法(pseudoanalytical method,简称PAM)提出一种新的基于伪解析法的TTI介质解耦波动方程,能在时间方向和空间方向都达到高精度纯P波模拟.数值模拟结果显示,本文构造的基于伪解析法的TTI介质解耦波动方程具有高精度且没有任何伪横波误差.     

3DTTI介质逆时偏移计算关键问题研究

TTI介质逆时偏移是复杂构造区地震成像的一种重要有效技术.然而,其3D数据实际应用仍然存在一些挑战.本文以TTI拟声波方程为基础,系统梳理了TTI逆时偏移计算中的几个关键问题并给出解决方案.首先,P波与SV波耦合问题.本文系统分析了各向异性参数匹配方法伪横波压制效果以及其对波传播的影响;其次,为了提高计算效率,本文基于CUDA平台实现GPU加速计算方案.主要解决了两个方面的问题,其一为扩展了随机边界在TTI介质波场传播中的应用,确定了TTI介质速度和其他四个参数对应的"均匀各向异性"随机边界;其二,本文将单GPU拓展到多GPU实现三维TTI介质逆时偏移,并且采用了一种异步流技术,大大提高了GPU计算的高并发策略,使双GPU实现的加速比与单GPU相比接近理想的2∶1.最后通过3D模型试验证明了方法的效率和准确性.     

TTI介质有限差分逆时偏移的稳定性探讨

在沉积学中,可假设在相同时期的沉积层具有相近的物理性质和演化过程.因此,沿层传播的地震波和垂直于地层传播的地震波具有各向异性的特点.在纵波资料的处理中,考虑各向异性对逆时偏移的影响,通常假设介质的横波速度为零,这样可以得到纵波在TTI介质中的传播方程,但是该方程在实际计算中仍存在数值稳定性问题.本文加入横波分量可有效解决数值稳定性问题,并选取适当的横波速度减小对纵波成像的影响,实现地震波在TTI介质中的逆时偏移.实际测算表明,P-SV波的方程中包含横波分量,若假设SV的速度为零,则会导致方程的差分格式不稳定;若加入SV波,选择合适的SV波速度可以使SV波的全区各向异性和反射系数达到极小,并可有效的抑制SV波对纵波勘探的影响.本文的方法是一种稳定的TTI介质中的逆时偏移方法.     

基于拟声波一阶应力-速度方程的TI介质叠前逆时偏移

随着计算机硬件技术发展和对野外数据高分辨率勘探需求的增加,叠前逆时偏移技术逐渐的运用到实际生产中;又鉴于实际地层中各向异性的普遍存在性,故实现各向异性介质的叠前逆时偏移成像很有必要.文中从各向异性介质中弹性波的基本理论出发,导出各向异性介质中的拟声波一阶速度-应力方程,并基于此方程提出一套稳健的叠前逆时偏移策略.对逆冲模型的数值模拟表明在旋转角变化的TTI介质(具有倾斜对称轴的横向各向同性介质)中该拟声波方程稳定性要好于基于频散关系推导出的拟声波方程;VTI_HESS模型的逆时偏移结果表明较常规声波方程,新导出的拟声波方程成像更清晰,深处保幅效果更加明显.     

横向各向同性介质拟声波方程及其在逆时偏移中的应用

地下岩石的速度各向异性影响地震波的传播与成像.横向各向同性(TI)介质为最普遍的等效各向异性模型.引入TI介质拟声波方程可以避免复杂的弹性波方程求解以及各向异性介质波场分离,以满足对纵波成像的实际需要.本文从垂直横向各向同性(VTI)介质弹性波方程出发,推导出正应力表达的拟声波方程以及相应的纵波分量的表达式,进而分析从频散关系得到的拟声波方程的物理意义,而后将拟声波方程扩展到更一般的倾斜横向各向同性(TTI)介质中.波前快照与群速度平面的对比验证了拟声波方程可以很好地近似描述qP波的运动学特征.在此基础上,将拟声波方程应用在逆时偏移中并与其特例声波近似方程进行对比,讨论了计算效率、稳定性等实际问题.数值试验表明VTI介质情况下采用声波近似方程可以提高计算效率,而TTI介质qP-qSV波方程则在效率相当的情况下可以保证稳定性.SEG/HESS模型和逆冲模型逆时偏移试验验证了本文TI介质拟声波方程的实用性.