深反射地震剖面数据处理的主要技术方法

深反射地震剖面技术以其探测精度高的优势被作为岩石圈精细结构研究的先锋技术,并取得了突飞猛进的发展.但是,因为探测剖面长,地表地形复杂等原因,对深反射地震数据处理的改进一直是地球物理学家关注的关键内容.近十几年来,随着探测工作量的增加和处理技术的发展,涌现出不少新思路用以深地震反射数据处理,特别在静校正、速度分析、去噪和一致性处理等重要环节更是优化了原处理技术,提升了深反射地震数据处理的水平.本文阐述了深反射地震数据的基本特点,结合实例对当前深反射地震数据处理中使用的主要技术进行了系统的总结,对比了一些方法的优缺点,并对深反射地震数据处理技术的未来发展方向进行了展望.     

基于深度神经网络的地震强反射剥离方法

在储层预测工作中,储层弱反射信号淹没在强反射信号之中的情况非常常见,这不利于精确识别和描述储层结构.本文提出了一种基于深度神经网络的强反射剥离方法,用于辅助储层弱反射信号的检测工作.该方法在卷积模型的框架下将强反射预测问题分解为地震子波预测与强反射预测两个子优化问题,并采用AIDNN与U-Net两个深度神经网络分别求解.通过训练直接得到地震数据与强反射之间的映射关系,避免了经验性调参过程,计算速度快,适用于海量地震数据处理.模型数据和实际资料试算结果表明,本文方法能够预测并剥离地震数据中的强反射且保真保幅性好;在该方法的基础上进行的储层砂体展布预测工作取得了良好效果.     

斜缆数据频率域高分辨率Radon变换与虚反射消除方法（英文）

检波器虚反射是斜缆地震数据采集的一种主要干扰波。检波器虚反射与一次波形成的陷波点随缆深而变化,引起频率成份的缺失,严重影响了地震资料处理与成像效果。研究了一种适用于斜缆地震数据的频率域高分辨率Radon变换方法,拟合出地震记录中的虚反射,将其从原始记录中减去实现了虚反射的消除。通过高分辨率线性Radon变换的基本原理分析,结合斜缆数据采集中检波器虚反射与一次波的传播特性,得到了高保真的频率域斜缆地震数据Radon正变换的表达式,根据虚反射传播路径建立了虚反射Radon反变换公式。将该方法应用于模型数据和实际资料,准确地拟合了斜缆数据中的虚反射。利用商业软件的匹配相减模块,从原始数据中减去检波器虚反射,有效地消除了虚反射对一次波的影响,陷频点得到了补偿,拓宽了地震资料的有效频带宽度。     

基于S变换的软阈值滤波在深地震反射数据处理中的应用

深地震反射原始单炮数据是非平稳的弱能量反射信号,信噪比较低.如何提高信噪比一直是深地震反射数据前处理中的一大难题.S变换是一种适用于分析非平稳信号的时频变换方法.同其他分析时变信号的方法相比,S变换的基本小波不必满足小波在时间域均值为零的容许性条件,它的时频分辨率与分析信号的频率有关,且其在时间域的积分可以得到傅里叶频谱,其反变换也简单.因此,S变换容易表示深地震反射信号复杂的时频特性.本文在S变换的基础上,利用软阈值滤波方法对深地震反射数据进行处理,实验结果表明,该方法有效地提高了信噪比,压制了有效频带范围内的混频干扰,突出了弱反射信号,使得波组信息更加丰富,有利于连续追踪有效反射波组和识别薄地层,特别是提高了深部Moho界面反射层位的分辨率,为深地震反射剖面后续处理和准确解释奠定了基础.     

利用基于子波估计的频谱校正方法提高Q值估计精度(英文)

用Q值刻画的地震衰减在地震信号处理和解释中具有很广泛的应用。利用反射地震资料进行Q值估计需要解决地震子波和反射系数序列耦合的问题。从反射地震资料中去除反射系数序列的影响,这个过程称为频谱校正。本文提出了一种基于子波估计的求取Q值的方法,进而设计了一个反Q滤波器。该方法利用反射地震资料的高阶统计量进行子波估计,并利用所估计子波实现频谱校正。我们利用合成数据实验给出了质心频移法与频谱比法这两种常用的Q值估计方法在不同参数设置下的性能。人工合成数据和实际数据处理表明,利用本文提出的方法进行频谱校正后,可以得到可靠的Q值估计。经过反Q滤波,地震数据的高频部分得到了有效地恢复。     

利用深反射地震数据构建的多阶面波频散曲线反演近地表横波速度结构——以跨班公湖—怒江缝合带深反射地震资料为例

深反射地震剖面法为了获取深部结构特征常常采取大的偏移距采集数据.目前公开发表的相关资料中,鲜有利用深反射地震炮集数据获取近地表的结构特征.为此,本文通过正演测试了相关数据处理流程,即利用有限差分正演了起伏地表模型的大偏移距地震单炮弹性波场特征,通过共检波点域面波信号F-K频谱叠加构建新方法,从深反射地震数据集中提取了高品质的多阶面波频散曲线,再利用多阶面波联合反演获得了近地表的结构特征.在前述正演流程基础上,利用跨越班公湖—怒江缝合带的SinoProbe深反射地震剖面中的实际炮集数据,求取了基阶和一阶瑞利波频散曲线,联合反演后得到近地表横波速度结构.该结果与初至波走时反演获取的纵波速度结构具有较好的一致性,且在近地表的浅层分辨率较纵波速度结构特征更高,而更与已有地质认识相吻合.本文提供的相关数据处理流程表明利用深反射地震炮集数据,也能够获取近地表浅层的横波速度结构.     

再论地震数据偏移成像

利用地震波正向传播方程对属于波形线性反演问题近似求解方法的地震数据偏移成像进行重新推导,得到了适合散射地震数据的散射偏移成像方法和适合反射地震数据的反射偏移成像方法.以地震波传播的散射理论为出发点,首先根据描述一次散射波正向传播的线性方程研究建立散射地震数据的偏移成像方法理论;利用高频近似对产生散射波场的地下速度扰动函数的空间变化进行近似,推导出地下反射率函数,再由散射波传播方程推导出基于反射率函数的反射波传播方程,然后根据描述一次反射波正向传播的线性方程研究建立反射地震数据的偏移成像方法理论.本文指出和修正了Claerbout偏移成像方法中的不足,提出的地震数据偏移成像方法是对当前偏移成像方法理论的完善,使反射地震数据偏移成像具有了更坚实的数学物理理论基础,得到的偏移成像结果相位正确、位置准确、分辨率提高.     

海上地震资料子波提取鬼波压制技术

海洋地震勘探由于受到海面强反射界面的影响,记录数据中存在鬼波,致使数据频谱中产生陷波效应,限制数据的频带宽度,影响数据的分辨率,干扰地层反射的识别.近年来发展的鬼波压制方法,能够较好的压制地震数据中的鬼波,然而现有鬼波压制技术容易产生延续相位,且计算效率较低.为进一步提高鬼波压制效果和工作效率,本文对海上地震资料特性进行分析,研究认为,海上地震资料中的海底反射信号真实记录了海底有效反射、震源鬼波和检波点鬼波的形态,可以利用提取的海底反射子波对数据进行确定性子波反褶积处理,能够有效的压制震源鬼波和检波点鬼波.模型海上实际数据测试结果表明,该方法具有鬼波压制效果好、计算效率高的优势,处理之后的数据频带得到有效拓宽,低频和高频信息均得到加强.     

现代功率谱估算在反射地震剖面随机介质参数估算中的应用

由于地壳深部物质呈现出多尺度的非均质性,导致了地震剖面上反射能量弱、同相轴连续性差、多成带状分布.传统的反射地震解释方法难以有效应用到深反射地震资料解释中.统计分析的方法为深反射地震资解释提供了新的工具,地下介质这种随机分布的非均质性可以借助随机介质模型来进行描述.深反射地震剖面继承了地下介质的非均质性,因此可以通过从反射地震剖面估算随机介质参数来研究地下介质非均质性的空间分布.由于地震剖面长度有限,随机介质参数估算中自相关函数计算误差较大,导致了估算的不准确.现代功率谱估算克服了短数据自相关函数难以准确求取的问题,提高了参数估算准确性.模型计算和实际数据测算都证明了该方法能有效提高自相关函数计算精度和随机介质参数剖面的分辨率.     

用深反射大炮对大巴山-秦岭结合部位的地壳下部和上地幔成像

根据近垂直反射原理对秦岭和大巴山结合部位的深地震反射测线上的5个大药量的反射大炮进行静校正、去噪、速度求取、CDP分选等处理获得沿深地震反射测线的大炮单次覆盖剖面.单次剖面和原始单炮数据结果初步显示:深地震反射测线最南端TWT 16 s处存在地幔反射,该地幔反射界面向北可以追踪到凤凰山下TWT 18.5 s.地幔反射在原始大炮数据表现为信噪比较高的连续震相.莫霍面反射在凤凰山南10km,TWT在16 s,从该处向南莫霍面震相难以识别,不能追踪;从该处向北,莫霍面反射可连续追踪,界面变化平缓、南深北浅,在测线最北端莫霍面反射TWT 14.2 s.初步认为测线最南端的地幔反射为扬子板块向北俯冲证据.     

用波速层析成像方法提高储层横向预测精度

以数值模拟方法为基础，考察了反射地震层析成像技术对油气储层波速横向检测的效果，然后用渤海油田的地震资料进行了实用试验．对数值模拟结果和实际数据的试验结果进行了分析与评价，认为油气勘探中的反射地震资料具有地震射线近似垂直分布的特点，适合用反射地震层析成像技术进行目标层横向非均匀性检测．     

汶川地震发震构造带过WSFD-2钻井深反射地震数据3-D处理技术研究

2008年5月12日汶川发生里氏8.0级地震的发震机制与龙门山断裂带的构造特征紧密相关,应用反射地震探测方法精细探测发震构造内部的结构特征,对研究汶川地震的发震机理有重要意义.5.12地震后在地表错断最剧烈的区域之一虹口乡白庙村穿过发震构造和深钻科研井WSFD-1和WSFD-2布置了兼顾浅中深层信息的反射地震探测剖面,该区域断裂带硬岩出露,地形地貌复杂多变,起伏剧烈,只能采取弯线地震测线布设方式采集二维反射地震数据,而且弯曲度大,造成CDP点分散严重.区域内构造复杂,褶皱逆推构造发育,地层和构造倾角大,采用常规二维地震数据处理方法进行叠加处理时,易将不同地层的反射信息叠加到同一反射层.弯线叠加剖面上侧面波,混波干扰严重,会显著地降低原本低信噪比数据的分辨率,或剖面可靠性低,容易在二维剖面中造成解释陷阱.本文利用弯线地震采集的三维特性,用拟三维地震叠加技术处理汶川地震科学钻井附近横跨北川-映秀断裂带的二维弯线地震数据,弥补常规二维弯线地震数据处理技术的不足,将不同地层的反射信息分离归位到不同的三维叠加剖面上,解决二维弯线地震数据处理时混波干扰严重的难题.通过理论分析和实验选取合适的共中心点面元,获得了高分辨率的三维叠加剖面.相比于二维弯线叠加剖面,三维叠加剖面切片成像更真实,剔除了不同地层反射信息混叠的影响,能得到更准确的断点信息,并可以获得沿断层走向横向的信息,显著提高了构造解释的可靠性和精度.应用拟三维地震数据处理方法处理龙门山断裂带的二维弯线地震数据,获取高分辨率的构造信息,有利于断裂带内汶川地震发震机制和龙门山隆起机制的解释.     

地壳测深的地震数据的收集与处理方法的最新进展

地壳测深的地震数据的收集与处理方法的最新进展５．１宽角地震数据方法与结果５．１．１用宽角地震数据进行速度和界面的同时成像由于技术的进步，海洋宽角地震数据现在可以用类似于常规反射测量的炮点密度进行采集。数据的大量增加超过了现有的正向模拟方法的能力。英国...     

基于反射地震数据的时频域包络反演

包络信号含有丰富的低频分量,即使在地震数据缺失低频条件下,包络目标函数也能有效缓解全波形反演的周期跳跃现象.但是,当初始速度模型较为平滑时,观测数据中的反射地震事件在模拟数据中没有与之相对应的波形,导致包络反演初期无法很好地利用反射波信号进行速度建模.本文提出基于反射地震数据的时频域包络反演方法,通过结合反射波全波形反演理论,构建反射波时频域包络目标函数,来提高包络反演的速度建模精度.本文首先利用Gabor变换获取时频域地震数据,并提取振幅信息,即为时频域包络信号.然后,推导反射波时频域包络反演的伴随震源和梯度算子.Marmousi模型数据测试结果表明,基于反射地震数据的时频域包络反演方法可以为全波形反演提供一个较好的初始速度模型.     

基于反射地震数据的时频域包络反演

包络信号含有丰富的低频分量,即使在地震数据缺失低频条件下,包络目标函数也能有效缓解全波形反演的周期跳跃现象.但是,当初始速度模型较为平滑时,观测数据中的反射地震事件在模拟数据中没有与之相对应的波形,导致包络反演初期无法很好地利用反射波信号进行速度建模.本文提出基于反射地震数据的时频域包络反演方法,通过结合反射波全波形反演理论,构建反射波时频域包络目标函数,来提高包络反演的速度建模精度.本文首先利用Gabor变换获取时频域地震数据,并提取振幅信息,即为时频域包络信号.然后,推导反射波时频域包络反演的伴随震源和梯度算子.Marmousi模型数据测试结果表明,基于反射地震数据的时频域包络反演方法可以为全波形反演提供一个较好的初始速度模型.     

大陆及其边缘地震探测国际讨论会综述

1988年7月1—8日于澳大利亚堪培拉举行了大陆及其边缘地震探测国际讨论会,会上主要交流了世界各国近年来用地震反射方法探测岩石圈的工作经验及成果。全球3万公里的深反射地震剖面及三维人工地震观测,给出了大陆及其边缘地区的岩石圈细结构图像,为研究地壳演化及深部过程提供了重要依据。本文综合介绍了这次国际讨论会的主要内容,包括:全球深地震反射剖面的主要成果及展望、岩石圈细结构及反射性质的研究、人工地震的地壳三维成像,及数据搜集和处理方法的进展。     

海洋中尺度涡与内波的地震图像

海洋反射地震通常用于调查、研究海底地质构造,勘探油气与天然气水合物资源.近期研究表明多道反射地震方法也可以对水柱的热盐细结构成像.中尺度涡与内波是重要的物理海洋现象,但是常规的物理海洋调查是在间隔若干公里的离散测站上进行的,水平分辨率较低,因此对中尺度涡的结构与内波的横向分布了解较差.本文利用在大西洋东部、南海采集的地震数据给出了低频反射地震可以对中尺度涡与内波清晰成像的新的证据.反射地震方法较传统海洋观测手段,具有明显的优势,主要体现在高的水平分辨率和短时间内对整个海水剖面进行成像方面．从地震剖面上,能够清楚地观测到中尺度涡、内波造成的反射特征变化,从而有助于改进对能量在不同尺度的海水运动之间传递过程的认识.     

空间反射结构正则化多道反褶积（英文）

在高分辨率地震资料处理中通常采用稀疏脉冲反褶积方法来展宽地震数据的有效频带,增强地震数据刻画薄层结构的能力。但是,现有的稀疏脉冲反褶积方法采用逐道运行模式,未考虑地震信号在相邻道之间的空间约束关系。受高频噪声的影响,反褶积结果不仅具有较强的多解性,且降低了地震信号的空间连续性。为此,本文提出了一种空间反射结构正则化多道稀疏脉冲反褶积方法。该方法从地震数据本身提取空间反射结构表征算子,将表征算子引入到多道稀疏脉冲反演的正则化条件,增强了稀疏脉冲反褶积方法的抗噪性和稳定性。理论模型实验和实际资料处理证实了该方法的正确性和可行性。     

地震逆散射波场和算子的谱分解

本文对地震逆散射的研究,旨在于为抑制层间多次波和地震波场多重散射对一次反射干扰效应提供理论依据.这对薄互层地层滤波的高频恢复、保幅弹性反演、衍射地震勘探及海洋地震勘探中的干扰消除皆具重要意义.本文基于上下行波分解及弹性波互易定理,导出横向变速介质条件下线性预测算子的表达式和反射数据的广义谱分解方程. 文中先由上覆地层广义反射透射矩阵的元素定义线性预测算子,并将其表示成一系列单程波算子的线性组合,之后将横向变速介质条件下线性预测方程表达为反射数据与线性预测算子及其逆的乘积. 对该方程的求解可获得上覆地层的线性预测算子,从而可借以求出相应的反射透射算子. 本文先将水平层状介质条件下垂直入射的一维线性预测方程推广到斜入射的情况,以此为参照,导出横向非均匀介质条件下反射数据的地震逆散射广义谱分解方程.文中也揭示了单程波地震逆散射算子、反射透射算子的性态.本文还针对水平层状介质条件,给出斜入射的数值结果.     

南海东北部东沙海域海水层反射特征分析

本文基于对南海东北部东沙海域近期采集的多道反射地震资料进行重新处理获得新的地震海洋学数据,分析了该海域内孤立波/内孤立波包、沙丘上方和陡坎附近特殊反射结构特征,从而提供了新的海水层与海底相互作用依据.研究结果表明,除之前已发表文章中地震海洋学资料显示存在的第一模态内孤立波/波包和沙丘上方常见的反射样式-披毛状发射外,地震海洋学资料上还发现了第二模态内孤立波、陡坎上方的上抬型波动反射结构样式.在新的地震海洋学数据中,第一模态内孤立波振幅均小于50 m,宽度上都小于5 km,单个内孤立波的最大振幅约为45 m.内孤立波包的内孤立波振幅都相对较小,均小于40 m,并且与之前不同的是,彼此之间振幅相差不大,没有明显的排列规律.此次地震海洋学数据记录到的第二模态内孤立波,形态较为完整,上层和下层反射的振幅相差不大,在30 m左右;中间层大约在水深130 m位置处,垂向结构的整体大小大于200 m.沙丘上方反射结构普遍存在弱反射层,可能是湍流边界层,并且存在特殊反射样式-披毛状反射.但并不是沙丘上方都存在披毛状反射样式,本文分析它出现在地震海洋学资料上可能是受测线与沙丘走向之间夹角的影响.陡坎区域的水层反射结构则表现为上抬型波动,并常常伴随着同相轴连续性的变化.该波动的大小及反射同相轴的连续性可能取决于陡坎的高度/坡度及水层动力的强度,新数据中出现的一个上抬型波动,高度达20~30 m,它的附近水层还存在一个形态不完整的内孤立波.陡坎附近的水层反射也常常出现弱反射带和小的波动.