Waveform-based simulated annealing of crosshole transmission data: a semi-global method for estimating seismic anisotropy
- 15 October 2014
- journal article
- research article
- Published by Oxford University Press (OUP) in Geophysical Journal International
- Vol. 199 (3), 1586-1607
- https://doi.org/10.1093/gji/ggu307
Abstract
We successfully apply the semi-global inverse method of simulated annealing to determine the best-fitting 1-D anisotropy model for use in acoustic frequency domain waveform tomography. Our forward problem is based on a numerical solution of the frequency domain acoustic wave equation, and we minimize wavefield phase residuals through random perturbations to a 1-D vertically varying anisotropy profile. Both real and synthetic examples are presented in order to demonstrate and validate the approach. For the real data example, we processed and inverted a cross-borehole data set acquired by Vale Technology Development (Canada) Ltd. in the Eastern Deeps deposit, located in Voisey's Bay, Labrador, Canada. The inversion workflow comprises the full suite of acquisition, data processing, starting model building through traveltime tomography, simulated annealing and finally waveform tomography. Waveform tomography is a high resolution method that requires an accurate starting model. A cycle-skipping issue observed in our initial starting model was hypothesized to be due to an erroneous anisotropy model from traveltime tomography. This motivated the use of simulated annealing as a semi-global method for anisotropy estimation. We initially tested the simulated annealing approach on a synthetic data set based on the Voisey's Bay environment; these tests were successful and led to the application of the simulated annealing approach to the real data set. Similar behaviour was observed in the anisotropy models obtained through traveltime tomography in both the real and synthetic data sets, where simulated annealing produced an anisotropy model which solved the cycle-skipping issue. In the real data example, simulated annealing led to a final model that compares well with the velocities independently estimated from borehole logs. By comparing the calculated ray paths and wave paths, we attributed the failure of anisotropic traveltime tomography to the breakdown of the ray-theoretical approximation in the vicinity of strong velocity discontinuities.This publication has 51 references indexed in Scilit:
- Waveform tomography imaging of a megasplay fault system in the seismogenic Nankai subduction zoneEarth and Planetary Science Letters, 2012
- Enhancing base-metal exploration with seismic imagingThis article is one of a series of papers published in this Special Issue on the themeLithoprobe — parameters, processes, and the evolution of a continent.Canadian Journal of Earth Sciences, 2010
- Thematic Set: Full waveform inversion of a deep water ocean bottom seismometer datasetFirst Break, 2010
- Random Heterogeneous Model with Bimodal Velocity Distribution for Methane Hydrate ExplorationExploration Geophysics, 2005
- MONTE CARLO METHODS IN GEOPHYSICAL INVERSE PROBLEMSReviews of Geophysics, 2002
- P- and S-wave anisotropy and shear-wave splitting at pressure and temperature in possible mantle rocks and their relation to the rock fabricPhysics of the Earth and Planetary Interiors, 1993
- A SIMULATED ANNEALING APPROACH TO SEISMIC MODEL OPTIMIZATION WITH SPARSE PRIOR INFORMATION1Geophysical Prospecting, 1991
- Statistical mechanics of combinatorial optimizationPhysical Review A, 1988
- Optimization by Simulated AnnealingScience, 1983
- ACM Transactions on Mathematical Software, 1982