1887
Volume 17, Issue 2
  • ISSN: 1569-4445
  • E-ISSN: 1873-0604

Abstract

ABSTRACT

We implement two different inversion strategies for solving the full‐waveform inversion of spectral induced polarization data: the first is a local inversion based on the Levenberg–Marquardt algorithm, whereas the second is a global optimization that makes use of the Particle Swarm Optimization algorithm. In addition, the analysis of the residual function maps and the singular value decomposition of the inversion kernel are used to analyse the complexity of the objective function and the ill‐posedness of this inverse problem. We limit the attention to synthetic data with the aim to maintain the inversion at a simple level and to draw essential conclusions about the ill‐conditioning of the considered inverse problem and about the suitability of the two inversion approaches we use. We consider two differexnt double‐dispersion reference models that generate resistivity amplitude and phase spectra with different characteristics. Realistic noise is added to the synthetic data to better simulate a field data set. We also apply a dedicated processing sequence to increase the signal‐to‐noise ratio of the observed data. It turns out that the full‐waveform inversion of spectral induced polarization is a well‐posed problem in case of double‐peaked resistivity spectra, whereas it becomes hopelessly ill‐conditioned when the subsurface model generates single‐peaked spectra. In particular, the analysis of the residual function maps demonstrates that in case of double‐peaked spectra the objective function is characterized by a well‐defined single minimum. Conversely, elongated valleys with similar misfit values arise for single‐peaked spectra. In this case, the eigenvectors in model space demonstrate that the estimate of the and parameters is a non‐unique problem and for this reason it would be advisable to reparametrize the inverse problem by considering a single‐dispersion model. In addition, our tests demonstrated the superior exploitation capability of the linearized inversion with respect to the global optimization, that is the ability to converge towards the minimum in case of objective functions with low gradient values.

Loading

Article metrics loading...

/content/journals/10.1002/nsg.12034
2019-03-08
2024-04-18
Loading full text...

Full text loading...

References

  1. AleardiM. and TognarelliA.2016. The limits of narrow and wide angle ava inversion for high vp/vs ratios: an application to elastic seabed characterization. Journal of Applied Geophysics131, 54–68.
    [Google Scholar]
  2. AleardiM. and MazzottiA.2017. 1D elastic FWI and uncertainty estimation by means of a hybrid genetic algorithm‐gibbs sampler approach. Geophysical Prospecting65, 64–85.
    [Google Scholar]
  3. AleardiM., CiabarriF. and MazzottiA.2017. Probabilistic estimation of reservoir properties by means of wide‐angle AVA inversion and a petrophysical reformulation of the Zoeppritz equations. Journal of Applied Geophysics147, 28–41.
    [Google Scholar]
  4. AsterR.C., BorchersB. and ThurberC.H.2011. Parameter Estimation and Inverse Problems. Academic Press.
    [Google Scholar]
  5. BellmanR.E.1957. Dynamic Programming. Princeton University Press.
    [Google Scholar]
  6. ButtkusB.2000. Spectral analysis and filter theory in applied geophysics. Springer.
    [Google Scholar]
  7. ClercM.1999. The swarm and queen: towards a deterministic and adaptive particle swarm optimization. Proceedings of the IEEE transactions on evolutionary computation, 1951–1957.
  8. DahlinT., LerouxV. and NissenJ.2002. Measuring techniques in induced polarization imaging. Journal of Applied Geophysics50, 279–298.
    [Google Scholar]
  9. DahlinT., RosqvistH. and LerouxV.2010. Resistivity‐IP for landfill applications. First Break28, 101–105.
    [Google Scholar]
  10. De NicolaoA., DrufucaG. and RoccaF.1993. Eigenvalues and eigenvectors of linearized elastic inversion. Geophysics58, 670–679.
    [Google Scholar]
  11. EbbesenS., KiwitzP. and GuzzellaL.2012. A generic particle swarm optimization matlab function. American Control Conference Fairmont Queen Elizabeth. Montréal, Canada.
  12. EffersøF.2000. Resolution of Cole‐Cole parameters based on induced polarization data. In: Methods and Applications of Inversion (eds P.C.Hansen , B.H.Jacobsen and K.Mosegaard ), pp. 120–128. Springer.
    [Google Scholar]
  13. Fernández‐MartínezJ.L., Fernández‐MuñizM.Z. and TompkinsM.J.2012. On the topography of the cost functional in linear and nonlinear inverse problems. Geophysics77, W1–W15.
    [Google Scholar]
  14. FiandacaG., AukenE., ChristiansenA.V., and GazotyA.2012. Time‐domain‐induced polarization: full‐decay forward modelling 1D laterally constrained inversion of Cole‐Cole parameters. Geophysics77, E213–E225.
    [Google Scholar]
  15. GodioA. and SantilanoA.2018. On the optimization of electromagnetic geophysical data: application of the PSO algorithm. Journal of Applied Geophysics148, 163–174.
    [Google Scholar]
  16. GurinG.V, TarasovA.V., IlyinY.T. and TitovK.V.2013. Spectral induced polarization of disseminated ore. Seriya Geologiya i Geografiya, 14–30.
    [Google Scholar]
  17. HördtA., BlaschekR., KemnaA. and ZisserN.2007. Hydraulic conductivity estimation from induced polarization data at the field scale – The Krauthausen case history. Journal of Applied Geophysics66, 33–46.
    [Google Scholar]
  18. JosephS., InghamM. and GouwsG.2016. Spectral‐induced polarization measurements on sieved sands and the relationship to permeability. Water Resources Research52, 4226–4246.
    [Google Scholar]
  19. KennedyJ. and EberhartR.C.2001. Swarm Intelligence. Morgan Kaufmann Publishers.
    [Google Scholar]
  20. KingmanJ.E.E., DonohueJ.G. and RitchieT.J.2007. Distributed acquisition in electrical geophysical systems. Proceedings of Exploration 07: Fifth Decennial International Conference on Mineral Exploration, Decennial Mineral Exploration Conferences, 425–432.
  21. LerouxV., DahlinT. and RosqvistH.2010. Time‐domain IP and resistivity sections measured at four landfills with different contents. 16th European Meeting of Environmental and Engineering Geophysics (Near Surface), Expanded Abstracts, P09.
  22. MajorJ. and SilicJ.1981. Restrictions on the use of Cole‐Cole dispersion model in complex resistivity interpretation. Geophysics46, 916–931.
    [Google Scholar]
  23. MariniF. and WalczakB.2015. Particle Swarm optimization (PSO). A tutorial. Chemometrics and Intelligent Laboratory Systems149, 153–165.
    [Google Scholar]
  24. MenkeW.2012. Geophysical Data Analysis: Discrete Inverse Theory. Academic Press.
    [Google Scholar]
  25. OldenburgD.W. and LiY.1994. Inversion of induced polarization data. Geophysics59, 1327‐1341.
    [Google Scholar]
  26. OllsonP‐I, FinadacaG., LarsenJ.J., DahlinT. and AukenE.2016. Doubling the spectrum of time‐domain induced polarization: removal of non linear self‐potential drift, tapered gain and uncertainty estimation. Proceedings of the 4th International Workshop on Induced Polarization, 14–16.
  27. PaineJ. and CopelandA.2003. Reduction of noise in induced polarization data using full time‐series data. Exploration Geophysics34, 225–228.
    [Google Scholar]
  28. PedersenM.E.H.2010. Good parameters for particle swarm optimization. Hvass Laboratories Copenhagen, Denmark Technical Report HL1001.
  29. PeltonW.H., WardS.H., HallofP.G., SillW.R. and NelsonP.H.1978. Mineral discrimination and removal of inductive coupling with multifrequency IP. Geophysics43(3), 588–609.
    [Google Scholar]
  30. PeltonW.H., SillW.R. and SmithB.D.1983. Interpretation of complex resistivity and dielectric data part I. Geophysical Transactions29, 297–330.
    [Google Scholar]
  31. RobinsonJ. and Rahmat‐SamiiY.2004. Particle swarm optimization in electromagnetics. IEEE Transactions on Antennas and Propagation52, 397–407.
    [Google Scholar]
  32. RiedelM. and TheilenF.2001. AVO investigations of shallow marine sediments. Geophysical Prospecting49, 198–212.
    [Google Scholar]
  33. ReynoldsJ.M.2011. An Introduction to Applied and Environmental Geophysics. Wiley‐Blackwell.
    [Google Scholar]
  34. SajevaA., AleardiM., GaluzziB., StucchiE., SpadavecchiaE. and MazzottiA.2017a. Comparing the performances of four stochastic optimisation methods using analytic objective functions, 1D elastic full‐waveform inversion, and residual static computation. Geophysical Prospecting65, 322–346.
    [Google Scholar]
  35. SajevaA., AleardiM., MazzottiA., StucchiE. and GaluzziB.2014. Comparison of stochastic optimization methods on two analytic objective functions and on a 1D elastic FWI. 76th EAGE Conference and Exhibition 2014.
  36. SajevaA., AleardiM. and MazzottiA.2017b. Genetic algorithm full‐waveform inversion: uncertainty estimation and validation of the results. Bollettino di Geofisica Teorica ed Applicata58, 395–414.
    [Google Scholar]
  37. SambridgeM.1999. Geophysical inversion with a neighbourhood algorithm—II. Appraising the ensemble. Geophysical Journal International138, 727–746.
    [Google Scholar]
  38. SenM.K. and StoffaP.L.2013. Global Optimisation Methods in Geophysical Inversion. Cambridge University Press.
    [Google Scholar]
  39. ShawR. and SrivastavaS.2007. Particle swarm optimisation: a new tool to invert geophysical data. Geophysics72, F75–F83.
    [Google Scholar]
  40. SogadeJ.A, Scira‐ScappuzzoF., VichabianY., ShiW. and MorganF.D.2006. Induced polarization detection and mapping of contaminant plumes. Geophysics71, B75–B84.
    [Google Scholar]
  41. VinciguerraA., CostantiniP. and AleardiM.2018. Inversion of FW‐IP data ‐ sensitivity analysis and synthetic inversions through local and global optimization strategies. 24th European Meeting of Environmental and Engineering Geophysics.
  42. WeigandM., and KemnaA.2016. Debye decomposition of time‐lapse spectral induced polarization data. Computers & Geosciences86, 34–45.
    [Google Scholar]
  43. YilmazÖ.2001. Seismic Data Analysis. Society of Exploration Geophysicists.
    [Google Scholar]
  44. ZimmermannE., KemnaA., BerwixJ., GlassW., MünchH.M. and HuismanJ.A.2008. A high‐accuracy impedance spectrometer for measuring sediments with low polarizability. Measurement Science and Technology19, 105603.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1002/nsg.12034
Loading
/content/journals/10.1002/nsg.12034
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): Inversion; IP; Optimization

Most Cited This Month Most Cited RSS feed

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