基于迭代再加权最小二乘的地震资料稀疏反演方法.pdf

第 卷 第 期 地 球 物 理 学 进 展 ， ２８ ５ Ｖｏｌ．２８ Ｎｏ．５ 　 　 　 　 　 　 年 月（页码： ） ， ２０１３ １０ ２５３６２５４６ ＰＲＯＧＲＥＳＳ ＩＮ ＧＥＯＰＨＹＳＩＣＳ Ｏｃｔ． ２０１３ 　 　 袁义明，孙 晨，杨长春 基于迭代再加权最小二乘的地震资料稀疏反演方法 地球物理学进展， ， （）： ， ： ． ． ２０１３２８５ ２５３６２５４６ｄｏｉ 　 ／ １０．６０３８ ｐｇ ， ， ＹＵＡＮ Ｙｉｍｉｎ ＳＵＮ Ｃｈｅｎ ＹＡＮＧＣｈａｎｃｈｕｎ．Ｓｅｉｓｍｉｃｄａｔａｓａｒｓｅｉｎｖｅｒｓｉｏｎｍｅｔｈｏｄｂａｓｅｄｏｎｉｔｅｒａｔｉｖｅｒｅｗｅｉｈｔｅｄ ｇ ｇ ｐ ｇ （ ）， ， （）： ， ： ／ ｔｅｃｈｎｉｕｅ．犘狉狅狉犲狊狊犻狀犌犲狅犺狊． ｉｎＣｈｉｎｅｓｅ ２０１３ ２８５ ２５３６２５４６ｄｏｉ１０．６０３８ ｑ 犵 狆狔 ｐｇ 基于迭代再加权最小二乘的地震资料稀疏反演方法 ， ， １３ ２３ ２ 袁义明 ， 孙 晨 ， 杨长春 　 　 　 （ 中国科学院地质与地球物理研究所兰州油气资源研究中心，兰州 ； １． ７３００００ ２．油气资源研究重点实验室，中国科学院地质与地球物理研究所，北京 １００２９； ３．中国科学院研究生院，北京 １０００４９） 摘 要 为了由地震数据得到高分辨率的反射系数序列，本文提出了使用迭代再加权最小二乘的寻优方法；该方 　 　 法通过再加权的方法将稀疏优化问题转化为一系列权不断变化的线性方程组，在保证精度的前提下增加了求解的 速度 对寻优过程中的线性方程组求解采用针对性的预条件化，并在权重的计算中引入伸缩变换，优化了整个寻优 ． 过程的稳定性和求解速度 经过大量的数值试验，得出了合理的参数选择范围，通过与真值的比较以及实际资料处 ． 理验证了本方法的实用性． 关键词 波阻抗反演，迭代再加权，预条件化，伸缩变换 　 ： ／ 中图分类号 文献标识码 ｄｏｉ１０．６０３８ Ｐ６３１ Ａ 　　　　　 　 　　　　 　 ｐｇ 犛犲犻狊犿犻犮犱犪狋犪狊犪狉狊犲犻狀狏犲狉狊犻狅狀犿犲狋犺狅犱