基于ISTA算法的稀疏约束反演谱分解及应用

doi:10.7657/XJPG20200413

新疆石油地质 ›› 2020, Vol. 41 ›› Issue (4): 471-476.doi: 10.7657/XJPG20200413

基于ISTA算法的稀疏约束反演谱分解及应用

高秋菊^1a(), 张云银^1b, 曲志鹏^1a, 徐彦凯²(), 王宗家^1a, 王千军^1b

1.中国石化胜利油田分公司 a.物探研究院;b.勘探开发研究院,山东东营 257000;
2.中国石油大学（北京）信息科学与工程学院,北京 102249

收稿日期:2019-04-15 修回日期:2019-10-31 出版日期:2020-08-01 发布日期:2020-08-05
通讯作者: 徐彦凯 E-mail:527266669@qq.com;xuyk163@163.com
作者简介:高秋菊（1972-）,女,辽宁北镇人,教授级高级工程师,油气勘探,（Tel）0546-8791583（E-mail）527266669@qq.com
基金资助:
国家科技重大专项(2017ZX05072)

Decomposition and Application of Constrained Sparse Inversion Spectrum Based on ISTA Algorithm

GAO Qiuju^1a(), ZHANG Yunyin^1b, QU Zhipeng^1a, XU Yankai²(), WANG Zongjia^1a, WANG Qianjun^1b

1. Sinopec Shengli Oilfield Company, a.Institute of Geophysics; b.Research Institute of Exploration and Development, Dongying, Shandong 257000, China;
2. School of Geosciences and Information Engineering, China University of Petroleum, Beijing 102249, China

Received:2019-04-15 Revised:2019-10-31 Online:2020-08-01 Published:2020-08-05
Contact: XU Yankai E-mail:527266669@qq.com;xuyk163@163.com

摘要/Abstract

摘要：

常规谱分解的分辨率难以满足地震解释的需求,稀疏约束反演谱分解可以很好地解决该问题。稀疏约束反演谱分解是把地震信号看成已知的子波矩阵库和伪反射系数的褶积,从而将谱分解问题转化为一个反问题来求解,其核心内容是如何快速获得最优解。采用L₁范数正则化的L₂范数作为稀疏约束反演谱分解的目标函数,并通过迭代阈值算法求解稀疏反问题。为了进一步提高计算速度,基于Ricker子波构建一种新的算子,并用ISTA算法进行计算。在此基础上,将稀疏约束反演谱分解用于模拟信号的数值试验,并与常规谱分解结果进行比较。结果表明,稀疏约束反演谱分解具有更好的时频聚集性和更高的时频分辨率;进一步应用于济阳坳陷沾化凹陷渤南洼陷义176井区地震资料,处理结果对油气响应十分敏感,可以较好地识别油气储集层。

关键词: ISTA算法, 稀疏约束反演, 谱分解, 连续小波变换, 油气检测, L₁范数正则化, Ricker子波, 低频阴影

Abstract:

The resolution of conventional spectral decomposition can’t meet the needs of seismic interpretation, which can be well solved by the spectral decomposition of constrained sparse inversion. In this method, seismic signal is regarded as the convolution of known wavelet matrix library and pseudo-reflection coefficient, then the spectral decomposition is transformed into an inverse problem which focuses on how to obtain an optimum solution. The L₂ norm regularized by L₁ norm is used as the objective function for the constrained sparse inversion spectrum decomposition, then the iterative threshold algorithm is used to obtain the solution of the inverse problem. In order to improve the calculation speed of spectral decomposition, a new kind of operator is established on the basis of Ricker wavelet and the optimal solution can be obtained by ISTA algorithm. Based on which, the spectral decomposition of constrained sparse inversion is applied into numerical modeling and the results are compared with those of conventional spectral decomposition. The results show that the spectral decomposition of constrained sparse inversion has high time-frequency resolution. The actual application of the algorithm in the Wellblock Yi176 of Bonan sag in Jiyang depression indicates that the processing result is sensitive to oil and gas, which can be used to identify oil and gas reservoirs.

Key words: ISTA algorithm, constrained sparse inversion, spectral decomposition, continuous wavelet transform, hydrocarbon detection, L₁ norm regularization, Ricker wavelet, low-frequency shadow

中图分类号:

P631

高秋菊, 张云银, 曲志鹏, 徐彦凯, 王宗家, 王千军. 基于ISTA算法的稀疏约束反演谱分解及应用[J]. 新疆石油地质, 2020, 41(4): 471-476.

GAO Qiuju, ZHANG Yunyin, QU Zhipeng, XU Yankai, WANG Zongjia, WANG Qianjun. Decomposition and Application of Constrained Sparse Inversion Spectrum Based on ISTA Algorithm[J]. Xinjiang Petroleum Geology, 2020, 41(4): 471-476.

图/表 5

参考文献 19

[1]	杨子川, 陶磊, 李海英, 等. 塔河油田S48井区流体多波属性综合预测[J]. 新疆石油地质, 2018,39(6):742-747.
	YANG Zichuan, TAO Lei, LI Haiying, et al. Comprehensive fluid prediction with multi-wave attribute seismic technology in Wellblock S48, Tahe oilfield[J]. Xinjiang Petroleum Geology, 2018,39(6):742-747.
[2]	彭军, 周家雄, 马光克, 等. 高精度三参数小波谱分析法在薄储集层预测中的应用[J]. 新疆石油地质, 2018,39(2):230-235.
	PENG Jun, ZHOU Jiaxiong, MA Guangke, et al. Application of high-precision three-parameter wavelet spectrum analysis in thin reservoir prediction[J]. Xinjiang Petroleum Geology, 2018,39(2):230-235.
[3]	MARFURT K J, KIRLIN R L. Narrow-band spectral analysis and thin-bed tuning[J]. Geophysics, 2001,66(4):1 274-1 283.
[4]	SUN S, CASTAGNA J P, SIEGFRIED R W. Examples of enhanced spectral processing in direct hydrocarbon detection[R]. Houston,Texas,USA：72nd Annual International Meeting,SEG,Expanded Abstracts, 2002: 457-460.
[5]	CASTAGNA J P, SUN S, SIEGFRIED R W. Instantaneous spectral analysis:detection of low-frequency shadows associated with hydrocarbons[J]. The Leading Edge, 2003,22(1):120-127.
[6]	EBROM D. The low-frequency gas shadow on seismic sections[J]. The Leading Edge, 2004,23(8):772-779.
[7]	GABOR D. Theory of communication[J]. Journal of the IEEE, 1946,93(26):429-457.
[8]	SINHA S, ROUTH P, ANNO P, et al. Spectral decomposition of seismic data with continuous-wavelet transform[J]. Geophysics, 2005,70(6):19-25.
[9]	PINNEGAR C R, MANSINHA L. The S-transform with windows of arbitrary and varying shape[J]. Geophysics, 2003,68(10):381-385.
[10]	WU X Y, LIU T Y. Seismic spectral decomposition and analysis based on Wigner-Ville distribution for sandstone reservoir characterization in west Sichuan depression[J]. Geophysics, 2010,75(2):126-134.
[11]	MALLAT S G, ZHANG Z F. Matching pursuits with time frequency dictionaries[J]. IEEE Transactions on Signal Processing, 1993,41(12):3 397-3 415.
[12]	ZHANG X W, HAN L G, WANG Y, et al. Seismic spectral decomposition fast matching pursuit algorithm and its application[J]. Geophysical Prospecting for Petroleum, 2010,49(1):1-6.
[13]	PORTNIAGUINE O, CASTAGNA J. Inverse spectral decomposition[R]. Denver,Colorado,USA：74th Annual International Meeting,SEG,Expanded Abstracts, 2004: 1 786-1 789.
[14]	BONAR D C, SACCHI M D. Complex spectral decomposition via inversion strategies[R]. Denver,Colorado,USA：80th Annual International Meeting,SEG,Expanded Abstracts, 2010: 1 408-1 412.
[15]	张生强, 韩立国, 李才, 等. 基于高分辨率反演谱分解的储层流体流度计算方法研究[J]. 石油物探, 2015,54(2):142-149.
	ZHANG Shengqiang, HAN Liguo, LI Cai, et al. Computation method for reservoir fluid mobility based on high-resolution inversion spectral decomposition[J]. Geophysical Prospecting for Petroleum, 2015,54(2):142-149.
[16]	黄广谭, 李景叶, 陈小宏, 等. 基于稀疏约束反演谱分解的频变AVO储层识别方法[J]. 地球物理学报, 2017,60(10):3 984-3 994.
	HUANG Guangtan, LI Jingye, CHEN Xiaohong, et al. Frequency-dependent AVO reservoir identification method based on sparse constrained inversion spectral decomposition[J]. Chinese Journal of Geophysics, 2017,60(10):3 984-3 994.
[17]	李倩, 狄帮让, 魏建新. 基于稀疏约束反演谱分解的缝洞储层叠后数据去噪应用效果分析[J]. 石油物探, 2017,56(5):684-693.
	LI Qian, DI Bangrang, WEI Jianxin. Application of denoising method based on sparse constrained inverse spectral decomposition in poststack seismic data of cave-fractured reservoirs[J]. Geophysical Prospecting for Petroleum, 2017,56(5):684-693.
[18]	DONOHO D. For most large underdetermined systems of linear equations,the minimalell-1 norm solution is also the sparsest solution[J]. Communications on Pure and Applied Mathematics, 2006,59(7):907-934.
[19]	DAUBECHIES I, DEFRISE M, MOL C. An iterative thresholding algorithm for linear inverse problems with a sparsity constraint[J]. Communications on Pure and Applied Mathematics, 2004,57(11):1 413-1 457.

基于ISTA算法的稀疏约束反演谱分解及应用

Decomposition and Application of Constrained Sparse Inversion Spectrum Based on ISTA Algorithm

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表 5

参考文献 19

相关文章 11

编辑推荐

Metrics

本文评价

[1]	丁传奇, 周路, 钟克修, 吴勇, 钟斐艳, 刘怡. 川东大猫坪地区长兴组生物礁储集层预测[J]. 新疆石油地质, 2020, 41(2): 164-171.
[2]	崔震，曹思远，刘伟，马媛媛. 基于变窗函数谱重排的谱分解技术[J]. , 2015, 36(5): 1-1.
[3]	兰锋，李明，李艳东，田瀚. 地震低频阴影的数值模拟分析[J]. , 2014, 35(4): 1-1.
[4]	熊冉，赵继龙，厚刚福. 用频谱分解和地震峰值属性分析预测薄砂岩储集层[J]. , 2013, 34(2): 1-1.
[5]	刘伟, 贺振华, 陈学华, 翁斌, 姜秀娣. 浊积体的低频伴影油气检测[J]. 新疆石油地质, 2011, 32(6): 609-612.
[6]	胡太平, 吉云刚, 潘杨勇, 董瑞霞, 吴梅莲. 傅立叶频谱分解技术在塔中45 井区油气勘探中的应用[J]. 新疆石油地质, 2011, 32(3): 308-310.
[7]	董文波, 吴雨韩, 吴采西, 毛丹凤. 应用地震解释技术识别湖底扇——以克拉玛依油田为例[J]. 新疆石油地质, 2011, 32(2): 183-184.
[8]	齐宇, 刘震, 魏建新, 狄帮让, 李娟. 基于小波变换的谱分解技术在地震模型解释中的应用[J]. 新疆石油地质, 2010, 31(4): 417-418.
[9]	刘兰锋, 曹思远, 夏竹, 熊伟. 广义 S 变换及其在储集层频谱成像中的应用[J]. 新疆石油地质, 2009, 30(5): 629-632.
[10]	周怀来, 李录明. 广义S 变换在地震信号特征信息提取中的应用[J]. 新疆石油地质, 2008, 29(6): 758-760.
[11]	胡水清, 韩大匡, 刘文岭, 王玉学. 应用谱分解技术预测河道砂体[J]. 新疆石油地质, 2008, 29(3): 370-372.