CN107144880A - 一种地震波波场分离方法 - Google Patents

一种地震波波场分离方法 Download PDF

Info

Publication number
CN107144880A
CN107144880A CN201710335916.6A CN201710335916A CN107144880A CN 107144880 A CN107144880 A CN 107144880A CN 201710335916 A CN201710335916 A CN 201710335916A CN 107144880 A CN107144880 A CN 107144880A
Authority
CN
China
Prior art keywords
wave
seismic
signal
frequency
field separation
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
Application number
CN201710335916.6A
Other languages
English (en)
Other versions
CN107144880B (zh
Inventor
柴贺军
黄河
阎宗岭
贾学明
毛友平
张小松
徐建强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Merchants Chongqing Communications Research and Design Institute Co Ltd
Original Assignee
China Merchants Chongqing Communications Research and Design Institute Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Merchants Chongqing Communications Research and Design Institute Co Ltd filed Critical China Merchants Chongqing Communications Research and Design Institute Co Ltd
Priority to CN201710335916.6A priority Critical patent/CN107144880B/zh
Publication of CN107144880A publication Critical patent/CN107144880A/zh
Application granted granted Critical
Publication of CN107144880B publication Critical patent/CN107144880B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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. analysis, for interpretation, for correction
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • 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. analysis, for interpretation, for correction
    • 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. analysis, for interpretation, for correction
    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/30—Noise handling
    • G01V2210/32—Noise reduction
    • G01V2210/322—Trace stacking
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/30—Noise handling
    • G01V2210/32—Noise reduction
    • G01V2210/324—Filtering
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/40—Transforming data representation
    • G01V2210/44—F-k domain
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/40—Transforming data representation
    • G01V2210/46—Radon transform
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/40—Transforming data representation
    • G01V2210/47—Slowness, e.g. tau-pi

Abstract

本发明涉及一种地震波波场分离方法，采用检波器、小波滤波器和工程地震仪数据预处理软件获取并分析地震波；包括如下步骤：采用检波器获取初始地震波信号，采用小波滤波器对初始地震波信号进行降噪；将降噪后的地震波信号进行F‑K频域变换和反F‑K频域变换，以及线性Radon变换和反线性Radon变换，分离出地震波信号中的直达波、反射波、折射波和面波信号。有效提高波场分离的精度和可靠性。

Description

一种地震波波场分离方法
技术领域
本发明涉及多波地震勘探数据处理方法，尤其是能对勘探过程中提取到的地震波进行波场分离的方法。
背景技术
地震勘探在能源矿产勘探中的应用得到广泛的认可。随着技术的进步,多波地震勘探技术得到很大的发展。多波地震勘探具有携带丰富地下信息等特点。但是，不同波之间的相互干扰降低了成像的分辨率,从而严重影响其对地质解释的准确性。目前的分离方法主要是依靠地震波的动力学特征(偏振特性)和运动学特征(视速度)对波场进行分离，目前只能分离出横波和纵波。例如CN201310553067.3公开的基于出射角的多分量波场分离方法，同时考虑了地震波的运动学和动力学特征，利用地层速度计算地震波的出射角，通过矢量旋转实现多分量的波场分离，能有效分离折射纵波和横波。
现有的地震波波场分离方法在实际运用中往往不能很好地实现波场分离，主要原因在于：实际的地震记录较为复杂，不仅含有各种噪音，而且实际的地层介质非常复杂，如非均匀性、各向异性、含有裂缝及充填物等。而合成记录不含有噪音，且所用的物理模型比较简单；所用的合成记录的合成与分解往往基于同样的方法原理和简化假设，难以应用于实际研究工作。现代地质勘探需要更多的波场分量，不同的波场分量包含有不同的地址信息，反应更丰富的地下信息，然而现有技术尚未解。
发明内容
针对现有技术中存在的上述不足，本发明的目的是提供一种地震波波场分离方法，解决现有技术不能满足多波场分离的需求问题，将地震波信号分离为直达波、面波、反射波、折射波，满足实际运用中的对地震波进行波场分离的需要。
实现上述目的，本发明采用如下技术方案：一种地震波波场分离方法，采用检波器和小波滤波器获取和分析地震波；具体包括如下步骤：
步骤1：采用检波器获取初始地震波信号，采用第一小波滤波器对初始地震波信号进行降噪；
步骤2：将降噪后的地震波信号进行F-K频域变换，在F-K域进行波场分离，利用面波低频特性，设计相应的低通滤波器，滤出地震波信号中的面波；剩下含有反射波，直达波与折射波的地震波信号；
步骤3：将步骤2得到包含反射波、直达波与折射波的地震波信号经第二小波滤波器对该地震波信号进行降噪；
步骤4：将步骤3降噪后的地震波信号进行F-K频域变换，在F-K域进行波场分离，利用频率特性，利用直达波速度一定，在F-K域分布固定分离出了地震波信号中的直达波，并得到包含反射波与折射波的地震波信号；
步骤5：将步骤4得到包含反射波与折射波的地震波信号采用第三小波滤波器器对该地震波信号进降噪；
步骤6：将步骤5降噪后的地震波信号进行线性Radon变换，τ-p变换将反射波与折射波进行切分，分离出地震波信号中的反射波和折射波信号。
进一步，所述步骤2和步骤在4F-K域进行波场分离，采用频率波数域的二维视速度方法进行波场分离。所述步骤2切除地震波信号中的面波是以面波频率为4Hz左右，而其他有效地震波信号远远高于该频率，设计相应的低通滤波器，可滤出面波。
进一步，所述步骤4切除掉直达波频率以直达波速度为基准；由于直达波速度固定，在F-K域中能量分布在一定的区域上，而散射波和折射波随着信号源的不同来向而变化，在F-K域中是散点状，不在固定的分布区域；直达波信号，在地震波数据显示中为一条直线，可算出直达波速度V，以直达波速度为基准，设定滤波范围。滤除直达波，剩下含有反射波和折射波的地震波信号；再进行F-K反变换。
相比现有技术，本发明具有如下有益效果：
1、本发明首先降低地震波信号的噪声信号，然后将多波地震勘探包含众多信息进行分离，分离出包括直达波、面波、反射波以及折射波,满足现代地质勘探需要更多的波场分量，反应更丰富的地下信息。再进行后续的数据处理，有效提高波场分离的精度和可靠性。
2、本发明方法具有巨大的优势。在现有的技术中，波长分离技术主要分离出的是纵波与横波，在本发明将分离精度进一步提高，最终分离出面波、直达波、反射波、折射波，使之后的数据分析有更多的地震资料，能够更准确地分析出地层结构。
3、该波场分离技术易于实现，可广泛应用于其他领域。
附图说明
图1为本发明地震波波场分离方法的流程方框图。
图2为τ-p变换图。
图3-图7为地震波原始波形图、分离出反射波、直达波、折射波和面波图。
具体实施方式
下面结合附图和具体实施方式对本发明的技术方案作进一步的说明。
如图1所述，一种地震波波场分离方法，采用检波器和小波滤波器获取和分析地震波；具体包括如下步骤：
步骤1：采用检波器获取初始地震波信号，采用第一小波滤波器对初始地震波信号进行降噪；
步骤2：将降噪后的地震波信号进行F-K频域变换，在F-K域进行波场分离，利用频率特性，分离是利用频率特性，采用落锤法获得震源信号(地震波源信号从不同角度进入地下)，分出地震波信号中的面波，并得到包含反射波、直达波与折射波的地震波信号；
步骤3：将步骤2得到包含反射波、直达波与折射波的地震波信号经第二小波滤波器对该地震波信号进行降噪；
步骤4：将步骤3降噪后的地震波信号进行F-K频域变换，在F-K域进行波场分离，利用频率特性，分离出了地震波信号中的直达波，并得到包含反射波与折射波的地震波信号；
步骤5：将步骤4得到包含反射波与折射波的地震波信号采用第三小波滤波器器对该地震波信号进降噪；
步骤6：将步骤5降噪后的地震波信号进行线性Radon变换，τ-p变换将反射波与折射波进行切分，分离出地震波信号中的反射波和折射波信号。
τ-p变换过程中，在做倾斜叠加时会产生大量的假频，亦称作周波跳跃，是一种很强的干扰，严重影响τ-p变换的质量。假频是指假的频率成分，它是由于采样不足，使得信号中高于生采样频率的频率成分折叠到低频区域中而形成的低频成分。因此，再一次用小波滤波的方法消除假频现象。通过小波滤波的方法对含噪信号进行处理，对其进行小波分解，选择小波确定层数N，对信号进行N层分解。对分层信号进行阈值量化，最后根据小波分解的第N层低频系数和第一到第N层的高频系数进行一维重构。通过选择合适的小波函数和阈值实现小波降噪。
参见图2，t-p变换原理：线性Radon变换(或τ|p变换)，是将数据沿着一系列倾斜线t＝τ+px进行叠加，τ为截距时间，p为射线参数。经过正变换后，x-t域的线性同相轴、双曲线形态的反射同相轴分别被映射为τ|p域的点和椭圆。根据面波、反射波、折射波的不同特性，将不同波的传播特性经过Radon变换到τ-p域内，进行分离。在时域交叉在一起的面波，反射波，折射波，在τ-p发生了分离。将x-t域数据变换到τ|p域的更利于地震波场分离。折射波在t-p域为固定点分布，反射波为椭圆分布，并且没有交集，设计τ-p域滤波器τ＝t(t0～t1)。
经过τ-p变换后的直线聚焦成为一个“点”。设计合适的滤波器，将在T-P域的反射波进行分离，在进行t-p反变换，剩下的地震波就是折射波分量。
第一、第二、第三小波滤波器采用了空域相关滤波法。
原理：由于信号与噪声在小波域具有不同的表现，空域相关滤波算法通过取相邻尺度的小波系数相乘进行相关运算，使得在保留信号细节的同时而抑制噪声。空域相关滤波算法是针对信号与噪声具有不同的Lipschitz常数，在小波域具有不同的表现。在不同的尺度上，信号的小波系数有较强的相关性，而噪声则没有这种相关性，其小波系数会随着尺度的增加而迅速减小。因此可以取相邻尺度的小波系数相乘进行相关运算，通过相乘能够锐化信号的细节同时抑制噪声。
采用本发明方法，能将如图3所示的地震波原始波形分离为图4-图7为分离出反射波、直达波、折射波和面波。为后续研究和分析提供参考依据。
面波传播的特征：近震源处发育、震幅大、传播速度低。面波频率较低，一般为4Hz左右，而其他波分量的有效波成分主要集中在60-100HZ。所以f-k变换后，面波分量能够得到充分剥离。使用低通滤波器，即可分离出面波。
直达波是地震波中的一种干扰波，地震波信号进行f-k变换后。在频域内设计扇形滤波器，能够达到去除直达波的目的。
最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (4)

1.地震波波场分离方法，采用检波器和小波滤波器获取和分析地震波；具体包括如下步骤：
步骤1：采用检波器获取初始地震波信号，采用第一小波滤波器对初始地震波信号进行降噪；
步骤2：将降噪后的地震波信号进行F-K频域变换，在F-K域进行波场分离，利用面波低频特性，滤出地震波信号中的面波；剩下含有反射波，直达波与折射波的地震波信号；
步骤3：将步骤2得到包含反射波、直达波与折射波的地震波信号经第二小波滤波器对该地震波信号进行降噪；
步骤4：将步骤3降噪后的地震波信号进行F-K频域变换，在F-K域进行波场分离，利用频率特性，利用直达波速度一定，在F-K域分布固定分离出了地震波信号中的直达波，并得到包含反射波与折射波的地震波信号；
步骤5：将步骤4得到包含反射波与折射波的地震波信号采用第三小波滤波器器对该地震波信号进降噪；
步骤6：将步骤5降噪后的地震波信号进行线性Radon变换，τ-p变换将反射波与折射波进行切分，分离出地震波信号中的反射波和折射波信号。
2.根据权利要求1所述地震波波场分离方法，其特征在于：所述步骤2和步骤在4F-K域进行波场分离，采用频率波数域的二维视速度方法进行波场分离。
3.根据权利要求1所述地震波波场分离方法，其特征在于：所述步骤2切除地震波信号中的面波是以面波频率为4Hz左右，而其他有效地震波信号远远高于该频率，设计相应的低通滤波器，滤出面波。
4.根据权利要求1所述地震波波场分离方法，其特征在于：步骤4切除掉直达波频率以直达波速度为基准；由于直达波速度固定，在F-K域中能量分布在一定的区域上，而散射波和折射波随着信号源的不同来向而变化，在F-K域中是散点状，不在固定的分布区域；直达波信号，在地震波数据显示中为一条直线，可算出直达波速度V，以直达波速度为基准，设定滤波范围。滤除直达波，剩下含有反射波和折射波的地震波信号；再进行F-K反变换。
CN201710335916.6A 2017-05-12 2017-05-12 一种地震波波场分离方法 Active CN107144880B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710335916.6A CN107144880B (zh) 2017-05-12 2017-05-12 一种地震波波场分离方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710335916.6A CN107144880B (zh) 2017-05-12 2017-05-12 一种地震波波场分离方法

Publications (2)

Publication Number Publication Date
CN107144880A true CN107144880A (zh) 2017-09-08
CN107144880B CN107144880B (zh) 2019-07-05

Family

ID=59777414

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710335916.6A Active CN107144880B (zh) 2017-05-12 2017-05-12 一种地震波波场分离方法

Country Status (1)

Country Link
CN (1) CN107144880B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107894613A (zh) * 2017-10-26 2018-04-10 中国石油天然气集团公司 弹性波矢量成像方法、装置、存储介质及设备
CN108363095A (zh) * 2018-05-14 2018-08-03 招商局重庆交通科研设计院有限公司 一种面波法快速检测装置
CN109521471A (zh) * 2017-09-18 2019-03-26 中国石油化工股份有限公司 用于分离纵横波场的数据变换方法及系统
CN110716231A (zh) * 2019-09-17 2020-01-21 中国地质大学（武汉） 基于共聚焦域稀疏反演的海上多震源波场分离方法及系统
CN110780346A (zh) * 2019-11-20 2020-02-11 李志勇 一种隧道超前探测复杂地震波场的分离方法
CN111856551A (zh) * 2020-06-22 2020-10-30 山东电力工程咨询院有限公司 浅层横向高分辨率瑞雷波勘探方法及系统
CN112147690A (zh) * 2019-06-28 2020-12-29 中国石油化工股份有限公司 提高首波带中浅层弱反射信噪比的方法及系统
CN111897001B (zh) * 2020-06-19 2021-04-13 成都理工大学 波场确定方法与装置
CN112698400A (zh) * 2020-12-04 2021-04-23 中国科学院深圳先进技术研究院 反演方法、反演装置、计算机设备和计算机可读存储介质
CN114199489A (zh) * 2021-12-10 2022-03-18 江南大学 基于数字滤波调整天然地震波拟合反应谱的方法和装置
CN117872478A (zh) * 2024-03-13 2024-04-12 中国石油大学(华东) 一种针对阵列声波测井数据的弱反射信号提取方法及装置

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1404582A (zh) * 2000-01-21 2003-03-19 施鲁博格控股有限公司 用于地震波场分离的系统和方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1404582A (zh) * 2000-01-21 2003-03-19 施鲁博格控股有限公司 用于地震波场分离的系统和方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
宁刚: "基于F-K滤波的TSP负反射提取方法研究", 《工程地球物理学报》 *
戴亦军: "多波反射地震勘探数据采集及波场分离技术研究与应用", 《中国优秀博硕士学位论文全文数据库（硕士）·基础科学辑》 *
李彩芹 等: "小波变换与F-K联合滤波在面波分离中的应用", 《中国煤田地质》 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109521471A (zh) * 2017-09-18 2019-03-26 中国石油化工股份有限公司 用于分离纵横波场的数据变换方法及系统
CN107894613B (zh) * 2017-10-26 2019-07-26 中国石油天然气集团公司 弹性波矢量成像方法、装置、存储介质及设备
CN107894613A (zh) * 2017-10-26 2018-04-10 中国石油天然气集团公司 弹性波矢量成像方法、装置、存储介质及设备
CN108363095B (zh) * 2018-05-14 2023-07-18 招商局重庆交通科研设计院有限公司 一种面波法快速检测装置
CN108363095A (zh) * 2018-05-14 2018-08-03 招商局重庆交通科研设计院有限公司 一种面波法快速检测装置
CN112147690A (zh) * 2019-06-28 2020-12-29 中国石油化工股份有限公司 提高首波带中浅层弱反射信噪比的方法及系统
CN110716231A (zh) * 2019-09-17 2020-01-21 中国地质大学（武汉） 基于共聚焦域稀疏反演的海上多震源波场分离方法及系统
CN110780346A (zh) * 2019-11-20 2020-02-11 李志勇 一种隧道超前探测复杂地震波场的分离方法
CN111897001B (zh) * 2020-06-19 2021-04-13 成都理工大学 波场确定方法与装置
CN111856551A (zh) * 2020-06-22 2020-10-30 山东电力工程咨询院有限公司 浅层横向高分辨率瑞雷波勘探方法及系统
CN112698400B (zh) * 2020-12-04 2023-06-23 中国科学院深圳先进技术研究院 反演方法、反演装置、计算机设备和计算机可读存储介质
CN112698400A (zh) * 2020-12-04 2021-04-23 中国科学院深圳先进技术研究院 反演方法、反演装置、计算机设备和计算机可读存储介质
CN114199489A (zh) * 2021-12-10 2022-03-18 江南大学 基于数字滤波调整天然地震波拟合反应谱的方法和装置
CN114199489B (zh) * 2021-12-10 2023-05-09 江南大学 基于数字滤波调整天然地震波拟合反应谱的方法和装置
CN117872478A (zh) * 2024-03-13 2024-04-12 中国石油大学(华东) 一种针对阵列声波测井数据的弱反射信号提取方法及装置

Also Published As

Publication number Publication date
CN107144880B (zh) 2019-07-05

Similar Documents

Publication Publication Date Title
CN107144880B (zh) 一种地震波波场分离方法
Mousavi et al. Automatic noise-removal/signal-removal based on general cross-validation thresholding in synchrosqueezed domain and its application on earthquake data
Akram et al. A review and appraisal of arrival-time picking methods for downhole microseismic data
Neelamani et al. Coherent and random noise attenuation using the curvelet transform
Yu et al. Attenuation of noise and simultaneous source interference using wavelet denoising
Bekara et al. Random and coherent noise attenuation by empirical mode decomposition
US7953556B2 (en) Geophone noise attenuation and wavefield separation using a multi-dimensional decomposition technique
CN106932824B (zh) 陆地地震勘探资料的降维自适应层间多次波压制方法
CN102939546A (zh) 用于在地震处理中的局部属性匹配的系统和方法
CN104007469A (zh) 一种基于曲波变换的弱地震信号重构方法
CN106896409B (zh) 一种基于波动方程边值反演的变深度缆鬼波压制方法
CN104020492A (zh) 一种三维地震资料的保边滤波方法
CN105652322A (zh) 多分量地震数据的t-f-k域极化滤波方法
CN102590859A (zh) Vti介质准p波方程各向异性逆时偏移方法
CN105445801B (zh) 一种消除二维地震资料随机噪音的处理方法
CN107884829A (zh) 一种联合压制浅海obc地震资料多次波的方法
Chopra et al. High-frequency restoration of surface seismic data
CN106526678A (zh) 一种反射声波测井的波场分离方法及装置
CN105182417A (zh) 一种基于形态成分分析的面波分离方法及系统
CN107390270A (zh) 一种基于弹性波逆时偏移ADCIGs的AVA分析方法
Liu et al. A review of variational mode decomposition in seismic data analysis
CN106950600A (zh) 一种近地表散射面波的去除方法
Huang et al. P/S-wave separation of multicomponent seismic data at the land surface based on deep learning
CN103558636A (zh) 一种从叠后地震数据采集脚印衰减的方法
CN104570114A (zh) 一种基于波场分解的逆时偏移噪音压制方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant