CN103576193B - 一种消除逆时偏移低频假像的方法 - Google Patents
一种消除逆时偏移低频假像的方法 Download PDFInfo
- Publication number
- CN103576193B CN103576193B CN201210273299.9A CN201210273299A CN103576193B CN 103576193 B CN103576193 B CN 103576193B CN 201210273299 A CN201210273299 A CN 201210273299A CN 103576193 B CN103576193 B CN 103576193B
- Authority
- CN
- China
- Prior art keywords
- big gun
- wave field
- reverse
- time migration
- single big
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明在于提供一种在进行滤波的同时,尽可能地保持地震波动力学特征的消除逆时偏移低频假像的方法。采用深度域层速度模型和正演炮集数据,首先计算补偿后的震源子波,在深度域进行单炮波动方程正演,存储所有时刻单炮源波场;然后用单炮记录作为初始边界,沿时间轴逆向进行波动方程正演得到接收波场,在波场逆推时读取相应时刻单炮源波场,计算补偿后的单炮逆时偏移成像,得到低频假像消除后的单炮逆时偏移成像。在构造成像和波形相对保持、沿层地震属性空变稳定性方面,本发明提出的逆时偏移低频假像的消除方法要好于Laplace算子直接滤波方法。
Description
技术领域
本发明涉及勘探和油藏地球物理技术,具体是一种消除地震逆时偏移成像中低频假像的方法。
背景技术
在勘探和油藏地球物理中,地震偏移成像是其中最关键的一步,成像的好坏直接影响到后续的地震和地质解释,从而影响油气储层的识别。地震偏移成像可以分为几大类:(1)基于射线理论的Kirchhoff积分法偏移;(2)基于单程波动理论的F-K偏移和有限差分偏移;(3)基于双程波动理论的逆时偏移。(1)和(2)属于常规地震偏移,对横向速度剧烈变化、高陡倾角等复杂地质构造,成像能力有限。而逆时偏移是基于双程波动理论,理论上它能解决常规地震偏移所面临的问题,当然它也会带来很多新的问题,其中首要问题是成像时引入了低频假像干扰。这是由于互相关成像时地质界面上存在反射波而导致的,主要分布在波的传播路径上,并且在波阻抗越大的地方,这种假像越强,严重干扰了成像效果(Claerbout,1971;Biondi和Shan,2002)。
为了有效消除和压制这种低频假像,Youn和Zhou(2001)等给出了Laplace算子直接滤波来消除逆时偏移假像,这种方法适应复杂介质,但会破坏有效信号的特征。Yoon等(2004)引入Poynting矢量方法来消除逆时偏移假像,这种方法需要判断波的传播方向,对于复杂地质构造而言,往往比较困难。Fletcher等(2006)提出用压制层间反射的方法来消除逆时偏移假像,这种方法破坏了波场传播特性。Bulcao等(2007)采用波场分离方法来消除逆时偏移假像,这种方法忽略了回转波对复杂构造的成像贡献。Guitton等(2007)提出最小二乘滤波方法来消除逆时偏移假像,由于反演的多解性,影响了滤波的效果。逆时偏移低频假像的消除好坏,直接关系到后续的储层解释和描述等。在以上这些方法中,Laplace算子滤波是一种能够适应复杂介质、简单可行的滤波方法,在实际当中比较常用。但直接滤波会破坏地震波的动力学特征,影响成像效果。
发明内容
本发明目的在于提供一种在进行滤波的同时,尽可能地保持地震波动力学特征的消除逆时偏移低频假像的方法。
本发明通过以下步骤实现:
1)采集地震勘探数据,处理得到深度域层速度模型和正演炮集数据;
2)利用公式计算补偿后的震源子波;(创新技术)
其中S°(t)为补偿后的震源子波,S(t)为震源子波,*为褶积,u(t)为阶跃函数,t为时间;
3)利用补偿后的震源子波作为震源,利用下式在深度域进行单炮波动方程正演:
并存储所有时刻补偿后的单炮源波场;(创新技术)
其中:v为地震波速度,Ps°(X,t)为补偿后的单炮源波场,Δ为Laplace算子,δ(X-Xs)为单位脉冲函数,X=(x,z)为波场所在二维空间位置坐标,Xs=(xs,zs)为震源所在二维空间位置坐标,x和xs分别为波场和震源横坐标,z和zs分别为波场和震源纵坐标;
4)利用单炮地震记录作为初始边界条件,用下式从最大记录时刻开始,沿时间轴逆向进行波动方程正演,得到接收波场:
其中:Pg(X,t)为接收波场,G(Xg,t)为单炮地震记录,Xg=(xg,zg)为检波器所在二维空间位置坐标,xg为检波器横坐标,zg为检波器纵坐标;
5)在接收波场逆推时,读取相应时刻补偿后的单炮源波场,利用以下公式计算补偿后的单炮逆时偏移成像:
I°(X)=∫Ps°(X,t)Pg(X,t)dt;(创新技术)
6)利用公式ISF(X)=S·F{I°(X)}进行逆时偏移低频假像消除,得到低频假像消除后的单炮逆时偏移成像;(创新技术)
其中:ISF(X)为消除低频假像的逆时偏移成像,S为满足二阶导数连续的平滑算子以保证成像波形的连续性(所述的平滑算子采用三次样条函数插值平滑算子),F=v2Δ为滤波算子;
7)将所有单炮逆时偏移成像进行叠加,得到消除低频假像的叠加逆时偏移成像。
本发明试验结果表明:在构造成像和波形相对保持、沿层地震属性空变稳定性方面,本发明提出的逆时偏移低频假像的消除方法要好于Laplace算子直接滤波方法。
附图说明
图1是理论模型与单炮地震记录。(a)为逆掩断层理论模型;(b)为其中一炮地震记录。
图2是滤波前后的逆时偏移成像。(a)未滤波;(b)Laplace算子直接滤波方法;(c)本发明提出的滤波方法;(d)为(b)滤掉的噪音部分;(e)为(c)滤掉的噪音部分;(f)为(a)的局部放大;(g)为(b)的局部放大;(h)为(c)的局部放大。
图3是滤波前后的逆时偏移成像波形(散射点P2)。(a)未滤波;(b)Laplace算子直接滤波方法;(c)本发明提出的滤波方法;(d)为(b)滤掉的噪音波形残差;(e)为(c)滤掉的噪音波形残差。
图4是沿层L2地震属性分析。(a)沿层L2层拉平剖面(Laplace算子直接滤波方法);(b)沿层L2层拉平剖面(本发明提出的滤波方法);(c)沿层L2瞬时振幅;(d)沿层L2瞬时频率;(e)沿层L2瞬时相位;(上)Laplace算子直接滤波方法;(下)本发明提出的滤波方法;(f)瞬时频率色标;(g)瞬时相位色标。注:由于瞬时振幅尺度范围不同,此处只关心相对空间变化特性,故省略瞬时振幅色标。
具体实施方式
以下结合附图和实施例详细说明本发明。
本发明通过以下步骤实现:
1)采集地震勘探数据,处理得到深度域层速度模型(图1(a),主要由一组单斜构造、一个逆掩断层和三个散射点组成)和正演炮集数据(图1(b),为其中一炮地震记录);
2)利用公式计算补偿后的震源子波;
其中S°(t)为补偿后的震源子波,S(t)为震源子波,*为褶积,u(t)为阶跃函数,t为时间;
3)利用补偿后的震源子波作为震源,利用下式在深度域进行单炮波动方程正演:
并存储所有时刻补偿后的单炮源波场;
其中:v为地震波速度,Ps°(X,t)为补偿后的单炮源波场,Δ为Laplace算子,δ(X-Xs)为单位脉冲函数,X=(x,z)为波场所在二维空间位置坐标,Xs=(xs,zs)为震源所在二维空间位置坐标,x和xs分别为波场和震源横坐标,z和zs分别为波场和震源纵坐标;
4)利用单炮地震记录作为初始边界条件,用下式从最大记录时刻开始,沿时间轴逆向进行波动方程正演,得到接收波场:
其中:Pg(X,t)为接收波场,G(Xg,t)为单炮地震记录,Xg=(xg,zg)为检波器所在二维空间位置坐标,xg为检波器横坐标,zg为检波器纵坐标;
5)在接收波场逆推时,读取相应时刻补偿后的单炮源波场,利用以下公式计算补偿后的单炮逆时偏移成像:
I°(X)=∫Ps°(X,t)Pg(X,t)dt;
6)利用公式ISF(X)=S·F{I°(X)}进行逆时偏移低频假像消除,得到低频假像消除后的单炮逆时偏移成像;
其中:ISF(X)为消除低频假像的逆时偏移成像,S为满足二阶导数连续的平滑算子以保证成像波形的连续性(所述的平滑算子采用三次样条函数插值平滑算子),F=v2Δ为滤波算子;
7)将所有单炮逆时偏移成像进行叠加,得到消除低频假像的叠加逆时偏移成像(图2(c))。
以上步骤2)到7)是放在一个程序包里实现的,直接运行即可得到最终的结果。从构造成像(图2(b)、(g)与(c)、(h))和滤掉的噪音(图2(d)、(e))可以看出,在逆掩断层及其下盘有效信号保持方面,本发明提出的滤波方法要好于Laplace算子直接滤波方法。从成像波形(图3(a)、(b)、(c))和波形残差(图3(d)、(e))可以看出,在波形相对保持上,本发明提出的滤波方法要好于Laplace算子直接滤波方法。从沿层L2地震瞬时振幅(图4(c))、瞬时频率(图4(d))和瞬时相位(图4(e))沿层空变稳定性上,本发明提出的滤波方法要好于Laplace直接滤波方法。
综上所述,在构造成像、波形相对保持以及沿层地震属性(瞬时振幅、频率和相位)空变稳定性上,本发明提出的滤波方法要好于Laplace算子直接滤波方法。
Claims (1)
1.一种消除逆时偏移低频假像的方法,特点是采用以下步骤实现:
1)采集地震勘探数据,处理得到深度域层速度模型和正演炮集数据;
2)利用公式计算补偿后的震源子波;
其中S°(t)为补偿后的震源子波,S(t)为震源子波,*为褶积,u(t)为阶跃函数,t为时间;
3)利用补偿后的震源子波作为震源,利用下式在深度域进行单炮波动方程正演:
并存储所有时刻补偿后的单炮源波场;
其中:v为地震波速度,为补偿后的单炮源波场,Δ为Laplace算子,δ(X-Xs)为单位脉冲函数,X=(x,z)为波场所在二维空间位置坐标,Xs=(xs,zs)为震源所在二维空间位置坐标,x和xs分别为波场和震源横坐标,z和zs分别为波场和震源纵坐标;
4)利用单炮地震记录作为初始边界条件,用下式从最大记录时刻开始,沿时间轴逆向进行波动方程正演,得到接收波场:
其中:Pg(X,t)为接收波场,G(Xg,t)为单炮地震记录,Xg=(xg,zg)为检波器所在二维空间位置坐标,xg为检波器横坐标,zg为检波器纵坐标;
5)在接收波场逆推时,读取相应时刻补偿后的单炮源波场,利用以下公式计算补偿后的单炮逆时偏移成像:
其中:为补偿后的单炮源波场,Pg(X,t)为接收波场;
6)利用公式进行逆时偏移低频假像消除,得到低频假像消除后的单炮逆时偏移成像;
其中:S为满足二阶导数连续的平滑算子以保证成像波形的连续性,所述的平滑算子采用三次样条函数插值平滑算子,F=v2Δ为滤波算子;
7)将所有单炮逆时偏移成像进行叠加,得到消除低频假像的叠加逆时偏移成像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210273299.9A CN103576193B (zh) | 2012-08-02 | 2012-08-02 | 一种消除逆时偏移低频假像的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210273299.9A CN103576193B (zh) | 2012-08-02 | 2012-08-02 | 一种消除逆时偏移低频假像的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103576193A CN103576193A (zh) | 2014-02-12 |
CN103576193B true CN103576193B (zh) | 2016-08-10 |
Family
ID=50048359
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210273299.9A Active CN103576193B (zh) | 2012-08-02 | 2012-08-02 | 一种消除逆时偏移低频假像的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103576193B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105510973B (zh) * | 2014-09-23 | 2018-06-15 | 中国石油化工股份有限公司 | 一种用于压制逆时偏移成像噪音的拉普拉斯滤波方法 |
US10962669B2 (en) * | 2018-12-28 | 2021-03-30 | China Petroleum & Chemical Corporation | Method and device for enhanced seismic imaging based on one-way wave equation |
CN109975873B (zh) * | 2019-04-23 | 2020-06-30 | 中国石油大学(华东) | 一种逆时偏移成像去除低频噪音的方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102116869A (zh) * | 2011-02-12 | 2011-07-06 | 中国石油大学(华东) | 高精度叠前域最小二乘偏移地震成像技术 |
CN102269820A (zh) * | 2010-06-01 | 2011-12-07 | 潜能恒信能源技术股份有限公司 | 一种基于gpu小存储量交错网格三维地震叠前逆时偏移成像方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101092668B1 (ko) * | 2009-06-17 | 2011-12-13 | 서울대학교산학협력단 | 파형 역산을 이용한 지하 구조의 영상화 장치와 방법 |
-
2012
- 2012-08-02 CN CN201210273299.9A patent/CN103576193B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102269820A (zh) * | 2010-06-01 | 2011-12-07 | 潜能恒信能源技术股份有限公司 | 一种基于gpu小存储量交错网格三维地震叠前逆时偏移成像方法 |
CN102116869A (zh) * | 2011-02-12 | 2011-07-06 | 中国石油大学(华东) | 高精度叠前域最小二乘偏移地震成像技术 |
Non-Patent Citations (3)
Title |
---|
Imaging conditions for prestack reverse-time migration;Sandip Chattopadhyay 等;《Geophysics》;20080502;第73卷(第3期);第S81-S89页 * |
逆时偏移成像技术研究进展;丁亮 等;《地球物理学进展》;20110630;第26卷(第3期);第1085-1100页 * |
逆时偏移成像条件研究;王娟 等;《地球物理学进展》;20120630;第27卷(第3期);第1173-1182页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103576193A (zh) | 2014-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Borisov et al. | Application of 2D full-waveform inversion on exploration land data | |
CN106526677B (zh) | 一种海上自适应压制鬼波的宽频逆时偏移成像方法 | |
CN102156296A (zh) | 地震多分量联合弹性逆时偏移成像方法 | |
CN103149585B (zh) | 一种弹性偏移地震波场构建方法及装置 | |
CA2910781A1 (en) | Deghosting with adaptive operators | |
US9234976B2 (en) | Device and method for removal of multiples from seismic data vintages | |
CN103576193B (zh) | 一种消除逆时偏移低频假像的方法 | |
CN103018777B (zh) | 一种能够消除逆时偏移噪声的成像方法及装置 | |
CN105652322A (zh) | 多分量地震数据的t-f-k域极化滤波方法 | |
CN104570124A (zh) | 一种适合井间地震大角度反射条件的延拓成像方法 | |
CN107884828A (zh) | 一种在空间‑频率域基于格林理论的鬼波压制方法 | |
CN103412325B (zh) | 一种地震纵、横波波场分离与去噪的方法 | |
CN104570114B (zh) | 一种基于波场分解的逆时偏移噪音压制方法 | |
CN105277983B (zh) | 压制可控震源滑动扫描谐波干扰的方法 | |
Shragge et al. | Wave-equation migration from topography | |
CN103576201A (zh) | 基于均炮源波场照明的逆时偏移空间振幅补偿方法 | |
CN104459792B (zh) | 一种构造约束下的边缘保护滤波方法 | |
Chiu | 3D attenuation of aliased ground roll on randomly undersampled data | |
AU2011248910B2 (en) | System and method for compensating time and offset varying near-surface effects in seismic data background | |
Wang et al. | Ground roll wave suppression based on wavelet frequency division and radial trace transform | |
Poletto et al. | Multidimensional deconvolution of seismic-interferometry Arctic data | |
D'Agosto et al. | Modeling and removal of ground roll from horizontal component of C-waves | |
Li* et al. | Simultaneous de-noising and regularization method for 3D field data based on the MSSA | |
Naidu et al. | Surface related multiple elimination: A case study from East Coast India | |
Abma et al. | Sparse inversion of simultaneous sources |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |