1887
image of Acoustic full‐waveform inversion of surface seismic data using the Gauss‐Newton method with active constraint balancing

Abstract

ABSTRACT

We propose a full‐waveform inversion algorithm using the Gauss‐Newton inversion method with active constraint balancing that uses the spatially variant damping factor and source normalized wavefield approach for surface seismic data in the frequency domain. The active constraint balancing technique automatically determines the optimum distribution of damping factors that control the stability and resolution in Gauss‐Newton inversion by using a parameter resolution matrix and spread function analysis. Through numerical experiments, we present that the active constraint balancing scheme provides stable inversion results without a severe loss of resolution compared with the conventional Gauss‐Newton method. The reconstructed image using the active constraint balancing method more closely resembles the true image for the region with low sensitivity. Also, the estimated value converges faster to the smaller RMS error level than those estimated by the conventional Gauss‐Newton method. We also implement the normalized wavefield method to overcome the lack of precise knowledge on the source. The source normalized wavefield approach effectively removes the potential inversion errors from source estimation because the source spectrum is eliminated during the normalization procedure. Our inversion algorithm, using the source normalization scheme, provides excellent inversion results even though the data are generated by two slightly different source wavelets. We present that the frequency selection scheme proposed by Sirgue and Pratt, which is based on the average amplitude of the whole received data, provides a useful guideline for selecting the proper frequencies for our inversion. Our novel inversion algorithm successfully reconstructs the velocity model within 10–30 iterations despite its starting from a homogeneous or linearly increasing velocity model. In addition, for testing the performance of our inversion algorithm on a more complicated structure, we apply the algorithm to the SEG/EAGE overthrust model. Successful inversion is achieved as the reconstructed image approaches the true model with the consistently converging RMS error even though insufficient data are used.

Loading

Article metrics loading...

/content/journals/10.1111/j.1365-2478.2012.01112.x
2012-10-18
2020-03-29
Loading full text...

Full text loading...

References

  1. AkiK. and RichardsP.G.1980. Quantitative Seismology: Theory and Methods , Vol. 1 . W.H. Freeman and Company, San Francisco, California .
    [Google Scholar]
  2. BérengerJ.P.1994. A perfectly matched layer for the absorption of electromagnetic waves. Journal of Computational Physics 114, 185–200.
    [Google Scholar]
  3. BrossierR., OpertoS. and VirieuxJ.2009. Seismic imaging of complex onshore structures by 2D elastic frequency‐domain full‐waveform inversion. Geophysics 74, 105–118.
    [Google Scholar]
  4. ChewW.C. and WeedonW.H.1994. A 3D perfectly matched medium from modified Maxwell's equations with stretched coordinates. Microwave and Optical Technology Letters 43, 599–604.
    [Google Scholar]
  5. DavisT.A. and DuffI.S.1997. An unsymmetric‐pattern multifrontal method for sparse LU factorization. SIAM Journal on Matrix Analysis and Applications , 18(1), 140–158.
    [Google Scholar]
  6. ForguesE., ScalaE. and PrattR.G.1998. High resolution velocity model estimation from refraction and reflection data. 68th Annual International Meeting, SEG, ExpandedAbstracts, 1211–1214.
  7. GauthierO., VirieuxJ. and TarantolaA.1986. Two‐dimensional non‐linear inversion of seismic waveforms: Numerical results. Geophysics 51, 1387–1403
    [Google Scholar]
  8. HuW., AbubakarA. and HabashyT.M.2009. Simultaneous multifrequency inversion of full‐waveform seismic data. Geophysics 74, P.R1–R14.
    [Google Scholar]
  9. JacksonD.D.1972. Interpretation of inaccurate, insufficient, and inconsistent data. Geophysical Journal of the Royal Astronomical Society 28, 97–109
    [Google Scholar]
  10. JuppD.L. and VozoffK.1975. Stable iterative methods for the inversion of geophysical data. Geophysical Journal of the Royal Astronomical Society 42, 67–72.
    [Google Scholar]
  11. KimH.J., SongY. and LeeK.H.1999. Inequality constraint in least‐squares inversion of geophysical data. Earth Planets Space 51, 255–259.
    [Google Scholar]
  12. KolbP., CollinoF. and LaillyP.1986. Pre‐stack inversion of a 1‐D medium. Proceedings of the IEEE 74
    [Google Scholar]
  13. LaillyP.1983. The seismic inverse problem as a sequence of before stack migrations. In: Conference on Inverse Scattering: Theory and Application , (eds J.B.Bednar , R.Redner , E.Robinson and  Weglein ). Society for Industrial and Applied Mathematics, Philadelphia , PA .
    [Google Scholar]
  14. LeeL.H.2003. 3‐D full waveform inversion of seismic data; Part I. Theory. Lawrence Berkeley National Laboratory, LBNL Paper, LBNL‐52665.
  15. LeeL.H. and KimH.J.2003. Source‐independent full‐waveform inversion of seismic data. Geophysics 68, 2010–2015.
    [Google Scholar]
  16. LevenbergK.1944. A Method for the Solution of Certain Non‐Linear Problems in Least Squares. The Quarterly of Applied Mathematics 2, 164–168.
    [Google Scholar]
  17. LiaoQ. and McMechanG.A.1996. Multifrequency viscoacoustic modeling and inversion. Geophysics 61, 1371–1378.
    [Google Scholar]
  18. LinesL.R. and TreitelS.1984. Tutorial: A review of least‐squares inversion and its application to geophysical problems. Geophysical Prospecting 32, 159–186.
    [Google Scholar]
  19. MarfurtK.J.1984. Accuracy of finite‐difference and finite‐element modeling of the scalar and elastic wave‐equations. Geophysics 49, 533–549
    [Google Scholar]
  20. MarquardtD.W.1963. An algorithm for least‐squares estimation of nonlinear parameters. Journal of the Society of Industrial and Applied Mathematics 11, 431–441.
    [Google Scholar]
  21. MenkeW.1984. Geophysical Data Analysis: Discrete Inverse Theory . Academic Press Inc.
    [Google Scholar]
  22. MoraP.R.1987. Nonlinear two‐dimensional elastic inversion of multi‐offset seismic data. Geophysics 52, 1211–1228.
    [Google Scholar]
  23. NemethT., NormakE. and QinF.1997. Dynamic smoothing in crosswell traveltime tomography. Geophysics 62, 168–176.
    [Google Scholar]
  24. OristaglioM.L. and WorthingtonM.H.1980. Inversion of surface and borehole electromagnetic data for two dimensional electric conductivity models. Geophysical Prospecting 28, 633–657.
    [Google Scholar]
  25. PrattR.G.1990. Inverse theory applied to multi‐source cross‐hole tomography. Part II : Elastic wave‐equation method. Geophysical Prospecting 38, 311–330.
    [Google Scholar]
  26. PrattR.G.1999. Seismic waveform inversion in the frequency domain, Part 1: Theory and verification in a physical scale model. Geophysics 64, 888–901.
    [Google Scholar]
  27. PrattR.G., ShinC.S. and HicksG.J.1998. Gauss‐Newton and full Newton methods in frequency‐space seismic waveform inversion. Geophysical Journal International 133, 341–362.
    [Google Scholar]
  28. PrattR.G. and WorthingtonM.H.1990. Inverse theory applied to multi‐source cross‐hole tomography. Part I : Acoustic wave‐equation method. Geophysical Prospecting 38, 287–310.
    [Google Scholar]
  29. RaoY. and WangY.2009. Fracture effects in seismic attenuation images reconstructed by waveform tomography. Geophysics 74, R25–R34
    [Google Scholar]
  30. SenM.K. and RoyI.G.2003. Computation of differential seismograms and iteration adaptive regularization in prestack waveform inversion. Geophysics 68, 2026–2039.
    [Google Scholar]
  31. SheenD.H., TuncayK., BaagC.E. and OrtolevaP.2006. Time domain Gauss‐Newton seismic waveform inversion in elastic media. Geophysical Journal International 167, 1373–1384.
    [Google Scholar]
  32. ShinC.1988. Nonlinear elastic wave inversion by blocky parameterization . PhD dissertation, The University of Tulsa.
    [Google Scholar]
  33. ShinC., KooN.H., ChaY.H. and ParkK.P.2010. Sequentially ordered single‐frequency 2‐D acoustic waveform inversion in the Laplace‐Fourier domain. Geophysical Journal International 181, 935–950.
    [Google Scholar]
  34. ShinC. and MinD.J.2006. Waveform inversion using a logarithmic wavefield. Geophysics 71, no. 3, R31–R42.
    [Google Scholar]
  35. SirgueL. and PrattR.G.2004. Efficient waveform inversion and imaging: A strategy for selecting temporal frequencies. Geophysics 69, 231–248.
    [Google Scholar]
  36. TarantolaA.1984. Inversion of seismic reflection data in the acoustic approximation. Geophysics 49, 1259–1266.
    [Google Scholar]
  37. TarantolaA.1987. Inverse Problem Theory: Methods for Data Fitting and Parameter Estimations . Elsevier Science Publishing Co., New York .
    [Google Scholar]
  38. WoodhouseJ.H. and DziewonskiA.M.1984. Mapping the upper mantle: Three‐dimensional modeling of earth structure by inversion of seismic waveforms. Journal of Geophysical Research. 89, 5953–5986.
    [Google Scholar]
  39. YiM.J., KimJ.H. and ChungS.H.2003. Enhancing the resolving power of least‐squares inversion with active constraint balancing. Geophysics 68, 931–941.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/j.1365-2478.2012.01112.x
Loading
/content/journals/10.1111/j.1365-2478.2012.01112.x
Loading

Data & Media loading...

  • Article Type: Research Article
Keywords: Inversion; Full waveform; Seismics
This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error