SUMMARY
We provide a two-stage approach to extract spectral induced polarization (SIP) information from time-domain IP data. In the first stage we invert dc data to recover the background conductivity. In the second, we solve a linear inverse problem and invert all time channels simultaneously to recover the IP parameters. The IP decay curves are represented by a stretched exponential (SE) rather than the traditional Cole–Cole model, and we find that defining the parameters in terms of their logarithmic values is advantageous. To demonstrate the capability of our simultaneous SIP inversion we use synthetic data simulating a porphyry mineral deposit. The challenge is to image a mineral body that is hosted within an alteration halo having the same chargeability but a different time constant. For a 2-D problem, we were able to distinguish the body using our simultaneous inversion but we were not successful in using a sequential (or conventional) SIP inversion approach. For the 3-D problem we recovered 3-D distributions of the SIP parameters and used those to construct a 3-D rock model having four rock units. Three chargeable units were distinguished. The compact mineralization zone, having a large time constant, was distinguished from the circular alteration halo that had a small time constant. Finally, to promote the use of the SIP technique, and to have further development of SIP inversion, all examples presented in this paper are available in our open source resources (https://github.com/simpeg-research/kang-2018-spectral-inducedpolarization).
