一种提高地震资料分辨率的方法技术

本发明专利技术公布一种提高地震资料分辨率的方法，涉及多次波发育地震资料的叠前处理，该方法首先基于SRME和聚焦变换理论推导得出利用多次波的多维子波反褶积模型；然后利用引入整形正则化的非稳态回归自适应匹配滤波方法将聚焦域内由多次波恢复的高分辨率数据分离出来，实现原始数据的高分辨率转换。本发明专利技术利用多次波提高地震资料分辨率，开辟了新的多次波应用领域。该方法实现基于波动理论，可以适用于任意复杂情况；不存在现有方法中的各种假设；每一道输出结果中所有炮记录都参与运算，从空间上加以约束，在提高纵向分辨率的同时可以很好的改善数据的横向分辨率。模型和实际数据的应用都证明了该方法对于提高地震资料分辨率非常实用和有效。

【技术实现步骤摘要】

本专利技术属于地球物理勘探领域，涉及多次波发育地震资料的叠前处理，具体是一 种提高地震资料分辨率的方法。
技术介绍
地震资料中的多次波一直是困扰地震勘探数据处理以及解释的棘手难题。为了消 除多次波对一次反射波成像的影响，传统数据处理中将其视为噪音，在偏移成像前需要对 其进行压制处理。 为了适应精细勘探的发展，深入挖掘及利用不同地震波场中蕴含的地下构造信息 成为一种新的发展趋势。由于多次波也是来自地下界面的真实反射，与一次波相比，多次波 具有地下传播路径长、覆盖区域广及反射角度小的特点。近年来，不同于传统处理中将多次 波视为噪音进行去除，国内外地球物理专家和学者将多次波作为有效信号，对多次波的有 效利用进行了深入的研宄。但是现有技术中，对多次波的有效利用仅局限于利用多次波成 像获取地下更丰富的地下构造信息，对多次波信息的有效利用尚存在较大的发展空间。 另外，野外地震资料受采集因素（如海上地震勘探中震源激发产生的气泡效应及 其他采集因素）影响，采集的资料分辨率较低，因此数据处理中需要进行反褶积来提高地 震资料的分辨率。然而，传统反褶积处理技术及其改进技术都是基于一定的子波相位和反 射系数假设条件，基于反褶积理论通过压缩子波实现地震数据分辨率的提高。该类方法存 在以下不足：1)假设条件在实际应用中很难满足，严重影响了传统反褶积方法在工业生产 中的应用效果；2)准确提取地震子波是传统反褶积实现的前提，但是子波提取方法多基于 特定假设且受地质复杂性等多种因素影响，提取的子波无法适应复杂地质条件下高精度处 理的需求；3)传统反褶积模型基于反褶积理论，仅考虑了地震数据的纵向分辨率，忽视了 对地震数据的空间约束，对复杂地质条件适应性差。
技术实现思路
针对上述技术难题，本专利技术的目的在于提供，目 的在于：1)深入挖掘地震数据中多次波有效信息；2)解决传统反褶积方法中由于无法满足 子波相位和反射系数假设、子波提取不准确对提高地震资料分辨率的影响。 依据本专利技术的技术方案，包括以下步骤： 1)对采集的地震资料采用拉东变换进行数据规则化，为一次波估计和聚焦变换准 备输入数据； 2)利用表层多次波去除方法（SRME)从地震数据中提取一次波估计，并利用估计 的一次波构建正聚焦变换算子F ; F = G-1= ΔΡ_1 (1) 式中，G是反聚焦变换算子，ΔP是一次波估计对应的Berhout数据矩阵（Berkhout 数据矩阵是指由Berkhout提出的用于描述地震波传播的频率域数据矩阵，以下简称数据 矩阵）； 3)利用正聚焦变换算子对地震数据进行正聚焦变换，多次波实现降阶并表现为原 始数据的多维子波反褶积； 按照下式进行正聚焦变换： Q = FP = ΔΡ_1Ρ = I+AP (2) 式中，P = ΔΡ+Μ，是地震数据对应的数据矩阵;M = (ΔΡΑ)Ρ，表示多次波数据矩 阵;A是地表相关加权算子；Q是正聚焦变换结果；I = Λ F1 Λ Ρ，表示聚焦域内位于原点周 围的一次波聚焦能量；多次波则恢复为原始数据的加权，表示为AP ;该过程实现了多次波 的降阶。 根据波场传播（WRW)模型，地表算子A表不如下： A=EST1IT (3) 式中，S+为下行的震源波场数据矩阵，每一列代表一个子波震源向量；IT为自由地 表反射系数矩阵，可将R_表示为负的单位矩阵-1，因此， A = - (4) 式中，算子A的每一列表示负的反子波频率域向量，在时间域内表示为子波反褶 积滤波器。 在计算过程中每一炮都参与运算，根据式（3)结合反褶积理论可知，式（2)中第二 项由多次波恢复的数据相当于对原始数据进行了多维子波反褶积。利用该多维子波反褶积 模型通过提取多次波恢复的数据即可达到提高地震数据分辨率的目的。 4)对聚焦域结果和地震数据分别进行τ -P正变换，将聚焦域内一次波和多次波 映射到不同的分布区域； τ -P域转换将一次波映射到截距时间τ = 〇处的直线上，多次波则表现为位于 τ >〇区域具有一定曲率的椭圆，有效避免了一次波和多次波能量在聚焦域内的干涉，利于 高分辨率数据的提取；步骤4)所述的τ -P正变换采用的是高分辨率τ -P正变换，该方法 能够在有效提高τ -P变换精度、更好的实现数据保幅的同时有效压制聚焦域内的随机噪 音。 5)利用自适应匹配滤波方法对聚焦域数据和地震数据的正τ -P变换结果进行自 适应匹配，分离出τ -P域内由多次波恢复的高分辨率数据； 步骤5)所述的自适应匹配滤波方法采用的是引入整形正则化的非稳态回归自适 应匹配滤波方法。相比传统的匹配滤波方法，该方法考虑了滤波系数的空间变化，不需要进 行分时窗处理，整形正则化避免了不稳定问题的产生，该方法对应的目标函数表示如下：(5) 式中，X是多维空间坐标，p' (X)是聚焦域响应，Pk(X)是由地震数据经不同空间 移动得到的序列，ak是滤波系数，S是正则化算子，λ是正则化系数。 6)将分离的高分辨率数据进行τ -P反变换得到利用多次波获取的高分辨率数 据，实现原始数据的高分辨率转换。 本专利技术的有益效果是：1)从利用多次波提高地震数据分辨率角度出发，通过深入 挖掘多次波中蕴含的有效信息开辟了新的多次波应用领域。2)相比传统反褶积模型，本发 明基于波动理论推导得出的多维子波反褶积模型不存在传统反褶积方法中的各种子波相 位和反射系数假设；不需要提取子波；基于波动理论推导得到，可以适用于任意复杂地质 情况；每一道输出结果中所有炮记录都参与了运算，从空间上加以约束，在提高纵向分辨率 的同时可以很好的改善数据的横向分辨率；因此，本专利技术相比现有技术能够更好的实现地 震资料分辨率的提高。【附图说明】 图1为本专利技术提高地震资料分辨率方法流程图 图2为应用本专利技术得到的模型处理结果 图2 (a)为原始单炮记录； 图2 (b)为利用SRME提取的一次波估计； 图2(c)为正聚焦变换结果； 图2 (d)为聚焦域结果的τ -P域映射； 图2 (e)为原始单炮记录的τ -P域映射； 图2 (f)为τ -P域内利用匹配滤波提取的多次波恢复响应； 图2 (g)为反τ -P变换得到的利用多次波构建的高分辨率数据； 图3为模型数据的多炮叠加结果对比图 图3 (a)为原始数据的叠加结果； 图3(b)为传统反褶积后数据的叠加结果； 图3(c)为利用多次波的多维子波反褶积后数据的叠加结果； 图4为模型数据的多炮叠加结果的振幅谱对比图 图4(a)为原始数据叠加结果的振幅谱； 图4(b)为传统反褶积后数据叠加结果的振幅谱； 图4(c)为利用多次波的多维子波反褶积后数据叠加结果的振幅谱； 图5为应用本专利技术得到的实际资料处理结果 图5 (a)为原始数据近偏移距剖面； 图5 (b)为应用本专利技术得到的高分辨率数据的近偏移距剖面。【具体实施方式】 为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例对 本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并 不用于限定本专利技术。 图1是本专利技术所述的一种提高地震资料分辨率方法的基本流程图，包括以下步 骤： SlOl :对采集的地震资料采用拉东变换进行数据规则化，为一次波估计和聚焦变 换准备输入数据； S102 :利用表层多本文档来自技高网...

【技术保护点】
一种提高地震资料分辨率的方法，其特征在于，包括以下步骤：1)对采集的地震资料采用拉东变换进行数据规则化，为一次波估计和聚焦变换准备输入数据；2)利用表层多次波去除方法(SRME)从地震数据中提取一次波估计，并利用估计的一次波构建正聚焦变换算子；3)利用正聚焦变换算子对地震数据进行正聚焦变换，多次波实现降阶并表现为原始数据的多维子波反褶积；4)对聚焦域结果和地震数据分别进行τ‑p正变换，将聚焦域内一次波和多次波映射到不同的分布区域；5)利用自适应匹配滤波方法对聚焦域数据和地震数据的正τ‑p变换结果进行自适应匹配，分离出τ‑p域内由多次波恢复的高分辨率数据；6)将分离的高分辨率数据进行τ‑p反变换得到利用多次波获取的高分辨率数据，实现原始数据的高分辨率转换。

【技术特征摘要】

【专利技术属性】
技术研发人员：李志娜，李振春，王鹏，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：山东;37