CN103645497A - 一种基于出射角的多分量波场分离方法 - Google Patents
一种基于出射角的多分量波场分离方法 Download PDFInfo
- Publication number
- CN103645497A CN103645497A CN201310553067.3A CN201310553067A CN103645497A CN 103645497 A CN103645497 A CN 103645497A CN 201310553067 A CN201310553067 A CN 201310553067A CN 103645497 A CN103645497 A CN 103645497A
- Authority
- CN
- China
- Prior art keywords
- component
- mean
- wave
- root
- velocity
- 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
- 238000000926 separation method Methods 0.000 title claims abstract description 32
- 238000005070 sampling Methods 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 17
- 230000010287 polarization Effects 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000001615 p wave Methods 0.000 description 2
- 229910052704 radon Inorganic materials 0.000 description 2
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种基于出射角的多分量波场分离方法,同时考虑了地震波的运动学和动力学特征,利用地层速度计算地震波的出射角,通过矢量旋转实现多分量的波场分离。本发明能有效分离折射纵波和横波,有效分解一次反射波的多分量波场,时间空间域实现波场矢量保真性保持好,与其它域实现方法相比无波形畸变,方法简单、易于实现。
Description
技术领域
本发明涉及多波多分量地震资料处理技术,具体地,涉及一种基于出射角的多分量波场分离方法。
背景技术
多波多分量地震勘探就是采用多分量激发和接收,综合利用纵波、横波和转换波等多种地震波信息,实现改善构造成像、进行岩性分析、检测储层裂缝以及直接预测油气等目的的地震技术。目前,在海上和陆上实施最多的多波多分量地震勘探方法是转换波地震勘探。转换波地震勘探是以纵波震源激发,地面三分量或海底四分量检波器接收的多分量勘探方法。
多分量地震资料记录的是一个矢量波场。在野外采集时,激发点与接收点总有一定的距离,这就使得在同一个接收点位置,纵波与转换波出射时总存在一定的出射角,这样水平分量上除接收到转换波外,还接收到部分纵波能量;同理,垂直分量接收到的也不全是纵波能量,也有少量转换波能量。在浅层或大炮检距时,接收的资料受出射角影响更大。这样,每个分量上都同时含有纵波和横波两种成份,特别是对于接收到的转换波,在达到一定的炮检距后才能接收到好的转换波资料,如果在资料处理中不能对波场进行分离,将会直接影响后续成像的质量。
目前用于波场分离的方法已有不少,主要有拉冬变换,F-K滤波法,偏振滤波法以及波动方程分解法,从文献资料看,各种分离方法用于合成记录都能得到令人满意的分离结果,但是实际应用的成功的例子并不多。拉冬变换,F-K滤波法是利用多波资料的运动学特征(即视速度特征)进行分离,偏振滤波法主要是利用多波资料的动力学特征(即偏振特性)进行波场分离,这些仅考虑了波的单一特征的分离方法,其应用范围均受到一定限制,也影响了分离效果,理论模型计算时,波动方程法的效果较好,但由于低速层的影响,实际资料计算时不能满足波动方程的条件,分离效果并不是很理想。需要采用新的方法解决多分量地震资料处理中波场分离的问题。
发明内容
本发明目的是提供一种基于出射角的多分量波场分离方法,同时考虑了地震波的运动学和动力学特征,解决因接收点出射角的影响,造成纵波和转换波波场混叠问题。
本发明通过以下技术方案实现:
一种基于出射角的多分量波场分离方法,具体步骤包括:
1)采集多分量地震数据,并进行预处理,
步骤1)所述的预处理包括加载观测系统、多分量地震道匹配、获得纵波近地表速度和地下均方根速度;
2)输入一道水平分量的径向分量R和垂直分量Z,
步骤2)所述的水平分量的径向分量R是指两个野外采集水平分量X和Y经坐标旋转,也就是检波器重定向得到的径向分量;
3)读取该地震道对应的均方根速度,进行均方根速度内插,并利用小时窗进行速度平滑;
步骤3)所述的进行均方根速度内插是指对输入的均方根速度进行插值,使得该道每个样点都有对应的均方根速度值;
4)由均方根速度和样点时间计算反射点深度,再由炮检距得该样点近地表近似入射角;
5)由近地表速度和第一个样点的均方根速度以及步骤4)计算的近似入射角,得该样点出射角的正弦和余弦,进行矢量旋转分离;
6)第一个样点至该道最后一个样点重复步骤4)和步骤5)即完成一道矢量旋转分离;
7)输出矢量旋转分离后的水平分量R和垂直分量Z。
本发明同时考虑了地震波的视速度、偏振特性及地层速度的变化,能有效分离纵波和横波,有效分解一次反射波的多分量波场,时间空间域实现波场矢量保真性保持好,与其它域实现方法相比无波形畸变,方法简单、易于实现。
附图说明
图1转换波传播路径示意图;
图2纵波和转换波波场分离原理示意图;
图3纵波和转换波波场分离前单炮剖面;
图4纵波和转换波波场分离后单炮剖面。
具体实施方式
本发明目的是提供一种基于出射角的多分量波场分离方法,利用表层速度与深层速度比值求取P波的出射角,通过矢量旋转实现波场分离,解决因接收点出射角的影响,造成纵波和转换波波场混叠问题。
本发明同时考虑了地震波的视速度、偏振特性及传播速度变化,根据描述这三者关系的波场分离矩阵,在时间空间域变换过程中进行分量分解,可以实现有效的波场分离。用近地表的P波速度确定旋转角度,旋转后的径向分量(R)上不再有纵波(PP)能量,垂直分量包括全部的P波能量,但仍有一些转换波(PS)能量;同样,用近地表的S波速度也可以改进PP波道集的质量。
本发明的实现原理如下:
本发明利用表层速度与深层速度比值求取P波与PS波的出射角,通过矢量旋转实现波场分离。
在各向均匀水平层状介质的条件下转换波的传播路径如图1所示,A为激发点,B为接收点,V0为近地表地层速度,V1为深层速度,假设h’/h≈0,则C’B’≈CB=x
再由斯奈尔定律可得
对于同一激发点和接收点,纵波和转换波的出射角分别为:
而Vps0=Vp0/γ0,Vps1=Vp1/γ1
所以
其中,βp和βps分别为近地表纵波和转换波的出射角,Vp0和Vps0分别为地表纵波和横波的速度,Vp1和Vps1分别为深层纵波和横波的均方根速度,xp和xps分别为纵波共中心点到共接收点的距离和转换波的共转换点到共接收点的距离,γ0和γ1分别为近地表和深层的纵横波速度比。
P波与PS波在接收点位置出射时,由于与地面存在一定的夹角,如图2所示,R分量上除接收到PS转换波外,还接收到部分P波能量,同理,Z分量接收到的也不全是P波能量,也有少量转换波能量,P波和PS波在Z分量和R分量上的能量投影关系为:
由上式可以得到:
由于 ,所以
这样,通过矢量旋转实现了波场分离。
本发明实现简便、适应大数据量处理。
本发明在转换波数据上采用以下步骤应用:
1)采集多分量地震数据,并进行预处理,完成观测系统加载,多分量地震道匹配等基础处理,获得纵波近地表速度和地下均方根速度;
2)输入一道水平分量的径向分量(R)和垂直分量(Z);
3)读取该地震道对应的均方根速度,进行均方根速度线性插值,使得该道每个样点都有对应的均方根速度值,并利用小时窗进行速度平滑,小时窗长度为50毫秒;
4)由均方根速度和样点时间计算反射点深度,再由炮检距,可得该样点近地表近似入射角;
5)近地表速度和高速顶层速度(第一个样点的均方根速度)以及步骤4)计算的近似入射角,可得该样点出射角的正弦和余弦,进行矢量旋转分离;
6)第一个样点至该道最后一个样点重复步骤4)和步骤5)即完成一道矢量旋转分离;
7)输出矢量旋转分离后的水平分量R和垂直分量Z。
图3是纵波和转换波波场分离前单炮剖面;图2纵波和转换波波场分离后单炮剖面。可看出应用本发明后,转换波的明显得到改善,特别初至部分。
Claims (1)
1.一种基于出射角的多分量波场分离方法,具体步骤包括:
1)采集多分量地震数据,并进行预处理,
步骤1)所述的预处理包括加载观测系统、多分量地震道匹配、获得纵波近地表速度和地下均方根速度;
2)输入一道水平分量的径向分量R和垂直分量Z,
步骤2)所述的水平分量的径向分量R是指两个野外采集水平分量X和Y经坐标旋转,也就是检波器重定向得到的径向分量;
3)读取该地震道对应的均方根速度,进行均方根速度内插,并利用小时窗进行速度平滑;
步骤3)所述的进行均方根速度内插是指对输入的均方根速度进行插值,使得该道每个样点都有对应的均方根速度值;
4)由均方根速度和样点时间计算反射点深度,再由炮检距得该样点近地表近似入射角;
5)由近地表速度和第一个样点的均方根速度以及步骤4)计算的近似入射角,得该样点出射角的正弦和余弦,进行矢量旋转分离;
6)第一个样点至该道最后一个样点重复步骤4)和步骤5)即完成一道矢量旋转分离;
7)输出矢量旋转分离后的水平分量R和垂直分量Z。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310553067.3A CN103645497B (zh) | 2013-11-08 | 2013-11-08 | 一种基于出射角的多分量波场分离方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310553067.3A CN103645497B (zh) | 2013-11-08 | 2013-11-08 | 一种基于出射角的多分量波场分离方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103645497A true CN103645497A (zh) | 2014-03-19 |
CN103645497B CN103645497B (zh) | 2015-04-22 |
Family
ID=50250745
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310553067.3A Expired - Fee Related CN103645497B (zh) | 2013-11-08 | 2013-11-08 | 一种基于出射角的多分量波场分离方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103645497B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104101897A (zh) * | 2014-07-11 | 2014-10-15 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种利用矢量合成纵波进行地震勘探的方法 |
CN105467440A (zh) * | 2015-10-28 | 2016-04-06 | 中国石油天然气股份有限公司 | 一种全向矢量地震数据处理方法及装置 |
CN109283582A (zh) * | 2018-10-18 | 2019-01-29 | 中国石油化工股份有限公司 | 一种静校正的补偿方法 |
CN109358361A (zh) * | 2018-10-18 | 2019-02-19 | 河海大学 | 一种提取地震射线矢量的方法 |
CN110398776A (zh) * | 2019-08-21 | 2019-11-01 | 西南石油大学 | 一种基于右手法则的波场分离方法 |
CN110780346A (zh) * | 2019-11-20 | 2020-02-11 | 李志勇 | 一种隧道超前探测复杂地震波场的分离方法 |
US10718875B2 (en) | 2015-10-28 | 2020-07-21 | Petrochina Company Limited | Omnidirectional vector electrostatic levitation geophone |
US10739477B2 (en) | 2015-10-28 | 2020-08-11 | Petrochina Company Limited | Omnidirectional vector geophone |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4766574A (en) * | 1987-03-31 | 1988-08-23 | Amoco Corporation | Method for depth imaging multicomponent seismic data |
CN102012522A (zh) * | 2010-11-19 | 2011-04-13 | 中国石油天然气股份有限公司 | 地震采集数据的评价方法及装置 |
CN102053267A (zh) * | 2010-10-22 | 2011-05-11 | 中国石油化工股份有限公司 | 一种地震剖面资料处理中基于参数反演的vsp波场分离方法 |
CN102667529A (zh) * | 2009-10-27 | 2012-09-12 | 雪佛龙美国公司 | 使用射束层析成像进行地震成像和地层建模的方法和系统 |
-
2013
- 2013-11-08 CN CN201310553067.3A patent/CN103645497B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4766574A (en) * | 1987-03-31 | 1988-08-23 | Amoco Corporation | Method for depth imaging multicomponent seismic data |
CN102667529A (zh) * | 2009-10-27 | 2012-09-12 | 雪佛龙美国公司 | 使用射束层析成像进行地震成像和地层建模的方法和系统 |
CN102053267A (zh) * | 2010-10-22 | 2011-05-11 | 中国石油化工股份有限公司 | 一种地震剖面资料处理中基于参数反演的vsp波场分离方法 |
CN102012522A (zh) * | 2010-11-19 | 2011-04-13 | 中国石油天然气股份有限公司 | 地震采集数据的评价方法及装置 |
Non-Patent Citations (1)
Title |
---|
鲁传恒等: "τ-p偏振滤波地震波场分离", 《石油地球物理勘探》, no. 05, 30 October 1993 (1993-10-30) * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104101897A (zh) * | 2014-07-11 | 2014-10-15 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种利用矢量合成纵波进行地震勘探的方法 |
CN105467440A (zh) * | 2015-10-28 | 2016-04-06 | 中国石油天然气股份有限公司 | 一种全向矢量地震数据处理方法及装置 |
US10578757B2 (en) | 2015-10-28 | 2020-03-03 | Petrochina Company Limited | Omnidirectional vector seismic data processing method and apparatus, computer-readable storage medium and device |
US10718875B2 (en) | 2015-10-28 | 2020-07-21 | Petrochina Company Limited | Omnidirectional vector electrostatic levitation geophone |
US10739477B2 (en) | 2015-10-28 | 2020-08-11 | Petrochina Company Limited | Omnidirectional vector geophone |
CN109283582A (zh) * | 2018-10-18 | 2019-01-29 | 中国石油化工股份有限公司 | 一种静校正的补偿方法 |
CN109358361A (zh) * | 2018-10-18 | 2019-02-19 | 河海大学 | 一种提取地震射线矢量的方法 |
CN110398776A (zh) * | 2019-08-21 | 2019-11-01 | 西南石油大学 | 一种基于右手法则的波场分离方法 |
CN110398776B (zh) * | 2019-08-21 | 2021-03-30 | 西南石油大学 | 一种基于右手法则的波场分离方法 |
CN110780346A (zh) * | 2019-11-20 | 2020-02-11 | 李志勇 | 一种隧道超前探测复杂地震波场的分离方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103645497B (zh) | 2015-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103645497B (zh) | 一种基于出射角的多分量波场分离方法 | |
Cheng et al. | Frequency–wavenumber (FK)-based data selection in high-frequency passive surface wave survey | |
CN102033242B (zh) | 一种深层倾斜裂缝储层地震振幅预测方法 | |
WO2017024523A1 (zh) | 一种射线弹性参数的反演方法 | |
CN104730579B (zh) | 一种基于表层横波速度反演的纵横波联合静校正方法 | |
CN105068117B (zh) | 用于裂缝性介质的avo反演方法、装置和设备 | |
She et al. | Shallow crustal structure of the middle‐lower Yangtze River region in eastern China from surface‐wave tomography of a large volume airgun‐shot experiment | |
CN103064115B (zh) | 一种射线参数域纵波与转换波匹配方法 | |
CN102636811B (zh) | 一种海上二维地震资料中多次波的消除方法 | |
CN101598803B (zh) | 一种直接得到转换波叠加剖面的方法 | |
CN101915939B (zh) | 一种面波压制方法 | |
CN105372705B (zh) | 一种基于多波资料的地层切片方法 | |
CN103424776A (zh) | 一种碳酸盐岩油气储层裂缝地震检测方法 | |
CN104820242B (zh) | 一种面向叠前反演的道集振幅分频补偿方法 | |
CN104316966B (zh) | 一种流体识别方法及系统 | |
CN104280777A (zh) | 一种压制陆上地震资料多次波干扰的方法 | |
CN103592683A (zh) | 一种用于转换波地震资料的地表一致性振幅补偿方法及装置 | |
CN104316965A (zh) | 一种裂缝方位和强度的预测方法及系统 | |
CN103675911A (zh) | 一种基于纵波和转换横波联合反演截距和梯度的方法 | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
CN104237938A (zh) | 利用动态时间调整进行纵波和转换波数据匹配的方法 | |
CN102854526A (zh) | 一种多分量地震资料的处理方法 | |
Özmen et al. | Microtremor array measurements for shallow S‐wave profiles at strong‐motion stations in Hatay and Kahramanmaras Provinces, Southern Turkey | |
CN103576197A (zh) | 一种转换波角道集抽取方法 | |
Plescia et al. | Teleseismic P‐wave coda autocorrelation imaging of crustal and basin structure, Bighorn Mountains Region, Wyoming, USA |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150422 |