基于源波场重构的子波自适应频移算子

地震信号中的多尺度信息对于分辨率、成像精度和反演结果有非常重要的意义,本文提出一种基于直达波模拟波场重建震源子波的地震数据频移算子,以期能应用于基于波动的地震信号多尺度分频.频移过程主要分为两步,第一步利用直达波反传构建震源子波,第二步借助震源子波和频移算子进行多尺度分频.与常规数字滤波器相比,频移算子突破了滤波造成信号波形特征改变和震源子波形态畸变的限制,频移地震数据与数值模拟地震数据完美匹配,同时频移算子具有理论子波自适应功能,更易于后续波动类应用展开.最后通过层状模型和实际资料进行测试,并与常规滤波器结果对比证明本方法的准确性和稳定性.     

利用直达波提取地震子波进行相关滤波

利用可控震源进行地震勘探时,通常都需要进行相关滤波,相关滤波结果的好坏取决于所选取的子波与有效反射波的相关性.海洋地震勘探常用的电火花震源激发能量通常可以调控,理论上,同一激发能量产生的震源子波应当是一致的,因此可以仿效可控震源资料的处理手段,对单道地震资料进行相关处理.但在复杂多变的海洋环境下,震源参数也要随勘探条件的变化而改变,采用理论计算或实验环境下获得的子波模型进行相关滤波往往难以达到预期的效果,需要寻找与实际勘探中激发的地震波相似度更高的地震子波.单道地震记录中直达波的能量较强,波形畸变较小,与震源子波较为接近.本文建立了一种利用直达波提取地震子波的方法,并应用到实际资料.实际资料的地震子波提取和相关滤波处理结果显示,利用直达波提取地震子波进行相关滤波能够有效地压制干扰波、突出有效波,提高了地震记录的信噪比.需要注意的是:1)提取子波的过程中要降低环境因素的影响,可利用统计平均的方法来实现;2)相关滤波时要根据实际处理效果的改善程度选择相关子波长度.     

基于特征能量梯度算子的时间域全波形反演

全波形反演中构建常规梯度算子过程中需要三步骤:震源子波的正向传波场传播,波场残差的反传波场和波场互相关构建梯度算子,其过程存在数据量大、效率低等缺点,为提高反演的效率,本文针对常规时间域梯度算子进行优化,提出了基于特征能量的梯度算法.在正传过程中计算每一个网格点上的正传波场的最大激发能量及其对应的时间步,保存一个子波时间长度的利用特征能量以构建梯度算子.在构建梯度算法中利用保存的子波长度的特征能量进行构建梯度算子.该算法无需保存震源的正传波场,可以减少运算过程中的磁盘读写,提高全波形反演的计算效率.在Mamousi模型梯度测试和实际资料的反演测试中表明:该算法在可以保证梯度算子的精度,具有数据读写量小的优点,效率高的优点.     

地震波频散效应与反Q滤波相位补偿

本文以Futterman提出的地震波振幅衰减和相速度频散表达式为基础,从井震匹配的角度出发,推导了不进行反Q滤波、仅反Q滤波振幅补偿、仅反Q滤波相位补偿,以及反Q滤波相位与振幅同时补偿四种情况下,地震波速度频散与相速度及地震记录振幅谱间的表达式,并从理论上说明了反Q滤波相位补偿的必要性.通过相同观测系统,炸药震源和可控震源分别激发采集的零偏移距VSP资料实例,验证了本文所推导的速度频散表达式的合理性;通过井震标定实例,进一步说明了反Q滤波相位补偿,可有效消除地震子波速度频散,提升地面地震资料与零偏移距VSP走廊叠加剖面的匹配度,最终提高地震资料成果的可信度.     

基于构造导向滤波的多震源最小二乘逆时偏移方法研究（英文）

多震源同时采集技术能够有效提高采集效率,得到来自多个震源的混合地震数据,该技术能缩短采集周期,降低采集成本,但对混合地震数据直接成像会在成像结果中引入串扰噪音,影响成像质量。因此,本文在实现多震源最小二乘逆时偏移的基础上,引入构造导向滤波算子作为多震源最小二乘逆时偏移的预条件算子,沿着构造走向应用非平稳滤波在有效保护构造信息的前提下压制串扰噪音,通过共轭梯度法使得反偏移数据与观测数据之间的误差达到最小,最终得到信噪比更高的成像结果。合成数据的数值实验表明,本文所提出的方法能够有效压制串扰噪音,提高成像精度。     

消除探地雷达数据的子波衰减和频散的反滤波方法

消除探地雷达数据的子波衰减和频散可以很好地提高探地雷达的勘探深度和勘探分辨率.常用的消除探地雷达数据的子波衰减和频散方法为反Q滤波方法.该方法需要利用地下介质的Q参数,但是正确求取地下介质的Q参数很困难.针对这一问题,本文提出了一种消除探地雷达数据的子波衰减和频散的反滤波方法.该方法以地下介质反射系数是随机数为前提,利用地下介质等效滤波器具有最小相位这个特性,通过求取等效滤波器的振幅谱来求取等效滤波器的反滤波器.最后,利用该反滤波器对探地雷达数据进行反滤波,实现消除探地雷达数据的子波衰减和频散.     

GI枪结构震源子波模拟研究

GI枪不同于常规空气枪,它虽然使子波峰-峰值减少,但是具有较好的压制气泡振荡的优点.现在GI枪已经成为海上地震勘探主要震源激发方式之一.因此研究GI枪子波特性和模拟子波过程就有理论与现实意义.GI枪震源子波模型是依据GI枪激发特性,基于空气枪震源子波模型理论,以准静态开放式热力学系统与气泡振荡方程为依托建立的.模型根据模拟经验参数模拟震源子波结果,并与GI枪实测子波资料进行了对比,模拟结果表明能较好的吻合.文中同时给出子波模拟经验参数,使模型标准化.     

反卷积在人工震源数据处理中的应用

震源到接收台站之间的地层响应函数能够反映地下介质信息。对地震波传播过程中的卷积模型进行推导:记录信号是众多震源子波经过时移加权叠加的结果;通过反卷积方法可去除震源子波信息,提取震源到接收台站之间的地层响应函数;地层响应函数中第一个突跳值对应的时间即为P波走时。在河北赤城—张北地区进行人工震源实验,通过反卷积计算得到该地区地层响应函数剖面图,得出P波波速约6 km/s。利用人工震源系统还可以对地下介质波速变化进行长期动态监测,对地震预测具有一定意义。     

地震子波处理的二步法反褶积方法研究

针对玛湖斜坡区三块三维地震资料和赛汉塔拉凹陷二块三维地震资料连片处理中的特点,结合地质任务和处理目标要求,提出了地震数据连片处理中的地震子波处理的方法.该方法主要体现了两次反褶积,一次是采用地表一致性反褶积,将不同震源的频带拓宽到一个标准上;再一次采用相位校正反褶积,将不同震源的数据校正到相同相位上.为了保证提取的相位校正反褶积算子稳定,采用叠后地震道提取（主要考虑到叠后地震道信噪比高,算子稳定性强）,然后将该算子应用到叠前地震道,进行相位校正.     

逆时偏移震源子波影响分析

逆时偏移作为一种先进的地震偏移成像方法,其成像结果的好坏取决于很多因素,其中成像奈件是关键的一个因素.成像条件一般采用互相关成像条件,它是将源波场沿时间进行正推,接收波场沿时间进行逆推,然后将源波场和接收波场进行互相关,从而得到偏移成像结果.源波场正推和接收波场逆推都是采用同一个波动方程,这里采用波动方程一阶应力-速度形式,数值计算方法采用交错网格有限差分方法.采用振幅补偿拉普拉斯滤波方法压制逆时偏移成像中的低频噪声.采用三次样条插值方法解决逆时偏移成像中的波形不光滑问题.在接收波场逆推过程中,是将地震记录作为边界条件.而在源波场正推过程中,需要给定震源子波作为初值条件.理论上来讲,震源子波为脉冲时是最佳的,但在实际当中难以实现.在源波场正推时,一般采用子波函数.子波函数形式非常多,最常用的是雷克子波.不同主频震源子波产生的源波场频率也不同,从而导致逆时偏移互相关成像效果也不同.那么震源子波主频对逆时偏移成像影响到底有多大呢?为了简便,本文基于不同主频雷克子波,初步分析了相应的逆时偏移成像效果,同时对速度模型存在不同误差时的情形也进行了试验分析.结果 表明:在不同速度误差情况下,震源子波主频越低,成像同相轴连续性越好,散射点聚焦能量越强;震源子波主频越高,虽然分辨率提高,但成像同相轴连续性变差,散射点聚焦能量变弱,同时背景噪声增强.本文研究结果对逆时偏移成像技术在实践中的应用具有一定的参考价值.     

零偏移距VSP资料层Q反演及其应用研究(英文)

地层品质因子Q的可用于地震资料高分辨率处理,而从VSP资料下行直达波更容易获取准确的地层品质因子。通过对零偏移距VSP资料的监控子波和下行初至波的频谱进行综合分析,仿照Ricker子波频谱的表达式,本文提出了震源子波频谱新的表达式。在震源子波频谱新的表达式基础上,我们介绍了改进的频谱拟合法和改进的谱比法的层Q值反演方法及相应的处理流程。基于本文提出的层Q值反演方法,利用实际的零偏移距VSP资料的下行直达波,反演稳定的层Q值,并用于零偏移距VSP资料及井旁地面地震资料的反Q滤波振幅补偿处理,提高了地震资料的分辨率。     

衰减干涉假象的海洋一次波与多次波同时成像方法

海水与空气间的强波阻抗差使得海洋地震资料普遍发育自由表面相关多次波,如何利用好多次波所携带的有效信息已成为提高海洋地震资料成像品质的新突破点.基于面炮偏移的一次波与多次波同时成像方法能够避免多次波预测精度的影响,但是,正向传播的震源子波与反向延拓的自由表面相关多次波所产生的干涉假象严重制约了该技术的应用,本文提出了一种基于单程波偏移算子,可在成像域压制干涉假象的一次波与多次波同时成像方法.其中包含了三个步骤:第一,传统单程波偏移成像方法中的震源子波替换为一次波、多次波与震源子波,初始上行延拓波场为一次波与多次波,基于单程波算子的波场延拓与互相关成像条件的应用得到包含干涉假象的一次波与多次波同时成像;第二,以子波为震源,自由表面相关多次波为记录,按照传统单程波偏移成像方法得到干涉假象;第三,基于最小二乘匹配滤波算法,将第一步的成像结果与第二步的干涉假象进行匹配相减,得到干涉假象衰减后的一次波与多次波同时成像,避开了由于实际资料子波无法准确提取而造成一次波与多次波对成像能量级的不一致性.Sigsbee2B模型测试验证了本方法的有效性,并在我国某探区深海实际资料处理中得到了成功应用,深层基底得到了清晰刻画,并且照明均衡度明显改善.     

用于深度域地震剖面衰减与频散补偿的反Q滤波方法

根据地震波在非弹性介质中的传播规律，本文提出了一种在深度域地震剖面进行反Q滤波的新方法. 在深度域的Q滤波算子符合地震波在衰减介质中的传播规律，不仅考虑了介质吸收对地震波振幅的影响，而且还保证了所造成的波形畸变满足因果规律，即地震体波具有某种频散性质. 利用该方法处理了2条总长度约60km的某浅海地区叠后深度偏移地震剖面. 合成记录剖面的正演和反演计算结果表明，在对地震剖面进行深度域反Q滤波之后 ，不仅达到了压缩中、深部目的层子波宽度提高分辨率和信噪比的目的，而且保持了原始地震剖面的整体形态和能量分布.     

基于Gabor变换反褶积技术在渤海某工区的应用研究

使用常规的Wiener反褶积必须假设震源子波在地层旅行过程中是平稳的即一成不变的,这个前提条件与实际野外地震资料采集差别较大,而基于Gabor变换反褶积技术考虑到地震能量的衰减、子波的形变等非平稳性特征.地震道在Gabor域可因式分解成三项即震源子波、衰减函数和反射系数,该技术设计POU窗函数,并利用此函数在Gabor域对地震信号进行局部时频分解.Gabor域反褶积算法在Gabor域通过除以衰减函数和震源子波的乘积来估算地层反射系数,然后再做Gabor反变换可求得时间域的地层反射系数.理论模型的测试和实际地震资料的应用均表明,与Wiener反褶积相比较,基于Gabor变换反褶积可补偿中深层的能量衰减并因此拓宽有效频带和提高时间分辨率.     

基于矢量预测的SIMO系统混合相位地震子波提取方法研究

针对高阶统计量混合相位地震子波提取方法的局限性, 提出一种基于矢量预测的单输入多输出系统(SIMO)的混合相位地震子波提取方法. 该方法利用二阶循环平稳统计量包含的系统相位信息, 将CMP道集视为一个单输入多输出系统的输出, 利用道集中相邻两道或多道数据通过矢量预测来构建反子波计算的方程式, 进一步进行混合相位子波提取. 利用提取出的子波相位信息对CMP道集进行纯相位滤波, 能够替代相位校正技术, 并且利用获取的反褶积算子对CMP道集进行反褶积处理, 使提高分辨率后的不同道子波振幅、 频率、 波形相一致, 提高叠加的质量. 模型试算和实际资料处理结果表明, 文中方法适用于任意相位的子波提取及反褶积处理, 且处理精度较高, 具有较高的实际应用价值.      

同态解卷在地震测深中的应用

在地震测深资料的处理分析中,震相的识别是至关重要的一步,只有震相识别正确,后续的一系列反演工作才能顺利进行。而目前震相识别在国内外地震测深中仍是一个令人困惑的问题。 本文将同态解卷技术应用到地震测深中,从一个新的角度开展了工作。首先用功率谱滤波法对记录进行预处理,降低噪音的影响。而后,根据同态解卷对震源子波估计的有效性,提出了多道时域迭加法,使子波估计更加真实可靠。在此基础上,利用震相间子波的相似性进行震相识别。理论模拟及实际资料处理的结果表明:该法是行之有效的,它为震相识别提供了一条新的途径。     

基于误差准则和循环迭代的时移地震匹配滤波方法

分析了用于时移地震互均衡处理中的传统匹配滤波算法适用范围，结果表明该方法在滤波算子的反算子为最小相位时效果最好.对该方法难以解决的混合相位问题，推导出通用公式，并提出基于最小平方误差准则和循环迭代的求解方法.理论和实际数据都验证了该方法较已有方法的优越性.     

基于非稳态反演的气枪阵列子波方向性反褶积（英文）

由于气枪阵列所激发的子波具有能量强,气泡比高等优点,因此,海上地震勘探大多采用气枪阵列进行激发。但是气枪阵列具有一定的长度和宽度,加之海水虚反射的影响;致使气枪阵列子波具有明显的方向性效应,破坏了地震子波一致性,改变了反射振幅随入射角和方位角的函数关系。针对这一问题,本文提出了一种基于非稳态反演的方向性反褶积方法。该方法首先根据气枪阵列的空间配置和近场子波计算与方位角和出射角有关的远场子波。然后,基于速度模型计算不同时间地震反射在震源位置的出射角,构建方向性滤波算子。最后,采用非稳态反演方法将不同方向的地震子波整形为阵列正下方的远场子波,实现气枪阵列子波方向性反褶积处理。     

一种改进的广义稳定反Q滤波方法

由于地层的黏弹性和非均质性,地震波在传播过程中振幅衰减,相位畸变,导致地震资料分辨率和信噪比降低.反Q滤波是补偿地层衰减的一种常用方法.常规的反Q滤波方法的振幅补偿算子只与时间、频率有关,其振幅补偿频带范围、强度固定不变,不能够根据地震记录的自身特点灵活地进行振幅补偿.基于此,本文在广义稳定反Q滤波法的基础上,提出了一种改进的广义稳定反Q滤波方法.新方法提出了一种新的稳定因子,并通过参数p对振幅补偿算子进行了优化.其振幅补偿算子与地震资料信噪比、走时、频率、品质因子相适应,可以更加精确地对振幅进行补偿.理论数据测试和实际地震资料应用均表明,新方法在振幅补偿和噪声压制上优于广义稳定反Q滤波法,其处理后的地震资料振幅恢复保真度、分辨率和信噪比更高,对储层预测和精细解释工作具有重要意义.     

Ricker子波核最小二乘支持向量回归滤波方法的稳健性研究

除了信噪比、有效子波畸变等,稳健性(Robustness)也是度量滤波方法效果的一个重要的物理量,它刻画了滤波系统应对异常点值的能力.一般用影响函数作为评价稳健性的工具.支持向量机方法已较成功地应用于信号与图像的滤波中,尤其Ricker子波核方法更适于地震勘探信号处理.通过考察Ricker子波核最小二乘支持向量回归(L...