Seismic Ray Impedance Inversion

This thesis investigates a prestack seismic inversion scheme implemented in the ray parameter domain. Conventionally, most prestack seismic inversion methods are performed in the incidence angle domain. However, inversion using the concept of ray impedance, as it honours ray path variation following the elastic parameter variation according to Snell’s law, shows the capacity to discriminate different lithologies if compared to conventional elastic impedance inversion. The procedure starts with data transformation into the ray-parameter domain and then implements the ray impedance inversion along constant ray-parameter profiles. With different constant-ray-parameter profiles, mixed-phase wavelets are initially estimated based on the high-order statistics of the data and further refined after a proper well-to-seismic tie. With the estimated wavelets ready, a Cauchy inversion method is used to invert for seismic reflectivity sequences, aiming at recovering seismic reflectivity sequences for blocky impedance inversion. The impedance inversion from reflectivity sequences adopts a standard generalised linear inversion scheme, whose results are utilised to identify rock properties and facilitate quantitative interpretation. It has also been demonstrated that we can further invert elastic parameters from ray impedance values, without eliminating an extra density term or introducing a Gardner’s relation to absorb this term. Ray impedance inversion is extended to P-S converted waves by introducing the definition of converted-wave ray impedance. This quantity shows some advantages in connecting prestack converted wave data with well logs, if compared with the shearwave elastic impedance derived from the Aki and Richards approximation to the Zoeppritz equations. An analysis of P-P and P-S wave data under the framework of ray impedance is conducted through a real multicomponent dataset, which can reduce the uncertainty in lithology identification.Inversion is the key method in generating those examples throughout the entire thesis as we believe it can render robust solutions to geophysical problems. Apart from the reflectivity sequence, ray impedance and elastic parameter inversion mentioned above, inversion methods are also adopted in transforming the prestack data from the offset domain to the ray-parameter domain, mixed-phase wavelet estimation, as well as the registration of P-P and P-S waves for the joint analysis. The ray impedance inversion methods are successfully applied to different types of datasets. In each individual step to achieving the ray impedance inversion, advantages, disadvantages as well as limitations of the algorithms adopted are detailed. As a conclusion, the ray impedance related analyses demonstrated in this thesis are highly competent compared with the classical elastic impedance methods and the author would like to recommend it for a wider application

Frequency-dependent AVO attribute: theory and example

Fluid-saturated rocks generally have seismic velocities that depend upon frequency. Exploring this property may help us discriminate different fluids from seismic data. In this paper, we introduce a scheme to calculate a frequency-dependent AVO attribute in order to estimate seismic dispersion from pre-stack data, and apply it to North Sea data. The scheme essentially combines the two-term approximation of Smith and Gidlow (1987) with the method of spectral decomposition based on the Wigner-Ville distribution, which is used to achieve high resolution. The result suggests the potential of this method for detection of seismic dispersion due to fluid saturation

Numerical treatment of seismic accelerograms and of inelastic seismic structural responses using harmonic wavelets

The harmonic wavelet transform is employed to analyze various kinds of nonstationary signals common in aseismic design. The effectiveness of the harmonic wavelets for capturing the temporal evolution of the frequency content of strong ground motions is demonstrated. In this regard, a detailed study of important earthquake accelerograms is undertaken and smooth joint time-frequency spectra are provided for two near-field and two far-field records; inherent in this analysis is the concept of the mean instantaneous frequency. Furthermore, as a paradigm of usefulness for aseismic structural purposes, a similar analysis is conducted for the response of a 20-story steel frame benchmark building considering one of the four accelerograms scaled by appropriate factors as the excitation to simulate undamaged and severely damaged conditions for the structure. The resulting joint time-frequency representation of the response time histories captures the influence of nonlinearity on the variation of the effective natural frequencies of a structural system during the evolution of a seismic event. In this context, the potential of the harmonic wavelet transform as a detection tool for global structural damage is explored in conjunction with the concept of monitoring the mean instantaneous frequency of records of critical structural responses

A non-separable stochastic model for pulse-like ground motions

A phenomenological non-separable non-stationary stochastic model is proposed to represent near-fault pulse-like ground motions (PLGMs) by means of a parametrically defined evolutionary power spectrum (EPSD). Numerical data pertaining to ensembles of EPSD compatible realizations and considering statistical analysis of peak elastic and inelastic spectral ordinates demonstrate the applicability of the model to capture the salient effects of PLGMs to structural responses. To this aim, the model parameters are calibrated against a field recorded PLGM. Further numerical data considering stochastic processes compatible with the response spectrum of the European aseismic code (EC8) are furnished to demonstrate the potential of the proposed model for including near-fault effects to spectrum compatible representations of the seismic action. It is foreseen that this model can be a useful tool in accounting for the low-frequency content of PLGMs in both Monte Carlo simulation-based analyses and in statistical linearization based studies