CN103901466A - 一种三维地震数据插值方法 - Google Patents
一种三维地震数据插值方法 Download PDFInfo
- Publication number
- CN103901466A CN103901466A CN201210585013.0A CN201210585013A CN103901466A CN 103901466 A CN103901466 A CN 103901466A CN 201210585013 A CN201210585013 A CN 201210585013A CN 103901466 A CN103901466 A CN 103901466A
- Authority
- CN
- China
- Prior art keywords
- interpolation
- frequency
- data
- window
- even number
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明是节约采集成本的三维地震数据插值方法,沿时间方向变换得到频率-空间域数据,将奇数、偶数道地震数据分为两组,变换至频率-波数域得到奇数道、偶数道波数域数据,将奇数道、偶数道波数域数据相除得到预测算子,与偶数频率数据相乘,并将该乘积通过傅立叶反变换得到时间-空间域插值结果,对当前计算窗以及前一计算窗进行均值处理,线性插值或者三次函数插值处理,完成插值。本发明相对于二维频率-波数域插值,插值结果具有更好的一致性,插值结果同相轴清晰、无明显抖动,背景自然;插值结果在信号与背景上均优于二维频率-波数域插值。
Description
技术领域
本发明涉及地球物理勘探技术,属于地震数据处理中用于加密地震数据、提高空间采样率和信噪比、节约采集成本的一种三维地震数据插值方法。
背景技术
在地震数据的采集过程中,由于施工环境以及工作量等客观条件的限制,采集得到的地震记录常常存在某些方向上道间距较大的问题,使得空间采样率严重不足,对偏移成像处理等会产生严重影响。解决这个问题最直接的方法自然是在野外采集时通过增加测线数量来减小空间采样间隔;但随之而来的是采集成本的大幅度提高,而且有些工区由于实际地表条件的限制,很难提高野外作业的空间采样率。
地震数据道内插值是在不增加采集成本的前提下,减小道间距的理想方法。频率-空间域地震道插值技术利用平面波的可预测性,在频率-空间域利用最小二乘法通过求解埃尔米特线性方程组来求得预测算子进行插值。频率-波数域地震道插值技术将频率-空间域插值技术发展,不需要在频率-空间域求解预测算子,插值过程均在频率-波数域完成,从而可以大幅度提高插值效率。
现有频率-波数域插值技术为二维插值,即每次插值时仅处理单个剖面。对于三维数据体应用二维插值时,每次仅利用了单个剖面上的信息,插值过程仅考虑了单一方向上的数据变化,而忽略了地震数据的整体性,容易造成插值结果不稳定。
发明目的
本发明的目的在于提供一种插值结果稳定、同相轴清晰、背景自然的三维频率-波数域地震数据插值方法。
发明内容
本发明通过如下步骤实现:
1)采集三维地震数据,做分时间、空间窗处理;
2)在当前计算窗内沿时间方向做傅立叶变换,得到频率-空间域数据,并将奇数、偶数道地震数据分为两组;
3)将步骤2)中奇数道、偶数道地震数据变换至频率-波数域,得到奇数道、偶数道波数域数据;将频率-空间域数据变换至频率-波数域得到偶数频率数据;
4)在频率-波数域中将步骤3)中得到的奇数道、偶数道波数域数据相除得到预测算子;
5)将步骤4)中得到的预测算子与步骤3)中得到的偶数频率数据相乘,并将该乘积通过傅立叶反变换得到时间-空间域插值结果;
6)对当前计算窗以及前一计算窗进行均值处理,线性插值或者三次函数插值处理;
7)当前计算窗如果是最末计算窗则完成插值,输出结果;如果不是最末计算窗则跳至下一计算窗并返回步骤2)循环,直至最末计算窗。
本发明相对于二维频率-波数域插值,插值结果具有更好的一致性,插值结果同相轴清晰、无明显抖动,背景自然;插值结果在信号与背景上均优于二维频率-波数域插值。
附图说明
图1:本发明处理流程图。
图2:某地区地震数据剖面(300CMP)
图3:本发明对该数据插值后所得剖面(300CMP)。
图4:二维频率-波数域插值方法对该数据插值后所得剖面(800LINE)。
图5:本发明对该数据插值后所得剖面(800LINE)。
具体实施方案
本发明按如下步骤实施,根据处理流程图1:
1)步骤S1,采集三维地震数据,记为smn(t),m=0,1,2,…,M,n=1,2,…,N,N表示X方向线数,M表示Y方向道数,t表示时间。
4)步骤S4,将奇数道、偶数道地震数据和变换至频率-波数域,得到奇数道、偶数道波数域数据,分别记为和将频率-空间域数据变换至频率-波数域得到频率-波数域数据,记为再沿频率方向间隔取值得到偶数频率数据,记为kx,ky分别表示X方向和Y方向的波数。
线性插值,即,
或者三次函数插值,即,
这里α为比例系数。
8)步骤S8,当前计算窗第i,j窗如果是最末计算窗则完成插值,至步骤S9输出结果;如果不是最末计算窗则跳至下一计算窗第i+1,j窗并返回步骤S3循环,直至最末计算窗。
图2、图3、图4和图5为本发明具体实施例。通过图2、图3比较可以看出,本发明插值结果与原始数据具有高度一致性,同相轴清晰,背景自然。
通过图4、图5对比可以看出,相对于二维频率-波数域插值,本发明插值结果具有更好的一致性,插值结果同相轴清晰、无明显抖动,背景自然;本发明插值结果在信号与背景上均优于二维频率-波数域插值。
Claims (1)
1.一种三维地震数据插值方法,通过如下步骤实现:
1)采集三维地震数据,做分时间、空间窗处理;
2)在当前计算窗内沿时间方向做傅立叶变换,得到频率-空间域数据,并将奇数、偶数道地震数据分为两组;
3)将步骤2)中奇数道、偶数道地震数据变换至频率-波数域,得到奇数道、偶数道波数域数据;将频率-空间域数据变换至频率-波数域得到偶数频率数据;
4)在频率-波数域中将步骤3)中得到的奇数道、偶数道波数域数据相除得到预测算子;
5)将步骤4)中得到的预测算子与步骤3)中得到的偶数频率数据相乘,并将该乘积通过傅立叶反变换得到时间-空间域插值结果;
6)对当前计算窗以及前一计算窗进行均值处理,线性插值或者三次函数插值处理;
7)当前计算窗如果是最末计算窗则完成插值,输出结果;如果不是最末计算窗则跳至下一计算窗并返回步骤2)循环,直至最末计算窗。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210585013.0A CN103901466B (zh) | 2012-12-28 | 2012-12-28 | 一种三维地震数据插值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210585013.0A CN103901466B (zh) | 2012-12-28 | 2012-12-28 | 一种三维地震数据插值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103901466A true CN103901466A (zh) | 2014-07-02 |
CN103901466B CN103901466B (zh) | 2016-08-10 |
Family
ID=50992914
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210585013.0A Active CN103901466B (zh) | 2012-12-28 | 2012-12-28 | 一种三维地震数据插值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103901466B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104990774A (zh) * | 2015-06-24 | 2015-10-21 | 中国石油天然气集团公司 | 一种地震数据插值方法 |
CN105182407A (zh) * | 2015-08-20 | 2015-12-23 | 中国石油天然气集团公司 | 一种三维地震数据道内插方法 |
CN105242304A (zh) * | 2014-07-04 | 2016-01-13 | 中国石油化工股份有限公司 | 基于频率空间域波场延拓的地震资料缺口补偿方法 |
CN109270576A (zh) * | 2018-11-06 | 2019-01-25 | 中国石油天然气股份有限公司 | 多通道三维地震物理模拟数据采集方法及装置 |
CN110595516A (zh) * | 2019-09-18 | 2019-12-20 | 华中科技大学 | 一种fpi腔长解调方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4922465A (en) * | 1989-05-30 | 1990-05-01 | Geco A/S | Interpolation of severely aliased events |
CN101201406A (zh) * | 2006-12-12 | 2008-06-18 | 中国石油天然气集团公司 | 一种高效地表一致性反褶积的方法 |
CN102221708A (zh) * | 2011-06-03 | 2011-10-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于分数阶傅里叶变换的随机噪声压制方法 |
CN102298155A (zh) * | 2011-05-23 | 2011-12-28 | 中国海洋石油总公司 | 一种基于高维小波变换的地震资料不连续性检测方法 |
CN103325092A (zh) * | 2013-03-13 | 2013-09-25 | 中国科学院电子学研究所 | 一种生成二维相位解缠质量图的方法及装置 |
-
2012
- 2012-12-28 CN CN201210585013.0A patent/CN103901466B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4922465A (en) * | 1989-05-30 | 1990-05-01 | Geco A/S | Interpolation of severely aliased events |
CN101201406A (zh) * | 2006-12-12 | 2008-06-18 | 中国石油天然气集团公司 | 一种高效地表一致性反褶积的方法 |
CN102298155A (zh) * | 2011-05-23 | 2011-12-28 | 中国海洋石油总公司 | 一种基于高维小波变换的地震资料不连续性检测方法 |
CN102221708A (zh) * | 2011-06-03 | 2011-10-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于分数阶傅里叶变换的随机噪声压制方法 |
CN103325092A (zh) * | 2013-03-13 | 2013-09-25 | 中国科学院电子学研究所 | 一种生成二维相位解缠质量图的方法及装置 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105242304A (zh) * | 2014-07-04 | 2016-01-13 | 中国石油化工股份有限公司 | 基于频率空间域波场延拓的地震资料缺口补偿方法 |
CN104990774A (zh) * | 2015-06-24 | 2015-10-21 | 中国石油天然气集团公司 | 一种地震数据插值方法 |
CN104990774B (zh) * | 2015-06-24 | 2018-01-05 | 中国石油天然气集团公司 | 一种地震数据插值方法 |
CN105182407A (zh) * | 2015-08-20 | 2015-12-23 | 中国石油天然气集团公司 | 一种三维地震数据道内插方法 |
CN105182407B (zh) * | 2015-08-20 | 2017-07-07 | 中国石油天然气集团公司 | 一种三维地震数据道内插方法 |
CN109270576A (zh) * | 2018-11-06 | 2019-01-25 | 中国石油天然气股份有限公司 | 多通道三维地震物理模拟数据采集方法及装置 |
CN110595516A (zh) * | 2019-09-18 | 2019-12-20 | 华中科技大学 | 一种fpi腔长解调方法及*** |
CN110595516B (zh) * | 2019-09-18 | 2020-12-18 | 华中科技大学 | 一种fpi腔长解调方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN103901466B (zh) | 2016-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2015205510B2 (en) | Determining a component of a wave field | |
Hennenfent et al. | Nonequispaced curvelet transform for seismic data reconstruction: A sparsity-promoting approach | |
US8103453B2 (en) | Method of seismic data interpolation by projection on convex sets | |
CN102854533B (zh) | 一种基于波场分离原理提高地震资料信噪比的去噪方法 | |
CN103901466A (zh) | 一种三维地震数据插值方法 | |
CN107678062B (zh) | 双曲Radon域综合预测反褶积和反馈循环方法压制多次波模型构建方法 | |
CN103064115B (zh) | 一种射线参数域纵波与转换波匹配方法 | |
CN102269820A (zh) | 一种基于gpu小存储量交错网格三维地震叠前逆时偏移成像方法 | |
CN109001813B (zh) | 一种压制多次波的方法、装置及*** | |
CN102692644B (zh) | 生成深度域成像道集的方法 | |
Dutta | Sparse least-squares reverse time migration using seislets | |
CN102053276A (zh) | 一种地震数字信号的复数道集二维滤波方法 | |
CN109738950B (zh) | 基于三维稀疏聚焦域反演的噪声型数据一次波反演方法 | |
CN103645497A (zh) | 一种基于出射角的多分量波场分离方法 | |
CN103576198A (zh) | 一种快速二维海上地震资料自由表面多次波预测方法 | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
Wang et al. | High-resolution wave-equation AVA imaging: Algorithm and tests with a data set from the Western Canadian Sedimentary Basin | |
CN102830431B (zh) | 真地表射线追踪自适应插值方法 | |
Zhai et al. | A new fractal interpolation algorithm and its applications to self-affine signal reconstruction | |
Song et al. | An efficient wavefield inversion for transversely isotropic media with a vertical axis of symmetry | |
CN105319593A (zh) | 基于曲波变换和奇异值分解的联合去噪方法 | |
CN106950600A (zh) | 一种近地表散射面波的去除方法 | |
Herrmann et al. | A modified, sparsity-promoting, Gauss-Newton algorithm for seismic waveform inversion | |
CN103558636A (zh) | 一种从叠后地震数据采集脚印衰减的方法 | |
US10598807B2 (en) | Correction of sea surface state |
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 |