陈思远^1,,
曹思远^1,,,
孙耀光¹,
江雨濛¹,
高君²
1. 中国石油大学(北京), 北京 102249
2. 中国石化石油勘探开发研究院, 北京 100083

基金项目: 国家重点研发计划"面向E级计算的能源勘探高性能应用软件系统与示范项目"子课题"基于压缩感知的海量数据高效并行处理"（2017YFB0202900）资助


详细信息

作者简介: 陈思远, 男, 1995年生, 博士研究生, 主要从事信号处理、时频分析、稀疏反演等方向研究.E-mail: chensy_1362@outlook.com

通讯作者: 曹思远, 男, 1962年生, 教授, 博士生导师, 主要研究方向为地震信号处理、时频分析、储层预测等.E-mail: siyuan.caoff@gmail.com

中图分类号: P631

收稿日期:2020-06-05
修回日期:2021-07-29
上线日期:2021-10-10

Seismic data interpolation based on the frequency-structure dual weighted threshold

CHEN SiYuan^1,,
CAO SiYuan^1,,,
SUN YaoGuang¹,
JIANG YuMeng¹,
GAO Jun²
1. China University of Petroleum(Beijing), Beijing 102249, China
2. Sinopec Petroleum Exploration andProduction Research Institute, Beijing 100083, China


More Information

Corresponding author: CAO SiYuan,E-mail:siyuan.caoff@gmail.com

MSC: P631

--> Received Date: 05 June 2020
Revised Date: 29 July 2021
Available Online: 10 October 2021


摘要
摘要:地震数据的插值是地震数据处理的关键环节，插值效果的优劣关系到后续反演、去噪、反褶积、偏移等工作的进行.基于稀疏反演的常规频率-波数（FK）域插值方法往往使用广义阈值，而由于缺失地震数据仅沿波数方向呈稀疏态，在沿频率方向呈类子波频谱形态分布的非稀疏态，直接应用基于Lp范数的广义阈值会难以补全缺失数据的高低频信息，不利于后续依赖低频信息的反演及依赖高频的高分辨率处理工作的开展；基于该问题，本文设计沿频率方向类子波频谱形态变化的频率加权阈值，同时考虑完整地震数据在波数域的聚集分布，引入沿波数方向的结构加权阈值，提出基于频率-结构的双重加权阈值，在补全缺失地震道的同时，保证其高低频信息的有效恢复.模型和实际数据表明，相比于广义阈值，基于频率-结构的双重加权阈值可以减少插值误差，最大限度保证缺失地震道位置插值结果高低频信息的完整性.
关键词: 插值/
结构加权/
频率加权/
FK变换/
稀疏反演

Abstract:Seismic data interpolation is a crucial part of seismic data processing. The quality of the interpolation is related to the subsequent inversion, denoising, deconvolution, migration, and other work. Conventional FK interpolation methods based on sparse inversion often use generalized thresholds. However, the missing seismic data in these methods only shows a sparse state along the wavenumber direction and a non-sparse state with a wavelet-like spectrum distribution along the frequency direction. The direct application of the generalized threshold based on the Lp norm will make it difficult to complete the high- and low-frequency information of the missing data, which is not conducive to the subsequent inversion that relies on low-frequency information and the development of high-resolution processing that relies on high frequency. To solve this problem, we design a frequency-weighted threshold that changes along the frequency direction and introduces a structure-weighted threshold along the wavenumber direction considering the aggregation distribution of complete seismic data in the wavenumber domain. Combining the above two points, we propose a frequency-structure dual weighted threshold to ensure the effective recovery of high- and low-frequency information while completing missing traces. By comparing with the generalized threshold, the model and actual data show that the frequency-structure dual weighted threshold can reduce the interpolation error and ensure the integrity of the high- and low-frequency information to the largest extent.
Key words:Interpolation/
Structure weighted/
Frequency weighted/
FK transform/
Sparse inversion



PDF全文下载地址:

http://www.geophy.cn/data/article/export-pdf?id=be860884-6485-4bb2-bf1f-4bbed0455807