基于最小二乘反演的拖缆双检鬼波压制方法

基于波场传播原理,分别推导了海上拖缆双检地震采集中压力波场和速度波场鬼波算子,建立了检波器接收信号和上行波场线性关系,利用最小二乘反演和匹配滤波算法求解拖缆处上行波场,该方法无需地下介质信息,实现数据驱动鬼波压制,通过引入噪声影响因子,减少因水、陆检信号的信噪比差异对鬼波压制效果的影响。数值模拟资料和拖缆双检实际资料处理结果表明,基于最小二乘反演的拖缆双检鬼波压制方法可以有效压制鬼波能量,拓宽资料频带,有利于后续资料的处理和解释。     

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

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

基于波动方程转换的时间域多尺度全波形反演速度建模

为了避免全波形反演的周波跳跃现象,提出了基于波动方程转换的时间域多尺度全波形反演速度建模策略,在时间域实现了从低波数到高波数的多尺度全波形反演。首先从声波方程参数化模式出发,研究了阻抗-速度和速度-密度两种参数化模式下速度的辐射模式:在阻抗-速度参数化模式下,速度扰动主要引起大角度波场扰动;在速度-密度参数化模式下,速度扰动对各个角度的波场扰动贡献量完全相同。基于此,提出了先利用阻抗-速度方程构建低波数全波形反演速度模型,再以此作为初始模型,利用速度-密度方程构建高波数全波形反演速度模型的方法。该方法有效避免了混合域全波形反演中的数据转换问题以及频率域反演中的吉普斯现象,同时充分发挥了时间域全波形反演在波动方程数值模拟计算效率方面的优势,保留了时间域数据匹配易控制的特点。通过MarmousiⅡ模型数据测试,对比分析了两种参数化模式下的速度梯度特征,实现了从阻抗-速度方程的低波数全波形反演速度建模到速度-密度方程的高波数全波形反演速度建模,说明该方法能够在初始速度缺失低波数的条件下充分刻画出断层的形态和位置,使断面清晰,地层起伏与真实模型吻合。     

起伏地表逆时偏移在复杂山前带地震成像中的应用

针对山前带"复杂地表、复杂构造"双复杂地震地质条件造成的资料信噪比低、静校正问题严重、地震波场复杂等一系列地球物理难题,在"真地表"地震成像面的确定及高频静校正的基础上,以基于起伏地表的深度域速度分析与建模为重点,以起伏地表逆时叠前深度偏移为核心建立了一套高精度地震成像处理流程,将长波长静校正问题隐含在偏移成像过程中,直接从起伏地表进行波场延拓和偏移成像,以便更好地应对复杂山前带的地震成像问题。针对性模型试算和实际资料处理表明,该技术在应对双复杂地震资料偏移成像时具有更高的精度,是复杂山前带资料高精度地震成像更理想的技术手段。     

叠前逆时深度偏移成像的实现与应用

介绍了高阶有限差分求解双程波动方程和完全匹配层吸收边界条件的基本理论原理.根据逆时偏移实用化所面临的技术瓶颈,提出了GPU/CPU异构协同并行方案和震源波场重构存储策略,有效解决了逆时偏移对海量计算和存储的需求.从逆时偏移成像噪声的产生机理出发,根据噪声的传播方向一致特性、大角度特性和低频特性,采用组合压噪流程,建立针对复杂构造成像的实用化叠前逆时深度偏移技术,最大程度地压制偏移噪声、保留有效信号.理论模型测试和实际资料处理表明,该技术能够有效解决复杂地质条件下的地震成像问题,尤其是对盐丘边界和盐下等特殊构造的刻画均有显著效果.     

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

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

VTI介质自适应聚焦束偏移

自适应聚焦束是高斯束的改进,采用多次聚焦的思想,考虑局部速度场对波束宽度的影响,将地震波束有效能量控制在一个波长的范围内,与常规高斯束的不同点在于动态选择初始波束参数。将自适应聚焦束偏移方法应用到VTI介质,基于经典弹性参数表征下的各向异性射线追踪方程体系,结合自适应聚焦束处理强横向变速成像问题的优势,实现了VTI介质自适应聚焦束偏移。通过与传统高斯束形态的对比,突出了自适应聚焦束形态控制的优势。利用各向异性Sub-Sag模型验证了各向异性介质自适应聚焦束偏移方法的正确性和稳定性,各向异性SEG/Hess模型测试结果证明了该方法对复杂地质模型的适应能力;实际地震资料的处理结果表明,该方法成像效果优于传统各向异性高斯束偏移方法。     

波动方程保幅地震偏移成像方法

地震数据中饱含有丰富的旅行时信息和振幅信息。传统的波动方程地震偏移成像方法忽略了介质参数变化对地震波振幅的改造,不能准确描述地震渡传播的动力学特征。为此,基于严格的解耦理论,从全声波方程出发,进行单程波保幅分解,得到直观、高效率的的单程波动方程,直接面对地震波传播波场的压力分量进行延拓,从而将振幅误差补偿作为偏移的一部分在偏移过程中进行,实现地震偏移成像在获取正确位置的同时也获取真实振幅的信息。数值试算表明,该方法可使散射能量聚焦、归位,提高成像精度;可输出正确反映地下反射系数的振幅信息。     

波动方程保幅地震偏移成像方法

地震数据中饱含有丰富的旅行时信息和振幅信息。传统的波动方程地震偏移成像方法忽略了介质参数变化对地震波振幅的改造,不能准确描述地震渡传播的动力学特征。为此,基于严格的解耦理论,从全声波方程出发,进行单程波保幅分解,得到直观、高效率的的单程波动方程,直接面对地震波传播波场的压力分量进行延拓,从而将振幅误差补偿作为偏移的一部分在偏移过程中进行,实现地震偏移成像在获取正确位置的同时也获取真实振幅的信息。数值试算表明,该方法可使散射能量聚焦、归位,提高成像精度;可输出正确反映地下反射系数的振幅信息。     

各向同性介质弹性波交错网格有限差分正演模拟

在弹性波有限差分数值模拟方面，差分网格及边界条件是影响弹性波模拟成功与否的关键。从各向同性介质速度一应力方程出发，利用交错网格高阶有限差分和完全匹配层（PML）边界条件，提出了各向同性介质弹性波交错网格有限差分正演模拟方法。数值试验结果表明，该方法精度较高，数值频散较小，人为边界反射吸收较好，为后续的地震属性分析（如AVO／AVA）奠定基础。