CN105866838B - 一种地震数据低频信息补偿方法及装置 - Google Patents
一种地震数据低频信息补偿方法及装置 Download PDFInfo
- Publication number
- CN105866838B CN105866838B CN201610325643.2A CN201610325643A CN105866838B CN 105866838 B CN105866838 B CN 105866838B CN 201610325643 A CN201610325643 A CN 201610325643A CN 105866838 B CN105866838 B CN 105866838B
- Authority
- CN
- China
- Prior art keywords
- transfer function
- detector
- frequency
- low
- frequency compensation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000012546 transfer Methods 0.000 claims abstract description 191
- 238000012545 processing Methods 0.000 claims abstract description 17
- 230000009466 transformation Effects 0.000 claims description 26
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 abstract description 9
- 239000003208 petroleum Substances 0.000 abstract description 2
- 238000001914 filtration Methods 0.000 description 16
- 230000000694 effects Effects 0.000 description 11
- 238000006243 chemical reaction Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 5
- 238000013178 mathematical model Methods 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
- G01V1/366—Seismic filtering by correlation of seismic signals
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及石油勘探地震信号处理领域,尤其涉及一种地震数据低频信息补偿方法及装置。该方法包括:获取检波器传递函数通用表达式;根据所述检波器传递函数通用表达式获得二阶低频补偿滤波器的传递函数;根据所述二阶低频补偿滤波器的传递函数获得高阶低频补偿滤波器的传递函数;根据所述高阶低频补偿滤波器的传递函数对地震数据进行低频补偿处理,获得处理后的地震数据。本申请实施例的方法简单高效,实现了对实际地震数据低频信息的补偿,可以有效拓宽地震资料的低频段谱宽,为后续解释和反演提供更高质量的数据。
Description
技术领域
本发明涉及石油勘探地震信号处理领域,尤其涉及一种地震数据低频信息补偿方法及装置。
背景技术
低频信号由于其特有的抗屏蔽、抗吸收能力和稳定性,在深层目标和复杂构造成像、研究地质构造内部反射特性、提高叠前叠后反演阻抗的精确度、提高全波反演的稳定性及准确性、提高地震波传播速度分析精度等方面具有重要的意义。低频信息对于地震资料处理、解释及反演具有重要的作用。
目前,由于硬件成本等方面条件的限制,采用特殊设备的低频地震采集并未普及。采用常规震源和检波器采集的地震数据,由于采集过程中震源和接收器的限制,采集到的地震数据缺失低频成分,处理过程中面波衰减等过程也会导致低频分量进一步衰减,给后续的解释和反演工作带来不利的影响。因此,对于常规震源和普通检波器,如何补偿地震数据中的低频信息以满足低频地震勘探的需求是当前地震勘探中亟需解决的问题。
发明内容
本申请实施例提供了一种地震数据低频信息补偿方法及装置,以补偿地震数据中的低频信息。
为达到上述目的,一方面,本申请实施例提供了一种地震数据低频信息补偿方法,所述方法包括:
获取检波器传递函数通用表达式;
根据检波器传递函数的通用表达式获得二阶低频补偿滤波器的传递函数;
根据所述二阶低频补偿滤波器的传递函数获得高阶低频补偿滤波器的传递函数;
根据所述高阶低频补偿滤波器的传递函数对地震数据进行低频补偿处理,获得处理后的地震数据。
进一步地,所述二阶低频补偿滤波器的传递函数为:
其中s为拉普拉斯变换参数,Hc(s)为二阶低频补偿滤波器传递函数,ω0为检波器自然角频率,ωn为起始低频补偿角频率,0<ωn<ω0,ξ0为第一常系数,ξ为第二常系数,ξ=kξ0,k为第三常系数。
进一步地,所述高阶低频补偿滤波器的传递函数为:
其中Hcascade(s)为高阶低频补偿滤波器传递函数,Hc(s)为二阶低频补偿滤波器传递函数,m为级联级数,∏表示连乘。
进一步地,在根据检波器传递函数的通用表达式获得二阶低频补偿滤波器的传递函数之前,还包括:
计算所述检波器自然角频率的值和所述第一常系数的值。
进一步地,所述计算检波器自然角频率的值和所述第一常系数的值包括:
将检波器传递函数的通用表达式进行Z变换,获得检波器传递函数的Z变换表达式;
计算所述检波器传递函数的Z变换表达式中的自回归滑动平均参数的值;
根据所述自回归滑动平均参数的值计算检波器自然角频率的值和第一常系数的值。
另一方面,本申请实施例还提供了一种地震数据低频信息补偿装置,所述装置包括:
检波器传递函数通用表达式获取单元,用于获取检波器传递函数通用表达式;
二阶低频补偿滤波器传递函数获取单元,用于根据检波器传递函数的通用表达式获得二阶低频补偿滤波器的传递函数;
高阶低频补偿滤波器传递函数获取单元,用于根据所述二阶低频补偿滤波器的传递函数获得高阶低频补偿滤波器的传递函数;
低频补偿处理单元,用于根据所述高阶低频补偿滤波器的传递函数对地震数据进行低频补偿处理,获得处理后的地震数据。
进一步地,所述二阶低频补偿滤波器的传递函数为:
其中s为拉普拉斯变换参数,Hc(s)为二阶低频补偿滤波器传递函数,ω0为检波器自然角频率,ωn为起始低频补偿角频率,0<ωn<ω0,ξ0为第一常系数,ξ为第二常系数,ξ=kξ0,k为第三常系数。
进一步地,所述高阶低频补偿滤波器的传递函数为:
其中Hcascade(s)为高阶低频补偿滤波器传递函数,Hc(s)为二阶低频补偿滤波器传递函数,m为级联级数,∏表示连乘。
进一步地,所述装置还包括:
参数计算单元,用于计算所述检波器自然角频率的值和所述第一常系数的值。
进一步地,所述计算检波器自然角频率的值和所述第一常系数的值包括:
将检波器传递函数的通用表达式进行Z变换,获得检波器传递函数的Z变换表达式;
计算所述检波器传递函数的Z变换表达式中的自回归滑动平均参数的值;
根据所述自回归滑动平均参数的值计算检波器自然角频率的值和第一常系数的值。
本申请实施例根据低频补偿原理,利用检波器传递函数的通用表达式获得二阶低频滤波器的传递函数,再根据滤波器的级联特性,利用二阶低频滤波器的传递函数获得高阶低频滤波器的传递函数,再利用高阶低频滤波器的传递函数对地震数据进行低频补偿,从而实现了对地震数据低频信息的补偿。本申请实施例的方法简单高效,通过对实际地震数据进行低频补偿,可以有效拓宽地震资料的低频段谱宽,为后续解释和反演提供更高质量的数据。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请实施例的地震数据低频信息补偿方法流程图;
图2是本申请实施例的地震数据低频信息补偿装置结构图;
图3(a)~图3(c)为本申请一实施例的地震叠后记录对比图;
图4(a)~图4(c)为本申请一实施例的地震叠后记录的振幅谱对比图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
下面结合附图,对本申请实施例的具体实施方式作进一步的详细说明。
参考图1,本申请实施例的地震数据低频信息补偿方法,包括:
S1、获取检波器的传递函数通用表达式。
传递函数是指零初始条件下线性系统响应(即输出)量的拉普拉斯变换与激励(即输入)量的拉普拉斯变换之比。拉普拉斯变换是工程数学中常用的一种积分变换,又称为拉氏变换。拉氏变换是一个线性变换,可将一个有引数实数t(t≥0)的函数转换为一个引数为复数的函数。在地震勘探中,实际的检波器接收的信号是检波器对上传到地面的地震信号进行滤波处理后的信号,检波器的传递函数反映了检波器对上传到地面的地震信号的滤波作用,不同的检波器具有不同的检波器传递函数。本申请实施例中,根据地震勘探中常用检波器的传递函数的形式,可以获取检波器的传递函数通用表达式。本申请实施例中,所述检波器的传递函数通用表达式为:
其中s为拉普拉斯变换参数,H(s)为检波器的传递函数,c1和c2为拉普拉斯变换滤波器参数。
本申请实施例中所述检波器的传递函数表达式也可以写成如下形式:
其中ω0为检波器自然角频率,ξ0为第一常系数,对应的:
S2、根据所述检波器传递函数的通用表达式获得二阶低频补偿滤波器的传递函数。
在实际地震勘探中,由于检波器设备的限制,采集到的地震信号通常经过检波器上传后会缺失低频成分。检波器传递函数反映了检波器对上传到地面的地震信号的滤波作用,而低频补偿原理是指根据检波器传递函数对低频信号的滤波特性,设计出新的传递函数从而抵消检波器传递函数的低频滤波效应。
本申请实施例基于低频补偿原理,根据检波器传递函数获得新的传递函数,即二阶低频补偿滤波器的传递函数,用于抵消检波器传递函数的低频滤波效应。
具体地,本申请实施例为了抵消检波器传递函数对低频信号的滤波特性,需要将检波器传递函数的分母项尽可能的抵消,从而可以获得的二阶低频补偿滤波器的传递函数为:
其中s为拉普拉斯变换参数,Hc(s)为二阶低频补偿滤波器传递函数的拉普拉斯变换形式,c1′和c2′为二阶低频补偿滤波器参数。
本申请实施例中,所述的二阶低频补偿滤波器的传递函数也可以表示成如下形式:
其中ωn为起始低频补偿角频率,0<ωn<ω0,ξ为第二常系数,ξ=kξ0,k为第三常系数,对应的:
c1′=2ξωn
在步骤S2之前,本申请实施例还包括:
计算所述检波器的传递函数通用表达式中检波器自然角频率的值和第一常系数的值。
具体的,本申请实施例中所述计算检波器的传递函数表达式中检波器自然角频率的值和第一常系数的值,包括:
将所述检波器的传递函数通用表达式进行Z变换,获得检波器的传递函数Z变换表达式;
计算所述检波器的传递函数Z变换表达式中的自回归滑动平均参数的值;
根据所述自回归滑动平均参数的值计算检波器自然角频率的值和第一常系数的值。
Z变换(Z-transform)是将离散系统的时域数学模型转化为较简单的频域数学模型,以简化求解过程的一种数学工具。本申请实施例中,将所述检波器的传递函数通用表达式进行Z变换,获得检波器的传递函数Z变换表达式,具体包括:
采用双线性变换将所述检波器的传递函数表达式进行Z变换,获得检波器的传递函数Z变换表达式。
本申请实施例中,采用双线性变换其中z和T表示双线性变换参数,将所述检波器传递函数的表达式变换到Z域,获得检波器传递函数的Z变换表达式:
其中a1、a2、b0、b1、b2为相应的自回归滑动平均参数,并且
利用最小二乘法计算得到自回归滑动平均参数c1和c2的值,再根据公式:
从而可以获得检波器自然角频率ω0和第一常系数ξ0的值。
本申请实施例中,所述二阶低频补偿滤波器的传递函数中的起始低频补偿角频率ωn和第三常系数k可以进行预设,然后再通过交互分析获得最佳的起始低频补偿角频率和第三常系数的值,从而可以确定最终的二阶低频补偿滤波器的传递函数。
S3、根据所述二阶低频补偿滤波器的传递函数获得高阶低频补偿滤波器的传递函数。
本申请实施例中,通过滤波器的级联特性,利用二阶低频补偿滤波器的传递函数可以后的高阶低频补偿滤波器的传递函数,从而可以提高对地震数据的低频补偿效果。本申请实施例中,所述高阶低频补偿滤波器的传递函数为:
其中Hcascade(s)为高阶低频补偿滤波器传递函数,Hc(s)为二阶低频补偿滤波器传递函数,m为级联级数,∏表示连乘。
即
其中s为拉普拉斯变换参数,Hcascade(s)为高阶低频补偿滤波器传递函数,ω0为检波器自然角频率,ωn为起始低频补偿角频率,0<ωn<ω0,ξ为第二常系数,ξ=kξ0,k为第三常系数。
S4、根据所述高阶低频补偿滤波器的传递函数对地震数据进行低频补偿处理,获得处理后的地震数据。
本申请实施例中,根据所述高阶低频补偿滤波器的传递函数对地震数据进行低频补偿处理,获得处理后的地震数据,具体包括:
将地震数据的拉普拉斯变换X(s)与高阶低频补偿滤波器的传递函数Hcascade(s)相乘获得低频补偿数据的拉普拉斯变换X(s)Hcascade(s);
将所述地震信号的拉普拉斯变换X(s)与所述低频补偿数据的拉普拉斯变换X(s)Hcascade(s)相加获得低频补偿后的地震数据的拉普拉斯变换X(s)(1+Hcascade(s));
根据所述低频补偿后的地震数据的拉普拉斯变换X(s)(1+Hcascade(s))进行拉普拉斯反变换,从而获得处理后的地震数据。
本申请实施例根据低频补偿原理,利用检波器传递函数的通用表达式获得二阶低频滤波器的传递函数,再根据滤波器的级联特性,利用二阶低频滤波器的传递函数获得高阶低频滤波器的传递函数,从而可以利用高阶低频滤波器的传递函数对地震数据进行低频补偿,从而实现了对地震数据低频信息的补偿。本申请实施例的方法简单高效,通过对实际地震数据进行低频补偿,可以有效拓宽地震资料的低频段谱宽,为后续解释和反演提供更高质量的数据。
为了清楚的说明本申请实施例的有益效果,下面结合附图进行说明:
图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)所示含低频实际叠后记录相一致,验证了本发明低频补偿方法的有效性。
参考图2,本申请实施例的地震数据低频信息补偿装置,包括:
检波器传递函数通用表达式获取单元21,用于获取检波器的传递函数通用表达式。
传递函数是指零初始条件下线性系统响应(即输出)量的拉普拉斯变换与激励(即输入)量的拉普拉斯变换之比。拉普拉斯变换是工程数学中常用的一种积分变换,又称为拉氏变换。拉氏变换是一个线性变换,可将一个有引数实数t(t≥0)的函数转换为一个引数为复数的函数。在地震勘探中,实际的检波器接收的信号是检波器对上传到地面的地震信号进行滤波处理后的信号,检波器的传递函数反映了检波器对上传到地面的地震信号的滤波作用,不同的检波器具有不同的检波器传递函数。本申请实施例中,根据地震勘探中常用检波器的传递函数的形式,可以获取检波器的传递函数通用表达式。本申请实施例中,所述检波器的传递函数通用表达式为:
其中s为拉普拉斯变换参数,H(s)为检波器的传递函数,c1和c2为拉普拉斯变换滤波器参数。
本申请实施例中所述检波器的传递函数表达式也可以写成如下形式:
其中ω0为检波器自然角频率,ξ0为第一常系数,对应的:
二阶低频补偿滤波器传递函数获取单元22,用于根据所述检波器传递函数的通用表达式获得二阶低频补偿滤波器的传递函数。
在实际地震勘探中,由于检波器设备的限制,采集到的地震信号通常经过检波器上传后会缺失低频成分。检波器传递函数反映了检波器对上传到地面的地震信号的滤波作用,而低频补偿原理是指根据检波器传递函数对低频信号的滤波特性,设计出新的传递函数从而抵消检波器传递函数的低频滤波效应。
本申请实施例基于低频补偿原理,根据检波器传递函数获得新的传递函数,即二阶低频补偿滤波器的传递函数,用于抵消检波器传递函数的低频滤波效应。
具体地,本申请实施例为了抵消检波器传递函数对低频信号的滤波特性,需要将检波器传递函数的分母项尽可能的抵消,从而可以获得的二阶低频补偿滤波器的传递函数为:
其中s为拉普拉斯变换参数,Hc(s)为二阶低频补偿滤波器传递函数的拉普拉斯变换形式,c1′和c2′为二阶低频补偿滤波器参数。
本申请实施例中,所述的二阶低频补偿滤波器的传递函数也可以表示成如下形式:
其中ωn为起始低频补偿角频率,0<ωn<ω0,ξ为第二常系数,ξ=kξ0,k为第三常系数,对应的:
c1′=2ξωn
本申请实施例地震数据低频信息补偿装置还包括:
参数计算单元,用于计算所述检波器的传递函数通用表达式中检波器自然角频率的值和第一常系数的值。
具体的,本申请实施例中所述计算检波器的传递函数表达式中检波器自然角频率的值和第一常系数的值,包括:
将所述检波器的传递函数通用表达式进行Z变换,获得检波器的传递函数Z变换表达式;
计算所述检波器的传递函数Z变换表达式中的自回归滑动平均参数的值;
根据所述自回归滑动平均参数的值计算检波器自然角频率的值和第一常系数的值。
Z变换(Z-transform)是将离散系统的时域数学模型转化为较简单的频域数学模型,以简化求解过程的一种数学工具。本申请实施例中,将所述检波器的传递函数通用表达式进行Z变换,获得检波器的传递函数Z变换表达式,具体包括:
采用双线性变换将所述检波器的传递函数表达式进行Z变换,获得检波器的传递函数Z变换表达式。
本申请实施例中,采用双线性变换其中z和T表示双线性变换参数,将所述检波器传递函数的表达式变换到Z域,获得检波器传递函数的Z变换表达式:
其中a1、a2、b0、b1、b2为相应的自回归滑动平均参数,并且
利用最小二乘法计算得到自回归滑动平均参数c1和c2的值,再根据公式:
从而可以获得检波器自然角频率ω0和第一常系数ξ0的值。
本申请实施例中,所述二阶低频补偿滤波器的传递函数中的起始低频补偿角频率ωn和第三常系数k可以进行预设,然后再通过交互分析获得最佳的起始低频补偿角频率和第三常系数的值,从而可以确定最终的二阶低频补偿滤波器的传递函数。
高阶低频补偿滤波器传递函数获取单元23,用于根据所述二阶低频补偿滤波器的传递函数获得高阶低频补偿滤波器的传递函数。
本申请实施例中,通过滤波器的级联特性,利用二阶低频补偿滤波器的传递函数可以后的高阶低频补偿滤波器的传递函数,从而可以提高对地震数据的低频补偿效果。本申请实施例中,所述高阶低频补偿滤波器的传递函数为:
其中Hcascade(s)为高阶低频补偿滤波器传递函数,Hc(s)为二阶低频补偿滤波器传递函数,m为级联级数,∏表示连乘。
即
其中s为拉普拉斯变换参数,Hcascade(s)为高阶低频补偿滤波器传递函数,ω0为检波器自然角频率,ωn为起始低频补偿角频率,0<ωn<ω0,ξ为第二常系数,ξ=kξ0,k为第三常系数。
低频补偿处理单元24,用于根据所述高阶低频补偿滤波器的传递函数对地震数据进行低频补偿处理,获得处理后的地震数据。
本申请实施例中,根据所述高阶低频补偿滤波器的传递函数对地震数据进行低频补偿处理,获得处理后的地震数据,具体包括:
将地震数据的拉普拉斯变换X(s)与高阶低频补偿滤波器的传递函数Hcascade(s)相乘获得低频补偿数据的拉普拉斯变换X(s)Hcascade(s);
将所述地震信号的拉普拉斯变换X(s)与所述低频补偿数据的拉普拉斯变换X(s)Hcascade(s)相加获得低频补偿后的地震数据的拉普拉斯变换X(s)(1+Hcascade(s));
根据所述低频补偿后的地震数据的拉普拉斯变换X(s)(1+Hcascade(s))进行拉普拉斯反变换,从而获得处理后的地震数据。
本申请实施例根据低频补偿原理,利用检波器传递函数的通用表达式获得二阶低频滤波器的传递函数,再根据滤波器的级联特性,利用二阶低频滤波器的传递函数获得高阶低频滤波器的传递函数,从而可以利用高阶低频滤波器的传递函数对地震数据进行低频补偿,从而实现了对地震数据低频信息的补偿。本申请实施例的方法简单高效,通过对实际地震数据进行低频补偿,可以有效拓宽地震资料的低频段谱宽,为后续解释和反演提供更高质量的数据。
在一个或多个示例性的设计中,本申请实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现,这些功能可以存储与电脑可读的媒介上,或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如,这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置,或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。
以上所述的具体实施例,对本申请的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本申请实施例的具体实施例而已,并不用于限定本申请的保护范围,凡在本申请的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。
Claims (6)
1.一种地震数据低频信息补偿方法,其特征在于,所述方法包括:
获取检波器传递函数通用表达式;
计算所述检波器传递函数通用表达式中检波器自然角频率的值和第一常系数的值;
根据所述检波器传递函数通用表达式获得二阶低频补偿滤波器的传递函数;
根据所述二阶低频补偿滤波器的传递函数获得高阶低频补偿滤波器的传递函数;
根据所述高阶低频补偿滤波器的传递函数对地震数据进行低频补偿处理,获得处理后的地震数据;其中,
所述检波器传递函数通用表达式为:
其中,ω0为检波器自然角频率,ξ0为第一常系数,s为拉普拉斯变换参数,H(s)为检波器传递函数;其中,
计算检波器自然角频率的值和第一常系数的值包括:
将检波器传递函数通用表达式进行Z变换,获得检波器传递函数的Z变换表达式;
计算所述检波器传递函数的Z变换表达式中的自回归滑动平均参数的值;
根据所述自回归滑动平均参数的值计算检波器自然角频率的值和第一常系数的值。
2.如权利要求1所述的地震数据低频信息补偿方法,其特征在于,所述二阶低频补偿滤波器的传递函数为:
其中s为拉普拉斯变换参数,Hc(s)为二阶低频补偿滤波器的传递函数,ω0为检波器自然角频率,ωn为起始低频补偿角频率,0<ωn<ω0,ξ0为第一常系数,ξ为第二常系数,ξ=kξ0,k为第三常系数。
3.如权利要求2所述的地震数据低频信息补偿方法,其特征在于,所述高阶低频补偿滤波器的传递函数为:
其中Hcascade(s)为高阶低频补偿滤波器的传递函数,Hc(s)为二阶低频补偿滤波器的传递函数,m为级联级数,Π表示连乘。
4.一种地震数据低频信息补偿装置,其特征在于,所述装置包括:
检波器传递函数通用表达式获取单元,用于获取检波器传递函数通用表达式;
参数计算单元,用于计算所述检波器传递函数通用表达式中检波器自然角频率的值和第一常系数的值;
二阶低频补偿滤波器传递函数获取单元,用于根据检波器传递函数的通用表达式获得二阶低频补偿滤波器的传递函数;
高阶低频补偿滤波器传递函数获取单元,用于根据所述二阶低频补偿滤波器的传递函数获得高阶低频补偿滤波器的传递函数;
低频补偿处理单元,用于根据所述高阶低频补偿滤波器的传递函数对地震数据进行低频补偿处理,获得处理后的地震数据;其中,
所述检波器传递函数通用表达式为:
其中ω0为检波器自然角频率,ξ0为第一常系数,s为拉普拉斯变换参数,H(s)为检波器传递函数;其中,
计算检波器自然角频率的值和第一常系数的值包括:
将检波器传递函数通用表达式进行Z变换,获得检波器传递函数的Z变换表达式;
计算所述检波器传递函数的Z变换表达式中的自回归滑动平均参数的值;
根据所述自回归滑动平均参数的值计算检波器自然角频率的值和第一常系数的值。
5.如权利要求4所述的地震数据低频信息补偿装置,其特征在于,所述二阶低频补偿滤波器的传递函数为:
其中s为拉普拉斯变换参数,Hc(s)为二阶低频补偿滤波器的传递函数,ω0为检波器自然角频率,ωn为起始低频补偿角频率,0<ωn<ω0,ξ0为第一常系数,ξ为第二常系数,ξ=kξ0,k为第三常系数。
6.如权利要求5所述的地震数据低频信息补偿装置,其特征在于,所述高阶低频补偿滤波器的传递函数为:
其中Hcascade(s)为高阶低频补偿滤波器的传递函数,Hc(s)为二阶低频补偿滤波器的传递函数,m为级联级数,Π表示连乘。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610325643.2A CN105866838B (zh) | 2016-05-17 | 2016-05-17 | 一种地震数据低频信息补偿方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610325643.2A CN105866838B (zh) | 2016-05-17 | 2016-05-17 | 一种地震数据低频信息补偿方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105866838A CN105866838A (zh) | 2016-08-17 |
CN105866838B true CN105866838B (zh) | 2018-10-16 |
Family
ID=56635168
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610325643.2A Active CN105866838B (zh) | 2016-05-17 | 2016-05-17 | 一种地震数据低频信息补偿方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105866838B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109444958B (zh) * | 2018-11-01 | 2020-05-26 | 科吉思石油技术咨询(北京)有限公司 | 致密储层孔隙度定量预测的方法 |
CN111427088A (zh) * | 2020-03-11 | 2020-07-17 | 王仰华 | 用于识别薄互储层的地震数据低频补偿方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001029581A1 (en) * | 1999-10-21 | 2001-04-26 | Core Laboratories, Inc. | Transfer function method of seismic signal processing and exploration |
CN1643783A (zh) * | 2002-04-03 | 2005-07-20 | 安德鲁公司 | 使用并入了与频率有关的非线性的放大器模型的自适应数字预矫正 |
CN102812381A (zh) * | 2010-03-09 | 2012-12-05 | 格库技术有限公司 | 利用地震传感器传递函数来获得高保真地震成像 |
EP2685288A2 (en) * | 2012-07-10 | 2014-01-15 | PGS Geophysical AS | Methods and Systems for Reconstruction of Low Frequency Particle Velocity Wavefields and Deghosting of Seismic Streamer Data |
CN105093329A (zh) * | 2015-08-25 | 2015-11-25 | 中国石油天然气集团公司 | 地震数据的低频信息的自动补偿方法及装置 |
-
2016
- 2016-05-17 CN CN201610325643.2A patent/CN105866838B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001029581A1 (en) * | 1999-10-21 | 2001-04-26 | Core Laboratories, Inc. | Transfer function method of seismic signal processing and exploration |
CN1643783A (zh) * | 2002-04-03 | 2005-07-20 | 安德鲁公司 | 使用并入了与频率有关的非线性的放大器模型的自适应数字预矫正 |
CN102812381A (zh) * | 2010-03-09 | 2012-12-05 | 格库技术有限公司 | 利用地震传感器传递函数来获得高保真地震成像 |
EP2685288A2 (en) * | 2012-07-10 | 2014-01-15 | PGS Geophysical AS | Methods and Systems for Reconstruction of Low Frequency Particle Velocity Wavefields and Deghosting of Seismic Streamer Data |
CN105093329A (zh) * | 2015-08-25 | 2015-11-25 | 中国石油天然气集团公司 | 地震数据的低频信息的自动补偿方法及装置 |
Non-Patent Citations (1)
Title |
---|
磁电式振动速度传感器低频补偿算法;孙玉国;《仪器仪表学报》;20070430;第28卷(第4期);第444-445页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105866838A (zh) | 2016-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2599958A1 (en) | Removal of noise from seismic data using radon transformations | |
CN112817047B (zh) | 海洋地震自适应去鬼波方法、装置、电子设备及介质 | |
CN112835103B (zh) | 自适应去鬼波与宽频准零相位反褶积联合处理方法及系统 | |
Innocent Oboué et al. | Robust damped rank-reduction method for simultaneous denoising and reconstruction of 5D seismic data | |
CN103926623B (zh) | 一种压制逆时偏移低频噪音的方法 | |
CN112904426A (zh) | 一种解耦弹性波逆时偏移方法、系统及应用 | |
CN111458749B (zh) | 应用于被动源地震勘探的面波与体波的分离方法及系统 | |
CN110646841B (zh) | 时变稀疏反褶积方法及系统 | |
CN105866838B (zh) | 一种地震数据低频信息补偿方法及装置 | |
CN103784164A (zh) | 超声信号的预处理方法及系统 | |
CN105093328B (zh) | 一种滑动扫描谐波压制方法及装置 | |
CN109765608A (zh) | 一种基于联合字典的煤层巷道锚杆振动噪声压制方法 | |
CA2497296A1 (en) | Removal of noise from seismic data using improved radon transformations | |
CN104156925A (zh) | 对超声图像去除散斑和边界增强的处理方法及系统 | |
CN113009569B (zh) | 一种地震偏移成像方法及装置 | |
CN109212590B (zh) | 一种保幅的逆时偏移低频噪音压制方法及系统 | |
US10338250B2 (en) | Method of removing incoherent noise | |
CN113820742A (zh) | 一种粘声各向异性介质中的成像方法 | |
CN109459788B (zh) | 地层品质因子计算方法及系统 | |
CN117849875B (zh) | 一种地震信号分析方法、系统、装置及存储介质 | |
CN113158797B (zh) | 一种微震资料去噪方法、系统、电子设备及存储介质 | |
Yegnanarayana et al. | Signal reconstruction from partial data for sensor array imaging applications | |
CN112666609B (zh) | 基于稀疏包络的地震数据低频补偿方法及装置 | |
CN118210038A (zh) | 求取地震反射系数的方法、电子设备及介质 | |
CN112051612A (zh) | 一种消除地震数据随机噪声的方法和系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |