skip to main content
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Interparameter trade-off quantification and reduction in isotropic-elastic full-waveform inversion: synthetic experiments and Hussar land data set application

Abstract

The problem of inverting for multiple physical parameters in the subsurface using seismic full-waveform inversion (FWI) is complicated by interparameter trade-off arising from inherent ambiguities between different physical parameters. Parameter resolution is often characterized using scattering radiation patterns, but these neglect some important aspects of interparameter trade-off. More general analysis and mitigation of interparameter trade-off in isotropic-elastic FWI is possible through judiciously chosen multiparameter Hessian matrix-vector products. We show that products of multiparameter Hessian off-diagonal blocks with model perturbation vectors, referred to as interparameter contamination kernels, are central to the approach. We then apply the multiparameter Hessian to various vectors designed to provide information regarding the strengths and characteristics of interparameter contamination, both locally and within the whole volume. With numerical experiments, we observe that S-wave velocity perturbations introduce strong contaminations into density and phase-reversed contaminations into P-wave velocity, but themselves experience only limited contaminations from other parameters. Based on these findings, we introduce a novel strategy to mitigate the influence of interparameter trade-off with approximate contamination kernels. Furthermore, we recommend that the local spatial and interparameter trade-off of the inverted models be quantified using extended multiparameter point spread functions (EMPSFs) obtained with pre-conditioned conjugate-gradient algorithm. Compared to traditional pointmore » spread functions, the EMPSFs appear to provide more accurate measurements for resolution analysis, by de-blurring the estimations, scaling magnitudes and mitigating interparameter contamination. Approximate eigenvalue volumes constructed with stochastic probing approach are proposed to evaluate the resolution of the inverted models within the whole model. With a synthetic Marmousi model example and a land seismic field data set from Hussar, Alberta, Canada, we confirm that the new inversion strategy suppresses the interparameter contamination effectively and provides more reliable density estimations in isotropic-elastic FWI as compared to standard simultaneous inversion approach.« less

Authors:
ORCiD logo [1];  [2];  [2]
  1. Univ. of Calgary, AB (Canada). Dept. of Geoscience; Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
  2. Univ. of Calgary, AB (Canada). Dept. of Geoscience
Publication Date:
Research Org.:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE; Natural Sciences and Engineering Research Council of Canada (NSERC)
Contributing Org.:
Princeton Univ., NJ (United States). TIGER cluster; Compute Canada, Toronto (Canada). Lattice cluster
OSTI Identifier:
1473814
Report Number(s):
LA-UR-18-20731
Journal ID: ISSN 0956-540X
Grant/Contract Number:  
AC52-06NA25396; CRDPJ 461179-13
Resource Type:
Accepted Manuscript
Journal Name:
Geophysical Journal International
Additional Journal Information:
Journal Volume: 213; Journal Issue: 2; Journal ID: ISSN 0956-540X
Publisher:
Oxford University Press
Country of Publication:
United States
Language:
English
Subject:
58 GEOSCIENCES; 97 MATHEMATICS AND COMPUTING; Inverse theory; Numerical modelling; Waveform inversion

Citation Formats

Pan, Wenyong, Innanen, Kristopher A., and Geng, Yu. Interparameter trade-off quantification and reduction in isotropic-elastic full-waveform inversion: synthetic experiments and Hussar land data set application. United States: N. p., 2018. Web. doi:10.1093/gji/ggy037.
Pan, Wenyong, Innanen, Kristopher A., & Geng, Yu. Interparameter trade-off quantification and reduction in isotropic-elastic full-waveform inversion: synthetic experiments and Hussar land data set application. United States. doi:10.1093/gji/ggy037.
Pan, Wenyong, Innanen, Kristopher A., and Geng, Yu. Tue . "Interparameter trade-off quantification and reduction in isotropic-elastic full-waveform inversion: synthetic experiments and Hussar land data set application". United States. doi:10.1093/gji/ggy037. https://www.osti.gov/servlets/purl/1473814.
@article{osti_1473814,
title = {Interparameter trade-off quantification and reduction in isotropic-elastic full-waveform inversion: synthetic experiments and Hussar land data set application},
author = {Pan, Wenyong and Innanen, Kristopher A. and Geng, Yu},
abstractNote = {The problem of inverting for multiple physical parameters in the subsurface using seismic full-waveform inversion (FWI) is complicated by interparameter trade-off arising from inherent ambiguities between different physical parameters. Parameter resolution is often characterized using scattering radiation patterns, but these neglect some important aspects of interparameter trade-off. More general analysis and mitigation of interparameter trade-off in isotropic-elastic FWI is possible through judiciously chosen multiparameter Hessian matrix-vector products. We show that products of multiparameter Hessian off-diagonal blocks with model perturbation vectors, referred to as interparameter contamination kernels, are central to the approach. We then apply the multiparameter Hessian to various vectors designed to provide information regarding the strengths and characteristics of interparameter contamination, both locally and within the whole volume. With numerical experiments, we observe that S-wave velocity perturbations introduce strong contaminations into density and phase-reversed contaminations into P-wave velocity, but themselves experience only limited contaminations from other parameters. Based on these findings, we introduce a novel strategy to mitigate the influence of interparameter trade-off with approximate contamination kernels. Furthermore, we recommend that the local spatial and interparameter trade-off of the inverted models be quantified using extended multiparameter point spread functions (EMPSFs) obtained with pre-conditioned conjugate-gradient algorithm. Compared to traditional point spread functions, the EMPSFs appear to provide more accurate measurements for resolution analysis, by de-blurring the estimations, scaling magnitudes and mitigating interparameter contamination. Approximate eigenvalue volumes constructed with stochastic probing approach are proposed to evaluate the resolution of the inverted models within the whole model. With a synthetic Marmousi model example and a land seismic field data set from Hussar, Alberta, Canada, we confirm that the new inversion strategy suppresses the interparameter contamination effectively and provides more reliable density estimations in isotropic-elastic FWI as compared to standard simultaneous inversion approach.},
doi = {10.1093/gji/ggy037},
journal = {Geophysical Journal International},
number = 2,
volume = 213,
place = {United States},
year = {2018},
month = {5}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 9 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Resolution analysis by random probing: RESOLUTION ANALYSIS BY RANDOM PROBING
journal, August 2015

  • Fichtner, Andreas; Leeuwen, Tristan van
  • Journal of Geophysical Research: Solid Earth, Vol. 120, Issue 8
  • DOI: 10.1002/2015JB012106

Efficient stochastic estimation of the model resolution matrix diagonal and generalized cross–validation for large geophysical inverse problems
journal, January 2011

  • MacCarthy, J. K.; Borchers, B.; Aster, R. C.
  • Journal of Geophysical Research, Vol. 116, Issue B10
  • DOI: 10.1029/2011JB008234

Elastic full waveform inversion based on mode decomposition: the approach and mechanism
journal, February 2017

  • Wang, T. F.; Cheng, J. B.
  • Geophysical Journal International, Vol. 209, Issue 2
  • DOI: 10.1093/gji/ggx038

The adjoint method in seismology
journal, August 2006

  • Fichtner, A.; Bunge, H. -P.; Igel, H.
  • Physics of the Earth and Planetary Interiors, Vol. 157, Issue 1-2
  • DOI: 10.1016/j.pepi.2006.03.016

Optimal observables for multiparameter seismic tomography
journal, June 2014

  • Bernauer, Moritz; Fichtner, Andreas; Igel, Heiner
  • Geophysical Journal International, Vol. 198, Issue 2
  • DOI: 10.1093/gji/ggu204

The adjoint method in seismology—
journal, August 2006

  • Fichtner, A.; Bunge, H. -P.; Igel, H.
  • Physics of the Earth and Planetary Interiors, Vol. 157, Issue 1-2
  • DOI: 10.1016/j.pepi.2006.03.018

Anisotropic full-waveform inversion with tilt-angle recovery
journal, May 2017

  • Rusmanugroho, Herurisa; Modrak, Ryan; Tromp, Jeroen
  • GEOPHYSICS, Vol. 82, Issue 3
  • DOI: 10.1190/geo2016-0025.1

The Resolving Power of Gross Earth Data
journal, October 1968


A Semismooth Newton-CG Method for Constrained Parameter Identification in Seismic Tomography
journal, January 2015

  • Boehm, Christian; Ulbrich, Michael
  • SIAM Journal on Scientific Computing, Vol. 37, Issue 5
  • DOI: 10.1137/140968331

Subspace methods for large inverse problems with multiple parameter classes
journal, August 1988


A stochastic estimator of the trace of the influence matrix for laplacian smoothing splines
journal, January 1990

  • Hutchinson, M. F.
  • Communications in Statistics - Simulation and Computation, Vol. 19, Issue 2
  • DOI: 10.1080/03610919008812866

Wave equation least square imaging using the local angular Hessian for amplitude correction: Wave equation least square imaging
journal, March 2011


Seismic structure of the European upper mantle based on adjoint tomography
journal, February 2015

  • Zhu, Hejun; Bozdağ, Ebru; Tromp, Jeroen
  • Geophysical Journal International, Vol. 201, Issue 1
  • DOI: 10.1093/gji/ggu492

3D coupled acoustic-elastic migration with topography and bathymetry based on spectral-element and adjoint methods
journal, July 2013


Double-difference adjoint seismic tomography
journal, June 2016

  • Yuan, Yanhua O.; Simons, Frederik J.; Tromp, Jeroen
  • Geophysical Journal International, Vol. 206, Issue 3
  • DOI: 10.1093/gji/ggw233

Seismic AVO and the inverse Hessian in precritical reflection full waveform inversion
journal, August 2014

  • Innanen, Kristopher A.
  • Geophysical Journal International, Vol. 199, Issue 2
  • DOI: 10.1093/gji/ggu291

A guided tour of multiparameter full-waveform inversion with multicomponent data: From theory to practice
journal, September 2013


The point-spread function measure of resolution for the 3-D electrical resistivity experiment
journal, February 2009


Resolution tests revisited: the power of random numbers
journal, December 2012

  • Trampert, Jeannot; Fichtner, Andreas; Ritsema, Jeroen
  • Geophysical Journal International, Vol. 192, Issue 2
  • DOI: 10.1093/gji/ggs057

Accelerating Hessian-free Gauss-Newton full-waveform inversion via l-BFGS preconditioned conjugate-gradient algorithm
journal, March 2017

  • Pan, Wenyong; Innanen, Kristopher A.; Liao, Wenyuan
  • GEOPHYSICS, Vol. 82, Issue 2
  • DOI: 10.1190/geo2015-0595.1

Multiscale adjoint waveform-difference tomography using wavelets
journal, May 2014


A recipe for practical full-waveform inversion in anisotropic media: An analytical parameter resolution study
journal, May 2014


Elastic imaging and time-lapse migration based on adjoint methods
journal, November 2009

  • Zhu, Hejun; Luo, Yang; Nissen-Meyer, Tarje
  • GEOPHYSICS, Vol. 74, Issue 6
  • DOI: 10.1190/1.3261747

Inversion strategies for visco-acoustic waveform inversion
journal, June 2013

  • Kamei, R.; Pratt, R. G.
  • Geophysical Journal International, Vol. 194, Issue 2
  • DOI: 10.1093/gji/ggt109

Full waveform inversion and the truncated Newton method: quantitative imaging of complex subsurface structures: FWI and the truncated Newton method
journal, June 2014

  • Métivier, L.; Bretaudeau, F.; Brossier, R.
  • Geophysical Prospecting, Vol. 62, Issue 6
  • DOI: 10.1111/1365-2478.12136

Nonlinear two‐dimensional elastic inversion of multioffset seismic data
journal, September 1987


Target-oriented wave-equation least-squares migration/inversion with phase-encoded Hessian
journal, November 2009


Hessian kernels of seismic data functionals based upon adjoint techniques: Hessian kernels
journal, February 2011


Frequency-domain full-waveform inversion with non-linear descent directions
journal, January 2018

  • Geng, Yu; Pan, Wenyong; Innanen, Kristopher A.
  • Geophysical Journal International, Vol. 213, Issue 2
  • DOI: 10.1093/gji/ggy002

An improved scattering-integral approach for frequency-domain full waveform inversion
journal, July 2015

  • Liu, Yuzhu; Yang, Jizhong; Chi, Benxin
  • Geophysical Journal International, Vol. 202, Issue 3
  • DOI: 10.1093/gji/ggv254

Elastic full-waveform inversion and parametrization analysis applied to walk-away vertical seismic profile data for unconventional (heavy oil) reservoir characterization
journal, March 2018

  • Pan, Wenyong; Innanen, Kristopher A.; Geng, Yu
  • Geophysical Journal International, Vol. 213, Issue 3
  • DOI: 10.1093/gji/ggy087

Synthetic inversions for density using seismic and gravity data
journal, March 2017

  • Blom, Nienke; Boehm, Christian; Fichtner, Andreas
  • Geophysical Journal International, Vol. 209, Issue 2
  • DOI: 10.1093/gji/ggx076

Seismic imaging of complex onshore structures by 2D elastic frequency-domain full-waveform inversion
journal, November 2009

  • Brossier, Romain; Operto, Stéphane; Virieux, Jean
  • GEOPHYSICS, Vol. 74, Issue 6
  • DOI: 10.1190/1.3215771

On the influence of model parametrization in elastic full waveform tomography: On the influence of model parametrization
journal, August 2012


Resolution analysis in full waveform inversion: Resolution in full waveform inversion
journal, October 2011


Full Waveform Inversion and the Truncated Newton Method
journal, January 2013

  • Métivier, L.; Brossier, R.; Virieux, J.
  • SIAM Journal on Scientific Computing, Vol. 35, Issue 2
  • DOI: 10.1137/120877854

Seismic waveform inversion best practices: regional, global and exploration test cases
journal, June 2016

  • Modrak, Ryan; Tromp, Jeroen
  • Geophysical Journal International, Vol. 206, Issue 3
  • DOI: 10.1093/gji/ggw202

The Iceland–Jan Mayen plume system and its impact on mantle dynamics in the North Atlantic region: Evidence from full-waveform inversion
journal, April 2013

  • Rickers, Florian; Fichtner, Andreas; Trampert, Jeannot
  • Earth and Planetary Science Letters, Vol. 367
  • DOI: 10.1016/j.epsl.2013.02.022

Inversion of seismic reflection data in the acoustic approximation
journal, August 1984


Multiscale adjoint waveform tomography for surface and body waves
journal, September 2015

  • Yuan, Yanhua O.; Simons, Frederik J.; Bozdağ, Ebru
  • GEOPHYSICS, Vol. 80, Issue 5
  • DOI: 10.1190/geo2014-0461.1

Introduction to the spectral element method for three-dimensional seismic wave propagation
journal, December 1999


Improved amplitude preservation for prestack depth migration by inverse scattering theory
journal, September 2001


A review of the adjoint-state method for computing the gradient of a functional with geophysical applications
journal, November 2006


Target-oriented wave-equation inversion
journal, July 2006

  • Valenciano, Alejandro A.; Biondi, Biondo; Guitton, Antoine
  • GEOPHYSICS, Vol. 71, Issue 4
  • DOI: 10.1190/1.2213359

Scattering of homogeneous and inhomogeneous seismic waves in low-loss viscoelastic media
journal, July 2015

  • Moradi, Shahpoor; Innanen, Kristopher A.
  • Geophysical Journal International, Vol. 202, Issue 3
  • DOI: 10.1093/gji/ggv242

An overview of full-waveform inversion in exploration geophysics
journal, November 2009


Estimation of elastic constants for HTI media using Gauss-Newton and full-Newton multiparameter full-waveform inversion
journal, September 2016

  • Pan, Wenyong; Innanen, Kristopher A.; Margrave, Gary F.
  • GEOPHYSICS, Vol. 81, Issue 5
  • DOI: 10.1190/geo2015-0594.1

On the resolution of density anomalies in the Earth's mantle using spectral fitting of normal-mode data
journal, July 2002


Building good starting models for full-waveform inversion using adaptive matching filtering misfit
journal, September 2016


Coupling in amplitude variation with offset and the Wiggins approximation
journal, July 2013


Global adjoint tomography: first-generation model
journal, September 2016

  • Bozdağ, Ebru; Peter, Daniel; Lefebvre, Matthieu
  • Geophysical Journal International, Vol. 207, Issue 3
  • DOI: 10.1093/gji/ggw356

2D frequency-domain elastic full-waveform inversion using the block-diagonal pseudo-Hessian approximation
journal, September 2016


The scattering potential of partial derivative wavefields in 3-D elastic orthorhombic media: an inversion prospective
journal, July 2016

  • Oh, Ju-Won; Alkhalifah, Tariq
  • Geophysical Journal International, Vol. 206, Issue 3
  • DOI: 10.1093/gji/ggw238

Misfit functions for full waveform inversion based on instantaneous phase and envelope measurements: Misfit functions for full waveform inversion
journal, March 2011


Delay-time tomography of the upper mantle below Europe, the Mediterranean, and Asia Minor
journal, November 1991


Poststack migration deconvolution
journal, May 2001

  • Hu, Jianxing; Schuster, Gerard T.; Valasek, Paul A.
  • GEOPHYSICS, Vol. 66, Issue 3
  • DOI: 10.1190/1.1444984

Efficient pseudo-Gauss-Newton full-waveform inversion in the τ-p domain
journal, September 2015

  • Pan, Wenyong; Innanen, Kristopher A.; Margrave, Gary F.
  • GEOPHYSICS, Vol. 80, Issue 5
  • DOI: 10.1190/geo2014-0224.1

A strategy for nonlinear elastic inversion of seismic reflection data
journal, October 1986


The imprint of crustal density heterogeneities on regional seismic wave propagation
journal, January 2016

  • Płonka, Agnieszka; Blom, Nienke; Fichtner, Andreas
  • Solid Earth, Vol. 7, Issue 6
  • DOI: 10.5194/se-7-1591-2016

Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion
journal, May 1998


A Bayesian approach to estimate uncertainty for full-waveform inversion using a priori information from depth migration
journal, September 2016


On the use of sensitivity tests in seismic tomography
journal, March 2016

  • Rawlinson, N.; Spakman, W.
  • Geophysical Journal International, Vol. 205, Issue 2
  • DOI: 10.1093/gji/ggw084