改进的正交多项式变换在地震资料去噪中的应用

水平叠加虽然在很大程度上压制了噪声,提高了地震剖面的信噪比,但CMP遭集上还存在不少不是一次波的规则干扰和随机噪声,不利于叠前资料的岩性反演和叠后资料的波阻抗反演.本文提出了基于改进的正交多项式变换压制地震资料中随机噪声的方法,其优势在于:通过对不同时间信号的奇异值分解,确定有效信号正交多项式系数谱的阶数;再利用小波变换,改善有效信号和噪声在低阶上的混叠.文中给出了具体处理的过程,数据试验和实际资料的处理结果表明该方法不仅能有效地压制噪声,而且还能较好地保护地震数据中AVO变化特征.     

拉东投影法三维叠前深度偏移

对地下地质构造进行正确成像是地震勘探的最终目的,由于三维地震资料采集不可能都沿垂直构造走向的方向进行,为地震资料的三维处理带来了许多困难. 本文将三维叠后拉东投影偏移思想应用于三维叠前处理中,提出了三维叠前投影偏移算法. 利用拉东投影变换的原理,将整个三维叠前数据体投影到一系列各方向的径向线上,各方位角的构造都包含在其中某条或多条径向剖面上. 投影完成后,形成一系列的独立的二维叠前测线,可采用二维叠前深度偏移成像方法来实现各径向线的叠前偏移,当各径向剖面偏移完成后,在时间切片上进行反投影,从而最终形成三维叠前深度偏移结果. 实际应用表明,用本方法进行三维叠前深度偏移时,深度偏移剖面对横向分辨率有所提高,对陡地层和小断层的成像效果有所改善.     

基于贝叶斯反演的叠前数据同时插值和保幅去噪方法研究

近年来,地震勘探逐渐面向深层、复杂构造和岩性油气藏,与此同时,对地震资料处理也提出了更高的标准和要求.保幅降噪是地震资料处理中的重要环节.叠前地震资料的品质直接影响后续反演结果的精度,但道编辑后的实际数据除信噪比较低之外还存在缺道等现象,常规的保幅去噪方法无法同时进行道插值.本文基于贝叶斯反演框架提出了一种能够同时进行道插值的叠前地震资料保幅去噪方法.该方法充分考虑了现有方法的局限性,利用插值算子对缺失地震道进行插值处理,同时基于地震振幅随角度变化的连续性在差分域采用贝叶斯框架引入了先验约束,能够在保幅去噪的同时,进行道插值,极大改善了叠前地震资料的品质.合成和野外的例子均验证了该方法的有效性.     

叠前时间偏移技术在淮北某煤矿采区三维地震勘探中的应用与效果

淮北某煤矿采区三维地震勘探地质条件极为复杂,区内褶皱、断层构造发育,常规叠后时间偏移处理很难准确成像,给后期资料解释工作带来不利影响.针对这一问题,采用叠前时间偏移技术,结合煤田地震勘探特点,对影响叠前偏移处理质量的关键步骤（预处理、静校正、叠前去噪、时间域速度模型建立）进行了研究,取得了显著效果.利用叠前偏移数据体解释,最终圆满完成了地质任务.本文通过应用实例,说明了煤田地震勘探资料处理利用叠前时间偏移技术对复杂构造成像效果较好,信噪比提高,断点解释更可靠.     

广义Radon变换与叠前地震数据处理

在本文中,首先讨论了与几种地震层析成像对应的Radon变换公式,并导出了叠前地震记录的数学模型.在分析叠前地震记录与广义Radon变换的关系的基础上,讨论了速度分析、滤波、动校正、叠加等地震数据处理的数学物理意义.为展示广义Radon变换在地震数据处理中的应用,给出了用于滤波和消除多次波的方法及算例.     

起伏地表叠前逆时偏移的实现及应用

复杂地表复杂构造是近年来陆上地震勘探的重难点问题.其中,复杂构造成像问题可用叠前深度偏移技术如逆时偏移较好地解决,而地表起伏大、横向变速剧烈引起道间时差大问题则需要采用起伏地表偏移的思路.结合常规地震数据处理,本文提出一种基于起伏面的逆时偏移成像技术流程,该流程利用固定面与平滑起伏面高程差和替换速度将静校正后的叠前炮集反静校正到起伏面上,然后进行起伏面逆时偏移成像计算.将流程应用到模型数据和实际资料处理,结果均表明该流程在复杂地表复杂构造地震成像方面具有很好的效果.     

海上气田浅层气下覆复杂构造地区Q值补偿及逆时偏移成像——以巴西桑托斯盆地应用实例

针对海上资料特点及研究区上覆碳酸盐岩地层成像差、多次波发育、受浅层气影响的下层构造成像模糊,振幅异常及能量损失、地震反射信号弱、信噪比低的难点,以叠前综合去噪、预测反褶积、CMP道集高精度Radon变换衰减层间多次波为核心前期处理方法,结合基于黏弹性介质地震波衰减模型的Q值补偿技术,利用逆时偏移对复杂断层进行了精确成像,处理后的地震剖面较好地反映出了砂岩、碳酸盐岩储层的形态和构造特征,目的层断块、断点清楚,证明设计流程及方法可以较好地解决该地区目的层成像难的问题,并可应用于相似浅层气地区成像.     

叠前同步反演在古近系窄河道储层预测中的应用——以渤海A油田为例

在渤海海域,河流相储层较发育,一般具有单层厚度薄、纵向叠置关系复杂、岩性及物性横向变化快等特点.古近系窄河道砂岩储层预测还面临埋藏深,砂岩与泥岩纵波阻抗差异小;地震资料垂向分辨率不够等问题.常规方法进行储层预测及刻画难度大,严重制约着开发井网部署.本文以测井岩石物理分析为基础,发现密度能很好地区分储层与非储层.采用密点速度分析、拉东变换去多次和自适应道集拉平等处理技术,提高了叠加剖面的质量,使处理后的资料满足叠前研究需求.应用叠前AVO同步反演技术对A油田的三维地震资料进行叠前反演,反演出密度数据体,提取地震属性,清楚刻画出河道的平面展布形态,据此优化钻井方案.实钻开发井储层钻遇情况与钻前储层预测一致,规避开发井的风险,提高经济效益,该方法可供类似油田借鉴.     

叠前时间偏移技术在深层断陷区地震成像中的应用

Kirchhoff积分叠前时间偏移应用波动方程的Kirchhoff积分解实现地下反射层的偏移问题,该技术应用所有偏移距的地震资料,能适应纵横向速度变化较大的情况,是复杂地区地震资料成像的较理想的方法.叠前时间偏移是对偏移处理和偏移速度修改的一个迭代过程,偏移速度的精度直接影响到偏移效果.对大庆油田古龙断陷某工区的三维地震资料,应用叠前时间偏移技术进行成像处理,对于埋深较大的目的层,断陷结构复杂,介质的各向异性较为严重,所以应用考虑各向异性参数的速度修改公式,得到较为精确的偏移速度体,进而保证偏移CRP道集拉平,得到成像良好的偏移数据体.实际地震数据偏移处理结果表明叠前时间偏移技术可对复杂构造的地震数据准确成像,可在古龙地区断陷地层勘探中推广应用.     

复杂地表条件下共反射面元（CRS）叠加方法研究

在地表地形复杂的情况下,静校正不易做好,这是制约山地资料处理质量的一个很重要的因素.复杂地表共反射面元（CRS）叠加不需对叠前数据做静校正,而且在得到叠加剖面后可以利用叠加得到的波场参数剖面实现基准面重建.地震数据的试算表明,复杂地表CRS叠加得出的剖面与常规处理剖面相比有着较高的信噪比和同相轴连续性.与水平地表CRS叠加不同的是,在复杂地表CRS叠加的时距公式中,波场三参数耦合,难以通过简化CRS道集的方法将它们全部分离并逐个优化.引入模拟退火算法后,有效地解决了这一组合优化的难题.     

基于压缩感知的加权MCA地震数据重构方法

地震数据规则化重构是地震资料处理十分重要的基础性工作.压缩感知理论打破了香农采样定理的制约,利用信号在某个变换域的稀疏特性重构出完整的信号,在地震数据重构领域得到了很好的应用.深反射地震剖面大都布置在地质构造比较复杂的区段,复杂的地质构造使深反射地震剖面上的波阻特征复杂,采用单一稀疏变换不能最有效地表征数据的内部结构特征.MCA（形态成分分析）方法将信号分解为几种形态特征区别明显的分量来逼近数据的内部复杂结构,但是对各成分简单的叠加仍然无法有效地描述复杂构造数据的各种特征.结合两种方法的优点,本文提出了一种新的基于压缩感知的重构算法框架,在MCA方法的基础上对各稀疏字典进行加权,在迭代中不断更新各个稀疏字典的权值系数,对信号内部的各种特征进行最优描述,从而实现对信号的高质量重构.模型测试和实际资料处理结果表明：基于压缩感知的加权MCA方法不仅可以对地质构造复杂的地震数据进行高效的插值重建,而且可以很好的消除空间假频.

     

双稀疏字典和FISTA的地震数据去噪

地震数据的随机噪声去除是地震数据处理中的一项重要步骤,双稀疏字典提供了两层稀疏模型,比单层稀疏模型可以更好地去除噪声.该方法首先利用contourlet变换对地震数据进行稀疏表示,然后在contourlet域中使用快速迭代收缩阈值算法（fast iterative shrinkage-thresholding algorithm,FISTA）对初始字典系数进行更新,接着采用数据驱动紧标架（data-driven tight frame,DDTF）在contourlet域中得到DDTF字典并通过FISTA得到更新后的字典系数,最后通过DDTF字典和更新后的字典系数获得新的contourlet系数,并对新的contourlet系数进行硬阈值和contourlet反变换得到去噪后的数据.通过模拟数据和实际数据的实验证明：与固定基变换去噪方法相比,该方法可以自适应地对地震数据进行稀疏表示,在地震数据较为复杂时得到更高的信噪比;与字典学习去噪方法相比,该方法不仅拥有较快的去噪速度,而且克服了字典学习因为缺少先验约束造成瑕疵的缺点.

     

丘陵地带地震资料随机噪声压制新技术:高阶加权阈值函数的Shearlet变换

随着山地和丘陵地震勘探环境的复杂化,传统的消噪方法已经难以有效地压制地震记录中的随机噪声.Shearlet变换是一种新的多尺度多方向的时频分析方法,具有良好的稀疏表示特性,并且在每个尺度进行方向分解,非常适合用于地震信号随机噪声的压制.但是传统的Shearlet变换去噪方法采用的是硬阈值,在抑制随机噪声的同时也消除了很多有效信号,使得去噪之后的地震资料出现虚假的同相轴,为了解决这一问题我们提出高阶加权阈值函数.高阶加权阈值函数不但整体上连续性较好,而且克服了硬阈值函数存在剧烈的变化的缺点以及软阈值在处理较大Shearlet系数总存在恒定偏差的问题,同时保留了传统的软硬阈值函数的优点.实验结果表明这种基于高阶加权阈值函数的Shearlet变换去噪的方法,可以有效的消除模拟地震信号和实际丘陵地带地震信号中的随机噪声,同时很好的保留有效信号的幅度.

     

基于生成对抗网络的塔里木深层超深层叠前地震子波提取

深层超深层叠前地震子波的准确提取可有效提高全波形反演及叠前偏移成像等方法的准确性,对储层预测和油气分析具有重要意义.针对塔里木低信噪比叠前地震记录,本文提出一种基于生成对抗网络（GAN）的塔里木深层超深层叠前地震子波提取方法.结合柯东地区地质构造特点及噪声特性,采用有限差分法正演叠前地震记录,构建了塔里木柯东地区叠前地震记录及对应子波的训练样本集,利用地震子波、反射系数及噪声特征的区别训练生成对抗网络,实现了从叠前地震记录中有效分离地震子波.合成数据和实际地震资料处理结果表明,本文所提方法对塔里木深层超深层叠前地震子波提取具有有效性和较高的准确性,较传统方法更具优越性.

     

A new approach for high fidelity seismic data recovery by fractal interpolation

Recovering accurate data is important for both earthquake and exploration seismology studies, when data are sparsely sampled or partially missing. We present a method that allows for precise and accurate recovery of seismic data using a localized fractal recovery method. This method requires that the data are selfsimilar on local and global spatial scales. We present examples that show that the intrinsic structure associated with seismic data can be easily and accurately recovered by using this approach. This result, in turn, indicates that seismic data are indeed self-similar on local and global scales. This method is applicable not only for seismic studies, but also for any field studies that require accurate recovery of data from sparsely sampled datasets with partially missing data. Our ability to recover the missing data with high fidelity and accuracy will qualitatively improve the images of seismic tomography.     

预测多次波的逆散射级数方法与SRME方法及比较

逆散射级数方法和SRME方法是典型的两类基于波动方程的无需地下结构信息的自由表面多次波预测方法,本文详细讨论了这两类方法的基本思想、方法原理和实现思路,并对2维逆散射级数方法进行了降维简化,推导给出了1.5维逆散射级数衰减自由表面多次波的方法,减少了计算量,并降低了对3维规则观测系统的要求.利用经典SMAART模型数据测试比较了2维逆散射级数方法、1.5维逆散射级数方法及2维SRME方法预测自由表面多次波的效果及优缺点,比较了三者在方法实现、并行计算效率及对数据要求上的异同,并结合分析比较结果,给出了实际应用中方法选择的建议.     

基于波动方程的目标导向观测系统设计方法研究

本文提出了复杂构造地区的目标导向观测系统的设计方法.使用波动方程正演模拟来指导并在二维声波方程的一阶速度-应力方程中应用交错网格有限差分法实现.使用了四阶精度的差分算子和完全匹配层吸收边界条件.通过分析理论模型的模拟结果,展示了如何将地面地震响应与地下目标构造匹配.通过分析桥口地区实际地质模型的模拟结果,指出波动方程正演模拟在小断块、小背斜增生的复杂地区中相对于传统方法更精确,图像更清晰,更利于分析和指导观测系统设计.跟踪目标区的反射同相轴并得到其到达地面的接收范围来约束炮检距的范围,通过比较不同道距的模拟记录并综合考虑成本和任务目标来选取最佳道距.最终获得了实际效果达到设计要求的观测系统参数.     

小波域叠前地震资料衰减参数定性估计方法

叠前地震资料中,高频分量和低频分量随传播距离的衰减特性不同.本文给出了一种在小波域定性估计叠前地震资料衰减参数的方法.该方法利用连续小波变换提取共反射点道集的高、低频分量,以低频分量和高频分量之差定性反映地震波的衰减.通过累加不同偏移距的衰减,提高了估计的稳定性;采用幅度归一化方法,降低了信号幅值对衰减参数估计的影响.将本文提出方法与常用的基于叠后地震资料衰减估计方法用于某油田的地震资料处理,结果表明,本文方法得到的衰减估计结果能够更好地反映油气的空间展布.     

郯庐断裂带合肥段五河—合肥断裂构造特征

五河—合肥断裂是郯庐断裂带的西边界断裂,该断裂穿过合肥市城区,是1条规模较大、切割较深的隐伏活动断裂.为了研究该断裂的浅部结构特征、空间展布以及断裂活动性,我们利用2015年在合肥盆地完成的深地震反射剖面数据,采用初至波层析成像方法得到了郯庐断裂带合肥段的浅层P波速度结构和构造形态;考虑到仅根据速度结构剖面还难以确定断裂的准确位置、断层上断点埋深、断层的近地表构造组合样式等特征,研究中跨五河—合肥断裂还完成了2条高分辨率的浅层地震反射剖面.研究结果表明：郯庐断裂带合肥段是一个由多条主干断裂构成的复杂构造带,近地表速度结构表现为凹隆相间的构造特征,且沉积盖层厚度明显受到郯庐断裂带分支断裂的影响和控制.五河—合肥断裂在P波速度结构剖面表现为高速和低速区的分界,对断裂两侧的地层沉积具有重要的控制作用,该断裂向下错断了盆地基底,向上错断了埋深21~35 m的中更新统下部地层,其最新活动时代为中更新世早期.研究结果不仅为进一步认识五河—合肥断裂浅部构造形态提供了地震学依据,还可为该区断裂两侧的城镇规划和建设中避让活动断层提供基础资料.

     

维波场延拓复杂表层模型校正方法及应用

如何正确地消除复杂地表对地震波场的影响，提高地下构造成像的质量一直是中国西部复杂地区地震勘探中存在的难题.本文在三维复杂表层速度模型层析反演\[1\]的基础上，研究了关于复杂地表的静校正问题，提出用三维波动方程在炮集上对地震波场进行正、反向延拓，消除复杂地表对波场的影响，实现三维复杂表层模型校正.理论和实际应用证明，该方法已超越常规静校正的含义，属时变校正方法.用本方法处理复杂地表问题，不但能消除表层对不同深度反射波产生的不同时差影响，提高叠加剖面质量，而且能使校正后的地震波场保持波动特征不发生畸变，可为建立正确的深层速度模型和波动方程叠前深度偏移奠定良好的基础.