基于逆散射级数法虚反射压制方法.pdfVIP

·668· 中国地球物理2011 (18)储层地球物理 基于逆散射级数法的虚反射压制方法 Inverse for series scatteringdeghosting 王芳芳1+李景叶1刘华锋2 Li Liu Wang Huafen92 FangfangJingyel 中国石油大学(北京)cr,rvc物探重点实验室北京昌平102249 1．引言 在海上地震勘探中，虚反射(鬼波)是一种不可避免的干扰波。虚反射的波形、频率、视速度等都与 一次波相似，从而严重干扰～次反射波，降低地震剖面的分辨率，有时甚至会产生假的同相轴，给地震地 质解释造成困难Il】。因此消除虚反射的影响是海上地震资料处理的～个重要问题。当前压制虚反射的方法， 往往需要估计虚反射参数(例如虚反射时延等)，不适合于多维复杂介质。本文研究的逆散射级数法在频 率一波数一波数域内直接消除虚反射，该方法不需要估计虚反射参数，而且能适用于多维复杂介质情况。 2．方法原理 在散射理论中，散射场岱通过Lippmann-Schwinger方程给出 (1) G=G—Go—GoVG—GOVGo+GoVGoVGO+… 其中，G和G。分别表示实际波场和参考波场在某一瞬时频率W下的微分算子及格林函数，y为扰动算子。 对于海洋勘探，选择参考介质为半空间的水体，空气和水的分界位置为z=O处。对应的参考格林算子Go 包括两部分：Go=G0d+G∥。其中Go“表示从震源直接传播到接收点，酣是与自由界面相关的项。将散射 God是源点型虚反射，Go“矿G矿是接收点型虚反射，G∥yG∥是混合型虚反射。二阶、高阶散射类似。消 除虚反射就是(1)式左乘GoaGo～，右乘Go—Goal。2|。 在二维介质中，压制虚反射的具体表达式为： “’ 之娶冬』采热。Dt(缸，磊，恕，磊；们(2) 矽石习≯万面。∥％磊忆磊； 其中，D表示为频率一波数一波数域内原始散射场数据。D’表示消除虚反射后的散射场数据，qg、q。分别 是对应于k、k的垂直波数，z。|与磊分别是指震源、接收点相对于自由界面的深度。 基于逆散射理论的压制虚反射的具体步骤主要归纳为两步：第～，从测得的总场中减掉参考场，得到 只含有一次波和多次波的数据。第二，在频率一波数一波数域内计算压制虚反射后的数据，在二维介质中 作如(2)式处理，但是当分母项很小时，计算将不稳定，所以需要对分母项进行约束。该方法在三维数 据上的应用原理与二维情况基本一致。另外，基于逆散射级数法的虚反射压制方法是针对全波场数据进行 的，该方法要求炮间距与道间距相同，要求数据有近偏移距数据以及数据内没有缺失道。所以需要对数据 进行内插或外推等规则化处理得到规则化观测系统。 3．结论 经测试，该方法对合成数据和实际数据都有较好的结果。通过上述研究，得到以下结论：基于逆散射 级数法的虚反射压制方法，不需要估计虚反射参数，不需要任何地下结构信息，能适用于多维复杂介质， 计算速度快，在频率．波数．波数域可以较好的压制虚反射，能够有效地恢复一次波。但是该方法要求观测 系统规则，要求已知震源与检波点相对于自由界面的深度。 参考文献 m斫neseismic in mettIod ghostreflectionsf-k 999，l 【