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

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

二维频率空间域粘声波正演模拟研究

地球介质本身的非弹性导致大地吸收效应,黏滞性影响波场的所有频率成分且对高频成分的影响更大,从而降低了地震记录和地震成像的垂向分辨率.通常使用经验公式描述吸收衰减.在时间域模拟吸收衰减需大量的计算机内存、计算效率较低,模拟Q值非常困难;而频率空间域可以通过引入复速度和Q值,使用任意的经验公式而进行粘滞介质的数值模拟.本文...     

二维频率空间域25点优化系数差分格式弹性波数值模拟(英文)

频率空间域地震波数值模拟具有独特的优势:可以同时模拟多源的波传播、每个频率之间独立并行地计算、计算频带选择灵活、不存在累计误差、容易模拟粘弹性介质中地震波传播.但是该方法的最大瓶颈是对于计算机内存的巨大需求.我们使用压缩存储系数矩阵的方法,极大地减少了计算机内存的需求量.同时为了减少短筹分算子的数值频散,引用了频率空间域25点弹性波波动方程的差分格式,并使用了最小二乘意义下求出的优化差分系数.为了克服边界反射,采用了最佳匹配层吸收边界条件.数值模拟试验证明:用压缩存储系数矩阵及优化差分系数的频率空间域25点差分格式进行弹性波正演模拟,可以减少数值频散,提高计算精度.使用较大的网格间距,降低计算机内存需求,并保持较高的计算效率.该正演方法为后续弹性波偏移和弹性参数反演提供较好的基础.     

倾角分解共反射面元叠加方法

共反射面元（Common Reflection Surface）叠加是一种独立于宏观速度模型的零偏移距剖面成像方法,传统的CRS叠加实现是以数据驱动的方式对属性参数进行自动搜索并对其进行优化合成相应的CRS叠加算子,通过该算子进行叠加能够得到信噪比和连续性更高的零偏移距剖面.但是数据驱动的实现方式带来了不可避免的“倾角歧视现象”,它造成了弱有效反射信号损失和运动学特征失真的问题.本文提出的倾角分解CRS叠加方法成功解决了上述问题,使CRS叠加方法更具实用价值.     

窄方位地震数据双平方根方程偏移方法探讨

简单回顾双平方根（DSR）方程偏移一般理论之后,介绍了适应偏移速度分析的角度域成像快速算法.针对地震勘探中窄方位三维地震数据成像问题,探讨了限定输入数据体或限定波场传播方向的DSR方程偏移方法.数值试验表明,共方位角偏移等DSR方程叠前深度偏移方法在盐体或礁体、古潜山等强横向变速地区具有很好的应用前景.     

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

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

三维弹性波数值模拟中的吸收边界条件

在地震波传播数值模拟的过程中,需要使用吸收边界条件从而达到衰减人为边界 反射的目的. 本文基于傍轴近似法提出了计算三维弹性波方程的吸收边界条件公式,表示了 各边界面、边棱和角点处波场所满足的单程波方程,并在三维弹性波数值模拟中进行了应用 . 理论模型及三维盐丘地质模型波场切片快照试算结果表明,该吸收边界条件可以有效地吸 收人为边界反射,适用于较大入射角情况,从而消除了边界有效波信息的干扰. 由于采用四 阶近似方程,在保证计算精度的前提下,该方法具有节省计算工作量和易于实现的特点.     

输出道成像方式的共反射面元叠加方法I——理论

共反射面元（Common Reflection Surface）叠加是一种独立于宏观速度模型的零偏移距成像方法，该方法属于典型的克希霍夫型成像方法. 根据成像方式的不同，克希霍夫型成像方法可以分为两大类：输出道成像方式和输入道成像方式. 考察共反射面元叠加方法，它属于输入道成像方式. 本文基于理论模型数据，实现了输出道成缘方式的CRS叠加方法. 相比传统的输入道成像方式，它具有能够保证大偏移距反射信息的成像精度和计算效率较高的优点，而且更加容易推广到三维情形.     

三种阵列声波测井数据频散分析方法的应用与比较

针对正交偶极声波测井检测地层各向异性和确定构造应力中频散分析的必要性，本文利用合成和实测阵列声波测井数据分析和比较了三种频散分析方法：Prony方法、同态处理方法和谱域加权相似法. 这三种频率域处理方法提取得到了合理的、一致的慢度结果. 谱域加权相似法只能用于每一频率下单一波模式的慢度估计；同态处理方法还可以同时估计单一波模式的衰减、幅度和初相位；利用SVD_TLS算法确定波模式的阶数后，Prony方法可以有效地估计每一频率下多个模式的上述参数. 另外，通过谱域加权相似法对正交偶极声波测井阵列数据进行频散分析，观察到了弯曲波频散曲线的交叉现象，从而确定出了对应地层所受最大水平主应力的方向.     

海洋天然气水合物的地震识别方法研究(摘要)

天然气水合物作为21世纪新的自然能源将为人类的生存发展服务。20世纪60年代证实,俄罗斯西伯利亚的麦素亚哈气田为典型的天然气水合物形成的气田,70年代又在海底发现了固体天然气水合物岩样。1971年,R Stoll首先将地震剖面中的似海底反射层解释为海洋天然气水合物存在的标志,后来被深海钻探证