摘要
稀疏脉冲反褶积方法对提高地震资料分辨率有着重要作用,迭代阈值收缩算法(ISTA)是其核心算法,首先利用地震数据提取子波,再利用ISTA求解反射系数.当地震子波提取不准确时,反褶积效果不理想.为此,在ISTA基础上,结合循环神经网络(RNN)中反向传播(BPTT)的思想,研究形成了一种类RNN的改进ISTA稀疏脉冲反褶积方法.该算法首先使用常规手段从实际地震数据中提取地震子波,构建反褶积的子波字典;然后将构建的地震子波字典作为已知的初始条件,结合ISTA求取的反射系数;再根据BPTT算法思想,将求取的反射系数与子波褶积并与实际数据进行比较,反向修改地震子波;最终,经过多次迭代修改获得合理的地震子波字典,并利用该地震子波字典求解实际地震数据的反射系数序列.为验证算法的有效性,采用不同信噪比的理论地震记录,给定存在较大误差的初始子波,进行了反褶积计算.采用传统的ISTA和类RNN的改进ISTA进行对比处理,结果表明,改进ISTA具有较好的抗噪能力和子波自适应能力,可使实测地震资料的有效频带拓展约1.5倍,能够较好地适应实际地震资料的反褶积处理.
The sparse spike deconvolution plays an important role in improving the resolution of seismic data.A successful deconvolution requires accurate wavelet data for the calculation of the reflection coefficient,which is performed using the core algorithm of the Iterative Shrinkage Thresholding Algorithm (ISTA).An inaccurate extraction of wavelet data can affect the result of the deconvolution.In this study,an RNN-like ISTA algorithm is proposed,which combines the concept of Back-Propagation Through Time (BPTT) in Recurrent Neural Network (RNN) with the traditional ISTA algorithm.First,seismic wavelets are extracted from the actual seismic data to construct a wavelet dictionary,which is taken as the initial condition,and the reflection coefficient is calculated using ISTA algorithm.Subsequently,the reflection coefficient is convoluted with the wavelets to construct seismic data by means of the BPTT algorithm,and a seismic wavelet correction is performed by comparing the obtained seismic data with the actual data.Finally,a reasonable seismic wavelet dictionary is obtained after several iterations,which can be used to calculate the actual reflection coefficient series.Tests on theoretical seismic records with different signal-to-noise ratio showed that the improved algorithm has a better anti-noise and wavelet adaptive abilities than the conventional ISTA algorithm.Moreover,the frequency band of the actual seismic data can be expanded by about 1.5 times.These results demonstrate that the proposed method can be successfully applied to the deconvolution of actual seismic data.
作者
潘树林
闫柯
杨海飞
蒋从元
秦子雨
PAN Shulin;YAN Ke;YANG Haifei;JIANG Congyuan;QIN Ziyu(School of Earth Science and Technology,Southwest Petroleum University,Chengdu 610500,China;North Operation Branch of Changqing Oilfield,PetroChina,Yulin 719000,China;Department of Electrical and Electronic Engineering,Sichuan Vocational and Technical College,Suining 629000,China)
出处
《石油物探》
EI
CSCD
北大核心
2019年第4期533-540,共8页
Geophysical Prospecting For Petroleum
基金
国家自然基金项目(NSFC 41674095)
“油气藏地质及开发工程”国家重点实验室开放基金项目(PLN201733)
天然气地质四川省重点实验室开放基金项目(2015trqdz03)共同资助~~
关键词
稀疏脉冲反褶积
分辨率
ISTA
地震子波
信噪比
循环神经网络
反向传播
sparse spike deconvolution
resolution
Iterative Shrinkage Thresholding Algorithm (ISTA)
seismic wavelet
SNR
Recurrent Neural Network (RNN)
Back-Propagation Through Time (BPTT)