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

2.
各向异性介质qP波传播描述I:伪纯模式波动方程   总被引:6,自引:6,他引:0       下载免费PDF全文
程玖兵  康玮  王腾飞 《地球物理学报》2013,56(10):3474-3486
地球介质相对于地震波波长尺度的定向非均匀性会导致波速的各向异性,进而影响地震波场的运动学与动力学特征.各向异性弹性波动方程是描述该类介质波场传播的基本工具,在正演模拟、偏移成像与参数反演中起着关键作用.为了面向实际应用构建灵活、简便的各向异性波场传播算子,人们一直在寻求简化的各向异性波动方程.本文借鉴各向异性弹性波波型分离思想,通过对平面波形式的弹性波方程(即Christoffel方程)实施一种代表向波矢量方向投影的相似变换,推导出了一种适应任意各向异性介质、运动学上与原始弹性波方程完全等价,在动力学上突出qP波的新方程,即qP波伪纯模式波动方程.文中以横向各向同性(TI)介质为例,给出了相应的qP波伪纯模式波动方程及其声学与各向同性近似,并在此基础上开展了正演模拟和逆时偏移试验,展示了这种描述各向异性波场传播的新方程的特点与优势.  相似文献   

3.
借助Christoffel方程可求解出各向异性介质弹性波精确频散关系.利用近似方法进行处理,再通过傅里叶逆变换将频率波数域算子变换为时空域算子,可导出解耦的 qP波或 qS波波动方程.本文在 TTI介质弹性波精确频散关系的基础上,利用近似配方法推导了 qP波和 qSV波近似频散关系,通过傅里叶逆变换推导了 TTI介质 qP波和 qSV波解耦的波动方程.为了验证近似频散关系的有效性,利用两组模型参数对其进行数值计算,分析了相对误差在不同传播方向上的分布.随后使用有限差分方法分别对均匀、层状及复杂 TTI介质弹性波近似解耦波动方程进行数值模拟,结果显示 qP波和 qSV波完全解耦,并且在各向异性参数η<0 以及介质对称轴倾角变化较大的情况下,纯 qP波和纯 qSV波近似波动方程依然可以保持稳定.  相似文献   

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

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

6.
波动方程偏移成像阴影的照明补偿   总被引:6,自引:6,他引:6       下载免费PDF全文
受地下复杂构造和地震数据采集系统的影响,使地震波对地下目标的照明出现不均匀性,地震采集系统难以有效地获取地下某些目标的反射信息,进而使数据偏移成像在这些目标体上出现成像阴影. 根据波场和Green函数的窗口Fourier框架展开,利用角度域波动方程偏移成像和波动方程照明分析,并结合波动方程反演理论,提出一种角度域波动方程偏移成像阴影照明补偿方法. 这种补偿方法能同时考虑地震数据采集系统和波场传播路径对偏移成像的影响,消除复杂构造区的偏移成像阴影,改进波动方程叠前深度偏移成像在复杂构造区的成像效果.  相似文献   

7.
地震照明分析是地震数据采集和实际资料处理解释中的一种辅助手段,指在给定地下地质构造时,采用数值模拟的方法模拟出地震波在地下介质中的能量分布形态.自由表面相关多次波是海上地震数据中重要的波场信息,传统地震数据处理方法一直都将其视为干扰信息而进行去除,错失了很多有意义的构造信息.本文采用基于时间域有限差分法的双程波动方程进行波场数值模拟,将多次波能量考虑到波场照明中,通过单一透镜盐体模型试算表明,与传统一次波照明相比,多次波照明方法的波场信息更丰富、盐下构造照明更清晰、照明阴影区域更小.将多次波照明分析应用到Pluto速度模型,对比分析可以看出多次波照明分析方法更能适应复杂地质构造目标.通过多次波照明分析方法,我们能够对复杂地质构造的采集系统进行评价和优化,并预测地震偏移成像的质量.  相似文献   

8.
由于地层的各向异性与黏滞性,地震波在传播过程中表现出方向上的各向异性和振幅值的吸收衰减特征.若简单地将地下介质视为各向同性,实际地质结构的地震波响应可能会被曲解,甚至忽略了有效信息.为了得到深部储层的高精度成像,本文将标准线性固体模型(SLS)扩展到各向异性介质中,并结合各向异性拟声波方程,推导出时间空间域的各向异性介质黏声波方程组,该方程中用伪微分算子来表述地震波的衰减特征.在逆时偏移的过程中,为了解决高频不稳定问题,文中通过引入规则化算子来构建稳定的逆时传播算子.模型的正演模拟表明该方程既能准确的描述各向异性介质下的标量波的传播规律,又能体现地层的吸收衰减效应.  相似文献   

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

10.
逆时偏移是基于双程波波动方程进行波场逆时外推,再利用成像条件提取成像值,进而实现偏移成像的一种方法.目前,RTM已经由各向同性介质推广到了TI介质,包括VTI介质和TTI介质.考虑了各向异性因素的逆时偏移技术能更好地满足日益增长的勘探精度的要求,进一步提高和改善成像效果.在回顾各向异性介质逆时偏移发展历史的基础之上,本文针对于国内外TTI介质逆时偏移当前的研究现状进行了概述,并对目前TTI介质逆时偏移方法中所存在的特有的数值稳定性和qSV波人为干扰压制问题的解决方法进行了归纳总结.最后对逆时偏移未来的发展方向和研究重点进行了展望.  相似文献   

11.
尽管射线理论是基于波动方程高频假设下的一种近似理论,两者本应相互佐证共同发展完善,然而实际上确是平行独立的.鉴于此,本文提出了一种可以相互验证波场数值模拟或射线追踪数值模拟方法正确性和有效性的佐证方法.为了实现上述目的,我们首先将用于各向同性介质中多次波预测的分区多步不规则最短路径算法推广至各向异性TTI介质中,使其具备追踪各向异性TTI介质中多次波的能力.然后,通过伪谱法波场数值模拟和分区多步不规则最短路径算法预测得到的波场快照图、单炮地震记录剖面、以及合成理论地震图的分析对比,不但可以相互验证两种算法的正确性和有效性,同时对地震波在各向异性TTI介质中的传播规律和特点有了深刻的理解和认识.结果表明:本文提出的各向异性介质中多次波射线追踪算法是正确有效的.  相似文献   

12.
黏弹TTI介质中旋转交错网格高阶有限差分数值模拟   总被引:4,自引:2,他引:2       下载免费PDF全文
严红勇  刘洋 《地球物理学报》2012,55(4):1354-1365
以Carcione黏弹各向异性理论为基础,给出了适用于黏弹性具有任意倾斜对称轴横向各向同性介质(黏弹TTI介质)的二维三分量一阶速度-应力方程,采用旋转交错网格任意偶数阶精度有限差分格式求解该方程,并推导出了二维黏弹TTI介质完全匹配层(PML)吸收边界条件公式和相应的旋转交错网格任意偶数阶精度有限差分格式,实现了该类介质的地震波场数值模拟.数值模拟结果表明:该方法模拟精度高,边界吸收效果好,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征.  相似文献   

13.
TTI介质有限差分逆时偏移的稳定性探讨   总被引:6,自引:4,他引:2       下载免费PDF全文
李博  李敏  刘红伟  刘洪 《地球物理学报》2012,55(4):1366-1375
在沉积学中,可假设在相同时期的沉积层具有相近的物理性质和演化过程.因此,沿层传播的地震波和垂直于地层传播的地震波具有各向异性的特点.在纵波资料的处理中,考虑各向异性对逆时偏移的影响,通常假设介质的横波速度为零,这样可以得到纵波在TTI介质中的传播方程,但是该方程在实际计算中仍存在数值稳定性问题.本文加入横波分量可有效解决数值稳定性问题,并选取适当的横波速度减小对纵波成像的影响,实现地震波在TTI介质中的逆时偏移.实际测算表明,P-SV波的方程中包含横波分量,若假设SV的速度为零,则会导致方程的差分格式不稳定;若加入SV波,选择合适的SV波速度可以使SV波的全区各向异性和反射系数达到极小,并可有效的抑制SV波对纵波勘探的影响.本文的方法是一种稳定的TTI介质中的逆时偏移方法.  相似文献   

14.
TTI(Tilted Transversely Isotropic)各向异性是对地下岩石中广泛存在的规则发育的裂缝和层理的一种有效的弹性近似,基于TTI介质的地震波数值模拟技术是分析地震波在复杂各向异性介质中的传播机理的有效工具.同时,高精度的数值模拟算法也能为后续的逆时偏移技术提供重要的技术支撑.由于TTI介质中地震波方程的弹性参数众多且变化复杂,常规有限差分技术在解决TTI介质正演模拟问题时往往会产生严重的数值频散现象,降低了数值模拟精度.通量校正传输(FluxCorrected Transport,FCT)技术能够有效地压制由空间离散产生的数值频散.本文将FCT技术用于TTI介质中弹性波方程的交错网格高阶精度差分正演,在数值模拟过程中通过对波场进行漫射和反漫射校正实现了空间网格频散的压制.模型模拟结果表明,与常规有限差分算法相比,本文算法能够有效的压制大网格条件下的数值频散,提高模拟精度.  相似文献   

15.
三维TTI介质相速度和群速度   总被引:2,自引:0,他引:2       下载免费PDF全文
相速度和群速度是研究地震波传播规律和描述介质特性的重要参数,是弹性波传播理论中的核心内容,在理论研究和实际应用中有重要作用.本文根据VTI介质的刚度矩阵,利用Bond变换建立了TTI介质刚度矩阵.再利用TTI介质刚度矩阵,结合弹性动力学的本构方程、牛顿运动微分方程和几何方程,得到了三维TTI介质弹性波波动方程和Christoffel方程.通过本征值方法求解Christoffel方程,推导了三维TTI介质弹性波相速度的解析表达式.利用Berryman和Crampin推导各向异性介质群速度公式,根据三维TTI介质的相速度解析式推导了三维TTI介质群速度解析表达式.数值试例表明,随着各向异性介质参数改变,TI介质弹性波相速度变化较为平缓,群速度变化较为剧烈,qP波和SH波速度变化较为平缓,qSV波速度变化较为剧烈.  相似文献   

16.
以具有垂直对称轴的横向各向同性(Vertical Transversely Isotropic,简称VTI)介质中纯准P波方程为基础,通过投影变换推导出倾斜横向各向同性(Titled Transversely Isotropic,简称TTI)介质中的纯准P波方程,该方程从根本上消除伪横波的影响。方程形式简单并且在δε的情况下仍能稳定求解。文章还给出TTI介质中纯准P波方程数值求解的高阶有限差分格式并推导对应的稳定性条件,依据完全匹配层(Perfectly Matched Layer,简称PML)的分裂思路给出TTI介质中纯准P波方程的PML吸收边界条件。理论分析和数值计算结果均表明:本文推导的纯准P波方程能更好地描述TTI介质中的准P波,文章给出的正演算法也能够准确模拟TTI介质中准P波的传播过程。  相似文献   

17.
TI介质局部角度域射线追踪与叠前深度偏移成像   总被引:1,自引:1,他引:0       下载免费PDF全文
研究与实践表明,对于长偏移距、宽方位地震数据,忽略各向异性会明显降低成像质量,影响储层预测与描述的精度.针对典型的横向各向同性(TI)介质,本文面向深度域构造成像与偏移速度分析的需要,研究基于射线理论的局部角度域叠前深度偏移成像方法.它除了像传统Kirchhoff叠前深度偏移那样输出成像剖面和炮检距域的共成像点道集,还遵循地震波在成像点处的局部方向特征、基于扩展的脉冲响应叠加原理获得入射角度域和照明角度域的成像结果.为了方便快捷地实现TI介质射线走时与局部角度信息的计算,文中讨论和对比了两种改进的射线追踪方法:一种采用从经典各向异性介质射线方程演变而来的由相速度表征的简便形式;另一种采用由对称轴垂直的TI(即VTI)介质声学近似qP波波动方程推导出来的射线方程.文中通过坐标旋转将其扩展到了对称轴倾斜的TI(即TTI)介质.国际上通用的理论模型合成数据偏移试验表明,本文方法既适用于复杂构造成像,又可为TI介质深度域偏移速度分析与模型建立提供高效的偏移引擎.  相似文献   

18.
目前在地震勘探频带范围内通常假设品质因子Q与频率无关,且呈衰减各向同性.事实上,相比较速度各向异性,介质的衰减各向异性同样不可忽视.本文将衰减各向异性和速度各向异性二者与常Q模型相结合,建立了黏弹性衰减VTI介质模型,并基于分数阶时间导数理论,给出了对应的本构关系和波动方程.利用均匀平面波分析和Poynting定理,推导出准压缩波qP、准剪切波qSV和纯剪切波SH的复速度、相速度、能量速度以及品质因子的解析表达式.对模型的正确性进行了数值验证,并分析了qP,qSV和SH波在介质中的传播特性.数值试验结果表明:本模型能够实现理想的恒定Q行为,表现了品质因子和速度的各向异性特征,显示出黏弹性增强将导致能量速度和相速度的频散曲线变化剧烈;速度和衰减各向异性参数与传播角度之间的耦合效应对qP,qSV和SH波的速度和能量影响明显;qP,qSV和SH波的频散曲线和波前面随着衰减各向异性强度的改变发生显著变化,其中耦合在一起的qP和qSV波变化趋势相同,而SH波与它们呈现相反的变化规律.本研究为从常Q模型角度分析地震波在衰减各向异性黏弹性介质中的传播特征奠定了理论基础.  相似文献   

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

20.
VTI介质纯P波混合法正演模拟及稳定性分析   总被引:3,自引:3,他引:0       下载免费PDF全文
各向异性介质纯P波方程完全不受横波的干扰,在一定程度上可以减缓由于介质各向异性引起的数值不稳定,本文推导了具有垂直对称轴的横向各向同性(VTI)介质纯P波一阶速度-应力方程.由于纯P波方程存在一个分数形式的伪微分算子,无法直接采用有限差分法求解.针对该问题,本文采用伪谱法和高阶有限差分法联合求解波动方程,重点分析了混合法求解纯P波一阶速度-应力方程的稳定性问题,并给出了混合法求解纯P波方程的稳定性条件.数值模拟结果表明纯P波方程伪谱法和高阶有限差分混合法能够进行复杂介质的正演模拟,在强变速度、变密度的地球介质中仍然具有较好的稳定性.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号