基于双曲线拉东变换的散射噪声压制方法及装置制造方法及图纸

本发明专利技术提供了一种基于双曲线拉东变换的散射噪声压制方法及装置。该方法包括：将共炮点单炮记录按照中心位置分为左支记录及右支记录；对左支记录及右支记录进行双曲线拉东变换，并对双曲线拉东变换后的左支记录及右支记录进行散射波切除，得到分别对应于左支记录及右支记录的散射波；对散射波进行反拉东变换，得到时间域的散射波信号；左支记录及右支记录分别减去对应的时间域的散射波信号，得到左支去除散射波记录及右支去除散射波记录；将左支去除散射波记录及右支去除散射波记录合成，得到去除散射波的共炮点单炮记录。本发明专利技术通过双曲线拉东变换的方式压制散射噪声，并对散射噪声进行切除，在减少了计算量的同时，达到有效压制散射噪声的目的。压制散射噪声的目的。压制散射噪声的目的。

【技术实现步骤摘要】
基于双曲线拉东变换的散射噪声压制方法及装置


[0001]本专利技术涉及石油勘探地震资料处理领域，尤指一种基于双曲线拉东变换的散射噪声压制方法及装置。

技术介绍

[0002]长期以来，噪声对地震资料品质的影响一直是地震资料处理中的一个重要问题。在复杂地区，如西部山地，近地表非均匀、存在地质异常体及地形起伏比较严重，在地震资料中会产生大量散射噪声，影响地震资料品质。根据散射波时距曲线形态，可以分为线性散射和双曲型散射。正向散射、反向散射和双向散射都属于线性散射。针对具有双曲规律的散射噪声来说，F-K滤波方法存在一定缺陷。
[0003]由于目前西部复杂地区的地震勘探存在严重的近地表散射干扰，地震资料品质差。Ernst等(1999)研究了近地表的散射噪声的衰减，这种方法认为散射过程是具有确定性特点的，因此可以对散射波进行成像，然后计算散射噪声，并在此基础之上去除地震资料中的散射噪声。该方法对于复杂的近地表比较有效，但是缺点是计算量很大。

技术实现思路

[0004]为了解决上述问题，本专利技术实施例提供一种基于双曲线拉东变换的散射噪声压制方法，所述方法包括：
[0005]将共炮点单炮记录按照中心位置分为左支记录及右支记录；
[0006]对所述左支记录及所述右支记录进行双曲线拉东变换，并对双曲线拉东变换后的所述左支记录及所述右支记录进行散射波切除，得到分别对应于所述左支记录及所述右支记录的散射波；
[0007]对所述散射波进行反拉东变换，得到时间域的散射波信号；
[0008]所述左支记录及所述右支记录分别减去对应的时间域的散射波信号，得到左支去除散射波记录及右支去除散射波记录；
[0009]将所述左支去除散射波记录及所述右支去除散射波记录合成，得到去除散射波的共炮点单炮记录。
[0010]可选的，在本专利技术一实施例中，所述对所述左支记录及所述右支记录进行双曲线拉东变换包括：利用炮检距及地震波速度，对所述左支记录及所述右支记录进行针对共炮点时距双曲线的正拉东变换。
[0011]可选的，在本专利技术一实施例中，所述利用炮检距及地震波速度，对所述左支记录及所述右支记录进行针对共炮点时距双曲线的正拉东变换包括：利用正拉东变换公式对所述左支记录及所述右支记录进行变换，所述正拉东变换公式为：
[0012][0013]其中，u(τ,q)为拉东变换的结果，d为拉东变换的过程，n＝1,2...N
x
为炮检距范
围，x为炮检距，q表示地震波速度变化，t与τ表示时间。
[0014]可选的，在本专利技术一实施例中，所述对所述散射波进行反拉东变换包括：利用炮检距及地震波速度，对所述散射波进行反拉东变换。
[0015]可选的，在本专利技术一实施例中，所述利用炮检距及地震波速度，对所述散射波进行反拉东变换包括：利用反拉东变换公式对所述散射波进行反变换，所述反拉东变换公式为：
[0016][0017]其中，d'(t,x)为反拉东变换的结果，u为反拉东变换的过程，i＝1,2...N
p
为地震波速度变化范围，x为炮检距，q表示地震波速度变化，t与τ表示时间。
[0018]本专利技术实施例还提供一种基于双曲线拉东变换的散射噪声压制装置，所述装置包括：
[0019]记录划分模块，用于将共炮点单炮记录按照中心位置分为左支记录及右支记录；
[0020]拉东变换模块，用于对所述左支记录及所述右支记录进行双曲线拉东变换，并对双曲线拉东变换后的所述左支记录及所述右支记录进行散射波切除，得到分别对应于所述左支记录及所述右支记录的散射波；
[0021]反拉东变换模块，用于对所述散射波进行反拉动变换，得到时间域的散射波信号；
[0022]散射波去除模块，用于所述左支记录及所述右支记录分别减去对应的时间域的散射波信号，得到左支去除散射波记录及右支去除散射波记录；
[0023]记录合成模块，用于将所述左支去除散射波记录及所述右支去除散射波记录合成，得到去除散射波的共炮点单炮记录。
[0024]可选的，在本专利技术一实施例中，所述拉东变换模块包括：拉东变换单元，用于利用炮检距及地震波速度，对所述左支记录及所述右支记录进行针对共炮点时距双曲线的正拉东变换。
[0025]可选的，在本专利技术一实施例中，所述反拉东变换模块包括：反拉东变换单元，用于利用炮检距及地震波速度，对所述散射波进行反拉东变换。
[0026]本专利技术实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现以下步骤：
[0027]将共炮点单炮记录按照中心位置分为左支记录及右支记录；
[0028]对所述左支记录及所述右支记录进行双曲线拉东变换，并对双曲线拉东变换后的所述左支记录及所述右支记录进行散射波切除，得到分别对应于所述左支记录及所述右支记录的散射波；
[0029]对所述散射波进行反拉东变换，得到时间域的散射波信号；
[0030]所述左支记录及所述右支记录分别减去对应的时间域的散射波信号，得到左支去除散射波记录及右支去除散射波记录；
[0031]将所述左支去除散射波记录及所述右支去除散射波记录合成，得到去除散射波的共炮点单炮记录。
[0032]本专利技术实施例还提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
[0033]将共炮点单炮记录按照中心位置分为左支记录及右支记录；
[0034]对所述左支记录及所述右支记录进行双曲线拉东变换，并对双曲线拉东变换后的所述左支记录及所述右支记录进行散射波切除，得到分别对应于所述左支记录及所述右支记录的散射波；
[0035]对所述散射波进行反拉东变换，得到时间域的散射波信号；
[0036]所述左支记录及所述右支记录分别减去对应的时间域的散射波信号，得到左支去除散射波记录及右支去除散射波记录；
[0037]将所述左支去除散射波记录及所述右支去除散射波记录合成，得到去除散射波的共炮点单炮记录。
[0038]本专利技术通过双曲线拉东变换的方式压制散射噪声，并对散射噪声进行切除，在减少了计算量的同时，达到有效压制散射噪声的目的。
附图说明
[0039]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0040]图1为本专利技术实施例一种基于双曲线拉东变换的散射噪声压制方法的流程图；
[0041]图2为本专利技术实施例中散射波理论模型示意图；
[0042]图3A-图3C为本专利技术实施例中用双曲拉东变换去除理论模型中左半支散射同相轴的散射波噪声的示意图；
[0043]图4A-图4C为本专利技术实施例中用双曲拉东变换本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于双曲线拉东变换的散射噪声压制方法，其特征在于，所述方法包括：将共炮点单炮记录按照中心位置分为左支记录及右支记录；对所述左支记录及所述右支记录进行双曲线拉东变换，并对双曲线拉东变换后的所述左支记录及所述右支记录进行散射波切除，得到分别对应于所述左支记录及所述右支记录的散射波；对所述散射波进行反拉东变换，得到时间域的散射波信号；所述左支记录及所述右支记录分别减去对应的时间域的散射波信号，得到左支去除散射波记录及右支去除散射波记录；将所述左支去除散射波记录及所述右支去除散射波记录合成，得到去除散射波的共炮点单炮记录。2.根据权利要求1所述的方法，其特征在于，所述对所述左支记录及所述右支记录进行双曲线拉东变换包括：利用炮检距及地震波速度，对所述左支记录及所述右支记录进行针对共炮点时距双曲线的正拉东变换。3.根据权利要求2所述的方法，其特征在于，所述利用炮检距及地震波速度，对所述左支记录及所述右支记录进行针对共炮点时距双曲线的正拉东变换包括：利用正拉东变换公式对所述左支记录及所述右支记录进行变换，所述正拉东变换公式为：其中，u(τ,q)为拉东变换的结果，d为拉东变换的过程，n＝1,2...N
x
为炮检距范围，x为炮检距，q表示地震波速度变化，t与τ表示时间。4.根据权利要求1所述的方法，其特征在于，所述对所述散射波进行反拉东变换包括：利用炮检距及地震波速度，对所述散射波进行反拉东变换。5.根据权利要求4所述的方法，其特征在于，所述利用炮检距及地震波速度，对所述散射波进行反拉东变换包括：利用反拉东变换公式对所述散射波进行反变换，所述反拉东变换公式为：其中，d'(t,x)为反拉东变换的结果，u为反拉东变换的过程，i＝1,2...N
p
为地震波速度变化范围，x为炮检距，q表示地震波速度变化，t与τ表示时间。6.一种基于双曲线拉东变换的散射噪声压制装置，其特征在于，所述装置包括：记录划分模块，用于将共炮点单炮记录按照中心位置分为左支记录及右支记录；拉东变换模块，用于对所述左支记录及所述右支记录进行双曲线拉东变换，并对双曲线拉东变换后的所述左支记录及所述右支记录进行散射波切除...

【专利技术属性】
技术研发人员：王万里，杨午阳，魏新建，禄娟，何欣，陈德武，李冬，王伟，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：