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

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

基于频率-慢度域叠前地震记录的合成与并行方法研究

正演模拟是叠前弹性波反演的基础.采用慢度法计算层状介质的叠前地震记录,分别对频率和慢度进行积分变换得到时-空域的地震道集,并对在慢度积分过程中产生的计算噪音提出了解决方案.为得到高精度合成地震记录,需将地层细分,但地层层数很多时,计算量较大;而对地层粗分虽然会大大加快运算速度,但合成记录会丢失很多信息,文中给出了地层的划分原则.该方法能够计算出包括转换波和多次反射在内的全地震响应.但在提高合成记录精度的同时,也导致计算量增大、计算效率降低,因此,本文对基于慢度法全波场模拟进行了并行算法设计,采用计算域分割、工作池并行技术,建立了慢度法全波场正演模拟的并行算法,使得弹性波正演问题求解更加高效,为充分利用叠前地震资料进行叠前反演提供了研究基础.     

面向目标的合成曲面波叠前深度偏移及其应用

面向目标的控制照明叠前深度偏移属于面炮偏移方法的一种特殊情况.通过对目标区域波场照射方向的控制,实现波场在目的层的最佳入射,提高目标区域成像质量.将目的层定义的目标合成算子外推到地表形成地表合成算子,应用地表合成算子进行波面源和记录的合成,最后通过地表合成算子的一系列旋转,实现旋转控制照明叠前深度偏移.本文系统阐述了面向目标的合成曲面波控制照明叠前深度偏移方法的基本原理和实现方法,通过数值模型试算验证了该方法的有效性,最后将其应用到某探区实际资料的处理中,并进行了实际资料的应用效果分析,总结了其在实际资料应用中的特点.     

基于控制照明的合成震源记录交互剩余偏移速度分析

提出了一种新的偏移速度分析方法——基于控制照明的合成震源记录交互剩余偏移速度分析方法．与其他类似偏移速度分析方法的不同点在于：（1） 叠前深度偏移采用基于波动理论的快速合成震源记录算法;（2）偏移方法采用平面波震源,与速度分析方法一致;（3）应用控制照明技术,避免了因横向变速而导致的平面波震源波场在传播过程中的畸变,从而减小了速度分析的误差;（4）实用的速度谱设计,使交互偏移速度分析可行且易于操作．模型和新疆实际资料的试算表明该方法是一种有效和实用的偏移速度分析方法．     

福建省南平台记录吉林7级地震震源深度分析

用南平台模拟地震记录和数字观测记录 1 999年 4月 8日吉林东部 7级地震的不同深震震相来测定其地震深度 ,结果为南平台用地核界面波测定的震源深度与中国地震台网所测定的深度接近 ,用震中附近反射转换波定震源深度与中国地震台网所测定的深度有一定的误差。对此 ,本文进行了讨论。     

基于Gabor框架合成小波束叠前深度偏移

利用小波束域波场分解和传播在空间及方向上的双重局域性,通过Fourier传播算子进行深度外推波场,完成小波束的传播,最终应用成像条件求取成像值.合成小波束和小束源记录,进行部分小束源叠前偏移成像的数值试验,获得的成像结果表明,小波束叠前深度成像是切实有效的.     

用短周期地震仪的记录图纸分析远震和确定震源深度

0 引言 众所周知,固有周期不同的地震仪监测不同周期的地震波。短周期地震仪主要用于监测地方震、近震以及远震P波初至,而中长周期和长周期地震仪监测震中距大于1000km的远震和极远震。目前,地震台站都配备有短周期地震仪,笔者在分析乌鲁木齐台短周期地震仪的记录时,发现对部分远震,不管是DD-1(T_0=1.0s)可见记录仪还是62型(T_0=2.0s)短周期光记录仪,其图纸上S波、pP波sP波、PP波、PcP波、ScP波、ScS波的起始一般都很清楚(图1)。特别对于中源地震(h>60km)和深源地震(h>300km),在短周期地震仪记录图上的震相形态清楚,比中长周期和长周期地震仪都具有优势,显示出独特的作用。在中长周期和长周期地震仪记录图上有些震级小于6.0的中深源地震的震相往往根本看不到,而在短周期地震仪的记录图上却清楚可见。 有些大地震,由于震级较大,加之传播路经上地壳结构等因素的影响,在中长周期仪上记录的纵波段比较发育,横波的起始因纵波的迭加而不太清楚。短周期仪对大周期成份有滤波作     

基于共聚焦点道集的叠前深度偏移

共聚焦点(CFP)偏移技术是一种基于等时原理,将Kirchhoff积分法的一步偏移分两步聚焦（即激发聚焦和检波聚焦）来完成的叠前地震成像方法.该方法借助于逆时聚集算子和共聚焦点道集来实现叠前偏移成像.基于共聚焦点道集的叠前深度偏移是把基于共炮集的深度偏移的算法引入到CFP技术上来,基于波场延拓的理论来实现偏移成像,该方法首先生成共聚焦点道集,然后基于面炮合成的理论合成聚焦震源,最终通过相关成像来实现叠前偏移成像.该方法选取较少的聚焦点就可以实现对于地下构造的偏移成像,和炮域波动方程偏移相比,其计算效率得到了提高.通过模型试算和实际资料的试处理,验证了该方法在实现叠前深度偏移成像上的正确性和有效性.     

高精度混合法叠前深度偏移及其并行实现

叠前深度偏移是复杂构造成像的重要手段.本文基于波场分裂理论,首先给出了波场延拓的一般方程,即一个上下行波的耦合方程组.通常所用的波场延拓方程就是该耦合方程组的特例.根据平方根算子的近似,推导了一种新的高精度混合法偏移方法,用分裂法即可进行计算.通过对一个横向强烈变速模型的叠后偏移及Marmousi复杂构造模型的叠前偏移计算,说明了方法的有效性,精度较高.采用MPI并行编程实现并行计算,提高了计算效率.该方法可用于对横向强烈变速的复杂构造的精确成像.     

双平方根方程三维叠前深度偏移

从双平方根（DSR）形式的波动方程出发，基于沉降观测概念和地震波扰动理论，介绍了深度域的DSR全偏移算子及共成像道集的生成方法. 根据三维地震数据的方位角特征，通过对全偏移算子的稳相近似，依次导出了适应于零方位角道集、Cross line共偏移距道集以及共偏移距矢量道集的偏移算子. 理论分析与合成数据的数值试验表明，DSR全偏移算子、共方位角偏移算子对介质速度变化的适应性很强，而其余两种偏移算子仅适用于缓变速情况.     

傅里叶有限差分法保幅叠前深度偏移方法

地震数据中饱含有丰富的走时信息和振幅信息. 为解决传统偏移方法中几何扩散和入射角变化引起的振幅误差问题,本文提出了一种实用的波动方程保幅地震偏移方法. 该方法从全声波方程出发进行单程波保幅分解,得到直观、高效率的直接面对地震波传播波场的压力分量进行延拓的保幅偏移单程波方程,进而推导出一个含有6项的傅里叶有限差分法保幅偏移的算子方程;修改边界条件和成像条件,使修改后的边界条件和成像方程中考虑振幅补偿,从而从三方面补偿几何扩散损失和入射角变化对振幅的影响. 脉冲响应测试、单炮记录的数值试验以及SEG/EAGE盐丘模型的叠前偏移结果表明,该方法不但可以使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,为后续的地震属性分析(如AVO/AVA)提供更真实的地震信息.     

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

本文提出了一种基于双平方根算子的共炮检距道集波动方程保幅叠前深度偏移方法,将振幅误差补偿作为偏移的一部分与“运动学偏移”一起在偏移过程中实现.其基本内容包括：(1)从保幅的单平方根算子方程出发,推导出由双平方根算子定义的保幅单程波方程;(2)根据地震波摄动理论把速度场分裂为层内常速背景和变速扰动,分别在频率－波数域和频率－空间域求得波场深度延拓的偏移时移量及振幅校正系数,从而得到最终的DSR保幅波场延拓算子;(3)在高频假设条件下,把DSR保幅波场延拓公式中的积分运算进行稳相近似,得到保幅波场延拓的相移公式.理论分析和模型数值试验表明,该方法不但可以使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,为后续的地震属性分析（如AVO/AVA）提供更真实的地震信息.     

三维最小二乘弹性高斯束叠前深度偏移

与声波高斯束成像相比,弹性波高斯束偏移更适用于复杂油气藏多波多分量地震勘探.但是由于观测系统的局限性和深部构造的复杂性,该方法同样存在成像分辨率低、照明不均衡等问题.本文结合最小二乘偏移和弹性高斯束偏移的优势,提出了一种通过弹性高斯束叠加构建Born正演（反偏移）算子和偏移算子的三维最小二乘叠前深度偏移方法.依据最小二乘反演理论,建立基于反偏移数据与实际观测数据残差的目标函数,采用共轭梯度算法迭代更新来建立地下真实的反射率.与传统弹性高斯束偏移方法相比,该方法不仅提高了成像分辨率,而且使复杂构造特别是陡倾角地层的成像照明也得到了补偿.理论模型测试结果证明了本文方法的可行性和有效性.

     

波动方程的高阶广义屏叠前深度偏移

不同于常规广义屏传播算子的推导中使用散射理论,本文利用单平方根算子的渐近展开,推导出了单程波方程广义屏传播算子的高阶表达式.高阶广义屏传播算子不仅可提高常规广义屏传播算子的计算精度,而且还能改善广义屏传播算子对速度强横向变化介质的适应性.把高阶广义屏传播算子应用于波动方程叠前深度偏移,可得到比常规广义屏传播算子更好的效果.高阶广义屏传播算子的阶数越高,计算精度越高,但计算量也越多.以SEG-EAGE二维盐丘模型数据的波动方程叠前深度偏移为例,二阶广义屏传播算子相对于常规（一阶）广义屏传播算子增加了30%的计算量.高阶广义屏传播算子是常规广义屏传播算子理论的发展和完善.     

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

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

盐丘模型 推覆体模型 波动方程 叠前深度偏移

三维波动方程叠前深度偏移是复杂介质中进行构造成像、弹性参数反演的重要环节.由于其技术实现不仅涉及波场延拓理论的创新,而且需要大规模计算,因而研究难度较大. 本文以实验效果的取得为目的,完整地实现了SEG/EAEG盐丘和推覆体模型的三维波动方程辛几何算法的叠前深度偏移成像计算. 文中详细考察了所研制的波动方程三维叠前深度偏移软件系统及其对复杂地质构造的成像能力. 具体包括：1）对于盐丘模型,文中讨论了成像参数的选择、地震子波对成像精度的影响、完成二维及三维叠前深度偏移的比较;2）对推覆体模型,文中进行了脉冲响应测试;3）由两个模型的成像结果可见本文的波动方程三维叠前深度偏移软件系统已具有适应强速度横向变化、复杂构造的成像能力.     

Kirchhoff prestack depth migration in 3-D simple models: comparison of triclinic anisotropy with simpler anisotropies

We use Kirchhoff prestack depth migration to calculate migrated sections in 3-D simple anisotropic homogeneous velocity models in order to demonstrate the impact of anisotropy on migrated images. The recorded wave field is generated in models composed of two homogeneous layers separated by one either non-inclined or inclined curved interface. The anisotropy in the upper layer is triclinic. We apply Kirchhoff prestack depth migration to velocity models with different types of anisotropy: a triclinic anisotropic medium, an isotropic medium, transversely isotropic media with a horizontal (HTI) and vertical (VTI) symmetry axis. We observe asymmetry in migration caused by triclinic anisotropy and we show the errors of the migrated interface caused by inaccurate velocity models used for migration. The study is limited to P-waves.