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

Seismic data sparse inversion method based on iterative reweighted technique

作　　者：袁义明[1,2] 孙晨[3,2] 杨长春[3] 

机构地区：[1]中国科学院地质与地球物理研究所兰州油气资源研究中心,兰州730000 [2]中国科学院研究生院,北京100049 [3]油气资源研究重点实验室,中国科学院地质与地球物理研究所,北京10029

出　　处：《地球物理学进展》2013年第5期2536-2546,共11页Progress in Geophysics

摘　　要：为了由地震数据得到高分辨率的反射系数序列,本文提出了使用迭代再加权最小二乘的寻优方法;该方法通过再加权的方法将稀疏优化问题转化为一系列权不断变化的线性方程组,在保证精度的前提下增加了求解的速度.对寻优过程中的线性方程组求解采用针对性的预条件化,并在权重的计算中引入伸缩变换,优化了整个寻优过程的稳定性和求解速度.经过大量的数值试验,得出了合理的参数选择范围,通过与真值的比较以及实际资料处理验证了本方法的实用性.The purpose of this paper is to invert a high resolution reflection coefficients series from seismic data using iterative reweighted technique. This technique transform the sparse optimization into a series of wighted linear function, speeds up the computing while keep signal-to-noise rate high. In order to stabilize and accelerate the solving, pre-eonditioning according to the linear function is used and shrink operator is introduced in computing the weighing. Through lots of numerical experimentations, the range of parameters is settled, and the practicability is conformed together with the real data processing.

关 键 词：波阻抗反演 迭代再加权 预条件化 伸缩变换 

分 类 号：P631[天文地球—地质矿产勘探]