A fast non-smooth regularisation method is proposed to solve the difficulties involved in the classical least-squares Kirchhoff migration. It not only accounts for the irregular and incomplete data sampling, but also compensates for the anomalous ray coverage and multipathing problem. Numerical experiments show that the method works well.

Amplification characteristics are critical in evaluating the reliability in seismic design for engineering and seismic source and crustal attenuation characteristics for seismology. Horizontal-to-vertical spectral ratios were analysed using six seismograms observed at four sites near Iedang Reservoir in Korea. Amplification characteristics were compared using the S-wave, coda wave and background noise of each seismogram.

In this study, the geophysical strata rating (GSR), which is an empirical measure of rock competency, is calculated from petrophysical data. The GSR is then extended to the whole South Pars gas field in the framework of 3D seismic data through an acoustic impedance poststack seismic inversion.

The projection onto convex sets (POCS) method is deduced using the iterative hard thresholding (IHT) algorithm and a projection operator with more detailed physical illustrations. The interpolation performances on noise-free and noisy data are explained in detail and the reasons behind these performances are fully discussed, which provide clues to further improve interpolation accuracy.

To maintain the initial resolution of the seismic image from seismic sparse acquisition, we propose several ways to sample data irregularly but periodically. At every processing step, we quantified the effect of interpolation by comparing the results with those from the fully sampled data. The result shows that using 60% of the available data is feasible.

Building on the concepts of cohesion degree and local relaxation, we propose an integrated hierarchical equilibrium parallel finite-element reverse time migration (HEP-FE-RTM) algorithm, which has the distinct advantage in that parallel efficiency does not decrease as the number of processors increases.

We image pre-existing fractures using the elastic reverse time migration with source-independent converted phase (ERTM SICP) imaging condition based on observed seismic signals during microseismic monitoring. However, the results showed linear spurious events which are caused by velocity differences between P- and S-waves. To suppress these events, we modified the imaging condition by adding a weighting function calculated from the Poynting vector of the P- and S-waves.

A new methodology for quantifying reservoir heterogeneities based on a Monte Carlo parameter estimation technique using fractal parameters derived from the von Karman autocorrelation function is proposed. The reservoir heterogeneity parameters computed from sonic logs in an oilfield show that these parameters can be discriminants for depositional environments.

The reverse time migration (RTM) equations for staggered-grid high order finite difference scheme incorporating a perfectly matched layer boundary for vertically transversely isotropic (VTI) media are proposed, and checkpoints and GPU accelerated techniques are utilised for data storage and computation efficiency. Hess 2D model tests demonstrate the effectiveness of the proposed algorithm.

In this paper, we propose a new finite difference (FD) scheme which uses different staggered grid FD operators for different first order spatial derivatives in the first order acoustic wave equation. The dispersion analysis and numerical simulation demonstrated the effectiveness of the proposed method.

We propose the non-equispaced fast discrete curvelet transform-based reconstruction method designed for 3D seismic data that are non-uniformly sampled along two spatial coordinates. The method we propose, which has a strong anti-aliasing and anti-noise ability, can be used to reconstruct the non-uniform sampled data to a specified uniform grid.

We use high-resolution aeromagnetic and 3D reflection seismics to delineate geological structures and bodies defined as loss-of-ground features in the Bushveld Complex, South Africa. These include faults, replacement pegmatoids and slump features. This is done for more efficient mine planning and risk reduction.

Airborne collection of electromagnetic and potential field data is a common strategy for extensive resource exploration and reconnaissance. Since these datasets contain information about different properties at different depths, they are normally considered complementary and are interpreted separately. Using airborne data acquired in Western Australia, we explore the viability of their joint inversion and show the advantages of their combined analysis and interpretation.

Multiple-point geostatistical simulation (MPS) was applied to develop 3D ore models matched to surrounding geological information accompanying aeromagnetic data using training image. The present study proposes a method for reducing the uncertainty of the 3D ore model, applying MPS to create probabilistic ore models and analysing the correlation between the models and geophysical data.