Forward modelling in cross-hole seismic tomography using reciprocity

Abstract

Seismic transmission tomography provides a means of direct estimation of compressional wave velocities using first arrival traveltimes from cross-hole data. Forward modelling in seismic-tomography involves computation of traveltimes at the receiver positions and estimation of ray travelpaths for a subsurface velocity model and source-receiver geometry. Conventional ray-tracing techniques face many difficulties when applied to real-life situations. These can be overcome to a great extent by using a ray-tracing technique based on the principle of reciprocity. In the present work, an algorithm is developed which makes use of the reciprocity principle and the Format's principle for estimation of raypaths. In this technique seismic traveltime calculation is based on a two-dimensional dynamic programming approach, involving systematic mapping of traveltimes over a grid of constant-velocity cells. The algorithm works for even highly complicated velocity models. First arrival seismic energy can travel either as transmitted waves, diffracted waves, or head waves, and this technique simulates all of them. The algorithm is tested for various subsurface velocity models, some of which are presented in this paper. The results are represented as ray diagrams, superimposed over the actual velocity model. The isotime lines depicting the wavefronts are also presented, which reveal the effect of velocity contrast on seismic wave propagation.