Kirchhoff synthetics for seismic reflections and diving waves in two dimensions and application to the D" layer

The Kirchhoff-Helmholtz (KH) integration has been used to model the reflected and the diving waves from an interface with a positive velocity gradient. The modelling is carried out for a spherical boundary and for a sinusoidal topography with a long-scale wavelength.
An artefact, which is a major problem in modelling the seismic response using the KH integration, has been reduced by introducing a Hilbert transform sign manipulation. Cleaner synthetic seismograms with correct amplitudes have been produced by this method. A discretization in larger surface elements has been made possible by introducing a smoothing factor that suppresses the noise that normally follows the constructed signal if a large element size is taken.     

An approach to reserve estimation enhanced with 3-D seismic data

Reserve estimation for hydrocarbon reservoirs can be improved by incorporating values extracted from three-dimensional (3-D) seismic data with those obtained from more conventional data sources of data, such as drill-core and well-log data. An example of this improved method is illustrated by an application to the QW pool located in the Buohaiwan Basin in eastern China. Parameter values extracted from 3-D seismic data extend the knowledge about the spatial distributions of such reservoir parameters as net thickness, porosity, and oil saturation. To assist in the extraction of these values, different pattern-recognition techniques can be applied. The results that are obtained by this method offer a more reliable and more credible approach to reserve estimation and can be applied at every stage of resource extraction from exploration to development.     

Sensitivities of seismic traveltimes and amplitudes in reflection tomography

Seismic traveltimes and amplitudes in reflection-seismic data show different dependences on the geometry of reflection interfaces, and on the variation of interval velocities. These dependences are revealed by eigenanalysis of the Hessian matrix, defined in terms of the Fréchet matrix and its adjoint associated with different norms chosen in the model space. The eigenvectors and eigenvalues of the Hessian clearly show that for reflection tomographic inversion, traveltime and amplitude data contain complementary information. Both for reflector-geometry and for interval-velocity variations, the traveltimes are sensitive to the model components with small wavenumbers, whereas the amplitudes are more sensitive to the components with high wavenumbers. The model resolution matrices, after the rejection of eigenvectors corresponding to small eigenvalues, give us some insight into how the addition of amplitude information could potentially contribute to the recovery of physical parameters.
In order to cooperatively invert seismic traveltimes and amplitudes simultaneously, we propose an empirical definition of the data covariance matrix which balances the relative sensitivities of different types of data. We investigate the cooperative use of both data types for, separately, interface-geometry and 2-D interval-velocity variations. In both cases we find that cooperative inversions can provide better solutions than those using traveltimes alone. The potential benefit of including amplitude-data constraints in seismic-reflection traveltime tomography is therefore that it may be possible to resolve the known ambiguity between the reflector-depth uncertainty and the interval-velocity uncertainty better.     

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.     

An efficient measure of compactness for two-dimensional shapes and its application in regionalization problems

A measure of shape compactness is a numerical quantity representing the degree to which a shape is compact. Ways to provide an accurate measure have been given great attention due to its application in a broad range of GIS problems, such as detecting clustering patterns from remote-sensing images, understanding urban sprawl, and redrawing electoral districts to avoid gerrymandering. In this article, we propose an effective and efficient approach to computing shape compactness based on the moment of inertia (MI), a well-known concept in physics. The mathematical framework and the computer implementation for both raster and vector models are discussed in detail. In addition to computing compactness for a single shape, we propose a computational method that is capable of calculating the variations in compactness as a shape grows or shrinks, which is a typical application found in regionalization problems. We conducted a number of experiments that demonstrate the superiority of the MI over the popular isoperimetric quotient approach in terms of (1) computational efficiency; (2) tolerance of positional uncertainty and irregular boundaries; (3) ability to handle shapes with holes and multiple parts; and (4) applicability and efficacy in districting/zonation/regionalization problems.