波动方程共炮检距道集叠前深度偏移

目前所用的几种波场外推算法（分步傅里叶、傅里叶有限差分、广义屏和频率--空间域有限差分）可以在共炮道集上获得较好的成像效果。在上述算法中,共炮道集为了兼顾各种波传播角,每炮偏移都要考虑相当多的边道,并且炮点和检波点要分别向下外推,因而计算量很大。针对上述问题,本文提出了基于双平方根算子的波动方程共炮检距道集的深度偏移方法,它适用于二维和三维叠前深度成像,其步骤为：①在高频假设条件下,把共炮检距道集波场延拓公式中的积分运算进行稳相近似,得到波场延拓的相移公式;②把速度场分裂为层内常速背景和变速扰动,可求得整个均匀层波场深度延拓的偏移时移量、各层的偏移时移量及振幅校正系数,进而可得到最终波场延拓值。脉冲响应测试和理论模型试算表明,该方法具有较高的计算效率和成像精度,可适用于复杂地质体成像。     

双平方根算子波场延拓道集速度分析

在中心点一半偏移距域，利用波动方程的叠前双平方根算子．实现炮点及检波点波场向下延拓。通过速度扫描，在波场延拓后的深度域共中心点道集和经分选后的共成像点道集上，分析速度变化对成像点道集的影响。利用分偏移距成像，修正速度场，以达到提高成像精度的目的。     

基于双平方根方程的保幅地震偏移

 本文提出了一种基于双平方根方程的保幅叠前深度偏移方法，首先从保幅的单平方根算子方程出发，将保幅的上、下行波场进行组合近似表示为全声波波场，然后推导出由双平方根算子定义的单程波方程；根据摄动理论把速度场分裂为层内常速背景和变速扰动，求得整个均匀层波场深度延拓的偏移时移量及振幅校正系数、各层的偏移时移量及振幅校正系数，得到最终的波场延拓值，实现了共炮域公式到CMP域公式的转换。理论模型试算和实际资料处理表明，该方法不但具有较高的计算效率和成像精度，而且可以输出正确反映地下反射系数的振幅信息，有利于AVO/AVA地震属性分析，提供精确的岩性识别与储层预测地震数据。     

单步波场延拓法单炮时间偏移

偏移算法可以分为波场和射线两种类型,其中叠前时间偏移主要使用射线法,而射线法难以达到最佳振幅响应特性,因此,波场法一直都是研究的目标。现有的波场法叠前时间偏移主要是共炮检距和共平面波两种F-K域实现方式。一般而言,单炮总是最适合的选择,它具有方程简单、操作方便和容易实现的特点。根据这个认识,本文通过相移加内插的方法实现单炮时间偏移,采用均方根速度进行单步波场延拓,从而取得良好的成像质量和振幅特性。     

波动方程共方位角三维叠前偏移方法

本文针对双平方根(DSR)单程波动方程全三维偏移方法在地震波成像实际应用中面临的计算量大,对覆盖次数很低的窄方位三维地震资料会产生很强的相干噪声,影响成像效果等难题,结合Crossline共炮检距偏移理论,首先推导出了基于双域传播算子的共方位角叠前深度偏移公式,然后通过坐标变换提出了沿双程垂直走时方向进行递归波场延拓的共方位角叠前偏移新方法。通过对SEG/EAGE盐丘模型合成数据偏移试验,证实了共方位角叠前偏移深度域和时间域成像方法的有效性。将该法和常用的Kirchhoff叠前偏移法同时用于实际地震资料偏移处理对比,展示出该法在成像精度与分辨率方面具有明显的优势。     

基于双平方根方程的双域保幅地震偏移方法

从全声波方程出发,基于严格的解耦理论进行单程波保幅分解,得到保幅的单程波动方程;然后按照叠前深度偏移中的波场双重向下延拓理论实现共炮域公式到CMP域公式的转换,推导出由双平方根方程定义的保幅单程波动方程;基于反问题求解中常用的摄动理论,把速度场分裂为层内常速背景和变速扰动,分别在频率波数域和频率空间域求得波场深度延拓的偏移时移量及振幅校正系数,得到最终的双域保幅波场延拓算子。理论模型试算和实际资料处理表明,该方法不但具有较高成像精度,而且可以为后续的地震属性分析提供更真实的振幅信息。     

保幅型裂步傅里叶叠前深度偏移方法探讨

利用地震资料进行岩性分析,要求叠前深度偏移成像结果不但包含走时信息,还要包含振幅信息。但是应用传统的单程波方程偏移算子只能得到地震波的运动学特征,无法得出其动力学特征。裂步傅里叶叠前深度偏移方法是基于单程波方程的一种快速稳定的叠前深度偏移方法,它遵循速度场分裂的思想,分别通过频率-波数域的相移处理和频率2空间域的时移处理来实现。但是传统单程波方程在保幅性方面的缺陷制约了该方法在岩性成像中的应用。文中引入了裂步傅里叶保幅偏移算法,分别通过频率2波数域的相移保幅校正算子和频率-空间域的保幅时移校正算子实现了保幅偏移。文中针对Marmousi 模型的炮点照明分析结果及其叠前深度偏移成像结果的对比表明了该方法的正确性和有效性。     

炮域波动方程叠前深度偏移的偏移距域和角度域共成像点道集计算

炮域波动方程叠前深度偏移是一种实用化的成像算法,特别在稀疏炮集宽方位角观测系统采集的数据成像中非常有效。本文将偏移距定义为叠前偏移中下行震源波场和上行接收点波场之间的位移,在偏移过程中提取偏移距域共成像点道集,利用最早应用于炮—检偏移的径向记录道映射方法,将偏移距域道集转换到角度域。计算过程包括多偏移距成像和炮—检点坐标系向中心点坐标系的转换。该方法能够在炮集偏移过程中生成随角度变化的反射波成像结果,但在稀疏炮集观测系统下,共成像道集的质量仍会受到炮点假频的影响。     

平面波偏移、分角度成像与AVA道集生成

基于波场延拓的叠前深度偏移是实现复杂构造地质体成像的可靠方法，但存在着计算量大、对观测系统适应性差等缺点。平面波偏移是利用波动方程实现精确叠前成像的另一类方法，其基本原理是：通过地表延迟放炮的方式生成平面波震源，利用下行波方程进行波场正向延拓得到下行波场；对地表采集的炮集记录，以组合延迟放炮的方式叠加，得到地表平面波记录，利用上行波方程进行波场反向延拓得到上行波场；二者互相关求和，实现平面波地下波场成像。分析表明：平面波成像技术的精度与单平方根算子的共炮点道集偏移和双平方根算子的共中心点道集偏移相当，但计算速度要快得多，且易于并行计算。二维Marmousi模型数值计算表明，射线参数的范围和间隔是影响平面波成像质量的主要因素；不同角度入射的平面波对最终成像结果的贡献是不同的，据此可以有针对性地选择射线参数进行平面波成像。     

高精度屏算子地震偏移成像方法研究

从波场延拓的非稳态相移公式出发,基于反问题求解中常用的摄动理论,利用单平方根算子的渐进展开,推导出了波动方程广义屏叠前深度偏移算子方程的高阶形式;针对散射波场计算项对于横向变速介质的不稳定性问题,通过数学近似提出了有效提高稳定性的策略并应用到波场递归外推过程中,从而得到一种稳定的高精度屏算子地震偏移成像方法。数值试验和实际资料处理表明,该方法具有更高的精度,对宽角度信息成像更好,是对传统广义屏算子的有效补充和发展。     

相移插值法偏移中一种新的变延拓步长方法

在叙述相移插值法波动方程偏移原理、对该偏移过程中的波场外推公式及最终成像方程进行适当的改进变形后,重新标定了向下外推算子,使得成像过程可借助于 FFT 来完成,并给出了一种不同于时间域内插的新的变延拓步长的方法。本方法与传统的相移插值法偏移中时间域内插值法变延拓步长相比,向下外推步长可远远大干地震波视波长的四分之一,从而满意地解决了相位插值偏移算法效率低的问题。文中给出的理论模型计算与实际煤田及陡倾角地震数据偏移处理的结果表明:本方法是一种快速经济的,适合于陡倾角、复杂地质构造归位的 f-k 域偏移成像技术。     

快速面炮记录叠前深度偏移

本文提出了一种允许部分相位值依赖频率，部分相位值不依赖频率的混合相位编码技术。将这一新的相位编码技术应用到目标照明的面炮记录叠前深度偏移上，即可使面炮记录偏移效率成倍提高，又可降低因相位编码产生的人为误差，从而得到了一种新的快速面炮记录叠前深度偏移方法。该方法基于递归波场延拓，具有较高的成像质量。由于既使用了比单炮记录偏移高效的面炮技术，又使用了相位编码技术，用相当于一个平面波震源的面炮记录叠前深度偏移的计算量，近似地得到多个面炮记录叠前深度偏移叠加结果，因此具有高效的特点。Marmousi模型数据和实际墨西歌湾海底电缆纵波数据试算结果表明了本方法的正确性、有效性和实用性。     

一种高精度的偏移速度分析方法——广义屏偏移加Born反演

在波动方程叠前深度偏移中，速度是最关键的成像参数。首先分别对各种求取速度的方法进行了回顾和分析，然后按照判断速度正确性的不同标准将求取速度的方法分为3类：第1类是以正演结果与实际数据的拟合最优作为判断标准；第2类是以成像质量最优作为判断标准；第3类是解析类型的方法。在此基础上提出，基于逆散射理论反演慢度扰动的方法和双域偏移成像方法在理论基础上有相当的一致性，所以在利用波形信息的高精度速度分析中可以将二者有机地结合起来。进而提出了一种将广义屏偏移与Born反演结合在一起实现速度分析的新思路。该方法保证了速度分析和偏移在理论基础上的一致性，能够充分发挥波动方程偏移在振幅保真方面的优势，是一种稳定的、高精度的速度分析方法。     

խ��λ����˫ƽ��������ƫ��Ӧ���о�

双平方根（DSR）方程为波动方程偏移提供了一种新的理论框架。基于“沉降观测”概念的DSR方程叠前深度偏移已经成为了一类重要的地震波成像方法。与单程波方程偏移方法相比，其计算效率要高很多，便于输出角度域的共成像道集，有利于开展波动方程偏移速度分析。DSR方程全偏移方法对计算机条件要求太高，且对地震勘探中大量存在的窄方位三维地震数据不太适用。文章结合DSR方程叠前深度偏移的波场传播算子与成像条件，讨论了它在窄方位三维地震数据成像中的应用问题，对比分析了几种稳相近似DSR方程偏移方法的特点与适用条件。数值试例表明，基于DSR方程的窄（共）方位角叠前深度偏移方法对盐丘（或礁体）、古潜山等复杂构造具有很强的成像能力，值得在石油天然气勘探中推广应用。     

共炮检距剖面相位移时间偏移

本文提出了一种新的二维频率一波数域共炮检距剖面叠前时间偏移算法。本文方法在常速条件下将零炮检距剖面的叠后相位移公式推广到非零共炮检距剖面，该公式是用数值拟合方法得到的，运用平均速度或叠加速度可使公式的使用范围适用纵向速度变化，理论上本文的公式只在常速下是精确的，但合成记录和野外资料试算表明，即使速度变化较大，该方法也能有效成像，而且更重要的是还能减少计算量。