Vertical profile of effective turbidity reconstructed from broadening of incoherent body-wave pulses—I. General approach and the inversion procedure

A method is developed for the reconstruction of a non-uniform distribution of scattering properties in the upper layers of the Earth using data on broadening of an incoherent body-wave group or pulse along a number of rays. The theoretical basis for this reconstruction is a linear integral formula after Bocharov (1985, 1988), which is employed to design a linear inversion procedure. The inversion is performed in terms of a single scalar parameter of effective turbidity. This parameter presents an adequate generalization of the common turbidity parameter used in the isotropic scattering case; it describes, simultaneously, scattering attenuation, pulse broadening and backscattering or coda formation. As a preliminary step, necessary conditions of applicability of the transport equation approach for the analysis of regional high-frequency seismic waves are verified. A new compact derivation of Bocharov's formula is then presented. A linear least-squares inversion procedure for recovering a layered turbidity structure is proposed and tested on synthetic data of onset-to-peak delays of incoherent body-wave pulses. A few practical aspects of the application of the general approach to seismological data are analysed, including the correctness of the low-angle approximation, the use of peak delay observations instead of pulse centroid, the effects of a realistic spatial spectrum of inhomogeneity field, the potential bias produced by intrinsic loss, and the distortions produced by a non-spherical (double dipole) source radiation pattern. The latter point is considered as critically important, as one can expect significant data contamination by nodal arrivals. An efficient robust estimation procedure is designed and tested that is capable of suppressing distortions from nodal and near-nodal data.     

Modelling and inversion of single-ended refraction data from the shot gathers of multifold deep seismic reflection profiling—an approach for deriving the shallow velocity structure

A multifold crustal-scale deep seismic near-vertical reflection profile generates a large number of single-ended shot gathers, which provide redundant data sets because of overlapping coverage of the shallow refractors. We present an approach for deriving the shallow velocity structure by modelling and inversion of single-ended seismic refraction first arrival traveltime data. We apply this method to a data set acquired with a 12-km long spread with 100 m spacing of shots and receivers, of the Neoproterozoic Marwar basin in the NW Indian shield. The approach is shown to be quite successful for delineating the shallow refractor depths, steep dips and velocities, even in the absence of regular reverse refraction profiles. The study reveals two-layered sedimentary formations, Malani volcanics and a complicated basement configuration of the Marwar basin, and provides a measure of resolution and uncertainty of the estimated model parameters. A seismic section of the near-trace gather is found to be qualitatively consistent with the derived structural features of the basin. The relative highs and lows, observed in the Bouguer gravity profile, further corroborate the derived velocity model. The present approach can be especially useful in offshore areas and elsewhere, where the single-ended multifold seismic profiles are the only available data sets.