一种地震数据低频信息补偿方法及装置的制造方法
【专利摘要】本发明涉及石油勘探地震信号处理领域，尤其涉及一种地震数据低频信息补偿方法及装置。该方法包括：获取检波器传递函数通用表达式；根据所述检波器传递函数通用表达式获得二阶低频补偿滤波器的传递函数；根据所述二阶低频补偿滤波器的传递函数获得高阶低频补偿滤波器的传递函数；根据所述高阶低频补偿滤波器的传递函数对地震数据进行低频补偿处理，获得处理后的地震数据。本申请实施例的方法简单高效，实现了对实际地震数据低频信息的补偿，可以有效拓宽地震资料的低频段谱宽，为后续解释和反演提供更高质量的数据。
【专利说明】
一种地震数据低频信息补偿方法及装置
技术领域
[0001] 本发明涉及石油勘探地震信号处理领域，尤其涉及一种地震数据低频信息补偿方 法及装置。
【背景技术】
[0002] 低频信号由于其特有的抗屏蔽、抗吸收能力和稳定性，在深层目标和复杂构造成 像、研究地质构造内部反射特性、提尚置如置后反演阻抗的精确度、提尚全波反演的稳定性 及准确性、提高地震波传播速度分析精度等方面具有重要的意义。低频信息对于地震资料 处理、解释及反演具有重要的作用。
[0003] 目前，由于硬件成本等方面条件的限制，采用特殊设备的低频地震采集并未普及。 采用常规震源和检波器采集的地震数据，由于采集过程中震源和接收器的限制，采集到的 地震数据缺失低频成分，处理过程中面波衰减等过程也会导致低频分量进一步衰减，给后 续的解释和反演工作带来不利的影响。因此，对于常规震源和普通检波器，如何补偿地震数 据中的低频信息以满足低频地震勘探的需求是当前地震勘探中亟需解决的问题。
【发明内容】
[0004] 本申请实施例提供了一种地震数据低频信息补偿方法及装置，以补偿地震数据中 的低频信息。
[0005] 为达到上述目的，一方面，本申请实施例提供了一种地震数据低频信息补偿方法， 所述方法包括：
[0006] 获取检波器传递函数通用表达式；
[0007] 根据检波器传递函数的通用表达式获得二阶低频补偿滤波器的传递函数；
[0008] 根据所述二阶低频补偿滤波器的传递函数获得高阶低频补偿滤波器的传递函数；
[0009] 根据所述高阶低频补偿滤波器的传递函数对地震数据进行低频补偿处理，获得处 理后的地震数据。
[0010] 进一步地，所述二阶低频补偿滤波器的传递函数为：
[0012] 其中s为拉普拉斯变换参数，Hc(s)为二阶低频补偿滤波器传递函数，〇^为检波器 自然角频率，ω n为起始低频补偿角频率，〇〈 ω η〈 ω ο，为第一常系数，ξ为第二常系数，ξ = k ξ〇Λ为第三常系数。
[0013] 进一步地，所述高阶低频补偿滤波器的传递函数为：
[0015]其中Hcascade(S)为高阶低频补偿滤波器传递函数，Hc(s)为二阶低频补偿滤波器传 递函数，m为级联级数，Π表示连乘。
[0016] 进一步地，在根据检波器传递函数的通用表达式获得二阶低频补偿滤波器的传递 函数之前，还包括：
[0017] 计算所述检波器自然角频率的值和所述第一常系数的值。
[0018] 进一步地，所述计算检波器自然角频率的值和所述第一常系数的值包括：
[0019] 将检波器传递函数的通用表达式进行Z变换，获得检波器传递函数的Z变换表达 式；
[0020] 计算所述检波器传递函数的Z变换表达式中的自回归滑动平均参数的值；
[0021] 根据所述自回归滑动平均参数的值计算检波器自然角频率的值和第一常系数的 值。
[0022] 另一方面，本申请实施例还提供了一种地震数据低频信息补偿装置，所述装置包 括：
[0023] 检波器传递函数通用表达式获取单元，用于获取检波器传递函数通用表达式；
[0024] 二阶低频补偿滤波器传递函数获取单元，用于根据检波器传递函数的通用表达式 获得二阶低频补偿滤波器的传递函数；
[0025] 高阶低频补偿滤波器传递函数获取单元，用于根据所述二阶低频补偿滤波器的传 递函数获得高阶低频补偿滤波器的传递函数；
[0026] 低频补偿处理单元，用于根据所述高阶低频补偿滤波器的传递函数对地震数据进 行低频补偿处理，获得处理后的地震数据。
[0027] 进一步地，所述二阶低频补偿滤波器的传递函数为：
[0029] 其中s为拉普拉斯变换参数，Hc(s)为二阶低频补偿滤波器传递函数，ω〇为检波器 自然角频率，ω η为起始低频补偿角频率，〇〈 ω η〈 ω Q，为第一常系数，ξ为第二常系数，ξ = k ξ〇Λ为第三常系数。
[0030] 进一步地，所述高阶低频补偿滤波器的传递函数为：
[0032]其中Hcasc；ade(S)为高阶低频补偿滤波器传递函数，H c(s)为二阶低频补偿滤波器传 递函数，m为级联级数，Π表示连乘。
[0033] 进一步地，所述装置还包括：
[0034] 参数计算单元，用于计算所述检波器自然角频率的值和所述第一常系数的值。
[0035] 进一步地，所述计算检波器自然角频率的值和所述第一常系数的值包括：
[0036] 将检波器传递函数的通用表达式进行Z变换，获得检波器传递函数的Z变换表达 式；
[0037] 计算所述检波器传递函数的Z变换表达式中的自回归滑动平均参数的值；
[0038] 根据所述自回归滑动平均参数的值计算检波器自然角频率的值和第一常系数的 值。
[0039]本申请实施例根据低频补偿原理，利用检波器传递函数的通用表达式获得二阶低 频滤波器的传递函数，再根据滤波器的级联特性，利用二阶低频滤波器的传递函数获得高 阶低频滤波器的传递函数，再利用高阶低频滤波器的传递函数对地震数据进行低频补偿， 从而实现了对地震数据低频信息的补偿。本申请实施例的方法简单高效，通过对实际地震 数据进行低频补偿，可以有效拓宽地震资料的低频段谱宽，为后续解释和反演提供更高质 量的数据。
【附图说明】
[0040]为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提 下，还可以根据这些附图获得其他的附图。
[0041 ]图1是本申请实施例的地震数据低频信息补偿方法流程图；
[0042] 图2是本申请实施例的地震数据低频信息补偿装置结构图；
[0043] 图3(a)~图3(c)为本申请一实施例的地震叠后记录对比图；
[0044] 图4(a)~图4(c)为本申请一实施例的地震叠后记录的振幅谱对比图。
【具体实施方式】
[0045] 为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实 施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施 例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通 技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护 的范围。
[0046] 下面结合附图，对本申请实施例的【具体实施方式】作进一步的详细说明。
[0047] 参考图1，本申请实施例的地震数据低频信息补偿方法，包括：
[0048] S1、获取检波器的传递函数通用表达式。
[0049] 传递函数是指零初始条件下线性系统响应（即输出)量的拉普拉斯变换与激励（BP 输入)量的拉普拉斯变换之比。拉普拉斯变换是工程数学中常用的一种积分变换，又称为拉 氏变换。拉氏变换是一个线性变换，可将一个有引数实数t(t 2 0)的函数转换为一个引数为 复数的函数。在地震勘探中，实际的检波器接收的信号是检波器对上传到地面的地震信号 进行滤波处理后的信号，检波器的传递函数反映了检波器对上传到地面的地震信号的滤波 作用，不同的检波器具有不同的检波器传递函数。本申请实施例中，根据地震勘探中常用检 波器的传递函数的形式，可以获取检波器的传递函数通用表达式。本申请实施例中，所述检 波器的传递函数通用表达式为：
[0051]其中S为拉普拉斯变换参数，H(s)为检波器的传递函数，cdPC2为拉普拉斯变换滤 波器参数。
[0052]本申请实施例中所述检波器的传递函数表达式也可以写成如下形式：
[0054]其中ω〇为检波器自然角频率，ξ〇为第一常系数，对应的：
[0057] S2、根据所述检波器传递函数的通用表达式获得二阶低频补偿滤波器的传递函 数。
[0058]在实际地震勘探中，由于检波器设备的限制，采集到的地震信号通常经过检波器 上传后会缺失低频成分。检波器传递函数反映了检波器对上传到地面的地震信号的滤波作 用，而低频补偿原理是指根据检波器传递函数对低频信号的滤波特性，设计出新的传递函 数从而抵消检波器传递函数的低频滤波效应。
[0059]本申请实施例基于低频补偿原理，根据检波器传递函数获得新的传递函数，即二 阶低频补偿滤波器的传递函数，用于抵消检波器传递函数的低频滤波效应。
[0060]具体地，本申请实施例为了抵消检波器传递函数对低频信号的滤波特性，需要将 检波器传递函数的分母项尽可能的抵消，从而可以获得的二阶低频补偿滤波器的传递函数 为：
[0062] 其中s为拉普拉斯变换参数，Hjs)为二阶低频补偿滤波器传递函数的拉普拉斯变 换形式，c/和c/为二阶低频补偿滤波器参数。
[0063] 本申请实施例中，所述的二阶低频补偿滤波器的传递函数也可以表示成如下形 式：
[0065] 其中ωη为起始低频补偿角频率，0〈〇^〈(〇()，|为第二常系数，| = 1^()，1^为第三常系 数，对应的：
[0066] ci7 =2ξ ωη
[0067]
[0068] 在步骤S2之前，本申请实施例还包括：
[0069] 计算所述检波器的传递函数通用表达式中检波器自然角频率的值和第一常系数 的值。
[0070] 具体的，本申请实施例中所述计算检波器的传递函数表达式中检波器自然角频率 的值和第一常系数的值，包括：
[0071] 将所述检波器的传递函数通用表达式进行Z变换，获得检波器的传递函数Z变换表 达式；
[0072] 计算所述检波器的传递函数Z变换表达式中的自回归滑动平均参数的值；
[0073] 根据所述自回归滑动平均参数的值计算检波器自然角频率的值和第一常系数的 值。
[0074] Z变换(Z-transform)是将离散系统的时域数学模型转化为较简单的频域数学模 型，以简化求解过程的一种数学工具。本申请实施例中，将所述检波器的传递函数通用表达 式进行Z变换，获得检波器的传递函数Z变换表达式，具体包括：
[0075] 采用双线性变换将所述检波器的传递函数表达式进行Z变换，获得检波器的传递 函数Z变换表达式。
[0076] 本申请实施例中：
，其中z和T表示双线性变换参数， 将所述检波器传递函数的表达式变换到Z域，获得检波器传递函数的Z变换表达式：
[0078] 其中31、82、13()、131、匕2为相应的自回归滑动平均参数，并且
[0081]利用最小二乘法计算得到自回归滑动平均参数CjPC2的值，再根据公式：
[0084] 从而可以获得检波器自然角频率ω〇和第一常系数ξ〇的值。
[0085] 本申请实施例中，所述二阶低频补偿滤波器的传递函数中的起始低频补偿角频率 ωη和第三常系数k可以进行预设，然后再通过交互分析获得最佳的起始低频补偿角频率和 第三常系数的值，从而可以确定最终的二阶低频补偿滤波器的传递函数。
[0086] S3、根据所述二阶低频补偿滤波器的传递函数获得高阶低频补偿滤波器的传递函 数。
[0087] 本申请实施例中，通过滤波器的级联特性，利用二阶低频补偿滤波器的传递函数 可以后的高阶低频补偿滤波器的传递函数，从而可以提高对地震数据的低频补偿效果。本 申请实施例中，所述高阶低频补偿滤波器的传递函数为：
[0089]其中Hcasc；ade(S)为高阶低频补偿滤波器传递函数，H c(s)为二阶低频补偿滤波器传 递函数，m为级联级数，Π表示连乘。
[0091] 其中s为拉普拉斯变换参数，Hc^c^dJs)为高阶低频补偿滤波器传递函数，ω〇为检 波器自然角频率，ω η为起始低频补偿角频率，〇〈 ω η〈 ω 0，ξ为第二常系数，ξ = ，k为第三 常系数。
[0092] S4、根据所述高阶低频补偿滤波器的传递函数对地震数据进行低频补偿处理，获 得处理后的地震数据。
[0093] 本申请实施例中，根据所述高阶低频补偿滤波器的传递函数对地震数据进行低频 补偿处理，获得处理后的地震数据，具体包括：
[0094] 将地震数据的拉普拉斯变换X(s)与高阶低频补偿滤波器的传递函数Hcasc；ade( S)相 乘获得低频补偿数据的拉普拉斯变换X(s)Hrasc；a(aS);
[0095] 将所述地震信号的拉普拉斯变换X(s)与所述低频补偿数据的拉普拉斯变换X(s) Hrasc；a(a S)相加获得低频补偿后的地震数据的拉普拉斯变换XbKl+Hc^c^dJs));
[0096] 根据所述低频补偿后的地震数据的拉普拉斯变换乂(8)(1+!1。_^(8))进行拉普拉 斯反变换，从而获得处理后的地震数据。
[0097]本申请实施例根据低频补偿原理，利用检波器传递函数的通用表达式获得二阶低 频滤波器的传递函数，再根据滤波器的级联特性，利用二阶低频滤波器的传递函数获得高 阶低频滤波器的传递函数，从而可以利用高阶低频滤波器的传递函数对地震数据进行低频 补偿，从而实现了对地震数据低频信息的补偿。本申请实施例的方法简单高效，通过对实际 地震数据进行低频补偿，可以有效拓宽地震资料的低频段谱宽，为后续解释和反演提供更 高质量的数据。
[0098]为了清楚的说明本申请实施例的有益效果，下面结合附图进行说明：
[0099] 图3(a)~图3(c)为地震叠后记录对比图，图3(a)为含低频实际叠后记录，图3(b) 为高通滤波叠后记录，图3(c)为本发明低频补偿后叠后记录；图4(a)~图4(c)为地震叠后 记录的振幅谱对比图，图4(a)为含低频实际叠后记录的振幅谱，图4(b)为高通滤波叠后记 录的振幅谱，图4(c)为本发明低频补偿后叠后记录的振幅谱。通过图4(b)与图4(a)的对比， 可以看出高通滤波叠后记录的有效低频信息缺失了，有效频带宽度也比含低频实际叠后记 录的有效频带宽度小了很多；图3(b)所示高通滤波叠后记录经过本发明低频补偿处理之后 得到图3(c)所示结果，由图3(c)及图4(c)可见高通滤波后缺失的低频信息得到了有效的恢 复，与高通滤波后缺失低频信息的记录图3(b)比较可见反射界面及大小断层都更为清晰， 与图3(a)所示含低频实际叠后记录相一致，验证了本发明低频补偿方法的有效性。
[0100] 参考图2,本申请实施例的地震数据低频信息补偿装置，包括：
[0101] 检波器传递函数通用表达式获取单元21，用于获取检波器的传递函数通用表达 式。
[0102] 传递函数是指零初始条件下线性系统响应（即输出)量的拉普拉斯变换与激励（BP 输入)量的拉普拉斯变换之比。拉普拉斯变换是工程数学中常用的一种积分变换，又称为拉 氏变换。拉氏变换是一个线性变换，可将一个有引数实数t(t 2 0)的函数转换为一个引数为 复数的函数。在地震勘探中，实际的检波器接收的信号是检波器对上传到地面的地震信号 进行滤波处理后的信号，检波器的传递函数反映了检波器对上传到地面的地震信号的滤波 作用，不同的检波器具有不同的检波器传递函数。本申请实施例中，根据地震勘探中常用检 波器的传递函数的形式，可以获取检波器的传递函数通用表达式。本申请实施例中，所述检 波器的传递函数通用表达式为：
[0104] 其中s为拉普拉斯变换参数，H(s)为检波器的传递函数，cdPc2为拉普拉斯变换滤 波器参数。
[0105] 本申请实施例中所述检波器的传递函数表达式也可以写成如下形式：
[0107]其中ω〇为检波器自然角频率，ξ〇为第一常系数，对应的：
[0110] 二阶低频补偿滤波器传递函数获取单元22,用于根据所述检波器传递函数的通用 表达式获得二阶低频补偿滤波器的传递函数。
[0111] 在实际地震勘探中，由于检波器设备的限制，采集到的地震信号通常经过检波器 上传后会缺失低频成分。检波器传递函数反映了检波器对上传到地面的地震信号的滤波作 用，而低频补偿原理是指根据检波器传递函数对低频信号的滤波特性，设计出新的传递函 数从而抵消检波器传递函数的低频滤波效应。
[0112]本申请实施例基于低频补偿原理，根据检波器传递函数获得新的传递函数，即二 阶低频补偿滤波器的传递函数，用于抵消检波器传递函数的低频滤波效应。
[0113]具体地，本申请实施例为了抵消检波器传递函数对低频信号的滤波特性，需要将 检波器传递函数的分母项尽可能的抵消，从而可以获得的二阶低频补偿滤波器的传递函数 为：
[0115] 其中s为拉普拉斯变换参数，Hjs)为二阶低频补偿滤波器传递函数的拉普拉斯变 换形式，c/和c/为二阶低频补偿滤波器参数。
[0116] 本申请实施例中，所述的二阶低频补偿滤波器的传递函数也可以表示成如下形 式：
[0118]其中ωη为起始低频补偿角频率，0〈〇^〈(〇()，|为第二常系数，| = 1^()，1^为第三常系 数，对应的：
[0119] ci7 =2ξωη
[0120] c：2 = ω2η
[0121] 本申请实施例地震数据低频信息补偿装置还包括：
[0122] 参数计算单元，用于计算所述检波器的传递函数通用表达式中检波器自然角频率 的值和第一常系数的值。
[0123] 具体的，本申请实施例中所述计算检波器的传递函数表达式中检波器自然角频率 的值和第一常系数的值，包括：
[0124] 将所述检波器的传递函数通用表达式进行Ζ变换，获得检波器的传递函数Ζ变换表 达式；
[0125] 计算所述检波器的传递函数Ζ变换表达式中的自回归滑动平均参数的值；
[0126] 根据所述自回归滑动平均参数的值计算检波器自然角频率的值和第一常系数的 值。
[0127] Z变换(Z-transform)是将离散系统的时域数学模型转化为较简单的频域数学模 型，以简化求解过程的一种数学工具。本申请实施例中，将所述检波器的传递函数通用表达 式进行Z变换，获得检波器的传递函数Z变换表达式，具体包括：
[0128] 采用双线性变换将所述检波器的传递函数表达式进行Z变换，获得检波器的传递 函数Z变换表达式。
[0129] 本申请实施例中：
，其中z和T表示双线性变换参数， 将所述检波器传递函数的表达式变换到Z域，获得检波器传递函数的Z变换表达式：
[0131] 其中&1、&2、130、131、132为相应的自回归滑动平均参数，并且
[0134]利用最小二乘法计算得到自回归滑动平均参数CjPC2的值，再根据公式：
[0137] 从而可以获得检波器自然角频率ω〇和第一常系数ξ〇的值。
[0138] 本申请实施例中，所述二阶低频补偿滤波器的传递函数中的起始低频补偿角频率 ωη和第三常系数k可以进行预设，然后再通过交互分析获得最佳的起始低频补偿角频率和 第三常系数的值，从而可以确定最终的二阶低频补偿滤波器的传递函数。
[0139] 高阶低频补偿滤波器传递函数获取单元23,用于根据所述二阶低频补偿滤波器的 传递函数获得高阶低频补偿滤波器的传递函数。
[0140] 本申请实施例中，通过滤波器的级联特性，利用二阶低频补偿滤波器的传递函数 可以后的高阶低频补偿滤波器的传递函数，从而可以提高对地震数据的低频补偿效果。本 申请实施例中，所述高阶低频补偿滤波器的传递函数为：
[0142]其中Hcasc；ade(S)为高阶低频补偿滤波器传递函数，H c(s)为二阶低频补偿滤波器传 递函数，m为级联级数，Π表示连乘。
[0144] 其中s为拉普拉斯变换参数，Hc^c^dJs)为高阶低频补偿滤波器传递函数，ω〇为检 波器自然角频率，ω η为起始低频补偿角频率，〇〈 ω η〈 ω 0，ξ为第二常系数，ξ = ，k为第三 常系数。
[0145] 低频补偿处理单元24,用于根据所述高阶低频补偿滤波器的传递函数对地震数据 进行低频补偿处理，获得处理后的地震数据。
[0146] 本申请实施例中，根据所述高阶低频补偿滤波器的传递函数对地震数据进行低频 补偿处理，获得处理后的地震数据，具体包括：
[0147] 将地震数据的拉普拉斯变换X(s)与高阶低频补偿滤波器的传递函数Hras。-(s)相 乘获得低频补偿数据的拉普拉斯变换X(s)H rasc；a(aS);
[0148] 将所述地震信号的拉普拉斯变换X(s)与所述低频补偿数据的拉普拉斯变换X(s) Hrasc；a(a S)相加获得低频补偿后的地震数据的拉普拉斯变换XbKl+Hc^c^dJs));
[0149] 根据所述低频补偿后的地震数据的拉普拉斯变换乂(8)(1+!1。_^(8))进行拉普拉 斯反变换，从而获得处理后的地震数据。
[0150] 本申请实施例根据低频补偿原理，利用检波器传递函数的通用表达式获得二阶低 频滤波器的传递函数，再根据滤波器的级联特性，利用二阶低频滤波器的传递函数获得高 阶低频滤波器的传递函数，从而可以利用高阶低频滤波器的传递函数对地震数据进行低频 补偿，从而实现了对地震数据低频信息的补偿。本申请实施例的方法简单高效，通过对实际 地震数据进行低频补偿，可以有效拓宽地震资料的低频段谱宽，为后续解释和反演提供更 高质量的数据。
[0151] 在一个或多个示例性的设计中，本申请实施例所描述的上述功能可以在硬件、软 件、固件或这三者的任意组合来实现。如果在软件中实现，这些功能可以存储与电脑可读的 媒介上，或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电 脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以 是任何通用或特殊电脑可以接入访问的可用媒体。例如，这样的电脑可读媒体可以包括但 不限于RAM、ROM、EEPROM、⑶-ROM或其它光盘存储、磁盘存储或其它磁性存储装置，或其它任 何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理 器读取形式的程序代码的媒介。
[0152] 以上所述的具体实施例，对本申请的目的、技术方案和有益效果进行了进一步详 细说明，所应理解的是，以上所述仅为本申请实施例的具体实施例而已，并不用于限定本申 请的保护范围，凡在本申请的精神和原则之内，所做的任何修改、等同替换、改进等，均应包 含在本申请的保护范围之内。
【主权项】
1. 一种地震数据低频信息补偿方法，其特征在于，所述方法包括： 获取检波器传递函数通用表达式； 根据所述检波器传递函数通用表达式获得二阶低频补偿滤波器的传递函数； 根据所述二阶低频补偿滤波器的传递函数获得高阶低频补偿滤波器的传递函数； 根据所述高阶低频补偿滤波器的传递函数对地震数据进行低频补偿处理，获得处理后 的地震数据。2. 如权利要求1所述的地震数据低频信息补偿方法，其特征在于，所述二阶低频补偿滤 波器的传递函数为：其中s为拉普拉斯变换参数，Hc(s)为二阶低频补偿滤波器的传递函数，ω〇为检波器自 然角频率，ω η为起始低频补偿角频率，〇〈 ω η〈 ω 0，为第一常系数，ξ为第二常系数，ξ = k|o， k为第三常系数。3. 如权利要求2所述的地震数据低频信息补偿方法，其特征在于，所述高阶低频补偿滤 波器的传递函数为：其中Hcasc；ade(S)为高阶低频补偿滤波器的传递函数，H c(s)为二阶低频补偿滤波器的传 递函数，m为级联级数，Π 表示连乘。4. 如权利要求3所述的地震数据低频信息补偿方法，其特征在于，在根据所述检波器传 递函数的通用表达式获得二阶低频补偿滤波器的传递函数之前，还包括： 计算所述检波器传递函数的通用表达式中检波器自然角频率的值和第一常系数的值。5. 如权利要求4所述的地震数据低频信息补偿方法，其特征在于，所述计算检波器自然 角频率的值和第一常系数的值包括： 将检波器传递函数的通用表达式进行Z变换，获得检波器传递函数的Z变换表达式； 计算所述检波器传递函数的Z变换表达式中的自回归滑动平均参数的值； 根据所述自回归滑动平均参数的值计算检波器自然角频率的值和第一常系数的值。6. -种地震数据低频信息补偿装置，其特征在于，所述装置包括： 检波器传递函数通用表达式获取单元，用于获取检波器传递函数通用表达式； 二阶低频补偿滤波器传递函数获取单元，用于根据检波器传递函数的通用表达式获得 二阶低频补偿滤波器的传递函数； 高阶低频补偿滤波器传递函数获取单元，用于根据所述二阶低频补偿滤波器的传递函 数获得高阶低频补偿滤波器的传递函数； 低频补偿处理单元，用于根据所述高阶低频补偿滤波器的传递函数对地震数据进行低 频补偿处理，获得处理后的地震数据。7. 如权利要求6所述的地震数据低频信息补偿装置，其特征在于，所述二阶低频补偿滤 波器的传递函数为：其中S为拉普拉斯变换参数，Hc(s)为二阶低频补偿滤波器的传递函数，ωο为检波器自 然角频率，ω η为起始低频补偿角频率，〇〈 ω η〈 ω Ο，为第一常系数，ξ为第二常系数，ξ = k|o， k为第三常系数。8. 如权利要求7所述的地震数据低频信息补偿装置，其特征在于，所述高阶低频补偿滤 波器的传递函数为：其中Hcasc；ade(S)为高阶低频补偿滤波器的传递函数，H c(s)为二阶低频补偿滤波器的传 递函数，m为级联级数，Π 表示连乘。9. 如权利要求8所述的地震数据低频信息补偿装置，其特征在于，所述装置还包括： 参数计算单元，用于计算所述检波器传递函数的通用表达式中检波器自然角频率的值 和第一常系数的值。10. 如权利要求9所述的地震数据低频信息补偿装置，其特征在于，所述计算检波器自 然角频率的值和第一常系数的值包括： 将检波器传递函数的通用表达式进行Z变换，获得检波器传递函数的Z变换表达式； 计算所述检波器传递函数的Z变换表达式中的自回归滑动平均参数的值； 根据所述自回归滑动平均参数的值计算检波器自然角频率的值和第一常系数的值。
【文档编号】G01V1/36GK105866838SQ201610325643
【公开日】2016年8月17日
【申请日】2016年5月17日
【发明人】李海山, 杨午阳, 高建虎, 张巧凤, 周春雷
【申请人】中国石油天然气股份有限公司