再论地震数据偏移成像

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

三维井间地震高斯射线束正演

使用高斯射线束方法对三维井间地质模型进行波场正演,解决三维井间地震复杂构造波场识别问题.正演模拟时将每一个检波点与炮点拟合出一个接收平面,使得每一个检波点处的能量都由其接收平面上检波点的有效半宽度范围内所有射线能量高斯加权.高斯束正演具有地震波动力学特征且能够有效解决复杂构造盲区问题,同时保持波场的相对振幅特征,因而适用于不同地质构造模型.比较三维井间地震逐段迭代正演方法,高斯束正演克服了复杂地质模型的盲区与动力学特征问题;并且高斯束正演时反射点的分布范围广,大大提高了反射波成像精度.应用三维井间地震高斯束正演方法对X地区进行波场正演,结果表明该方法可以有效分辨波场特征,为进一步地震数据处理提供准确依据.     

不同参数对地震干涉实现VSP向SWP波场重构的影响

地震干涉是近几年勘探地球物理领域研究的一个热点方向,它是波场重构的有力工具.然而,地震干涉往往引入虚假同相轴,影响波场重构的质量.为进一步分析虚假同相轴产生的原因并改善波场重构的质量,本文基于稳相分析,详细探讨子波主频、激发炮数、检波器埋深、检波器间距和地层倾角五个因素对波场重构的影响.模型结果表明,震检排列方式和地层倾角等因素通过改变有限震检范围内稳相点的位置来影响波场重构的质量.将垂直地震剖面VSP中的下行直达波与下行反射波进行地震干涉处理,可有效重构来自高陡反射面的下行反射波,将传统的VSP转化为单井地震剖面SWP.直接利用重构SWP波场进行成像,不仅扩大传统VSP的成像范围,而且避免常规勘探面临的静校正和近地表速度建模问题,为高陡构造成像提供一种新方法.     

三维斜井XSP-CDP叠加成像

为了解决三维斜井井间地震数据成像问题,借鉴二维地面地震弯线及宽线处理思路,将三维井间地震成像问题转换成二维成像问题进行研究.利用三维逐段迭代方法对三维斜井井间地质模型进行正演模拟,计算上、下行反射波走时及反射点位置,通过对比正演记录与实际资料的特征,调整地质模型参数获取准确速度场.然后对波场分离后的共炮集数据进行反射波归位,得到共反射点道集数据.最后根据反射点的分布选择一条反射点分布最集中的剖面线,以一定的宽度将剖面线左右的反射点包裹起来;根据叠加宽度划分成大小相等的矩形面元,将面元内的共反射点数据投影到剖面线上进行二维叠加成像.通过模型测试,三维XSP-CDP叠加方法能有效反演地下结构特征,指导三维井间地震采集观测系统的设计并解决不共面斜井叠加成像问题.应用该方法对M地区实际资料进行成像研究,结果表明该方法可有效反演地下三维地质结构特征,为进一步数据处理提供依据.     

共反射元弧叠加方法在火山岩成像中的应用

常规叠加方法假设叠加的地震道集的反射波来自地下同一个反射点,而地表观测到 的反射波场在空间上小于菲涅尔体范围是不可分辨的. 本文根据火山岩体特征,用元弧描述 火山岩的反射界面. 将来自菲涅尔体范围内的反射波沿着元弧界面求和,进而扩大叠加次数 ,达到增加反射信号能量提高地震反射波的叠加成像质量的目的. 文中给出了反射波走时参 数与地下元弧反射界面之间的联系和叠加实现的步骤. 计算结果表明,元弧叠加方法能够提 高剖面的信噪比,增加反射波同相轴的连续性.     

基于小振幅波方程的起伏海面构建及其对地震波场传播特性的影响

海面起伏状态对波场传播以及模型空间成像有着非常重要的影响,在常规海洋地震数据采集和处理中,都是将海面作为水平界面处理,实际海洋地震数据采集中,由于海水受风浪、洋流等的影响,海面呈现高低起伏状态,这样来自底层反射界面的反射波经过海水起伏界面反射后会引起波场的强烈散射,从而造成地震记录反射波同相轴抖动和畸变,最终影响成像效果.为了更加深入了解起伏海面存在对波场传播以及偏移成像的影响,本文首先基于流体力学小振幅波方程,从边界条件出发,给出包含海浪谱信息的动态边界条件,并在此基础上构建起伏海面速度模型;然后在构建的起伏海面速度模型上通过冻结不同时刻海面开展波场传播特征分析,研究起伏海面在不同冻结时刻下波场传播规律以及对反射地震记录和偏移成像的影响.研究结果表明,不管是对于简单的层状介质模型还是海底较为复杂的大水深模型,海面不同形态的起伏变化直接影响着波场记录以及波场快照形态,同时会对偏移成像结果造成较大影响,随着海水深度的增加,该影响对于简单模型影响逐渐变小,但对于复杂模型,影响依然很严重.     

高分辨率Radon变换及其在声反射成像测井中的应用

声波反射成像测井是利用阵列声波测井的反射波进行井旁构造成像的测井新方法,由于是从滑行波能量较强的阵列声波信号中来获得反射波,因此反射波的提取问题尤为重要,文中采用能使波场更加聚焦的高分辨率Radon变换来提取反射波.分析得出滑行波和反射波到时在阵列数据道情况下表现为线性或近似线性特征,因而采用线性Radon变换方法.首先分析了声反射成像测井不同数据道集的波场几何特征,得出滑行波和反射波在不同道集都表现有不同程度的视速度差异,为反射波提取提供了依据.传统的Radon变换采用的最小二乘法不能使同相轴在Radon域足够聚焦,结合最大熵原理和贝叶斯原理实现了线性高分辨率Radon变换,模型数据处理结果表明后者可以在较复杂条件下提取出期望的波场.采用高阶交错网格有限差分模拟得到软、硬地层声反射成像测井波场,处理结果表明,相对于最小二乘Radon变换,采用高分辨率Radon变换后可以使波场更加聚焦,能提取出高质量的反射波,并能适应不同性质地层,相比于目前应用较多的反射波提取方法,其结果更接近于理论波形.实际数据处理也表明了该方法的有效性和准确性.该方法在声反射成像测井的速度分析、数据道插值以及偶极反射波提取方面有一定潜力,需要进一步研究.     

粘弹各向异性反射体模型方位地震AVO模拟及分析（英文）

基于岩石物理模型和广义各向异性Zoeppritz方法在频率域计算裂缝型反射体模型反射波方位地震AVO响应。反射体模型为粘弹各向异性、有限厚度的地质体,其地震反射波形序列包含如下动力学信息,即分界面处介质的波阻抗和非弹性差异、反射体内部波的各向异性传播、在传播路径上的频散与衰减,以及来自顶底界面的反射波的调谐与干涉等。计算表明,速度频散和衰减增顶界面反射波大入射角反射时的振幅,而减弱底界面反射振幅。对于固定入射角的方位地震方位地震响应,PP波反射特征表现为随方位角的增加反射波形序列延续时间变长,而PSV和PSH转换类型反射波的方位各向异性变化特征稳定且受储层厚度影响较小,表现为PSV波反射振幅随方位角增加而增加,PSH波在0°和90°方位无反射能量,在45°方位反射振幅最强。     

煤层陷落柱散射波数值模拟与成像

煤层陷落柱是煤田勘探开发中常见的一种典型的非均匀地质体.由于来自陷落柱的反射信号少、反射能量弱,使得基于反射波原理的常规地震成像方法难以有效识别陷落柱.本文以散射波理论为基础,采用数值模拟方法,研究了陷落柱的散射波场特征,研究表明地面接收的波场中含有来自陷落柱陡倾角界面的散射波场.通过共散射点道集波场的模拟,可以清晰地识别散射波,获得地下散射点和非均匀地质体的信息,判断散射点的位置,从而勾画出不均匀地质体的形态.采用等效偏移距假设抽取共散射点道集,在此基础上进行叠前偏移,对陷落柱成像;模拟与实际数据成像结果对比表明此方法能够合理地提取散射点的散射波场信息,对陷落柱形态及内部结构准确成像,是一种有效的煤田陷落柱成像方法.     

伊宁凹陷煤系地层对深层地震反射波能量衰减机理分析

伊宁凹陷侏罗系以下油气勘探未获进展,其原因是侏罗系煤系地层与砂岩地层交互,煤层厚度大部分处在四分之一的地震波长,不仅层间多次波发育严重,而且目的层反射能量较弱.为此,本文基于薄煤系地层地震层间多次波反射和透射的频率特性,来定量分析煤系地层使深层地震反射能量弱的机理,并基于伊宁凹陷已钻探井及初勘地震剖面构制地质模型,采用双程波动方程照明分析方法,验证薄煤层内层间多次波对下伏地层反射能量的衰减,旨在为靶区资料处理阶段提高深层反射能量提供依据.     

τ-p变换在隧道反射地震超前预报波场分离中应用的数值模拟研究

本文基于二维有限元数值方法,模拟了复杂条件下隧道内的地震全波波场,分析了隧道地震超前预报中的反射记录中的多方向和多波反射波特征,证实了分离前方有效反射波和其它方向的干扰反射波,以及有效反射波的P波、S波分离在预报准确性中的重要性.基于多个数值模型,研究了τ-p变换进行波场分离的有效性和影响因素.研究结果表明:在隧道前方反射层倾向与隧道轴线交角变化较大范围内(90°～45°),τ-p变换都可以有效的分离多方向和多波反射事件,提取主要来自掌子面前方的有效反射波剖面.来自隧道项、底的缓倾界面的具有双曲时距曲线特征的干扰波和S波的高视速度会对P波的波场分离产生较大干扰;而由于S波的低速和强能量特征,其波场分离受具有双曲时距曲线特征的干扰波的视速度影响较小,且不受P波的影响.     

松辽盆地长岭断陷火山岩相地震特征及其响应关系

火山岩储层发育受岩性岩相的控制,为了提高岩性岩相地震识别的精度,选择松辽盆地南部长岭断陷营城组和火石岭组典型钻井火山岩进行井旁地震相分析,分别建立了酸性和中基性火山岩的岩相识别的模版.酸性与中基性火山岩最明显的差别为爆发相和喷溢相之间的差异,酸性岩爆发相具有席状,板状,平行—亚平行反射,连续性好、局部中等,强振幅,低频的特征;酸性岩喷溢相具有楔状、局部透镜,波形反射特征,连续性中—差,中弱振幅、见中强振幅,中高频特征.中基性岩爆发相具有板状,楔状,蠕虫形反射、偶见亚平行反射,中高振幅,中高频,连续性差、见连续性中—好的特征;中基性岩喷溢相具有板状,楔状,平行—亚平行反射、局部波形反射,连续性中等、局部较差,中强振幅,中频、局部高频的特征.根据酸性、中基性火山岩不同的电性特征,结合火山岩段表现的地震内部反射结构,几何形态等地震相特征.总结出不同火山岩之间地震相与测井相的对应关系,为新区识别火山岩相,预测火山岩储层,划分有利储集相带提供依据.     

大陆架科学钻探CSDP-2井的垂直地震剖面测量

南黄海海相地层的地震波场特征和层位标定一直是困扰地震勘探的重要问题.为了近距离、高精度和高分辨率地观测井周围构造特征和岩石性质引起的波场变化,为地震资料的采集、处理与解释提供地震波衰减规律、速度与层位标定等信息,对大陆架科学钻探CSDP-2井实施了近零偏移距垂直地震剖面（VSP）观测.针对海相地层顶部强反射界面地震波穿透难的问题,采用了大容量气枪震源并设计了气枪阵列组合方式,提高了激发地震波的能量,获得了强反射界面之下清晰的PP、PS下行波和上行波信号.采用了三分量偏振合成、组合滤波和波场分离等处理方法,对VSP观测数据进行处理,获得了海相三叠系—志留系的精细的纵波、横波速度结构和地层吸收因子等物性数据,建立了钻井地层、测井、VSP上行波和多道地震剖面对应关系,实现了不同尺度的地质和地球物理属性资料的有效衔接,标定了钻井地质剖面上各深度地质体的地震反射特性,厘定了过井地震剖面上反射同相轴的地质属性.此次观测取得的纵波、横波速度信息,成为建立南黄海海相地层速度模型主要的资料来源,也是地震资料的岩性反演处理不可缺少的信息.     

碳酸盐岩裂缝型储层最小二乘偏移成像方法研究

碳酸盐岩裂缝型储层是海相油气储层以及西部探区的重要储层类型之一,多表现为缝洞储集体形态,埋藏深度大,地震反射波信号复杂,波场特征分析困难,因而对其进行高分辨成像具有一定的难度.本文推导并实现了基于最小二乘理论的偏移成像(LSM)方法,并将LSM用于西部探区碳酸盐岩裂缝型储层等效地震地质模型的偏移成像工作,成像结果表明:在理论上,LSM方法能够较好地分辨出埋藏深度在3000 m以下的中深部20 m尺度碳酸盐岩裂缝型储层,在西部碳酸盐岩裂缝型储集体地震勘探中,具有较好的适应性.     

VSP上下行反射波联合成像方法研究

VSP资料上下行波场发育丰富.本文在分析VSP直达波、上行反射波、下行反射波传播路径及其照明范围的基础上,指出了常规VSP波动方程偏移方法缺陷,进而通过修改波场延拓方式,提出了上下行反射波联合成像方法,并在高频近似下分析了该方法的成像原理.该方法不需要进行VSP上下行反射波场分离,能够同时对VSP资料中的一次反射波、自由表面多次波、层间多次波进行成像,比常规成像剖面具有更宽的成像范围和更好的成像效果.该方法能够对下行一次反射波进行成像,从而可以实现常规偏移方法难以处理的高陡倾角构造成像.模拟资料和实际资料处理证明了本文方法的正确性.     

煤矿陷落柱的地震响应特征与识别

根据煤矿陷落柱的地质特征,设计了两类非均匀陷落柱模型,通过有限差分法模拟陷落柱的地震波场.通过对模拟波场叠加剖面和偏移剖面的对比,得出了大、小陷落柱的地震响应特征,并指导了实际地震剖面的陷落柱识别.结果表明:大的陷落柱的地震响应在时间剖面上特征表现为标准反射波的中断;而较小的陷落柱的地震响应在时间剖面上特征表现为弱反射能量异常.     

基于匹配滤波的多次波压制方法研究

常规滤波方法对浅海多次波压制不能取得理想效果.通过改变模型表面边界条件,求解波动方程实现多情景地震数值模拟.利用自由边界条件模拟时获得含海平面强反射多次波地震记录,而利用吸收边界条件时得到理想条件下无多次波地震记录.对比分析多次波对有效地震信号的影响得到海平面强反射多次波改变了有效地震信号的振幅、频率和相位.利用归一化方法计算得到压制多次波滤波器算子,并作用于含海平面强反射多次波地震记录.叠加地震剖面对比分析表明,归一化滤波器算子较好的压制了地震记录中的多次波,提高地震资料分辨率.最后,根据多次波产生条件,给出了利用实际区域简化地质模型的地震数值模拟结果计算归一化滤波器算子的方法,从而利用该方法压制实际地震数据中多次波.     

井间透射-反射波场分离方法及应用

近十几年来,井间地震在石油勘探领域发挥着重要作用.然而井间地震波极其复杂,给同相轴拾取、层析成像、偏移成像等工作带来了困难.因此,波场分离成为井间地震数据处理的一项重要内容.目前,不少学者研究了井间上下行波分离,纵横波分离等多种方法,但尚未出现相关文献研究透射波和反射波的分离方法.本文根据井间透射波和反射波在不同道集(共炮点集或共接收点集)中视速度极性的差异,首次提出了一种"透射.反射"波场分离方法.该方法首先在不同道集中利用τ-p变换进行视速度分离,得到反射上行波、反射下行波、透射上行波、透射下行波四个波场,然后选择对应的波场叠加,从而得到透射波和反射波.其中,透射波为井间层析成像提供了相对清晰的初至同相轴,尤其是横渡初至同相轴,以便进行到时提取;所得反射波为井中偏移成像提供了可靠的反射波数据,比人工切除透射波方法更为准确有效.本文通过理论模型试验和实际数据处理,证明了该方法的有效性和实用性.此外,在本文基础上,再进行纵横波分离处理可以得到更多的波场分量,为井间成像提供更有效的数据.     

云南思茅—中甸地震剖面的地壳结构

云南思茅—中甸宽角反射/折射地震剖面切割松潘—甘孜、扬子和华南三个构造单元的部分区域. 我们利用初至波和壳内反射波走时层析成像获得地壳纵波速度结构. 在获得新的地壳速度结构模型基础上，利用地震散射成像思想和低叠加次数的叠前深度偏移方法重建了研究区的地壳、上地幔反射结构. 综合分析研究区地壳P波速度模型和壳内地震反射剖面发现：沿测线从北至南地壳厚度从约50 km减薄至35 km左右，地壳厚度的减薄量主要体现在下地壳，剖面北段下地壳厚度约为30 km，剖面南段下地壳厚度仅为15 km左右；上地幔顶部局部位置P波速度值偏低，一般为76～78 km/s，反映出云南地区是典型的构造活动区的特点.剖面沿线地壳内地震反射发育，其中莫霍强反射出现在景云桥下方；在景云桥弧形断裂带8～10 km深处出现宽约50 km的强反射带.     

云南思茅—中甸地震剖面的地壳结构

云南思茅—中甸宽角反射/折射地震剖面切割松潘—甘孜、扬子和华南三个构造单元的部分区域. 我们利用初至波和壳内反射波走时层析成像获得地壳纵波速度结构. 在获得新的地壳速度结构模型基础上，利用地震散射成像思想和低叠加次数的叠前深度偏移方法重建了研究区的地壳、上地幔反射结构. 综合分析研究区地壳P波速度模型和壳内地震反射剖面发现：沿测线从北至南地壳厚度从约50 km减薄至35 km左右，地壳厚度的减薄量主要体现在下地壳，剖面北段下地壳厚度约为30 km，剖面南段下地壳厚度仅为15 km左右；上地幔顶部局部位置P波速度值偏低，一般为76～78 km/s，反映出云南地区是典型的构造活动区的特点.剖面沿线地壳内地震反射发育，其中莫霍强反射出现在景云桥下方；在景云桥弧形断裂带8～10 km深处出现宽约50 km的强反射带.