基于峰值因子最大化的鬼波压制方法

针对传统拖缆采集得到的海上地震资料受鬼波影响而导致频带宽度变窄,笔者提出一种去鬼波方法,可以压制鬼波效应,拓宽频带。首先在频率-空间域推导出镜像数据生成公式,并基于峰值因子最大化的参数搜索方法来估计最优的鬼波延迟时间,将最优的鬼波延迟时代入镜像记录生成公式,获得最优的镜像记录,利用联合反褶积算法得到最优的去鬼波结果。将本文中的方法应用于模拟的变深度采集数据和海上实测的变深度采集数据,证实该方法可以有效地压制鬼波,消除频谱中的陷波点,得到宽频数据。     

深拖平缆采集资料中鬼波衰减方法探讨

鬼波衰减是深拖平缆采集资料处理中非常重要的一步,也是实现宽频处理的关键步骤。为此,笔者分析了深拖平缆资料特点,推导了鬼波的时距曲线方程,在此基础上,介绍了线性Radon域最小平方反演方法衰减鬼波的基本原理,该方法首先通过最小平方方法求得无鬼波的Radon域数据,再应用反变换算子中的一次波算子进行反变换,从而实现了对鬼波的压制。采用该理论对模型数据与实际地震资料进行了试处理,结果表明,深拖平缆资料通过本文的鬼波衰减技术处理后,陷频点得到了补偿,有效频带得到了拓宽,达到了宽频处理的效果。     

各向异性 Ｒadon 变换在深水地震多次波压制中的应用

在深海多缆、长拖缆地震勘探的各种噪音影响中,多次波对深层的有效反射影响难以消除。笔者介绍了一种高分辨率各向异性 Ｒadon 变换方法,在 Ｒadon 变换积分路径中增加非椭圆率 eta 参数,将传统 Ｒadon 变换积分曲线由慢度、偏移距两参数改为慢度、偏移距及非椭圆率三参数,更好地刻画了长偏移距下反射非双曲时差。通过对某海域实际采集的多缆长偏移距地震数据进行压制多次波处理,与传统方法相比取得了较为理想的效果,为后续深层构造成像提供了基础数据。     

各向异性Radon变换在深水地震多次波压制中的应用

在深海多缆、长拖缆地震勘探的各种噪音影响中,多次波对深层的有效反射影响难以消除。笔者介绍了一种高分辨率各向异性Radon变换方法,在Radon变换积分路径中增加非椭圆率eta参数,将传统Radon变换积分曲线由慢度、偏移距两参数改为慢度、偏移距及非椭圆率三参数,更好地刻画了长偏移距下反射非双曲时差。通过对某海域实际采集的多缆长偏移距地震数据进行压制多次波处理,与传统方法相比取得了较为理想的效果,为后续深层构造成像提供了基础数据。     

异常振幅衰减法压制海上船噪声干扰CSCD

海上船噪声干扰是海洋地震数据中的主要干扰波之一,由于其反射同相轴符合地震时距曲线规律,且其频谱与有效信号基本重叠,用滤波等常规噪声压制技术很难对其进行有效地压制。由于海上船相对于地震采集拖缆的位置具有随机性,共炮点域中的规则船干扰在共偏移距域内变为随机强能量噪音,为此提出一种用异常振幅衰减法压制海上船噪声干扰的方法,其思路是先将炮集抽成共偏移距道集,然后用异常振幅衰减法对船噪声干扰进行压制。通过对实际海洋地震资料进行处理,得到了良好的应用效果。     

异常振幅衰减法压制海上船噪声干扰

海上船噪声干扰是海洋地震数据中的主要干扰波之一,由于其反射同相轴符合地震时距曲线规律,且其频谱与有效信号基本重叠,用滤波等常规噪声压制技术很难对其进行有效地压制。由于海上船相对于地震采集拖缆的位置具有随机性,共炮点域中的规则船干扰在共偏移距域内变为随机强能量噪音,为此提出一种用异常振幅衰减法压制海上船噪声干扰的方法,其思路是先将炮集抽成共偏移距道集,然后用异常振幅衰减法对船噪声干扰进行压制。通过对实际海洋地震资料进行处理,得到了良好的应用效果。     

西湖凹陷斜缆采集关键参数优选研究

西湖凹陷是东海最大的含油气凹陷,其主要储层位于中深部,埋藏较深。增加低频信息是解决西湖凹陷主要地质问题的核心之一。斜缆地震采集技术可以有效地压制鬼波,拓宽频带,尤其是低频成分,并且具有野外操作简单、成本低等优势,对解决西湖凹陷的地质问题有利。基于斜缆平均鬼波滤波特性,对斜缆缆型、沉放深度范围等关键采集参数进行了理论分析;通过二维野外试验,对不同斜缆采集模式的单炮记录、偏移剖面开展了进一步的对比研究。结论认为:对于西湖凹陷深部储层,相较于常规水平缆资料,斜缆资料的构造刻画更清楚,频带得到了一定的拓宽。常用的两种斜缆缆型中,直线型斜缆的陷波补偿效应更好,抛物线型斜缆的低频能量更强。在斜缆最浅沉放深度固定时,加大斜缆沉放深度差,有利于获取更多的低频信息。     

基于宽频资料的扩展弹性阻抗反演方法在陆丰22洼陷低勘探区古近系岩性预测中的应用

低勘探程度区作为后备储量增长点越来越受到重视,但面临着许多地球物理难题.陆丰22洼勘探程度较低,钻井数量少,地质资料缺乏,依靠传统测井数据或叠加速度构建反演低频模型难度大、精度低,另外海洋地震资料受鬼波影响,压制了低频和高频波场成分,减小了地震资料的频宽,降低反演结果的保真度和精度.基于此,本文首先通过宽频处理技术获取...     

海上地震勘探斜缆采集中鬼波产生机理及压制效果分析

在海洋油气勘探中,为了准确获取勘探目的 层高精度、高分辨率的反射信息,实现更加细腻的成像,海洋地震宽频采集及处理技术近年来得到了极大的创新和发展.本文研究了海洋地震勘探中鬼波的产生机理、鬼波的分类和特性及其在实际地震资料中的表现形式.采用频率域高精度斜缆Radon变换作为主要处理方法,对某直斜缆宽频数据进行了鬼波的衰减...     

基于复杂海况鬼波压制的工区匹配拼接技术

海洋拖缆地震资料中鬼波会影响地震子波的形态,不同工区之间由于炮点、检点沉放深度不一致引起的鬼波周期和形态不一致,使得不同工区子波特征存在较大差异;由于常规拼接匹配技术没有考虑鬼波差异,解决不了鬼波的周期不一致问题,拼接处往往会存在明显的拼接痕迹,达不到理想的拼接效果。这里提出基于复杂海况鬼波压制(CSDG)的不同工区间的匹配拼接技术,对不同工区分别进行源、检鬼波及气枪信号进行压制,有效地解决了不同工区间的鬼波差异的问题,并在A工区及B工区数据匹配拼接中进行了验证,实现了两个工区的叠前无缝拼接处理,改善连片处理的数据品质,为后续的地质解释和油气预测提供便利。     

基于共轭梯度法的全波形反演

全波形反演是一种全新的地震成像方法,主要利用全波形信息反演地下介质参数,通过非线性优化波场理论值和观测值的残差实现波形反演。基于时间域声波方程,建立了时间域波场残差目标函数,以分层模型为例,分别从波场精度、目标函数收敛性和运行时间3个方面,比较了共轭梯度（Conjugate Gradient,CG）算法和拟牛顿算法（Broyden-Fletcher-Goldfarb-Shanno,BFGS）反演的效果。同时,应用共轭梯度法对正、逆断层模型和Marmousi模型进行了速度结构反演。反演结果表明:共轭梯度法计算效率较高,反演得到的速度模型精度更高,反演效果较好,是一种有效的波形反演方法。      

宽频地震处理技术在南海某工区的应用

在海洋地震勘探中产生的鬼波会抑制低频分量,造成频域的陷波点,降低地震记录的分辨率,使地震记录的同相轴发生混叠甚至产生假的同相轴。笔者采用FK域水平电缆压制鬼波法,将地震数据进行FK变换,对维纳滤波求取反褶积算子投影到FK域,并对二者进行褶积,将结果 FK反变换回时间域。鬼波压制后的模型数据频带得到了有效拓宽,低频、高频信息均得到了加强。鬼波压制后的实际数据目的层频带从12~60 Hz拓宽到8~75 Hz,三亚组能量与合成记录匹配较原始数据更好。     

一种加权抛物Radon变换地震波场重建方法

基于带限正反抛物Radon变换的基本原理和计算方法,笔者给出了一种加权抛物Radon变换（PRT）算法,该方法可应用于叠前地震波场恢复和重建,如缺失偏移距地震数据的插值外推和反假频重采样,其权系数完全由数据域来确定,也就是说在迭代计算过程中,该权系数也在变化。对于均匀采样的地震数据,笔者给出了近偏移距外推和波场的反假频重采样的计算结果,试算表明笔者给出的加权PRT方法可有效地应用于地震波场重建中,并具有算法稳健、计算精度高和效率高的特点。     

共偏移距可分表示法叠前深度偏移

基于双平方根方程的共偏移距可分表示法叠前深度偏移用于复杂介质成像, 该方法在中点-偏移距坐标中同时向下延拓炮点和检波点波场, 实现方式采用正反傅立叶变换.构造的双平方根方程波场延拓算子能够使波数域变量与空间(速度) 域变量分离, 波数域内进行相移计算, 在空间域对因介质横向变速引起的时移作修正.地震数据偏移不需要逐炮计算, 具有较高的效率.对Marmousi模型数据的偏移成像结果显示, 该方法较好地成像强横向变速介质中的复杂构造.      

The Offset-Domain Prestack Depth Migration with Optimal Separable Approximation

The offset-domain prestack depth migration with optimal separable approximation, based on the double square root equation, is used to image complex media with large and rapid velocity variations. The method downward continues the source and the receiver wavefields simultaneously. The mixed domain algorithm with forward Fourier and inverse Fourier transform is used to construct the double square root equation wavefield extrapolation operator. This operator separates variables in the wave number domain and variables in the space domain. The phase operation is implemented in the wave number domain, whereas the time delay for lateral velocity variation is corrected in the space domain. The migration algorithm is efficient since the seismic data are not computed shot by shot. The data set test of the Marmousi model indicates that the offset-domain migration provides a satisfied seismic migration section on which complex geologic structures are imaged in media with large and rapid lateral velocity variations.     

起伏地表条件下频率-空间域声波介质正演模拟

频率-空间域正演模拟是频率域及Laplace-Fourier域全波形反演的基础,起伏地表条件下波形反演算法的关键是正演算法中考虑起伏地表的影响。基于带PML吸收边界的声波波动方程,在已有最优9点有限差分正演算法的基础上构建了起伏地表条件下频率-空间域正演算法。通过应用变网格技术,进一步提高算法的计算效率、降低内存开销,使得大规模起伏地表模型的频率域正反演问题成为可能。理论分析及数值测试表明：通过对近地表区域进行局部网格加密,可有效地压制由于矩形网格离散引起的角点散射;结合变网格技术可较易获得5倍以上计算效率的提高及内存占用的降低,且随着模型尺度的增加及地表起伏高程差的减小,倍数将显著增加;在细网格与粗网格交界处产生的虚假反射振幅幅值控制在原始波场的2%以内,满足地震波场正反演的需求。     

波动方程角度域共成像道集

针对基于波场深度递推的波动方程叠前深度偏移成像方法难以准确地给出如Kirchhoff积分叠前深度偏移成像方法的偏移距域共成像道集的不足,在分析研究现有的各种波动方程偏移成像共成像道集方法的基础上,利用波场外推的单平方根算子和波场的窗口Fourier框架展开与重构方法,提出一种局部角度域共成像道集方法——成像点处反射角共成像道集方法.把这种角度域共成像道集方法应用于国际标准的Marmousi模型数据和一条实际二维地震数据都取得了理想的结果.提出的波动方程反射角度域共成像道集方法可为进一步的叠前偏移数据振幅随角度变化分析和偏移速度分析工作提供基础.     

井间地震时频域联合速度和衰减系数同步层析反演

井间地震资料中蕴含丰富的地下储层岩性和物性信息,利用井间地震直达波初至信息,通过层析反演可以得到两井之间的速度剖面。通过对井间地震直达波振幅衰减信息的处理,可以实现井间的吸收衰减层析成像。基于井间地震直达波射线路径的一致性,提出了井间地震时域直达波走时层析反演和频域质心频率衰减层析同步反演方法,增加了约束条件,增强了抗干扰能力。对模型和实际资料的处理结果表明,速度和吸收衰减同步反演方法,提高了速度层析反演精度,得到了对油气更为敏感的衰减属性,增强了井间地震在储层描述和油气预测中的应用能力。