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

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

基于成像射线的偏移剖面和速度时深转换评述

深度域地震波速度模型是地震勘探中最重要的参数之一,是获得精确深度域偏移剖面的基础.本文基于成像射线的概念,阐述了时间偏移域和深度域之间的关系,明确了时间偏移后的像点是通过成像射线与深度域的散射点联系起来的.在有横向变速的情况下,时间偏移剖面需要经过时深转换才能真实反映深度域层位,同样,从时间偏移剖面上提取的速度也需要经过时深转换才能得到深度域的速度.与前者不同的是,在速度的时深转换中,不仅速度的位置需要改变,且其数值也需要进行校正.该校正因子有几种等价的描述形式：度规量、速度扩散因子和几何扩散因子.     

基于Walkaway VSP的Thomsen各向异性参数估计与应用(英文)

Thomsen各向异性参数的求取对于正确的时深转换和深度域偏移成像处理至关重要。相比其它各向异性参数估算方法,从VSP资料中更容易获取准确的各向异性参数用于地面地震偏移成像。本文分析研究了利用Walkaway VSP资料估算VTI介质Thomsen各向异性参数的方法,该方法基于VTI介质近偏移距动校正公式利用Walkaway VSP近偏移距初至信息求取各向异性参数δ;基于各向异性介质纵波速度Thomsen近似公式采用射线追踪时差扫描方法求取各向异性参数ε。数值模型正演表明利用该方法估算的各向异性参数误差较小。利用塔里木盆地8个方位的Walkaway VSP实际资料求取了该区深度域Thomsen各向异性参数ε和δ值,同时结合地面三维地震资料建立了较为准确的各向异性深度一速度模型用于叠前深度偏移成像,进一步提高了碳酸盐岩储层的成像精度,减小了目标地质体的深度误差。     

任意复杂介质中主能量法地震波走时计算

积分法叠前深度偏移及层析成像的核心是复杂介质情况下的地震波走时计算. 复杂构造的高精度地震成像需要有稳健的走时计算方法。本文把 Nichols提出的用地震波主能量计算走时的方法由二维推广到三维,并推导出三维波动方程Helmholtz形式在球坐标系下用因式分解法求解的差分表达式.三维SEG/EAGE盐丘模型的理论走时计算和积分法叠前深度偏移的实践都验证了本文方法的正确性.     

基于走时的保幅偏移方法

振幅随偏移距变化是描述储层特征的重要方法之一,保幅偏移方法就是使偏移剖面能够反映出振幅随偏移距的变化.本论文中的保幅偏移是以走时为基础,主要的方法是采用走时的双曲线展开法,通过走时的二阶空间导数来确定波前曲率.该方法通过建立在大网格上的走时表来确定插值系数,将大网格插值成为较为精细的网格,这样就节省了数据的存储空间.对于相同的网格密度,通过插值来计算走时表比采用程函方程有限差分法直接计算走时要节省5至6倍的时间.走时的插值系数还可以用来计算几何扩散因子、权函数,不仅提高了成像质量,还大大节省了计算时间.     

关于共反射面元叠加方法在实际应用中的一些思考

共反射面元（Common Reflection Surface=CRS）叠加是一种特殊的零偏移距成像方法,实践中它具有独立于宏观速度模型和完全数据驱动实现的鲜明特色,CRS叠加理论认为在得到高质量的零偏移距剖面的同时,还可以得到三个有用的波场属性参数剖面反演宏观速度模型,CRS叠加剖面之后的叠后深度偏移质量将超过叠前深度偏移.虽然CRS叠加倡导的成像方式和承诺的上述理想境界带来了全新的启示,但是实践中这些特色同样带来了令人困扰的问题,为此我们提出了倾角分解CRS叠加方法解决这些问题.本文即是作者通过上述实践之后对CRS叠加方法形成的一些思考和总结.     

双平方根单程波动方程叠前τ偏移方法

本文将常规双平方根（DSR）单程波动方程从深度域变换到双程垂直走时（τ）域，由此推导出可从数学上实现“沉降观测”的单程波DSR传播算子. 其递归波场延拓算法包含波数域针对常速背景的相移处理和空间域针对横向速度扰动的相位校正，可以应对上覆地层速度横向变化对构造成像的影响. 结合零炮检距、零时间成像条件，提出了在τ域进行波场延拓与成像的DSR方程叠前偏移新方法. 为了克服其全三维偏移算法在实际应用中可能面临的困难，本文采用稳相近似，在crossline常炮检距偏移理论基础上推导了实用的共方位角叠前τ偏移方法. 数值试验表明，DSR方程叠前τ偏移在强横向非均匀介质中的成像精度与分辨率优于传统的时间域成像技术.     

双平方根单程波动方程叠前τ偏移方法

本文将常规双平方根（DSR）单程波动方程从深度域变换到双程垂直走时（τ）域，由此推导出可从数学上实现“沉降观测”的单程波DSR传播算子. 其递归波场延拓算法包含波数域针对常速背景的相移处理和空间域针对横向速度扰动的相位校正，可以应对上覆地层速度横向变化对构造成像的影响. 结合零炮检距、零时间成像条件，提出了在τ域进行波场延拓与成像的DSR方程叠前偏移新方法. 为了克服其全三维偏移算法在实际应用中可能面临的困难，本文采用稳相近似，在crossline常炮检距偏移理论基础上推导了实用的共方位角叠前τ偏移方法. 数值试验表明，DSR方程叠前τ偏移在强横向非均匀介质中的成像精度与分辨率优于传统的时间域成像技术.     

复杂介质下保真振幅Kirchhoff深度偏移

本文详细讨论地下复杂介质地震成像问题.当速度模型含有强横向变化时,射线场会折叠形成多次走时,在这种情况下,文中证明在保真振幅叠前深度偏移时,应当考虑所有的到时才能得到定量的像.通过实例证明,在多次走时情况下,常规的共炮集、共检波器集都存在强的假像;由此提出共角度偏移的概念.最后,在多次走时下,讨论多次偏移算子(走时、相移、振幅等)的高效重构,给出了三维叠前深度偏移的快速算法.     

基于Chebyshev多项式的弯曲射线Kirchhoff叠前时间偏移

弯曲射线Kirchhoff叠前时间偏移被认为是一种综合了叠前时间偏移效率和叠前深度偏移精度的方法.本文以保精度地减少高阶Kirchhoff叠前时间偏移走时计算量为目标,在分析了Chebyshev正交多项式性质的基础上,建立了Chebyshev多项式约简系数表,进而用模拟退火法对转换系数进行分段优化,从而实现了在大炮检距...     

利用τ-p波场反演法确定东沙群岛的地壳速度模型

本文利用τ-p波场反演法直接对OBS原始数据进行转换,得到一维速度深度函数.该方法包含倾斜迭加和向下延拓两个线性转换:倾斜迭加将原始数据由时间-偏移距域(t-x)转换为时间截距-慢度域(τ-p);向下延拓将时间截距-慢度域(τ-p)波场转换为深度-慢度域(z-p)波场.然后根据五个OBS台站得到的一维速度深度函数,合成东沙群岛地壳的二维速度结构模型.为了检验该速度结构模型,我们用拾取走时和利用模型计算的理论走时进行拟合,得到了满意的结果,从而为下一步进行层析成像提供了可靠的初始速度结构模型.     

基于射线追踪的成像域走时层析反演及其应用

层析成像是现今地震勘探开发处理中的常用手段,针对传统走时层析反演中角度覆盖问题,利用正则化约束加以解决,利用角道集的拾取拟合来判断实验模型的准确性,同时,不再利用原始炮记录拾取剩余时差,利用深度残差与走时残差的线性关系提高走时残差的拾取精度。以三层洼陷模型和实际资料进行试算,结果表明:本方法可以较好地反演更新偏移速度场,得到更精确的偏移成像结果,对于小尺度地质构造可以精细刻画。      

非椭圆各向异性介质高阶优化走时计算方法及应用

随着地震勘探的深入,对勘探精度的要求越来越高,目前,提高勘探精度的主要方法之一就是提高地震波走时计算的精度,地震走时精度直接影响如动校正、静校正和偏移成像质量.根据地质模型假设条件的不同,有不同的走时计算方法,常用的有各向同性介质假设和各向异性介质假设条件下的走时计算方法,各向异性假设具有比各向同性假设更好的适用性,具有更高的走时计算精度,能够满足大偏移距和高陡构造情况下成像的要求.本文根据Alkhalifah的VTI介质中走时计算方法,以Fomel的TI介质中群速度的计算为基础,利用各向异性介质中Snell定律,提出了在非椭圆各向异性介质中的高阶优化走时计算方法,通过理论和实际数据的应用,表明该方法克服了常规各向异性走时方法中不能够计算大偏深比情况下地震波走时的缺点,具有走时计算精度高、便于实际应用等特点,对叠加和偏移成像有广阔的应用.     

任意介质中的动态规划法地震波三维走时计算

任意介质中的地震波三维走时计算是复杂介质情况下Kirchhoff积分法三维叠前深度偏移及走时层析成像的核心.走时算法的效率及精度决定了成像方法的应用范围及效果,对复杂地质构造区域的地震波成像时需要有稳健的走时计算方法.本文把Schneider等提出的用动态规划法计算二维任意复杂介质中走时的方法推广到三维.此方法的核心是构造从源点到当前计算点的平均慢度,基于Fermat原理,用球面波近似导出走时计算所用的公式,并用动态规划法搜索到达当前计算点的初至走时.它适用于任意复杂的介质情况,对速度差异没有限制,计算过程中考虑到各个可能的方向到达当前计算点的初至波.首波及回转波的初至走时也能正确地计算出来.各种理论速度模型上的走时计算及胜利油田某探区的三维叠前深度偏移的成功实践验证了方法的正确性.     

各向异性弹性波有限元叠前逆时偏移

利用有限元法和最小走时射线追踪的界面点法，实现了各向异性弹性波的叠前逆时偏移.理论模型资料的偏移结果清晰准确，证实了该系统的可靠性.通过对各向异性模型资料的各向同性偏移处理，发现常规偏移剖面存在较大误差，地质体的垂向深度和横向位置与实际模型有偏离.     

2D共炮时间域高斯波束偏移

针对传统射线方法在奇异区成像精度不高,而2D频率域高斯波束叠前深度偏移需要计算成像点处每个频率的格林函数,影响计算效率的问题,本文通过使用复走时代替实走时,改变频率域下成像公式的积分顺序,给出了在时间域下进行高斯波束偏移的方法和计算公式.本文使用复杂数值模型验证了2D时间域高斯波束叠前偏移方法的正确性,并同传统射线偏移成像结果做了对比.对比结果表明时间域高斯波束偏移在成像精度上优于传统射线偏移.     

先进的共中点射线追踪法

本文介绍的射线追踪法是以对共反射元叠加方法而得出的做工为基础的。假设二维，平滑的横向变化介质中，偏移射线和走时都是根据垂直入射（中央）射线计算的。该方法使用射线的二阶非对称近似和计算偏移走时的中央波前附加的倾斜面近似。为了求解在共中点周围位置对称的源－接收器对的共中点构形的两点射线追踪问题，扰动中央射线以满足上述非对称分布。尽管对远偏移距计算精度受到限制，但是对反射深度不同一量级的距离来说，它仍是     

人工源宽角地震实验的谱元法模拟和逆时偏移成像方法研究

人工源宽角反射(折射)地震资料具有偏移距较大、信噪比较低等特点,通常用于地震波走时反演重建地壳速度结构。逆时偏移成像方法作为勘探地震学领域获取地下构造形态的有效手段之一,可以有效弥补走时反演方法的不足。本文针对大偏移距宽角反射(折射)地震实验,利用四边形网格谱元法进行波场模拟,结合了有限元法的灵活性和谱方法的指数收敛性,高效且高精度获取模型合成地震记录,后采用逆时偏移成像方法将合成地震记录偏移归位,获取地壳几何结构,验证了逆时偏移成像方法在宽角地震资料处理及结果解释中的适用性,为后期实际地震资料的偏移成像提供了理论依据和支持。     

利用τ\|p波场反演法确定东沙群岛的地壳速度模型

本文利用τ\|p波场反演法直接对OBS原始数据进行转换,得到一维速度深度函数.该方法包含倾斜迭加和向下延拓两个线性转换:倾斜迭加将原始数据由时间\|偏移距域(t-x)转换为时间截距\|慢度域(τ\|p);向下延拓将时间截距\|慢度域(τ\|p)波场转换为深度\|慢度域(z\|p)波场. 然后根据五个OBS台站得到的一维速度深度函数,合成东沙群岛地壳的二维速度结构模型.为了检验该速度结构模型,我们用拾取走时和利用模型计算的理论走时进行拟合,得到了满意的结果,从而为下一步进行层析成像提供了可靠的初始速度结构模型.     

横向均匀介质中地震走时的正演和反演

横向均匀介质可由若干水平分层足够好地逼近,本文引入零厚度层来描述成层介质,从而建立了回折波、反射波、首波和直达波统一的时距关系以及地震走时随模型参数变化的公式。本文详细讨论了水平分层介质里地震走时曲线的特征,合成地震走时的试验表明了由地震走时,用阻尼最小二乘法可以成功地反演速度模型参数,所反演的模型含有速度梯度层,也可以含有低速层。