CN107144880B - 一种地震波波场分离方法 - Google Patents
一种地震波波场分离方法 Download PDFInfo
- Publication number
- CN107144880B CN107144880B CN201710335916.6A CN201710335916A CN107144880B CN 107144880 B CN107144880 B CN 107144880B CN 201710335916 A CN201710335916 A CN 201710335916A CN 107144880 B CN107144880 B CN 107144880B
- Authority
- CN
- China
- Prior art keywords
- wave
- seismic
- signal
- seismic signal
- domain
- 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
- 238000000926 separation method Methods 0.000 title claims abstract description 32
- 230000009467 reduction Effects 0.000 claims abstract description 20
- 230000009466 transformation Effects 0.000 claims abstract description 12
- 229910052704 radon Inorganic materials 0.000 claims abstract description 7
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims abstract description 7
- 238000000034 method Methods 0.000 claims description 16
- 238000001914 filtration Methods 0.000 claims description 7
- 238000000354 decomposition reaction Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000001615 p wave Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000003786 synthesis reaction Methods 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
-
- 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
-
- 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
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
本发明涉及一种地震波波场分离方法,采用检波器、小波滤波器和工程地震仪数据预处理软件获取并分析地震波;包括如下步骤:采用检波器获取初始地震波信号,采用小波滤波器对初始地震波信号进行降噪;将降噪后的地震波信号进行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和步骤4在F-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,τ-p变换原理:线性Radon变换(或τ-p变换),是将数据沿着一系列倾斜线t=τ+px进行叠加,τ为截距时间,p为射线参数。经过正变换后,x-t域的线性同相轴、双曲线形态的反射同相轴分别被映射为τ-p域的点和椭圆。根据面波、反射波、折射波的不同特性,将不同波的传播特性经过Radon变换到τ-p域内,进行分离。在时域交叉在一起的面波,反射波,折射波,在τ-p发生了分离。将x-t域数据变换到τ-p域的更利于地震波场分离。折射波在τ-p域为固定点分布,反射波为椭圆分布,并且没有交集,设计τ-p域滤波器τ=t(t0~t1)。
经过τ-p变换后的直线聚焦成为一个“点”。设计合适的滤波器,将在τ-p域的反射波进行分离,在进行τ-p反变换,剩下的地震波就是折射波分量。
第一、第二、第三小波滤波器采用了空域相关滤波法。
原理:由于信号与噪声在小波域具有不同的表现,空域相关滤波算法通过取相邻尺度的小波系数相乘进行相关运算,使得在保留信号细节的同时而抑制噪声。空域相关滤波算法是针对信号与噪声具有不同的Lipschitz常数,在小波域具有不同的表现。在不同的尺度上,信号的小波系数有较强的相关性,而噪声则没有这种相关性,其小波系数会随着尺度的增加而迅速减小。因此可以取相邻尺度的小波系数相乘进行相关运算,通过相乘能够锐化信号的细节同时抑制噪声。
采用本发明方法,能将如图3所示的地震波原始波形分离为图4-图7为分离出反射波、直达波、折射波和面波。为后续研究和分析提供参考依据。
面波传播的特征:近震源处发育、震幅大、传播速度低。面波频率较低,一般为4Hz左右,而其他波分量的有效波成分主要集中在60-100HZ。所以f-k变换后,面波分量能够得到充分剥离。使用低通滤波器,即可分离出面波。
直达波是地震波中的一种干扰波,地震波信号进行f-k变换后。在频域内设计扇形滤波器,能够达到去除直达波的目的。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (1)
1.一种地震波波场分离方法,采用检波器和小波滤波器获取和分析地震波;具体包括如下步骤:
步骤1:采用检波器获取初始地震波信号,采用第一小波滤波器对初始地震波信号进行降噪;
步骤2:将降噪后的地震波信号进行F-K频域变换,在F-K域进行波场分离,利用面波低频特性,设计相应的低通滤波器,滤出地震波信号中的面波;剩下含有反射波,直达波与折射波的地震波信号;切除地震波信号中的面波频率为4Hz;
步骤3:将步骤2得到包含反射波、直达波与折射波的地震波信号经第二小波滤波器对该地震波信号进行降噪;
步骤4:将步骤3降噪后的地震波信号进行F-K频域变换,在F-K域进行波场分离,利用频率特性,利用直达波速度一定,在F-K域分布固定分离出了地震波信号中的直达波,并得到包含反射波与折射波的地震波信号;
其中,切除掉直达波频率以直达波速度为基准,滤除直达波,剩下含有反射波和折射波的地震波信号;再进行F-K反变换;
步骤5:将步骤4得到包含反射波与折射波的地震波信号采用第三小波滤波器对该地震波信号进降噪;
步骤6:将步骤5降噪后的地震波信号进行线性Radon变换,变换将反射波与折射波进行切分,分离出地震波信号中的反射波和折射波信号;
所述步骤2和步骤4中,在F-K域进行波场分离,采用频率波数域的二维视速度方法进行波场分离;
其中,第一、第二、第三小波滤波器采用空域相关滤波法。
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 CN107144880A (zh) | 2017-09-08 |
CN107144880B true 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) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109521471A (zh) * | 2017-09-18 | 2019-03-26 | 中国石油化工股份有限公司 | 用于分离纵横波场的数据变换方法及系统 |
CN107894613B (zh) * | 2017-10-26 | 2019-07-26 | 中国石油天然气集团公司 | 弹性波矢量成像方法、装置、存储介质及设备 |
CN108363095B (zh) * | 2018-05-14 | 2023-07-18 | 招商局重庆交通科研设计院有限公司 | 一种面波法快速检测装置 |
CN112147690A (zh) * | 2019-06-28 | 2020-12-29 | 中国石油化工股份有限公司 | 提高首波带中浅层弱反射信噪比的方法及系统 |
CN110716231B (zh) * | 2019-09-17 | 2020-11-24 | 中国地质大学(武汉) | 基于共聚焦域稀疏反演的海上多震源波场分离方法及系统 |
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 | 中国科学院深圳先进技术研究院 | 反演方法、反演装置、计算机设备和计算机可读存储介质 |
CN114199489B (zh) * | 2021-12-10 | 2023-05-09 | 江南大学 | 基于数字滤波调整天然地震波拟合反应谱的方法和装置 |
CN117872478A (zh) * | 2024-03-13 | 2024-04-12 | 中国石油大学(华东) | 一种针对阵列声波测井数据的弱反射信号提取方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU779802B2 (en) * | 2000-01-21 | 2005-02-10 | Schlumberger Holdings Limited | System and method for seismic wavefield separation |
-
2017
- 2017-05-12 CN CN201710335916.6A patent/CN107144880B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107144880A (zh) | 2017-09-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107144880B (zh) | 一种地震波波场分离方法 | |
Xue et al. | Amplitude-preserving iterative deblending of simultaneous source seismic data using high-order radon transform | |
Yu et al. | Attenuation of noise and simultaneous source interference using wavelet denoising | |
AU2011223757B2 (en) | System and method for local attribute matching in seismic processing | |
Chen et al. | Random noise reduction using a hybrid method based on ensemble empirical mode decomposition | |
US7616524B1 (en) | Wavelet based intercept attribute for seismic exploration | |
CN106896409B (zh) | 一种基于波动方程边值反演的变深度缆鬼波压制方法 | |
CN104007469A (zh) | 一种基于曲波变换的弱地震信号重构方法 | |
CN105549076B (zh) | 一种基于交替方向法和全变分理论的地震数据处理方法 | |
CN105445801B (zh) | 一种消除二维地震资料随机噪音的处理方法 | |
CN105182417B (zh) | 一种基于形态成分分析的面波分离方法及系统 | |
CN105652322A (zh) | 多分量地震数据的t-f-k域极化滤波方法 | |
CN107884829A (zh) | 一种联合压制浅海obc地震资料多次波的方法 | |
Hyslop et al. | Imaging lateral heterogeneity using reflected surface waves | |
CN106950600A (zh) | 一种近地表散射面波的去除方法 | |
Tao et al. | Processing of array sonic logging data with multi-scale STC technique | |
CN108919345B (zh) | 一种海底电缆陆检噪声的衰减方法 | |
CN111694055B (zh) | 一种强反射屏蔽下avo属性提取和分析的方法和装置 | |
Wang et al. | Ground roll wave suppression based on wavelet frequency division and radial trace transform | |
Du et al. | Research and application of Rayleigh wave extraction method based on microtremors signal analysis | |
Chopra et al. | Simultaneous acquisition of 3D surface seismic and 3D VSP data-processing and integration | |
Drijkoningen et al. | Nongeometrically converted shear waves in marine streamer data | |
Park et al. | Phase identification analysis using directionality and rectilinearity functions in three component seismograms | |
Li et al. | An efficient deep learning method for VSP wavefield separation: A DAS-VSP case | |
Viony et al. | Full Waveform Modelling for Subsurface Characterization with Converted-Wave Seismic Reflection: Residual PP Removal on PS Component Using FK Filter |
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 |