基于Bregman迭代的不...值模型地震数据重构方法分析_赵虎.pdf
《基于Bregman迭代的不...值模型地震数据重构方法分析_赵虎.pdf》由会员分享,可在线阅读,更多相关《基于Bregman迭代的不...值模型地震数据重构方法分析_赵虎.pdf(10页珍藏版)》请在咨信网上搜索。
1、2023 年第 38 卷 第1期2023,38(1):0409-0418地球物理学进展Progress in Geophysicshttp:/wwwprogeophyscnISSN 1004-2903CN 11-2982/P赵虎,赵子涵,陈伟,等 2023 基于 Bregman 迭代的不同阈值模型地震数据重构方法分析 地球物理学进展,38(1):0409-0418,doi:10 6038/pg2023FF0310ZHAO Hu,ZHAO ZiHan,CHEN Wei,et al 2023 Analysis of seismic data reconstruction methods for d
2、ifferent threshold models based on Bregmaniteration Progress in Geophysics(in Chinese),38(1):0409-0418,doi:106038/pg2023FF0310基于 Bregman 迭代的不同阈值模型地震数据重构方法分析Analysis of seismic data reconstruction methods for different threshold modelsbased on Bregman iteration赵虎1,2,赵子涵2,陈伟3,邸志欣4,韩嵩5,胥良君6,杨暾2ZHAO Hu1
3、,2,ZHAO ZiHan2,CHEN Wei3,DI ZhiXin4,HAN Song5,XU LiangJun6,YANG Tun2收稿日期2021-10-19;修回日期2022-08-31投稿网址http:/www progeophys cn基金项目中国石油-西南石油大学创新联合体科技合作项目(2020CX010201)和国家自然基金青年科学基金(41704134)联合资助第一作者简介赵虎,男,1983 年生,博士,教授,现在西南石油大学从事地震采集方法与解释教学和科研工作E-mail:cumtzhaohu163 com1 西南石油大学天然气地质四川省重点实验室,成都6105002 西南
4、石油大学地球科学与技术学院,成都6105003 中石油西南油气田公司勘探事业部,成都6100414 中石化石油工程地球物理有限公司科技研发中心,南京2100005 中石油西南油气田公司勘探开发研究院,成都6100416 中石油西南油气田分公司重庆气矿,重庆4007071 Sichuan Province University Key Laboratory of Natural Gas Geology,SWPU,Chengdu 610500,China2 School of Geoscisence and Technology,Southwest Petroleum University,Che
5、ngdu 610500,China3 Petrochina Southwest Oil and Gas field Company Exploration Division,Chengdu 610041,China4 Sinopec Geophysical Corporation,esearch Center,Nanjing 210000,China5 Exploration and Development esearch Institute,PetroChina Southwest Oil and Gasfield Company,Chengdu 610041,China6 PetroChi
6、na Southwest Oil and Gasfield Company Chongqing Gas District,Chongqing 400707,China摘要传统的地震数据采样必须严格遵循 Nyquist采样定理,而野外实际数据的采集可能由于施工条件或者地表障碍物的限制,不一定能记录到完整的地震波场,所以地震资料处理中的数据重建是非常重要的问题 压缩感知理论最先来自信号处理领域,它所包括的问题类型有信号的稀疏表征和数学组合优化,它给地震数据重建这类问题指明了思考方向 而其中如何选择最优的迭代算法是数据重建中的关键问题 本文将地震数据插值问题归纳到约束最优化问题,选择能有效稀疏表征地
7、震波场的傅里叶变换,对于压缩感知理论框架下的混合范数反问题,再用 Bregman 迭代方法去求解,在地震数据的重建过程中,传统的阈值参数收敛慢,为了降低迭代次数并且提高地震数据恢复的AbstractTraditional seismic data sampling must strictlyfollow the Nyquist sampling theorem,and actual field dataacquisition may not be able to record the complete seismicwave field due to construction conditio
8、ns or obstacles onthe ground Therefore,data reconstruction in seismic dataprocessing is an important factor problem Compressedsensingtheoryoriginatesfromthefieldofsignalprocessing It includes sparse representation of signals andmathematicalcombinationoptimizationproblemsItprovides a new solution for
9、 the reconstruction of seismicdata And how to choose the optimal iterative algorithm isthekeyissueindatareconstructionThispapersummarizes the seismic data interpolation problem into aconstrainedoptimizationproblem,choosesFouriertransform,which can effectively represent the seismic wavefield sparsely
10、,and then applies the Bregman iterativemethod to solve the mixed norm inverse problem under the地球物理学进展www progeophys cn2023,38(1)精度,总结出改进型指数衰减规律的阈值参数,选择用硬阈值算子来重建恢复地震数据 通过对理论模型和实际地震资料的处理结果表明该方法可以快速、有效的恢复地震波场的缺失数据关键词地震数据插值;Bregman 迭代;阈值参数;压缩感知中图分类号P631文献标识码Adoi:10 6038/pg2023FF0310frameworkofcompresse
11、dsensingtheoryInthereconstructionprocess,forthedisadvantageofslowconvergenceoftraditionalthresholdparameters,thethreshold parameters of the improved exponential decay laware summarizedIn order to reduce the number ofiterations andimprovetheaccuracyofseismicdatareconstruction,the hard threshold opera
12、tor is selected toperform the two-dimensional seismic dataebuild andrestore The processing results of theoretical model andactual seismic data show that this method can quickly andeffectively restore the missing data of seismic wave fieldKeywordsSeismic data interpolation;Bregman iteration;Threshold
13、 parameter;Compressed sensing0引言随着我国油气勘探程度的提高,如何提高复杂勘探区的地震勘探效果成为下一步工作的关键 在野外地震数据采集的进程中,因为受限于采集环境和经济条件,导致某些位置地震数据缺失,使得整个地震剖面不完整 从而会对室内资料处理和解释造成不良影响,因此如何进行地震数据规则化重建显得尤为关键地震数据重建问题,最早由 Larner 等(1981)提出 他最开始对不完整地震道的恢复进行深入的研究,之后陆续出现了不同的地震数据重建方法;onen(1987)提出把缺失道视为零道再代入波动方程部分偏移的叠前地震数据重构方法;Spitz(1991)提出了线性同相
14、轴在 f-k 域是能够被预测的地震数据重建方法;随后 Porsani(1999)对这一方法提出了改进,得出 f-k 域半步长预测滤波的地震道插值方法;Naghizadeh 和 Sacchi(2007)提出了基于多步自回归预测滤波的不规则缺失道插值重建;Kaplan 等(2010)采用最小平方反演方法进行数据重建 不难发现以上方法基本都是针对均匀采样数据进行重建而近年来提出的压缩感知理论对于解决此类问题 指 明 了 新 的 方 向,压 缩 感 知 由 EmmanuelCandes、omberg 和 Donoho 等人在 2004 年正式提出 基于压缩感知的地震勘探技术,目前国内外已经做了一些研究
15、,Herrmann 等(2008)、Herrmann(2010)提出了一种基于 Curvelet 的多尺度非线性地震数据处理方法 Candes 等(2006)基于鲁棒不确定性原理对高度不完整频率信息的精确信号重建进行了研究 Hennenfent 和 Herrmann(2008)提出了基于随机抖动的欠采样方案,对相邻观测点之间的最大距离进行控制,便于实际施工 Wu 等(2009)用 dreamlet 方 法 在 压 缩 域 进 行 成 像Naghizadeh 和 Sacchi(2010)研究了采样函数和傅里叶重建方法 Mosher 等(2012b)首次提出并且完善了一种根据约束条件选择最佳炮点和
16、检波点位置的非均匀采样方法,用于观测系统设计,取得了显著的成效 Mosher 等(2012a)对压缩感知成像这一方向有一个更深层次的探究 Brown 等(2017)在2015 年将非规则排列观测系统应用在阿拉斯加地区的陆上冰雪覆盖区域,通过利用数据重建的方法,得到一个相比传统采集方法更好的高分辨率地震资料数据重构中最为关键的问题是求解策略,常用方法是基于某种变换,Xu 等(2005)提出了基于反泄露傅里叶变换的地震数据重建方法 Jin(2010)提出了基于阻尼最小范数 Fourier 反演,实现了五维地震数 据 的 重 建 Herrmann 在 2008 年 提 出 基 于Curvelet 数
17、据重建的稀疏促进反演方法 Liu 和Fomel(2010)提出稀疏多尺度变换域中的 OC-seislet变换重建缺失的地震数据并消除随机噪声 Yang 和Fomel(2015)阐述了用快速 Fourier 反演策略代替非均匀快速 Fourier 变换的观点,采用共轭梯度法来优化最优化问题,用来提高计算的效率 还有复值曲波变换(徐卫等,2016)和 Shearlet 稀疏变换基(王常波,2018)对于地震数据重建的迭代算法,主要有凸集投影算法,最小加权范数法,不准确 Uzawa 方法,两步法以及迭代阈值法(IST)等 但是,已有的迭代算法往往难以同时兼顾计算精度和效率 为此,本文拟引用 Breg
18、man 迭代算法,分析不同迭代算法的重建特性,通过提出改进型指数衰减规律的阈值参数来改善在地震数据恢复过程中,传统的阈值参数收敛得较为缓慢这一不佳的表现,运用理论模型和实际地震资料经过方法验证,有效解决地震数据缺失数据重建问题0142023,38(1)赵虎,等:基于 Bregman 迭代的不同阈值模型地震数据重构方法分析(www progeophys cn)1方法原理1.1理论基础压缩感知理论和传统的 Nyquist 采样不同,它是同时在进行着采样和压缩这两个步骤,不被Nyquist 采样定理所约束,而是受控于稀疏性和非相干性这两个基本准则 采用稀疏变换去构建数据空间和模型空间之间的联系 对地
19、震数据进行重构,就是把残缺的地震数据构建成完整的数据,模型如下:d=m,(1)其中,d 为采集到的地震数据,m 为插值重建后的地震数据,为稀疏变化算子 利用信号具有稀疏性的先验条件,通过一些数学变换(Curvelet 变换、Seislet 变换和 Fourier 变换等)来构建数据空间和模型空间的联系 由于公式(1)为数学 g 问题,要通过增加另外的正则化条件来求得合理的 m1.2Bregman 迭代重建算法根据上述分析可得出,压缩感知方法是利用最小化策略这一关键的技术来对前面的欠定问题进行求解的 与 L2准则不同,压缩感知理论是通过稀疏反变换算子来构建一个稀疏的数据空间 m(其中m=Tm),
20、用来解决下面的约束最优化问题:minmm1s tKTm=d(2)该模型是利用求解带约束的 L1范数最小化问题,来求取线性系统 KTm=d 的稀疏解 通常时候,考虑用到最小二乘法来求解此问题:m=argminmKTm d22(3)将公式(2)的带约束最小化问题转化为公式(3)的无约束问题后,需要在公式(3)的右边添加正则项 m1:minm m1+12KTm d22,(4)其中 为惩罚参数 但事实上,公式(4)与公式(2)并不等价,为了满足 KTm=d 的条件,需要取为较小值Bregman 迭代方法是 Osher 等(2005)于图像处理领域提出来的方法,该方法的迭代框架为:dk+1=dk+(d
21、KTmk)mk+1=argminm m1+12KTm dk+122(5)求解公式(5)里面的非约束问题能够采用迭代硬阈值方法,再结合稀疏正变换算子,用于解决数据插值重建问题的 Bregman 迭代方法为:dk+1=dk+(d KTmk)mk+1=S T1dk+1,(6)其中,m0=0,d0=0,S 为硬阈值算子 这样,最终插值结果 m 即可通过 m=Tmk+1来求得 应用Bregman 迭代方法时,关键问题在于公式(6)中阈值参数的选取,接下来,我们讨论该参数的选取方法图 1阈值模型迭代次数变化图Fig 1Threshold model iteration number change grap
22、h图 2四种阈值参数的评价函数 J 曲线Fig 2Evaluation function J curve of four threshold parameters图 3重构数据的评价函数 J 值1 为线性、2 为指数型、3 为反比例型、4 为改进型指数Fig 3Evaluation function J value of reconstructed data1 is linear,2 is exponential,3 is inverse proportional,4 is improved exponent114地球物理学进展www progeophys cn2023,38(1)图 4(a)
23、理论模型;(b)理论数据 FK 谱Fig 4(a)Theoretical model;(b)The FK spectrum of theoretical data图 5(a)50%随机欠采样数据;(b)欠采样数据 FK 谱Fig 5(a)50%random undersampled data;(b)The FK spectrum of undersampled data图 6(a)线性阈值参数重建结果;(b)重建结果 FK 谱Fig 6(a)Linear threshold parameter reconstruction results;(b)The FK spectrum of recon
24、struction result2142023,38(1)赵虎,等:基于 Bregman 迭代的不同阈值模型地震数据重构方法分析(www progeophys cn)图 7(a)指数型阈值参数重建结果;(b)重建结果 FK 谱Fig 7(a)Exponential threshold parameter reconstruction results;(b)The FK spectrum of reconstruction result图 8(a)反比例型阈值参数重建结果;(b)重建结果 FK 谱Fig 8(a)econstruction results of inverse proportio
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 Bregman 模型 地震 数据 方法 分析 赵虎
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。