1887
Volume 67, Issue 8
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

Least‐squares reverse time migration is often formulated as an iterative updating process, where estimating the gradient of the misfit function is necessary. Traditional time domain shot‐profile least‐squares reverse time migration is computationally expensive because computing the gradient involves solving the two‐way wave equation several times in every iteration. To reduce the computational cost of least‐squares reverse time migration, we propose a double‐plane‐wave least‐squares reverse time migration method based on a misfit function for frequency‐domain double‐plane‐wave data. In double‐plane‐wave least‐squares reverse time migration, the gradient is computed by multiplying frequency‐domain plane‐wave Green's functions with the corresponding double‐plane‐wave data residual. Because the number of plane‐wave Green's functions used for migration is relatively small, they can be pre‐computed and stored in a computer's discs or memory. We can use the pre‐computed plane‐wave Green's functions to obtain the gradient without solving the two‐way wave equation in each iteration. Therefore, the migration efficiency is significantly improved. In addition, we study the effects of using sparse frequency sampling and sparse plane‐wave sampling on the proposed method. We can achieve images with correct reflector amplitudes and reasonable resolution using relatively sparse frequency sampling and plane‐wave sampling, which are larger than that determined by the Nyquist theorem. The well‐known wrap‐around artefacts and linear artefacts generated due to under‐sampling frequency and plane wave can be suppressed during iterations in cases where the sampling rates are not excessively large. Moreover, implementing the proposed method with sparse frequency sampling and sparse plane‐wave sampling further improves the computational efficiency. We test the proposed double‐plane‐wave least‐squares reverse time migration on synthetic models to show the practicality of the method.

Loading

Article metrics loading...

/content/journals/10.1111/1365-2478.12803
2019-06-14
2020-07-11
Loading full text...

Full text loading...

References

  1. AkbarF.E., SenM.K. and StoffaP.L.1996. Prestack plane‐wave kirchhoff migration in laterally varying media. Geophysics61, 1068–1079.
    [Google Scholar]
  2. AminzadehF., BracJ. and KunzT.1997. 3‐D salt and overthrust models: SEG/EAGE 3‐D modeling series No.1, SEG.
  3. AokiN. and SchusterG.T.2009. Fast least‐squares migration with a deblurring filter. Geophysics74, WCA83–WCA93.
    [Google Scholar]
  4. ClaerboutJ.F.1985. Imaging the Earth's Interior. Blackwell Scientific Publications, Inc.
    [Google Scholar]
  5. DaiW., FowlerP. and SchusterG.T.2012. Multi‐source least‐squares reverse time migration. Geophysical Prospecting60, 681–695.
    [Google Scholar]
  6. DaiW. and SchusterG.T.2013. Plane‐wave least‐squares reverse‐time migration. Geophysics78, S165–S177.
    [Google Scholar]
  7. DaiW., WangX. and SchusterG.T.2011. Least‐squares migration of multisource data with a deblurring filter. Geophysics76, R135–R146.
    [Google Scholar]
  8. DongS., CaiJ., GuoM., SuhS., ZhangZ., WangB.et al. 2012. Least‐squares reverse time migration: Towards true amplitude imaging and improving the resolution. 83rd Annual International Meeting, SEG, Expanded Abstracts, 1425–1429.
  9. DuanY., SavaP. and GuittonA.2016. Elastic least‐squares reverse time migration. 86th Annual International Meeting, SEG, Expanded Abstracts, 4152–4157.
    [Google Scholar]
  10. DuttaG. and SchusterG.T.2014. Attenuation compensation for least‐squares reverse time migration using the viscoacoustic‐wave equation. Geophysics79, S251–S262.
    [Google Scholar]
  11. EtgenJ.T.2005. How many angles do we really need for delayed‐shot migration? 75th Annual International Meeting, SEG, Expanded Abstracts, 1985–1988.
  12. FengZ. and SchusterG.T.2017. Elastic least‐squares reverse time migration. Geophysics82, S143–S157.
    [Google Scholar]
  13. FokkemaJ.T. and BergP.M.1992. Reflector imaging. Geophysical Journal International110, 191–200.
    [Google Scholar]
  14. GrayS.H.2013. Spatial sampling, migration aliasing, and migrated amplitudes. Geophysics78, S157–S164.
    [Google Scholar]
  15. HerrmannJ.F. and LiX.2012. Efficient least‐squares imaging with sparsity promotion and compressive sensing. Geophysical Prospecting60, 696–712.
    [Google Scholar]
  16. HouJ. and SymesW.W.2016. Accelerating extended least‐squares migration with weighted conjugate gradient iteration. Geophysics81, S165–S179.
    [Google Scholar]
  17. HuangY. and SchusterG.T.2012. Multisource least‐squares migration of marine streamer and land data with frequency‐division encoding. Geophysical Prospecting60, 663–680.
    [Google Scholar]
  18. HustedtB., OpertoS. and VirieuxJ.2004. Mixed‐grid and staggered‐grid finite‐difference methods for frequency‐domain acoustic wave modelling. Geophysical Journal International157, 1269–1296.
    [Google Scholar]
  19. KaplanS.T., RouthP.S. and SacchiM.D.2010. Derivation of forward and adjoint operators for least‐squares shot‐profile split‐step migration. Geophysics75, S225–S235.
    [Google Scholar]
  20. KwonT., SeolS.J. and ByunJ.2015. Efficient full‐waveform inversion with normalized plane‐wave data. Geophysical Journal International201, 53–60.
    [Google Scholar]
  21. LiuF., DouglasW.H., DanN.W., RichardS.D. and RobertH.S.2006. Toward a unified analysis for source plane‐wave migration. Geophysics71, S129–S139.
    [Google Scholar]
  22. LiS., FomelS. and ZhuH.2014. Frequency‐domain reverse‐time migration with a sparse‐frequency sampling. 84th Annual International Meeting, SEG, Expanded Abstracts, 3898–3904.
  23. LiX., AravkinY.A., van LeeuwenT. and HerrmannJ.F.2012. Fast randomized full‐waveform inversion with compressive sensing. Geophysics77, A13–A17.
    [Google Scholar]
  24. LiY., MetivierL., BrossierR., HanB. and VirieuxJ.2015. 2D and 3D frequency‐domain elastic wave modeling in complex media with a parallel iterative solver. Geophysics80, T101–T118.
    [Google Scholar]
  25. LiaoQ. and McMechanG.A.1996. Multifrequency viscoacoustic modeling and inversion. Geophysics61, 1371–1378.
    [Google Scholar]
  26. LuoS. and HaleD.2014. Least‐squares migration in the presence of velocity errors. Geophysics79, S153–S161.
    [Google Scholar]
  27. MarfurtK.J.1984. Accuracy of finite‐difference and finite‐element modeling of the scalar and elastic wave equations. Geophysics49, 533–549.
    [Google Scholar]
  28. MulderW.A. and PlessixR.E.2004. How to choose a subset of frequencies in frequency‐domain finite‐difference migration. Geophysical Journal International158, 801–812.
    [Google Scholar]
  29. NemethT., WuC. and SchusterG.T.1999. Least‐squares migration of incomplete reflection data. Geophysics64, 208–221.
    [Google Scholar]
  30. OpertoS., VirieuxJ., AmestoyP., L'ExcellentJ.‐Y., GiraudL. and AliH.B.H.2007. 3D finite‐difference frequency‐domain modeling of visco‐acoustic wave propagation using a massively parallel direct solver: A feasibility study. Geophysics72, SM195–SM211.
    [Google Scholar]
  31. PestanaC.R. and StoffaL.P.2010. Time evolution of the wave equation using rapid expansion method. Geophysics75, T121–T131.
    [Google Scholar]
  32. PlessixR.E.2006. A review of the adjoint‐state method for computing the gradient of a functional with geophysical applications. Geophysical Journal International167, 495–503.
    [Google Scholar]
  33. PlessixR.E. and MulderW.A.2004. Frequency‐domain finite‐difference amplitude‐preserving migration. Geophysical Journal International157, 975–987.
    [Google Scholar]
  34. PlessixR.E., MulderW.A. and PrattR.2001. Frequency‐domain finite‐difference migration with only few frequencies? 71st Annual International Meeting, SEG, Expanded Abstracts, 1103–1106.
  35. PrattR.G.1999. Seismic waveform inversion in the frequency domain; Part 1, theory and verification in a physical scale model. Geophysics64, 888–901.
    [Google Scholar]
  36. PrattR.G., ShinC. and HickG.J.1998. Gauss–Newton and full newton methods in frequency–space seismic waveform inversion. Geophysical Journal International133, 341–362.
    [Google Scholar]
  37. RenZ., LiuY. and SenM.K.2017. Least‐squares reverse time migration in elastic media. Geophysical Journal International208, 1103–1125.
    [Google Scholar]
  38. RenH., WangH. and ChenS.2013. Least‐squares reverse time migration in frequency domain using the adjoint‐state method. Journal of Geophysics and Engineering10, 035002.
    [Google Scholar]
  39. RiyantiC.D., KononovA., ErlanggaY.A., VuikC., OosterleeC.W., PlessixR.‐E.et al. 2007. A parallel multigrid‐based preconditioner for the 3D heterogeneous high‐frequnecy Helmholtz equation. Journal of Computational Physics224, 431–448.
    [Google Scholar]
  40. SchusterG.T.1993. Least‐squares cross‐well migration. 63rd Annual International Meeting, SEG, Expanded Abstracts, 110–113.
  41. SirgueL. and PrattR.G.2004. Efficient waveform inversion and imaging: A strategy for selecting temporal frequencies. Geophysics69, 231–248.
    [Google Scholar]
  42. StantonA. and SacchiM.D.2015. Least squares wave equation migration of elastic data. 77th EAGE Conference and Exhibition 2015.
  43. StantonA. and SacchiM.D.2017. Elastic least‐squares one‐way wave‐equation migration. Geophysics82, S293–S305.
    [Google Scholar]
  44. ŠteklI. and PrattR.G.1998, Accurate viscoelastic modeling by frequency‐domain finite differences using rotated operators. Geophysics63, 1779–1794.
    [Google Scholar]
  45. StoffaP.L., SenM.K., SeifoullaevR.K., PestanaR.C. and FokkemaJ.T.2006. Plane‐wave depth migration. Geophysics71, S261–S272.
    [Google Scholar]
  46. StorkC. and KapoorJ.2004. How many P values do you want to migrate for delayed shot wave equation migration? 74th Annual International Meeting, SEG, Expanded Abstracts, 1041–1044.
    [Google Scholar]
  47. SunJ., FomelS. and ZhuT.2015. Preconditioning least‐squares RTM in viscoacoustic media by Q‐compensated RTM. 85th Annual International Meeting, SEG, Expanded, 3959–3965.
  48. TangY.2009. Target‐oriented wave‐equation least‐squares migration/inversion with phase‐encoded Hessian. Geophysics74, WCA95–WCA107.
    [Google Scholar]
  49. TangY. and BiondiB.2009. Least‐squares migration/inversion of blended data. 79th Annual International Meeting, SEG, Expanded Abstracts, 2859–2863.
  50. TarantolaA.1984. Inversion of seismic reflection data in the acoustic approximation. Geophysics49, 1259–1266.
    [Google Scholar]
  51. TatalovicR., DillenM.W.P. and FokkemaJ.T.1991. Prestack imaging in the double transformed radon domain. 61st Annual International Meeting, SEG, Expanded Abstracts, 1285–1288.
  52. TrompJ., TapeC. and LiuQ.2004. Seismic tomography, adjoint methods, time reversal and banana‐doughnut kernels. Geophysical Journal International160, 195–216.
    [Google Scholar]
  53. ValencianoA.A., BiondiB. and GuittonA.2006. Target‐oriented wave‐equation inversion. Geophysics71, A35–A38.
    [Google Scholar]
  54. ValencianoA.A., CheminguiN., WhitmoreD. and Brandsberg‐DahlS.2011. Wave equation migration with attenuation and anisotropy compensation. 81st Annual International Meeting, SEG, Expanded Abstracts, 232–236.
  55. WangJ. and SacchiM.D.2007. High‐resolution wave‐equation amplitude‐variation‐with‐ray‐parameter (AVP) imaging with sparseness constraints. Geophysics72, S11–S18.
    [Google Scholar]
  56. XueZ., ChenY., FomelS. and SunJ.2016. Seismic imaging of incomplete data and simultaneous‐source data using least‐squares reverse time migration with shaping regularization. Geophysics81, S11–S20.
    [Google Scholar]
  57. ZhangY., DuanL. and XieY.2015. A stable and practical implementation of least‐squares reverse time migration. Geophysics80, V23–V31.
    [Google Scholar]
  58. ZhangY., SunJ., NotforsC., GrayS.H., ChernisL. and YoungJ.2005. Delayed‐shot 3D depth migration. Geophysics70, E21–E28.
    [Google Scholar]
  59. ZhaoZ., SenM.K. and StoffaP.L.2014. Double plane wave reverse time migration in the frequency domain. 84th Annual International Meeting, SEG, Expanded Abstracts, 1033–1039.
  60. ZhaoZ., SenM.K. and StoffaP.L.2016. Double‐plane‐wave reverse time migration in the frequency domain. Geophysics81, S367–S382.
    [Google Scholar]
  61. ZhaoZ., SenM.K. and StoffaP.L.2017a. Double plane‐wave reverse‐time migration. Geophysical Prospecting65, 1541–1558.
    [Google Scholar]
  62. ZhaoZ., SenM.K. and StoffaP.L.2017b. Reciprocity and double plane‐wave migration. Geophysics82, S453–S466.
    [Google Scholar]
  63. ZhaoZ., SenM.K., StoffaP.L. and ZhuH.2015. Double plane wave least squares reverse time migration. 85th Annual International Meeting, SEG, Expanded Abstracts, 4170–4174.
http://instance.metastore.ingenta.com/content/journals/10.1111/1365-2478.12803
Loading
/content/journals/10.1111/1365-2478.12803
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): Computing aspects , Imaging , Inversion and Numerical study
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