白云凹陷陆架—陆坡区多次波压制技术

白云凹陷位于中国南海北部陆架—陆坡区,由于水深变化剧烈,目标区地震资料中发育的表面多次波周期变化大,较难压制。传统的预测反褶积、SRME、Radon技术及其组合,压制这种水深变化剧烈的多次波效果不佳。为此,对近年发展的一种浅水表面多次波压制技术即确定性水层多次波压制技术DWD进行了研究,提出将其与SRME和Radon传统方法联合的多次波压制思路。DWD方法是先从水速动校后的最小炮检距的自相关剖面上拾取水底双程反射时间,然后在近道外推后将t-x域地震数据变换到τ-p域,再在τ-p域通过波场外推预测出多次波模型,最后与原始数据进行自适应匹配相减压制多次波。处理结果表明,DWD+SRME+Radon方法组合能有效地压制白云凹陷目标区多次波。与原处理剖面相比,此次处理后的成果剖面品质明显提高,降低了解释风险;同时,用此次处理得到的偏移道集进行叠前波阻抗反演,结果较好地刻画了目标砂体的特征和边界。     

渤海海域水平拖缆数据宽频处理关键技术

近年来海上浅水区宽频地震勘探技术进展缓慢,尤其像渤海这样平均水深不超过30m的浅水海域,宽频采集难以实施,只能通过常规采集、宽频处理才有可能实现这一目的,其中鬼波及浅水多次波是制约浅水资料宽频处理的两个主要因素。针对浅水海域水平拖缆鬼波时空变剧烈难以压制的问题,研究了τ-p域自适应参数估计的水平缆鬼波压制技术;针对渤海海域海底模型难以建立、浅水多次波压制困难的问题,研究了τ-p域自适应海底模型估计的浅水多次波压制技术。渤海M工区实际数据的应用结果表明,这两项关键技术能有效拓宽地震资料的频带宽度,实现渤海海域常规水平拖缆地震资料的宽频处理。     

多方位Walk-away VSP处理方法

从观测系统的定义、抽道集、速度分析、波场分离、纵波和转换波探度域成像等方面对多方位Walk-away VSP资料处理方法进行了讨论。首先利用数据分选技术从三维数据体中抽取不同方位的Walk-away VSP数据；然后通过综合三分量定位处理和中值滤波、f-k波波、τ-p变换、偏振分析等波场分离方法，将上行反射纵波（P-P波）和转换波（P-SV波）分离出来；最后根据速度分析结果建立速度模型，进行尝试域偏移成像，得到纵波和转换波的成像剖面。在垦71井区利用多方位Walk-away VSP资料处理方位法获得了很好的成像效果，偏移剖面与高精度三维地震剖面的品质相当。     

基于稀疏表示的非零井源距VSP波场分离方法

随着垂直地震剖面(VSP)技术的不断发展,高保真的波场分离成为VSP资料处理与应用的关键之一。引入欠定盲源分离思想,利用VSP波场的极化特征及多分量资料在时频域的稀疏性,提出一种基于稀疏表示的波场分离方法。该方法分为两步:第一步将地震信号变换到时频域,并估计极化矩阵;第二步在时频域分离各单一波场,再反变换回时间域。在第二步波场分离过程中,针对时频域波场稀疏性不足的情况,提出一种改进的最短路径分解法,在一定程度上解决了时频域局部信号混叠的问题。模型试算与实际资料处理结果表明,与传统方法相比该方法能够更准确地分离上、下行P-P波、P-SV波,分离后的波场保幅性较好,不存在混波、假频以及边界问题。     

多波多分量地震波场数值模拟及分析

以多波多分量地表资料处理和解释为目的，利用波动方程数值模拟方法对多波多分量地震波场进行了分析和研究。通过单界面和双界面模型正演，对反射纵波（PP波）和转换横波（P-SV波）的识别及波场响应特征进行了研究：①P-SV波速度低，频率低，能量随偏移距的增加而增加，零偏移距处能量为零；②界面反射系数为正时PSV波与PP波极性相反，界面反射系数为负时P-SV波与PP波极性一致；③Z分量和X分量地震记录都是PP波与PSV波的混合信息；④X分量的PP波和PSV波都是由两个极性相反的分支组成的。通过多界面模型正演，分析了转换波勘探的多解性，即地质上的同一个岩性界面有可能对应地震剖面上的两个甚至更多的同相轴。通过理论、模型和实际资料分析，探讨了多波多分量勘探中水平分量旋转处理存在的问题。即通过水平分量旋转处理获得的三分量记录仍然包含了全波场信息，指出通过极化分析，进行三分量同时旋转，可以实现纵波波场和横波波场的完全分离。最后讨论了PP波和PSV波的分辨率，认为P-SV波的纵、横向分辨率均低于PP波。     

λ-f域三维抛物Radon变换多次波压制方法

通常情况下,应用Radon变换压制多次波,一般是对正变换的模型空间进行处理,再反变换回时空域,而Radon域能量团的收敛性直接影响多次波去除效果,且求解过程中需反复进行矩阵的逆运算,计算成本偏高。为此,基于二维抛物Radon变换,借鉴Abbad等的方法,首次提出λ-f(λ为曲率q与频率f的乘积)域三维抛物Radon变换多次波压制方法,将常规三维抛物Radon变换的f-q_x-q_y(q_x、q_y分别为横、纵向离散曲率)域转换到一个全新的λ_x-λ_y-f(λ_x=q_xf、λ_y=q_yf)域,继承了二维λ-f域Radon变换思路。通过引入新的变量λ_x和λ_y,消除常规Radon算子对频率的依赖,减少了矩阵运算次数,显著提高了计算效率;根据一次波和多次波能量在λ-f域空间分布特征,设计三维椎体滤波器,更有效地分离一次波和多次波,降低空间截断效应引起的误差。理论模型及实际数据测试表明：1所提方法降低了空间截断效应的影响,消除了变换算子对频率的依赖,有效减少了矩阵求逆运算次数,较常规三维Radon变换方法的计算效率提高了约8倍以上;2所提方法仍然存在Radon类方法的弊端,即当一次波和多次波速度差异较小时,在近炮检距位置无法得到较理想的多次波压制效果。     

τ—P偏振滤波地震波场分离

视速度和偏振是地震波的两个重要特性。综合考虑这两个特性，可以对多分量地震波场有效地实现波场主能量分离。为改善和加强波场分离的效果，根据质点运动的直线性。在τ－Ｐ域设计一个二维非线性偏振滤波器，作用于主能量波场，便可压制和消除其它类型的波和澡声，使主能量波场更加突出。在频率－空间域实现τ－Ｐ变换，不需要对原始记录进行重量采样和插值，因而此算法快而精确；当参数满足一定要求时，还能有效地克服假频和端点效     

三维三分量Walk-away VSP处理方法及效果分析

三维三分量Walk-away VSP数据处理是一项系统工程,由于其观测系统的特殊性,对数据处理既不能借用常规的零与非零井源距多分量VSP处理方法,也不能采用目前地面地震多分量数据的处理方法。为此,针对三维三分量Walk-away VSP数据的特点,在观测系统定义、速度分析、波场分离、纵波和转换波的一维深度域成像等多方面进行了研究,初步形成了一套处理方法。通过对垦71井区VSP数据的处理,利用三分量处理和线性滤波、F-K滤波及基于τ-P变换和偏振分析的联合波场分离方法,有效地分离出上行纵波和转换波;根据速度分析和一维速度模型,在深度域对上行纵波和转换波进行成像,分别得到两个成像剖面。其中转换波成像显示出更高的分辨率和信噪比,具有很好的应用前景。     

三维三分量Walk-away VSP处理方法及效果分析

针对三维三分量Walk-awayVSP数据的特点,从观测系统定义、数据提取方法、抽道集、波场分离、纵波和转换波的一维深度域成像等方面进行了研究,初步形成了一套间接的三维VSP处理方法。在对垦71井区VSP数据的处理中,利用三分量处理和线性滤波、F—K滤波及基于τ—P变换和偏振分析的联合波场分离方法有效地分离出上行纵波和转换波;根据速度分析和一维速度模型在深度域对上行纵波和转换波进行成像,分别得到了两个成像剖面,其中转换波成像具有更高的分辨率和信噪比,显示出其良好的应用前景。     

基于Low-rank一步法波场延拓的黏声各向异性介质纯qP波正演模拟

各向异性介质纯qP波正演模拟及逆时偏移近年受到广泛关注,但它虽考虑了地下介质的各向异性特征,却忽略了黏滞性特征,使得最终偏移结果中噪声增加、分辨率降低。常规拟声波方程存在伪横波干扰、受模型参数限制（ε≥δ）、传播不稳定等因素影响,极大地限制了其应用。为此,引入一步法波场延拓方法,推导了黏声介质方程在空间—波数域的表达形式;结合空间—波数域各向异性介质延拓算子,构建一种适用于黏声各向异性介质的空间—波数域纯qP波波场延拓算子;引入Low-rank分解算法,实现基于Low-rank一步法波场延拓的黏声各向异性介质纯qP波正演模拟。数值模拟结果表明：①地震波场能同时表现出各向异性特征和黏滞性特征,更符合实际地下介质情况;②该方法克服了拟声波方程的局限性,消除了伪横波干扰,不受模型参数限制且地震波场能稳定传播;③在适当增大时间步长情形下无数值频散现象,所提算法能同时兼顾计算效率和计算精度,是一种稳定、高效的正演模拟方法,为基于Q补偿的各向异性介质逆时偏移提供了理论依据。     

矢量Radon变换波场分离方法研究现状与展望

保持地震波场矢量特征的波场分离是多分量地震数据处理的核心问题之一。为提高多分量波场分离精度,基于传统Radon变换的优势,探讨了矢量Radon变换保持波场特征的可行性。首先,在分析传统Radon变换基础上,阐述了基于组稀疏的矢量Radon变换实现保持矢量波场振幅特性的纵、横波分离方法,并将其应用于模型数据的纵、横波分离;其次,介绍了复矢量Radon变换压制多分量地震数据面波干扰的处理方法,并展望了复矢量Radon变换在多分量地震数据处理中纵、横波分离和快、慢横波分离的应用前景。模型数据应用结果表明,基于组稀疏的矢量Radon变换波场分离方法能够保持多分量地震数据的矢量特性,并能获得好于非矢量Radon变换的纵、横波分离结果。     

VSP三分量资料处理

在 VSP 观测中采用三分量检波器接收,每个分量所记录的信息都是空间波场振动在该分量上作矢量投影的合振动,因此,VSP 资料处理的关键,就是进行波场分离。为了获得好的波场分离效果,必须确定两水平分量的定向和求取偏振角。两水平分量的定向可以借助于直达 P 波的偏振方向在水平面内的投影加以确定。偏振角可通过矢端曲线和瞬时相位法求取。两水平分量定向后,上、下行 SH 波的水平分量被分离到 H_T(t)分量中,上、下行 P 波和转换 SV波的水平分量被分离到 H_S(t)中。将 Z(t)分量和 H_S(t)分量按在垂直平面内估算的下行 P 波偏振角进行坐标转换,即可实现下行波场的初步分离。然后,再联合使用波场分离和反转投影方法对上行波场进行分离。但是,反转投影只能实现上行 P 波和 SV 波的主能量分离,不能完成全波场的分离。文中展示了在PDP-11/24计算机上使用人机交互三分量处理软件的处理效果。     

基于直达波偏振分析的三分量检波器定向方法

本文提出一种基于三分量地震记录直达波偏振分析的三分量检波器定向方法。该方法首先利用三分量地震记录的直达波振幅构造协方差矩阵,然后通过求解该协方差矩阵的主特征向量得到直达波的偏振方向,再根据炮点和接收点的位置计算直达波的传播方向,在此基础上,利用直达波的偏振方向和传播方向的差异最终确定三分量检波器与设计方向的水平和垂直偏角。利用估算的三分量检波器的水平和垂直偏角,通过三分量旋转技术实现三分量地震记录的重定向,进而消除由于三分量检波器放置方向的误差对波场带来的影响。模型和实际数据试算结果表明,该方法具有很高的计算精度和很好的应用效果。     

勘探地震中的六分量观测

地震旋转分量在天然地震和工程地震等领域都有广泛的应用。在回顾地震旋转运动的基础上,首先介绍了六分量地震数值模拟研究,论证了不同震源类型与地震波型在六分量波场上存在能量与极性的差异。随后,回顾了利用模拟数据与实际观测数据验证差分法与行波法两种旋转分量间接换算方法的可行性与精度。通过数值模拟,分析了勘探地震中广泛应用的面波,在六分量记录上呈现的不同或相似的特征,以及旋转分量记录波场信息具有的不同于平动分量的频散特征。近年来,不断发展的地震旋转分量研究,对波型识别、波场分离、横波成像、微震的高精度预测、探索海底电磁和地震观测的关联等领域均具有重要意义与深入研究价值。     

采用τ-q变换法进行纵、横波波场分离

多波多分量地震勘探与传统的单一纵波地震勘探相比具有许多优势，同时采集纵、横波资料可以分别形成叠加剖面，有利于构造成像，更可贵的是，它可以用来提取更多的约束岩性，储层，裂隙分布的地球物理参数以及应力场信息，波场分离是多波多分量地震资料处理的重要环节之一，采用常规的τ－q变换法进行纵、横波波场分离时，由于双曲线经τ－q变换后为椭圆（其能量分布是渐变的）、实际的分离效果产不好，因此，根据江家屯地区多波多分量地震资料特点，提出了τ－q变换法，该方法将来自不同层位的反射波经τ－q变换后聚焦在不同的点上，更有利于波场分离，因此切实可行。     

变速视慢度波场分离

本文根据地震波视速度和偏振特性,导出波场分离矩阵是视慢度p和波速V_p、V_s的函数。由于波速是空间坐标函数,因此,在带有空间坐标的数域与带有视慢度参数户的数域(如时—空域与τ-p域、频—空域与f-k域)相互变換,可实现变速视慢度波场分离。文中还给出了使用τ-p变换法、f-k变换法和DRT法实现变速视慢度波场分离的三种具体方法。采用变速视慢度波场分离方法时,速度函数可采用零井源距直达波的速度或声波测井资料换算的速度。在进行纵、横波分离的同时,可对视慢度p限定不同范围,实现上、下行波的波场分离及f-k滤波、τ-p滤波等。     

一种简便的海洋P-SV转换波叠前偏移成像处理方法

在莺歌海盆地,地震模糊区的成像一直是勘探工作中的一大难题,1998年在该地区采集了多波地震资料,通过对P-SV波的处理,在模糊带成像方面取得了较好的处理效果[1]。在对常规P-SV波处理中的偏移技术等存在的特殊问题进行分析的基础上,提出了一种P-SV波叠前偏移处理的方法,该方法将输入道分选成共转换散射点道集(CCSP),在CCSP道集中,双程旅行时和等效炮检距之间满足双曲线关系,故而使得常规处理P-P波的Kirchhof叠前偏移在此也可以运用,至于速度分析它也可以运用P-P波的速度分析方式来进行,因此采用该方法简化了转换波的处理,同时也提高了P-SV波的成像效果。     

基于共转换点道集的二参数转换波速度分析及动校正方法

在不引入介质的各向异性参数并认为共反射点处同时产生波场的反射与转换的状况下,推导P波与P-SV波之间的运动学关系,进而得到转换波时-距曲线方程--转换横波叠加速度分析公式。借助地震全波场数值模拟技术,进行多波波场响应特征的分析研究及多波多分量地震模拟记录的生成,使用计算机C语言编程技术和SU编程技术实现了该方法的软件化。通过模拟数据和实际资料的验证与测试,结果表明所提出的转换波二参数方程精度高、经济实用、效率高,解决了转换波地震资料处理中叠加速度场建立和叠加成像的问题,对多波用于油气勘探具有指导作用。     

双检数据上下行波场分离技术研究进展

随着海洋油气勘探开发技术的发展，水、陆双检海底接收的应用日趋广泛，上下行波场分离作为双检数据处理的关键技术，决定了资料处理品质及应用效果。在双检数据上下行波场分离技术的发展过程中，国内外学者提出了许多切实可行的技术方案。本文通过系统调研国内外相关文献，对上下行波场分离方法和技术进行了归纳和总结。首先基于波动理论，把双检数据与上下行波场联系起来，建立了双检数据标定和上下行波场分离的理论基础：水、陆检数据接收的波场分别为压力波场和质点垂直速度波场，可以分解为上行压力波场和下行压力波场及上行质点垂直速度波场和下行质点垂直速度波场；下行压力波场与下行质点速度波场振幅成正比，且极性相同，上行压力波场与上行质点速度波场振幅成正比，且极性相反。然后针对海底电缆/海底节点（OBC/OBN）水、陆检数据，建立了基于消除海水鸣震的双检数据处理技术。最后，总结出7种双检数据上下行波分离方法，包括：常数标定因子分离、频率—波数域分离、镜像分离、去虚反射分离、最优去虚反射分离、检波器脉冲响应分离和τ-p域分离等方法。双检数据上下行波场分离之前，首先要进行检波器脉冲响应校正；然后进行双检数据匹配和标定处理等，使水检数据中上下行压力波场与陆检数据中上下行质点垂直速度波场振幅和频率趋于一致，且由于下行压力和质点垂直速度波场极性相同、上行压力和质点垂直速度波场极性相反，陷波问题得到补偿，检波点端虚反射得到合理压制；最后实现最佳的上下行波场分离处理。     

多波多分量地震勘探在松辽盆地的初步应用

根据松辽盆地陆相薄互层岩性勘探所遇到的困难，近年来在油田上开展了多波多分量地震勘探试验，探讨了一套多波多分量地震资料采集，处理和解释方法。采集时考虑纵，横波入射临界角，求出最大，最小炮检距，决定野外观测系统；处理时，采用τ－ｐ域各向异性介质的平面波理论与波的偏振投影，实现Ｐ波，ＳＶ波分离，对转换波（ＳＶ）处理考虑到射线路径的不对称性，故用共转换点（ＣＣＰ）道集求其时距曲线，从而得到最终转换波叠加剖