世界地质 2017, 36(2) 588-594 DOI:   10.3969/j.issn.1004-5589.2017.02.025  ISSN: 1004-5589 CN: 22-1111/P

本期目录 | 下期目录 | 过刊浏览 | 高级检索                                                            [打印本页]   [关闭]
地球物理
扩展功能
本文信息
Supporting info
PDF(3156KB)
[HTML全文]
参考文献[PDF]
参考文献
服务与反馈
把本文推荐给朋友
加入我的书架
加入引用管理器
引用本文
Email Alert
文章反馈
浏览反馈信息
本文关键词相关文章
全波形反演
互相关
最小二乘
本文作者相关文章
PubMed
互相关与最小二乘加权目标函数全波形反演
梁煌, 韩立国, 许卓, 胡勇, 邹佳儒
吉林大学地球探测科学与技术学院, 长春 130026
摘要

全波形反演是一个高度非线性的优化问题,当地震数据中缺少低频成分而初始速度与真实速度相差较远时,反演容易陷入局部极小值。笔者提出一种新的目标函数,将模拟地震记录和观测记录的归一化互相关与最小二乘结合。互相关侧重相位匹配,具有更强的线性,能减弱“跳周”现象。通过设置权重因子,在反演前期利用互相关先恢复低波数的背景速度模型,再加入最小二乘约束恢复高波数的模型细节。数值模拟试验结果表明,基于该目标函数的全波形反演不依赖精确的初始模型和低频信息,向全局极小值迅速收敛,能有效改善反演的稳定性,并获得比基于常规目标函数的全波形反演更精确的结果。

关键词 全波形反演   互相关   最小二乘  
Full waveform inversion based on weighted cross-correlation and least squares objective function
LIANG Huang, HAN Li-guo, XU Zhuo, HU Yong, ZOU Jia-ru
College of Geo-exploration Science and Technology, Jilin University, Changchun 130026, China
Abstract:

Full waveform inversion (FWI) is a highly non-linear optimization problem. When there is a lack of low-frequency components in the seismic data and the initial velocity is far from the true velocity, it is easily trapped into local minima. The authors propose a new objective function that combines the normalized cross-correlation of modeled and observed data with the least squares criterion. Cross-correlation emphasizes on phase matching and behaves in a more linear way, thus can mitigate the cycle-skipping issue. By setting a weighting factor, the authors use the cross-correlation norm to update the low-wave number components of the velocity model at the early stage of the inversion, while using the least squares criterion to update high-wave number details of the model. Numerical examples show that FWI based on the new objective function converges quickly toward global minima and does not rely on accurate initial velocity model as well as low-frequency information, it can improve the robustness of the inversion and generate more accurate inversion results than the conventional approach.

Keywords: full waveform inversion   cross-correlation   least squares  
收稿日期 2016-09-20 修回日期 2016-12-02 网络版发布日期  
DOI: 10.3969/j.issn.1004-5589.2017.02.025
基金项目:

国家自然科学基金项目(41304086)与国家高技术研究发展规划(863计划)重大课题项目联合资助.

通讯作者: 韩立国(1961),男,教授,博士生导师,主要从事地震数据处理解释工作.E-mail:hanliguo@jlu.edu.cn
作者简介:
作者Email: hanliguo@jlu.edu.cn

参考文献:

[1] Lailly P. The seismic inverse problem as a sequence of before stack migrations [C]//Conference on inverse scattering: theory and application. Society for Industrial and Applied Mathematics, 1983: 206-220.
[2] Tarantola A. Inversion of seismic reflection data in the acoustic approximation[J]. Geophysics, 1984, 49(8): 1259-1266.
[3] Tarantola A. A strategy for nonlinear elastic inversion of seismic reflection data[J]. Geophysics, 1986, 51(10): 1893-1903.
[4] Shin C, Cha Y H. Waveform inversion in the Laplace domain[J]. Geophysical Journal International, 2008, 173: 922-931.
[5] Warner M, Guasch L. Adaptive waveform inversion-FWI without cycle skipping-theory[C]//76th EAGE Conference and Exhibition, Expanded Abstracts, 2014.
[6] 白璐,韩立国,张盼,等. 最小平方滤波法时间域全波形反演[J]. 石油地球物理勘探,2016, 51(4): 721-729. BAI Lu, HAN Li-guo, ZHANG Pan, et al. Time-domain full waveform inversion based on least square filter [J]. Oil Geophysical Prospecting, 2016, 51(4): 721-729.
[7] Hu Y, Han L G, Zhang P, et al. Multi-step full waveform inversion based on waveform mode decomposition [C]//87th Annual International Meeting, SEG Expanded Abstracts, 2016: 1501-1505.
[8] Luo Y, Schuster G T. Wave-equation traveltime inversion[J]. Geophysics, 1991, 56: 645-653.
[9] Leeuwen T V, Mulder W A. A correlation-based misfit criterion for wave-equation traveltime tomography [J]. Geophysical Journal International, 2010, 182: 1383-1394.
[10] Luo S, Sava P. A deconvolution-based objective function for wave-equation inversion[C]//81st Annual International Meeting, SEG Expanded Abstracts, 2011:2788-2792.
[11] Choi Y, Alkhalifah T. Application of multi-source waveform inversion to marine streamer data using the global correlation norm [J]. Geophysical Prospecting, 2012, 60: 748-758.
[12] 杨贺龙,韩立国,陈雪. 波场相位相关时移全波形反演[J]. 石油物探,2016,55(4):568-575. YANG He-long, HAN Li-guo, CHEN Xue. Wavefield phase correlation shifting full waveform inversion [J]. Geophysical Prospecting for Petroleum, 2016, 55(4): 568-575.
[13] 胡光辉,贾春梅,夏洪瑞,等. 3D声波全波形反演的实现与验证[J]. 石油物探,2013,52(4): 417-425. HU Guang-hui, JIA Chun-mei, XIA Hong-rui, et al. Realization and validation of 3D acoustic full waveform inversion [J]. Geophysical Prospecting for Petroleum, 2013, 52(4): 417-425.
[14] 孟旭,刘四新,吴俊军,等. 时间域跨孔雷达全波形反演及实际应用[J]. 世界地质,2016,35(1): 256-263. MENG Xu, LIU Si-xin, WU Jun-jun, et al. Full waveform inversion of time-domain cross-hole radar and its field application [J]. Global Geology, 2016, 35(1): 256-263.
[15] 孙慧秋,韩立国,徐杨杨,等. 基于L-BFGS算法的全波形反演Q值方法研究[J]. 世界地质,2015,34(4): 1106-1112. SUN Hui-qiu, HAN Li-guo, XU Yang-yang, et al. Method of Q through full waveform inversion based on L-BFGS algorithm [J]. Global Geology, 2015, 34(4): 1106-1112.
[16] Romero L A, Ghiglia D C, Ober C C, et al. Phase encoding of shot records in prestack migration [J]. Geophysics, 2000, 65: 426-436.
[17] Krebs J R, Anderson J E, Hinkley D, et al. Fast full-wavefield seismic inversion using encoded sources [J]. Geophysics, 2009, 74(6):177-188.

本刊中的类似文章
1. 张元元, 杨国东, 王凤艳.丹麦法稳健估计在隧道椭圆拟合中的应用[J]. 世界地质, 2012,31(1): 199-203
2. 孟旭, 刘四新, 吴俊军, 傅磊, 邓莉.时间域跨孔雷达全波形反演及实际应用[J]. 世界地质, 2016,35(1): 256-263
3.邢慧婷, 刘财, 李炳秀, 刘洋, 魏忠宇.双曲Radon变换参数分析[J]. 世界地质, 2017,36(2): 579-587
4.张威, 韩立国, 李洪建, 叶林, 张齐.拉东域基于最小二乘反演的斜缆数据鬼波压制方法[J]. 世界地质, 2017,36(2): 595-601
5. 朱恒, 王德利.抛物 Radon 域地震干涉技术[J]. 世界地质, 2015,34(2): 484-490
6. 孙慧秋, 韩立国, 徐杨杨, 高晗, 周岩, 张盼.基于 L-BFGS 算法的全波形反演 Q 值方法研究[J]. 世界地质, 2015,34(4): 1106-1112
7.赫建伟, 周家雄, 王宇, 张盼, 王新领.崎岖海底OBC地震数据全波形反演策略[J]. 世界地质, 2017,36(1): 274-282

Copyright by 世界地质