解飞^1,2,,
黄建平^1,2,,,
李振春^1,2
1. 中国石油大学(华东)地球科学与技术学院地球物理系, 青岛 266580
2. 海洋国家实验室海洋矿产资源评价与探测技术功能实验室, 青岛 266071

基金项目: 泰山****工程专项经费，国家973子课题（2014CB239006），国家油气重大专项（2016ZX05014-001-008HZ，2016ZX05002-005-07HZ），中央科研经费（17CX05011）及青岛市自主创新（16-5-1-40-jch）联合资助


详细信息

作者简介: 解飞, 男, 1993年生, 硕士研究生, 主要从事频率域全波形反演工作.E-mail:xiefeilbq@foxmail.com

通讯作者: 黄建平, 1982年生, 博士, 教授, 主要从事地震波正演及偏移成像工作.E-mail:jphuang@upc.edu.cn

中图分类号: P631

收稿日期:2017-04-13
修回日期:2017-05-23
上线日期:2018-08-05

Full waveform inversion in the frequency domain with an efficient iterative solver

XIE Fei^1,2,,
HUANG JianPing^1,2,,,
LI ZhenChun^1,2
1. School of Geosciences, China University of Petroleum(East China), Qingdao 266580, China
2. Laboratory of Marine Mineral Resources Evaluation and Exploration Technology, Qingdao 266071, China


More Information

Corresponding author: HUANG JianPing,E-mail:jphuang@upc.edu.cn

MSC: P631

--> Received Date: 13 April 2017
Revised Date: 23 May 2017
Available Online: 05 August 2018


摘要
摘要:巨大的计算量是制约全波形反演（FWI）生产实用化的难题之一.为此，本文提出了一种高效的波场迭代解法，将其应用于频率域常密度声波方程FWI，并给出了详细的反演流程.通过建立用于波场迭代的目标函数，推导相应梯度、步长公式，新方法将反演中波场正传和残差波场反传过程转化为无约束优化问题，从理论上分析了新方法的计算效率显著高于常规FWI.在数值试验中，本文方法通过几次迭代便能获得高精度的正传、残差反传波场，收敛速度明显高于未经预处理的GMRES方法.进一步引入高效编码策略，新方法的计算时间约为常规编码FWI的1/8，与理论分析结果吻合（波场迭代次数为8，模型未知量个数约为7万），且波场迭代次数为6时，反演效果已与常规编码FWI相近.
关键词: 全波形反演/
计算量/
迭代解法/
编码

Abstract:The huge computational complexity is one of the major challenges in the application of full waveform inversion (FWI) technology. The most computationally expensive part of FWI is a forward modeling engine that numerically simulates wavefields. In this paper, an efficient iterative solver is proposed for wavefield continuation and applied to frequency-domain acoustic full-waveform inversion with constant density. We reformulate the forward propagation of the source wavefield and the back propagation of the residual wavefield as unconstrained optimization problems. We also derive calculating formula for gradient and step-length. Theoreticall analysis indicates that the computational efficiency of the new method is significantly higher than that of the conventional FWI. Numerical examples prove that our method can obtain high-precision wavefields with several iterations and improve convergence rate compared with original GMRES method. With high-efficiency source encoding technology, the calculation time of the new method is only about 1/8 of that of the conventional FWI, which is consistent with the theoretical analysis (The number of wavefield iteration is 8, and the number of unknown parameters is about 70 thousands). And when the number of wavefield iteration is 6, the inversion results of the new method is basically the same as that of the conventional FWI.
Key words:Full waveform inversion/
Computation/
Iterative solver/
Encoding



PDF全文下载地址:

http://www.geophy.cn/data/article/export-pdf?id=dqwlxb_14636