CN102338886A - 一种有效衰减三分量地震记录中面波的极化滤波方法 - Google Patents
一种有效衰减三分量地震记录中面波的极化滤波方法 Download PDFInfo
- Publication number
- CN102338886A CN102338886A CN2010102314787A CN201010231478A CN102338886A CN 102338886 A CN102338886 A CN 102338886A CN 2010102314787 A CN2010102314787 A CN 2010102314787A CN 201010231478 A CN201010231478 A CN 201010231478A CN 102338886 A CN102338886 A CN 102338886A
- Authority
- CN
- China
- Prior art keywords
- ground roll
- component
- data
- dominant frequency
- component data
- 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.)
- Granted
Links
Images
Abstract
本发明是地震资料处理有效衰减三分量地震记录中面波的极化滤波方法。先确定面波的分布范围,对三分量资料低通滤波;利用垂直分量与水平分量存在90度的相位差进行三分量资料的希尔伯特变换,消除垂直分量与水平分量的相位差;再通过希尔伯特变换构成复数道,复数道奇异值分解得到主本征图像,然后利用面波与有效波在频带和能量上的差异约束初始模型,得到可靠的面波模型;最后从原始记录中自适应减去,实现了面波的有效衰减。本发明方法实现简便、稳定且适应性强。
Description
技术领域
本发明涉及地震资料处理技术,是一种有效衰减三分量地震记录中面波的极化滤波方法。
背景技术
陆上多分量地震勘探是单点三分量数字检波器接收,能真实的记录地下全波场信息,充分的提供地震全波场资料。由于陆上地震勘探近地表环境变化复杂,数字检波器接收动态范围大、单点接收野外无组合和真实的记录低频和高频的波场信息,野外采集的多分量资料一般信噪比较低。面波是陆上地震勘探无法回避的一种强能量低频干扰,在三分量地震勘探资料中的垂直分量(Z)和两个相互垂直的水平分量(X、Y)上随时间变化,呈椭圆偏振特征,而三分量记录中的纵波和转换波在三分量记录上却呈线性偏振。面波在三分量地震记录上与常规的勘探(纵波勘探)记录上一样,与有效信号相比具有低频、低速和能量强的特征,只是三分量勘探资料完整的记录了面波的波场信息,在地震记录上呈“扫帚”状分布,有效反射信号淹没其中,特别是对于水平分量记录的转换波,其视速度较反射纵波低,大量转换波落在面波区。如果不对面波加以有效压制,则对转换波的速度分析以及静校正的计算造成很大的影响。即使在三分量的垂直分量上,面波的能量也较常规勘探野外组合记录的面波能量强。
常规压制的面波方法,比如:区域滤波、叠前相干滤波、局域FK等,应用于三分量数据处理上是有一定局限的。高通滤波能一定程度的压制面波,但同时也损失了一部分有效信号,特别是转换波,其主频较低。F-K滤波理论上能有效的衰减视速度与有效信号差异的线性干扰,但面波在单个分量记录上能量变化大,频散特征明显,视速度变化范围大,直接在单个分量上应用F-K滤波,效果往往不甚理想。在常规的单个分量处理中,为了在压制面波干扰的同时,最大程度的保护有效信号。区域滤波法、自适应面波能量衰减法和小波变换面波压制法等已取得很好的效果,但它们都是真对常规的纵波资料的,应用到三分量资料上仍有一定的局限性。
发明内容
本发明提供了一种有效衰减三分量记录中面波的极化滤波方法,该方法只是在单点三分量地震资料上处理,不受空间假频的影响,能有效衰减三分量地震资料中每个分量上的面波干扰。
本发明通过以下技术方案实现,具体步骤包括:
1)在野外采集三分量地震数据,经预处理,定义面波区,确定面波主频和终止频率,有效波主频和有效波带宽,滤波门槛值,极化分析时窗长度,极化分析时窗重叠百分比;
步骤1)所述的预处理是将野外采集的三分量地震数据加载观测系统。
步骤1)所述的定义面波区是通过面波最大视速度定义随偏移距变化的面波起始时间,或直接定义随偏移距变化的面波起始时间。
步骤1)所述的定确定面波主频,有效波主频是指利用频谱分析确定三分量中的面波主频和终止频率,有效波主频,滤波门槛值是利用频谱分析确定三分量中面波主频带和有效波主频带内的能量比。
2)利用面波主频和面波终止频率确定面波滤波器,利用有效波主频和有效波带宽确定带通滤波器;
步骤2)所述的利用面波主频和面波终止频率确定滤波器的面波滤波器是低通滤波器也可是带通滤波器。
3)输入一单点三分量地震数据,对三分量地震数据应用面波滤波器滤波,得到面波三分量数据,对三分量地震数据应用带通滤波器滤波得到有效波三分量数据;
4)对步骤3)的面波三分量数据进行90度相移,得到相移后的面波三分量数据;
步骤2)所述的进行90度相移是指对三分量数据中的垂直分量进行90度相移,
上述的垂直分量进行90度相移是用希尔伯特一次变换实现,水平分量不变;或者,
对三分量数据中的水平分量进行-90度相移,应用希尔伯特一次变换实现,垂直分量不变。
5)对相移后的面波三分量数据,进行希尔伯特变换,得到面波复数三分量数据;
6)对面波复数三分量数据,自面波起始时间,按时窗长度,重叠的滑动时窗进行复数奇异值分解,每个时窗取奇异值分解的主本征图像得到三分量模型数据,时窗重叠部分线性插值;
步骤6)所述的每个时窗取奇异值分解的本征图像是取可由复数奇异值分解得到的第一本征图像,或第一本征图像与第二本征图像之和。
7)对三分量模型数据进行与步骤4)相对应的反相移,得到面波模型三分量数据;
8)利用以下公式计算面波区每个时间样点的面波模型与有效信号的能量比值PATIO(t);
式中,MX、MY和MZ是面波模型三分量数据,RX、RY和RZ是有效波三分量数据,符号H表示希尔伯特变换,GF是面波主频,EF是面波终止频率,RF是有效波主频,RB是有效波带宽;
9)若时间样点的能量比值小于滤波门槛值,将面波模型三分量数据的样点的值置为0,得到最终的面波模型三分量数据;
10)从原始三分量地震数据中减去最终的面波模型三分量数据,实现了面波衰减。
步骤10)从原始三分量地震数据中减去是指按最小平方准则自适应减去。
本发明是根据面波的椭圆极化特性,通过构建复数道矩阵和对该矩阵SVD分解,求取特征值与特征向量,得到本征图像,适于计算椭圆极化的面波,计算结果稳定,利于识别和剔除面波,提高资料的信噪比;选择合适的滑动重叠时窗,可以大大的减少计算量,并保持了分析精度。
本发明利用面波椭圆偏振特性,去除面波的方法不受空间采样间隔的限制,避免了空间假频的影响,既适用于高密度的单点接收的多分量记录也适用于大道距的多分量记录。针对面波区处理,有效滤除面波后,矢量保真性好,面波区有效信号得到恢复,振幅保持好,更有利于后续的成像和反演。
附图说明
图1是利用极化滤波衰减面波前三分量记录;
图2是利用极化滤波衰减面波后三分量记录。
具体实施方式
本发明提供了一种利用极化滤波有效衰减陆上多分量勘探资料中面波的方法,综合利用了面波相对有效信号低频、低速、强能量和空间椭圆偏振的特点。
先确定面波的分布范围,对三分量资料低通滤波;恰当利用垂直分量与水平分量存在90度的相位差的特征,进行三分量资料的希尔伯特变换,消除垂直分量与水平分量存在的相位差;再通过希尔伯特变换构成复数道,复数道奇异值分解得到主本征图像,即面波初始模型;然后利用面波与有效波在频带和能量上的差异约束初始模型,得到可靠的面波模型;最后从原始记录中自适应减去就实现了面波的有效衰减。
本发明的实现原理如下:
三分量数据的复矩阵奇异值分解
对正交的三分量D(t)=(x(t),y(t),z(t))数据,其希尔波特变换HD(t)=(hx(t),hy(t),hz(t)),构成复数道S(t)=D(t)+jHD(t),进行奇异值分解:
式中,U是地震数据在记录坐标轴上的投影旋度,地震数据的量度由特征值控制;W是一个对角矩阵,元素σi是特征值,且σ1≥σ2≥σ3;V是极化主平面的矢量vi;Ei称为本征图像,E1,E2和E3分别代表了地震数据在极化主平面上的分量,由于面波能量强,E1基本代表了面波的主能量分量,这样在面波主频带内奇异值分解,利用E1建立面波模型,就能有效的衰减三分量记录中的面波。
本发明在采集的多分量资料上采用以下步骤应用:
1)在野外采集三分量地震数据(X、Y和Z),经预处理,采用面波最大视速度1200m/s定义面波区,确定面波主频为10Hz和终止频率是15Hz,有效波主频为30Hz和带宽是10Hz,滤波门槛值是1.5,分析时窗长度是0.1s,时窗重叠百分比为25;
2)利用面波主频和终止频率确定滤波器,利用有效波主频和带宽确定带通滤波器;
3)输入一单点三分量地震数据,对三分量地震数据应用面波滤波器滤波,得到面波三分量数据,对三分量地震数据应用带通滤波器滤波得到有效波三分量数据;
4)对步骤3)应用面波滤波器滤波,得到的面波三分量数据中的垂直分量进行90度相移,得到相移后的面波三分量数据;
5)对相移后的面波三分量数据,进行希尔伯特变换,得到面波复数三分量数据;
6)对面波复数三分量数据,自面波起始时间,按时窗长度,重叠的滑动时窗进行复数奇异值分解,每个时窗取奇异值分解的主本征图像得到三分量模型数据,时窗重叠部分线性插值;
步骤6)所述的每个时窗取奇异值分解的本征图像是指取可由复数奇异值分解得到的第一本征图像。
7)对三分量模型数据进行与步骤4)相对应的反相移,得到面波模型三分量数据;
8)利用公式(1)计算面波区每个时间样点的面波模型与有效信号的能量比值;
步骤8)中的公式(1):
能量比值:
式中,MX、MY和MZ是面波模型三分量数据,RX、RY和RZ是有效波三分量数据,符号H表示希尔伯特变换,GF是面波主频,EF是面波终止频率,RF是有效波主频,RB是有效波带宽,t是样点时间。
9)若时间样点的能量比值小于滤波门槛值,将面波模型三分量数据的该样点的值置为0,这样通过能量约束,保护了有效信号,得到最终的面波模型三分量数据;
10)从原始三分量地震数据中减去最终的面波模型三分量数据,即实现了面波的有效衰减。
图1从左至右分别是未进行面波衰减的X分量,Y分量和Z分量记录,图2是按本方法有效衰减面波后的三分量记录。三个分量上的面波都得到了很好的衰减,特别是面波区界定的明显,没有影响到面波区以外的有效信号,面波区的有效信号恢复出来,并确保三分量记录的矢量保真性。
本发明在单点三分量上应用,不受空间假频的影响,不用准备专门的道集,在多分量处理流程上也不受限制,方法实现简便、稳定且适应性强。
Claims (10)
1.一种有效衰减三分量记录中面波的极化滤波方法,特征是通过以下技术方案实现,具体步骤包括:
1)在野外采集三分量地震数据,经预处理后定义面波区,确定面波主频和终止频率,有效波主频和有效波带宽,滤波门槛值,极化分析时窗长度,极化分析时窗重叠百分比;
2)利用面波主频和面波终止频率确定面波滤波器,利用有效波主频和有效波带宽确定带通滤波器;
3)输入一单点三分量地震数据,对三分量地震数据应用面波滤波器滤波,得到面波三分量数据,对三分量地震数据应用带通滤波器滤波得到有效波三分量数据;
4)对步骤3)的面波三分量数据进行90度相移,得到相移后的面波三分量数据;
5)对相移后的面波三分量数据,进行希尔伯特变换,得到面波复数三分量数据;
6)对面波复数三分量数据,自面波起始时间,按时窗长度,重叠的滑动时窗进行复数奇异值分解,每个时窗取奇异值分解的主本征图像得到三分量模型数据,时窗重叠部分线性插值;
7)对三分量模型数据进行与步骤4)相对应的反相移,得到面波模型三分量数据;
8)利用以下公式计算面波区每个时间样点的面波模型与有效信号的能量比值PATIO(t);
式中,MX、MY和MZ是面波模型三分量数据,RX、RY和RZ是有效波三分量数据,符号H表示希尔伯特变换,GF是面波主频,EF是面波终止频率,RF是有效波主频,RB是有效波带宽;
9)若时间样点的能量比值小于滤波门槛值,将面波模型三分量数据的样点的值置为0,得到最终的面波模型三分量数据;
10)从原始三分量地震数据中减去最终的面波模型三分量数据,实现面波衰减。
2.根据权利要求1所述的方法,特征是步骤1)所述的预处理是将野外采集的三分量地震数据加载观测系统。
3.根据权利要求1所述的方法,特征是步骤1)所述的定义面波区是通过面波最大视速度定义随偏移距变化的面波起始时间,或直接定义随偏移距变化的面波起始时间。
4.根据权利要求1所述的方法,特征是步骤1)所述的确定面波主频,有效波主频是指利用频谱分析确定三分量中的面波主频和终止频率,有效波主频,滤波门槛值是利用频谱分析确定三分量中面波主频带和有效波主频带内的能量比。
5.根据权利要求1所述的方法,特征是步骤2)所述的利用面波主频和面波终止频率确定滤波器的面波滤波器是低通滤波器也可是带通滤波器。
6.根据权利要求1所述的方法,特征是步骤4)所述的进行90度相移是指对三分量数据中的垂直分量进行90度相移,
7.根据权利要求1或6所述的方法,特征是上述的垂直分量进行90度相移是用希尔伯特一次变换实现,水平分量不变;或者,
8.根据权利要求1所述的方法,特征是对三分量数据中的水平分量进行-90度相移,应用希尔伯特一次变换实现,垂直分量不变。
9.根据权利要求1所述的方法,特征是步骤6)所述的每个时窗取奇异值分解的本征图像是取可由复数奇异值分解得到的第一本征图像,或第一本征图像与第二本征图像之和。
10.根据权利要求1所述的方法,特征是步骤10)从原始三分量地震数据中减去是指按最小平方准则自适应减去。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010231478 CN102338886B (zh) | 2010-07-14 | 2010-07-14 | 一种有效衰减三分量地震记录中面波的极化滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010231478 CN102338886B (zh) | 2010-07-14 | 2010-07-14 | 一种有效衰减三分量地震记录中面波的极化滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102338886A true CN102338886A (zh) | 2012-02-01 |
CN102338886B CN102338886B (zh) | 2013-09-04 |
Family
ID=45514713
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010231478 Active CN102338886B (zh) | 2010-07-14 | 2010-07-14 | 一种有效衰减三分量地震记录中面波的极化滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102338886B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102944895A (zh) * | 2012-10-26 | 2013-02-27 | 中国石油大学(北京) | 一种地震记录面波干扰压制方法及装置 |
CN103869357A (zh) * | 2012-12-13 | 2014-06-18 | 中国石油集团东方地球物理勘探有限责任公司 | 一种地震数据频谱奇异性的检测方法与装置 |
CN104199104A (zh) * | 2014-08-15 | 2014-12-10 | 中国石油天然气集团公司 | 三维地震数据面波衰减方法及装置 |
CN105652322A (zh) * | 2016-01-07 | 2016-06-08 | 中国科学院地球化学研究所 | 多分量地震数据的t-f-k域极化滤波方法 |
CN111856555A (zh) * | 2020-06-19 | 2020-10-30 | 同济大学 | 一种基于面波多尺度窗分析的地下探测方法 |
US11686872B2 (en) | 2019-12-23 | 2023-06-27 | Saudi Arabian Oil Company | Attenuation of guided waves using polarization filtering |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6278949B1 (en) * | 1998-11-25 | 2001-08-21 | M. Aftab Alam | Method for multi-attribute identification of structure and stratigraphy in a volume of seismic data |
US20040098199A1 (en) * | 2002-11-19 | 2004-05-20 | Yi Luo | Seismic data processing method to enhance fault and channel display |
CN1981212A (zh) * | 2004-06-30 | 2007-06-13 | Pgs美洲公司 | 利用复变道分集滤波器来衰减地震数据中的噪声的方法 |
US7463552B1 (en) * | 2003-11-25 | 2008-12-09 | Michael John Padgett | Method for deriving 3D output volumes using filters derived from flat spot direction vectors |
CN101609160A (zh) * | 2008-06-19 | 2009-12-23 | 中国石油天然气集团公司 | 一种抗背景噪音强的地震数据初至波自动拾取方法 |
CN101644782A (zh) * | 2009-08-25 | 2010-02-10 | 中国石化集团胜利石油管理局 | 基于极化滤波的多波多分量地震资料的去噪方法 |
-
2010
- 2010-07-14 CN CN 201010231478 patent/CN102338886B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6278949B1 (en) * | 1998-11-25 | 2001-08-21 | M. Aftab Alam | Method for multi-attribute identification of structure and stratigraphy in a volume of seismic data |
US20040098199A1 (en) * | 2002-11-19 | 2004-05-20 | Yi Luo | Seismic data processing method to enhance fault and channel display |
US7463552B1 (en) * | 2003-11-25 | 2008-12-09 | Michael John Padgett | Method for deriving 3D output volumes using filters derived from flat spot direction vectors |
CN1981212A (zh) * | 2004-06-30 | 2007-06-13 | Pgs美洲公司 | 利用复变道分集滤波器来衰减地震数据中的噪声的方法 |
CN101609160A (zh) * | 2008-06-19 | 2009-12-23 | 中国石油天然气集团公司 | 一种抗背景噪音强的地震数据初至波自动拾取方法 |
CN101644782A (zh) * | 2009-08-25 | 2010-02-10 | 中国石化集团胜利石油管理局 | 基于极化滤波的多波多分量地震资料的去噪方法 |
Non-Patent Citations (2)
Title |
---|
贺振华: "用二维付氏变换及希氏变换作倾斜滤波", 《物探化探计算技术》, vol. 7, no. 03, 30 September 1985 (1985-09-30), pages 234 - 243 * |
陈?S等: "基于加窗Hilbert变换的复偏振分析方法及其应用", 《地球物理学报》, vol. 48, no. 04, 17 July 2005 (2005-07-17), pages 889 - 895 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102944895A (zh) * | 2012-10-26 | 2013-02-27 | 中国石油大学(北京) | 一种地震记录面波干扰压制方法及装置 |
CN102944895B (zh) * | 2012-10-26 | 2015-04-15 | 中国石油大学(北京) | 一种地震记录面波干扰压制方法及装置 |
CN103869357A (zh) * | 2012-12-13 | 2014-06-18 | 中国石油集团东方地球物理勘探有限责任公司 | 一种地震数据频谱奇异性的检测方法与装置 |
CN104199104A (zh) * | 2014-08-15 | 2014-12-10 | 中国石油天然气集团公司 | 三维地震数据面波衰减方法及装置 |
CN104199104B (zh) * | 2014-08-15 | 2017-03-08 | 中国石油天然气集团公司 | 三维地震数据面波衰减方法及装置 |
CN105652322A (zh) * | 2016-01-07 | 2016-06-08 | 中国科学院地球化学研究所 | 多分量地震数据的t-f-k域极化滤波方法 |
US11686872B2 (en) | 2019-12-23 | 2023-06-27 | Saudi Arabian Oil Company | Attenuation of guided waves using polarization filtering |
CN111856555A (zh) * | 2020-06-19 | 2020-10-30 | 同济大学 | 一种基于面波多尺度窗分析的地下探测方法 |
CN111856555B (zh) * | 2020-06-19 | 2021-09-03 | 同济大学 | 一种基于面波多尺度窗分析的地下探测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102338886B (zh) | 2013-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102338886B (zh) | 一种有效衰减三分量地震记录中面波的极化滤波方法 | |
Langston | Wave gradiometry in two dimensions | |
Voisin et al. | Seismic noise monitoring of the water table in a deep-seated, slow-moving landslide | |
CN102681014B (zh) | 基于多项式拟合的规则线性干扰压制方法 | |
CN101915939B (zh) | 一种面波压制方法 | |
CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
Cheng et al. | Automated data selection in the Tau–p domain: Application to passive surface wave imaging | |
CN104614769B (zh) | 一种压制地震面波的聚束滤波方法 | |
CN104808245A (zh) | 道集优化处理方法及其装置 | |
CN106547020B (zh) | 一种地震数据的保幅处理方法 | |
CN103698807A (zh) | 梯状二维宽线观测系统设计方法 | |
CN109884709B (zh) | 一种基于面波旅行时层析的转换波静校正方法 | |
CN103645497A (zh) | 一种基于出射角的多分量波场分离方法 | |
CN106707334A (zh) | 一种提高地震资料分辨率的方法 | |
Meng et al. | Analysis of seismic signals generated by vehicle traffic with application to derivation of subsurface Q‐values | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
CN104635264B (zh) | 叠前地震数据的处理方法及设备 | |
CN106950600B (zh) | 一种近地表散射面波的去除方法 | |
Zhang et al. | Retrieval of shallow S-wave profiles from seismic reflection surveying and traffic-induced noise | |
CN105510975A (zh) | 提高地震数据信噪比的方法及装置 | |
CN103412325A (zh) | 一种地震纵、横波波场分离与去噪的方法 | |
Askari et al. | Estimation of S-wave static corrections using CMP cross-correlation of surface waves | |
Rosa-Cintas et al. | Polarization analysis in the stationary wavelet packet domain: Application to HVSR method | |
CN103558636A (zh) | 一种从叠后地震数据采集脚印衰减的方法 | |
Ker et al. | Multiscale seismic attributes: a wavelet-based method and its application to high-resolution seismic and ground truth data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |