叠前逆时深度偏移中的激发时间成像条件

与其他偏移方法相比,逆时偏移基于精确的波动方程而不是对其近似,用时间外推来代替深度外推.因此,它具有良好的精度,不受地下构造倾角和介质横向速度变化的限制.激发时间成像条件的求取是叠前逆时偏移的难点之一,本文采用求解程函方程的方法得到地下各点的初至波走时,以此作为叠前逆时偏移的成像条件.基于任意矩形网格和局部平面波前近似的有限差分初至波走时计算方法精度较高并适用于强纵横向变速的复杂介质.试算结果表明,在复杂介质模型中利用叠前逆时深度偏移收到了很好的成像效果.     

VTI介质中准P波方程叠前逆时深度偏移

根据具有垂直对称轴的横向各向同性（VTI）介质中的一阶准P波方程,导出了该方程在交错网格中逆时延拓的高阶有限差分格式,给出了其稳定性条件,采用完全匹配层吸收边界条件解决边界反射问题,分别应用下行波最大能量法和归一化互相关成像条件, 实现了VTI介质中准P波方程的叠前逆时深度偏移．各向异性Marmousi模型的试算结果表明,VTI介质准P波方程叠前逆时深度偏移算法不受地下构造倾角和介质横向速度变化的限制,对复杂模型具有良好的成像能力;应用归一化互相关成像条件能得到更好的成像效果．对比该模型的各向异性和各向同性逆时偏移剖面表明,在各向异性地区采集的纵波数据用各向异性偏移算法理论上能得到更好的成像结果．      

起伏地表条件下分块法插值叠前深度偏移

为实现复杂起伏地表条件下基于单程波动方程的高效叠前深度偏移,本文采用"逐步一累加"思想,应用傅里叶有限差分法(FFD)算法实现复杂起伏地表务件下的波动方程叠前深度偏移.为提高效率,本文提出在整个成像实现过程中,将整个成像区域分两块区别对待,在地表高程最低点所在的水平面到基准面的区域内,应用常规精细延拓步长进行波场延拓及...     

分层坐标变换法起伏自由地表弹性波叠前逆时偏移

传统有限差分方法在处理起伏地表时存在一些困难, 而坐标变换法可将起伏地表映射为水平地表以克服此缺点.但同时, 地下构造被变换得更加复杂, 导致了波传播和成像的不准确.本文提出了一种分层的坐标变换方法, 并应用到了弹性波逆时偏移中, 此方法既可以克服起伏地表的影响, 又可以不破坏地下构造.波场正向延拓、逆时延拓和分离是在辅助坐标系下完成的, 而成像是在笛卡尔坐标系下完成的.通过对简单起伏模型和中原起伏模型的试算证明了本文提出方法的准确性.同时, 对两种极端起伏地层高程不准确的情况进行测试可以看出:分层坐标变换逆时偏移方法的成像效果远好于传统坐标变换方法.     

起伏地表条件下基于复Pade逼近的叠前深度偏移

叠前深度偏移是解决复杂地表和复杂构造地震成像的有效技术,而波场“直接下延”法实现了复杂地表条件下的地震成像.基于上述成果,结合高精度的波场延拓算子,本文提出了一种新的叠前深度偏移方法,这种方法是在波场延拓时,对声波方程中的平方根项进行复Pade逼近,通过推导得到基于复Pade逼近的傅里叶有限差分算子,结合波场“直接下延”法,实现了起伏地表条件下的叠前深度偏移,该算法减少了偏移噪音,从而得到准确、稳定的偏移成像结果.通过理论模型试算和实际资料试处理,验证了该方法的有效性.     

起伏地表叠前逆时偏移的实现及应用

复杂地表复杂构造是近年来陆上地震勘探的重难点问题.其中,复杂构造成像问题可用叠前深度偏移技术如逆时偏移较好地解决,而地表起伏大、横向变速剧烈引起道间时差大问题则需要采用起伏地表偏移的思路.结合常规地震数据处理,本文提出一种基于起伏面的逆时偏移成像技术流程,该流程利用固定面与平滑起伏面高程差和替换速度将静校正后的叠前炮集反静校正到起伏面上,然后进行起伏面逆时偏移成像计算.将流程应用到模型数据和实际资料处理,结果均表明该流程在复杂地表复杂构造地震成像方面具有很好的效果.     

任意广角波动方程叠前逆时深度偏移

从任意广角声波方程出发,在深度与时间方向采用二阶差分代替微分,在水平方向上采用高阶差分代替微分,得到任意广角声波方程逆时延拓的有限差分格式. 通过差分求解程函方程得到叠前逆时深度偏移的成像条件,在此基础上实现了二维各向同性介质中任意广角方程的叠前逆时深度偏移. 对理论模型的偏移处理表明,该方法能够有效地消除波场下行传播过程中的层间反射,减小偏移噪声. 同时,该方法还具有保真性高、 适用于大倾角地层的偏移等优点.     

地表旋转控制照明叠前深度偏移

基于炮集数据的常规叠前深度偏移波场外推的计算量非常大,控制照明叠前深度偏移是基于快速的组合炮技术,通过在目的层上控制震源波场,高效实现对地质目标的高质量成像. 组合炮所需的合成算子,由在目的层定义的震源波场反向外推到地表建立. 本文提出一种更快速的地表旋转控制照明方法,直接在地表对震源波场进行面向目标的照明控制,实现目标区域高质量成像. 该方法只需对某一核合成算子进行简单的几何旋转建立合成算子函数,避免了复杂的波场外推. 这种新的控制照明方法与快速的一般组合炮偏移效率相同,但对复杂目标区域的成像能力显著增强. 应用本文方法对Marmousi模型进行计算,得到深层复杂构造非常理想的成像结果.     

非零偏VSP弹性波叠前逆时深度偏移技术探讨

非零偏VSP地震资料是一种多分量资料,处理非零偏VSP资料,弹性波叠前逆时深度偏移技术无疑是最适合的处理技术.本文从二维各向同性介质的弹性波波动方程出发,研究了对非零偏VSP资料进行叠前逆时深度偏移的偏移算法,讨论了逆时传播过程中的边值问题和数值频散问题及其相应的解决方案;采用求解程函方程计算得到地下各点的地震波初至时间作为成像时间,实现了非零偏VSP资料的叠前逆时深度偏移.最后进行了模型试算和非零偏VSP地震资料的试处理,结果表明该方法不受地层倾角限制,较适用于高陡构造地区或介质横向速度变化较大地区的非零偏VSP地震资料处理.     

三维各向异性介质中的波动方程叠前深度偏移方法

基于三维VTI各向异性介质的频散关系,构建波数项和空间项分离的单程波算子表达式,以优化算法,确定算子的待定系数,实现广角逼近三维VTI介质的广义相移算子,发展了可灵活处理强或弱各向异性介质的波动方程叠前深度偏移方法.文中同时也针对其工业应用建议了三维VTI各向异性介质中可提高计算效率的频率相关变步长波场深度延拓算法及稀疏采样情况下可实现陡倾角构造正确成像的反假频波场延拓算法.SEG二维Hess VTI各向异性介质理论模型数据及野外数据集观测系统下的三维脉冲响应计算表明,本文提出了一种具有工业应用潜力的三维各向异性波动方程叠前深度方法.     

Prestack reverse-time depth migration of arbitrarily wide-angle wave equations

Based on arbitrarily wide-angle wave equations,a reverse-time propagation scheme is developed by substituting the partial derivatives of depth and time with central differences. The partial derivative of horizontal direction is replaced with high order difference. The imaging condition is computed by solving the eikonal equations. On the basis of above techniques,a prestack reverse-time depth migration algorithm is developed. The processing exam-ples of synthetic data show that the method can remove unwanted internal reflections and decrease the migration noise. The method also has the advantage of fidelity and is applicable of dip angle reflector imaging.     

谱元法叠前逆时偏移研究

谱元法是基于弹性力学方程的弱形式,运用有限元的思想在单元上进行谱展开,理论上该方法具备有限元适应任意复杂介质模型的韧性和伪谱法的精度.这里借助谱元法这一思想,在正演算法的基础上,进行偏移算法的改编,并在各向异性介质的二个模型中进行叠前逆时偏移.结果表明:这二个模型上进行的偏移结果与模型吻合较好.     

非稳态相移法叠前深度偏移

介绍一种能够适应介质速度横向变化的非稳态相移算子及其叠前深度偏移方法．为了克服常规相移偏移算法中要求速度横向不变的缺点,出现了基于非稳态滤波器理论的非稳态相移算子,即PSPI算子、NSPS算子和SNPS算子,其中SNPS算子是将前二者结合起来的一种对称的非稳态相移算子,它比前二者具有更高的精度和稳定性．为了提高运算速度,基于非稳态相移算子的叠前深度偏移算法采取了分片均匀近似的策略,Marmousi模型的叠前深度偏移结果证明了该算法的可行性和有效性。     

起伏地表下的直接叠前时间偏移

提出了一种新的叠前时间偏移方法和流程,可不必应用野外静校正,直接对起伏地表采集的地震数据进行叠前时间偏移.本文采用输入道成像方式,通过基于稳相点原理给出单道数据的走时和振幅计算方法,发展了一个表驱动的叠前时间偏移算法.偏移方法可依据同相轴是否被拉平确定叠加速度和修正近地表速度模型,也可依据拟成像的构造倾角,自适应地确定偏移孔径;后者既减少了偏移计算量,也压制了偏移噪声.文中用二维起伏地表的断陷盆地模型的理论数据验证了所发展方法的成像效果.     

结合基准面重建的叠前时间偏移方法

提出了一种结合虚拟界面、瑞利积分和相移法的混合的基准面重建方法.通过与叠前时间偏移方法结合,形成了针对起伏地表采集数据的叠前时间偏移方法和新流程.该方法能正确考虑波在近地表传播的实际路径,克服了高速层出露时静校正方法的误差;它也能自己确定虚拟层速度,避免了现行基于波场延拓的基准面重建方法需要准确近地表速度的困难.文中分别用近地表存在明显低速层和近地表有高速层出露这两类模型的理论数据,验证了所发展方法和流程的有效性.     

双复杂条件下的波动方程叠前深度偏移

复杂地表条件下的地震勘探越来越被人们所关注.双复杂条件下的叠前深度偏移方法是解决复杂地表条件和复杂地质构造成像的有效手段."波场上延"法能实现由非水平观测界面开始的偏移过程,解决复杂地表对地下构造成像的影响.复杂理论模型的试算以及实际资料处理表明,"波场上延"方法较好地克服了起伏地形对地下构造成像的影响,取得了令人满意的效果,实现了波动方程基准面校正和深度成像的有机结合.     

基于非结构化网格的最小二乘逆时偏移

相对于传统的逆时偏移,最小二乘逆时偏移具有更高的成像质量,这种改善是通过迭代反演来获得的,另外其精度与效率依赖于求解波动方程算法的精度与效率.本文给出了基于非结构化网格的最小二乘逆时偏移的方法,该方法能够充分结合最小二乘算法与非结构化网格精细刻画地下界面以及随速度自适应剖分的优点;并采用带补偿的拉普拉斯滤波算法,来消除梯度计算中的低频噪声,从而加速目标函数的收敛速度.通过简单倾斜模型以及复杂Marmousi模型测试,显示了该方法的有效性和潜力.

     

基于虚拟偏移距方法的各向异性转换波保幅叠前时间偏移

In this paper, we use the method of pseudo-offset migration （POM） to complete converted wave pre-stack time migration with amplitude-preservation in an anisotropic medium. The method maps the original traces into common conversion scatter point （CCSP） gathers directly by POM, which simplifies the conventional processing procedure for converted waves. The POM gather fold and SNR are high, which is favorable for velocity analysis and especially suitable for seismic data with low SNR. We used equivalent anisotropic theory to compute anisotropic parameters. Based on the scattering wave traveltime equation in a VTI medium, the POM pseudo-offset migration in anisotropic media was deduced. By amplitude-preserving POM gather mapping, velocity analysis, stack processing, and so on, the anisotropic migration results were acquired. The forward modeling computation and actual data processing demonstrate the validity of converted wave pre-stack time migration with amplitude-preservation using the anisotropic POM method.