Full waveform inversion with an auxiliary bump functional

被引:46
作者
Bharadwaj, Pawan [1 ]
Mulder, Wim [1 ,2 ]
Drijkoningen, Guy [1 ]
机构
[1] Delft Univ Technol, Dept Geosci & Engn, Delft, Netherlands
[2] Shell Global Solut Int BV, Rijswijk, Netherlands
关键词
Numerical solutions; Inverse theory; Tomography; Non-linear differential equations; Seismic tomography; GENETIC ALGORITHMS; FREQUENCY-DOMAIN; MODEL; REFLECTION; MIGRATION; TOMOGRAPHY; VELOCITY; FIELD;
D O I
10.1093/gji/ggw129
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
Least-squares inversion of seismic arrivals can provide remarkably detailed models of the Earth's subsurface. However, cycle skipping associated with these oscillatory arrivals is the main cause for local minima in the least-squares objective function. Therefore, it is often difficult for descent methods to converge to the solution without an accurate initial large-scale velocity estimate. The low frequencies in the arrivals, needed to update the large-scale components in the velocity model, are usually unreliable or absent. To overcome this difficulty, we propose a multi-objective inversion scheme that uses the conventional least-squares functional along with an auxiliary data-domain objective. As the auxiliary objective effectively replaces the seismic arrivals by bumps, we call it the bump functional. The bump functional minimization can be made far less sensitive to cycle skipping and can deal with multiple arrivals in the data. However, it can only be used as an auxiliary objective since it usually does not provide a unique model after minimization even when the regularized-least-squares functional has a unique global minimum and hence a unique solution. The role of the bump functional during the multi-objective inversion is to guide the optimization towards the global minimum by pulling the trapped solution out of the local minima associated with the least-squares functional whenever necessary. The computational complexity of the bump functional is equivalent to that of the least-squares functional. In this paper, we describe various characteristics of the bump functional using simple and illustrative numerical examples. We also demonstrate the effectiveness of the proposed multi-objective inversion scheme by considering more realistic examples. These include synthetic and field data from a cross-well experiment, surface-seismic synthetic data with reflections and synthetic data with refracted arrivals at long offsets.
引用
收藏
页码:1076 / 1092
页数:17
相关论文
共 39 条
[1]   Misfit functions for full waveform inversion based on instantaneous phase and envelope measurements [J].
Bozdag, Ebru ;
Trampert, Jeannot ;
Tromp, Jeroen .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2011, 185 (02) :845-870
[2]   Full waveform tomography for lithospheric imaging: results from a blind test in a realistic crustal model [J].
Brenders, A. J. ;
Pratt, R. G. .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2007, 168 (01) :133-151
[3]   MULTISCALE SEISMIC WAVE-FORM INVERSION [J].
BUNKS, C ;
SALECK, FM ;
ZALESKI, S ;
CHAVENT, G .
GEOPHYSICS, 1995, 60 (05) :1457-1473
[4]  
Chauris H., 2012, 74 EAGE C EXH INC SP
[5]   Full waveform inversion method using envelope objective function without low frequency data [J].
Chi, Benxin ;
Dong, Liangguo ;
Liu, Yuzhu .
JOURNAL OF APPLIED GEOPHYSICS, 2014, 109 :36-46
[6]  
Datta D., 2015, 77 EAGE C EXH 2015 E
[7]  
Donno D., 2013, 75 EAGE C EXH INC SP
[8]  
Engquist Bjorn, 2013, ARXIV13114581
[9]  
Fei T.W., 2012, SEG TECHN PROGR EXP
[10]  
Fichtner A, 2011, ADV GEOPHYS ENV MECH, P1, DOI 10.1007/978-3-642-15807-0_1