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

Application of wavelet analysis in tool wear evaluation using image processing method

Tool wear plays a significant role for proper planning and control of machining parameters to maintain the product quality. However, existing tool wear monitoring methods using sensor signals still have limitations. Since the cutting tool operates directly on the work-piece during machining process, the machined surface provides valuable information about the cutting tool condition. Therefore, the objective of present study is to evaluate the tool wear based on the workpiece profile signature by using wavelet analysis. The effect of wavelet families, scale of wavelet and statistical features of the continuous wavelet coefficient on the tool wear is studied. The surface profile of workpiece was captured using a DSLR camera. Invariant moment method was applied to extract the surface profile up to sub-pixel accuracy. The extracted surface profile was analyzed by using continuous wavelet transform (CWT) written in MATLAB. The re-sults showed that average, RMS and peak to valley of CWT coefficients at all scale increased with tool wear. Peak to valley at higher scale is more sensitive to tool wear. Haar was found to be more effective and significant to correlate with tool wear with highest R2 which is 0.9301

Waterborne GPR survey for estimating bottom-sediment variability: A survey on the Po River, Turin, Italy

We conducted an integrated geophysical survey on a stretch of the river Po in order to check the GPR ability to discriminate the variability of riverbed sediments through an analysis of the bottom reflection amplitudes. We conducted continuous profiles with a 200-MHzGPR system and a handheld broadband EM sensor.Aconductivity meter and a TDR provided punctual measurements of water conductivity, permittivity, and temperature. The processing and interpretation of the GEM-2 and GPR data were enhanced by reciprocal results and by integration with the punctual measurements of the EM properties of the water. We used a processing flow that improved the radargram images and preserved the amplitude ratios among the different profiles and the frequency content at the bottom reflection signal.We derived the water attenuation coefficient both from the punctual measurements using the Maxwell formulas and from the interpretation of the GPR data, finding an optimal matching between the two values. The GPR measurements provided maps of the bathymetry and of the bottom reflection amplitude. The high reflectivity of the riverbed, derived from the GPR interpretation, agreed with the results of the direct sampling campaign that followed the geophysical survey. The variability of the bottom-reflection-amplitudes map, which was not confirmed by the direct sampling, could also have been caused by scattering phenomena due to the riverbed clasts which are dimensionally comparable to the wavelength of the radar pulse

The origin of fine scale acoustic stratigraphy in deep-sea carbonates

In this paper we investigate the origin and geologic significance of the closely spaced high-frequency subbottom acoustic reflectors characteristic of pelagic carbonates. A detailed survey was conducted of a small area in the equatorial Pacific with the Marine Physical Laboratory\u27s Deep-Tow instrument package, providing high-resolution 4-kHz profiles and precise positioning of core samples. The cores were sampled at closely spaced intervals for sound velocity and saturated bulk density. Acoustic impedances were calculated, and a reflection coefficient log determined for the upper 10 m of the sediment column. The reflection coefficient log revealed no interfaces with large reflection coefficients that correlated with the reflectors seen on the Deep-Tow 4-kHz seismic profile. The calculated reflection coefficients were very low (typically 10−3–10−5) and varied about a wavelength that was on the order of the wavelength of the 4-kHz pulse, implying that interference plays a role in the composition of the seismic record. Convolving the outgoing 4-kHz pulse with the reflection coefficient log generated a synthetic seismogram that very closely resembled the 4-kHz reflection profile. Varying the frequency of the outgoing pulse changed the amplitude and position of the reflectors seen on the synthetic seismograms. Thus we conclude that the reflectors seen on the 4-kHz seismic profile were not caused by discrete geologic horizons but rather are the result of the interference of many small layers

Retrieving shallow shear-wave velocity profiles from 2D seismic-reflection data with severely aliased surface waves

The inversion of surface-wave phase-velocity dispersion curves provides a reliable method to derive near-surface shear-wave velocity profiles. In this work, we invert phase-velocity dispersion curves estimated from 2D seismic-reflection data. These data cannot be used to image the first 50 m with seismic-reflection processing techniques due to the presence of indistinct first breaks and significant NMO-stretching of the shallow reflections. A surface-wave analysis was proposed to derive information about the near surface in order to complement the seismic-reflection stacked sections, which are satisfactory for depths between 50 and 700 m. In order to perform the analysis, we had to overcome some problems, such as the short acquisition time and the large receiver spacing, which resulted in severe spatial aliasing. The analysis consists of spatial partitioning of each line in segments, picking of the phase-velocity dispersion curves for each segment in the f-k domain, and inversion of the picked curves using the neighborhood algorithm. The spatial aliasing is successfully circumvented by continuously tracking the surface-wave modal curves in the f-k domain. This enables us to sample the curves up to a frequency of 40 Hz, even though most components beyond 10 Hz are spatially aliased. The inverted 2D VS sections feature smooth horizontal layers, and a sensitivity analysis yields a penetration depth of 20–25 m. The results suggest that long profiles may be more efficiently surveyed by using a large receiver separation and dealing with the spatial aliasing in the described way, rather than ensuring that no spatially aliased surface waves are acquired.Fil: Onnis, Luciano Emanuel. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; Argentina. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Física; ArgentinaFil: Osella, Ana Maria. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; Argentina. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Física; ArgentinaFil: Carcione, Jose M.. Istituto Nazionale di Oceanografia e di Geofisica Sperimentale; Itali

Basis Pursuit Receiver Function

Receiver functions (RFs) are derived by deconvolution of the horizontal (radial or transverse) component of ground motion from the vertical component, which segregates the PS phases. Many methods have been proposed to employ deconvolution in frequency as well as in time domain. These methods vary in their approaches to impose regularization that addresses the stability problem. Here, we present application of a new time-domain deconvolution technique called basis pursuit deconvolution (BPD) that has recently been applied to seismic exploration data. Unlike conventional deconvolution methods, the BPD uses an L1 norm constraint on model reflectivity to impose sparsity. In addition, it uses an overcomplete wedge dictionary based on a dipole reflectivity series to define model constraints, which can achieve higher resolution than that obtained by the traditional methods. We demonstrate successful application of BPD based RF estimation from synthetic data for a crustal model with a near-surface thin layer of thickness 5, 7, 10, and 15 km. The BPD can resolve these thin layers better with much improved signal-to-noise ratio than the conventional methods. Finally, we demonstrate application of the BPD receiver function (BPRF) method to a field dataset from Kutch, India, where near-surface sedimentary layers are known to be present. The BPRFs are able to resolve reflections from these layers very well.Jackson Chair funds at the Jackson School of Geosciences, University of Texas, AustinCouncil of Scientific and Industrial Research twelfth five year plan project at the Council of Scientific and Industrial Research National Geophysical Research Institute (CSIR-NGRI), HyderabadInstitute for Geophysic

Geo-neutrinos and Earth Models

We present the current status of geo-neutrino measurements and their implications for radiogenic heating in the mantle. Earth models predict different levels of radiogenic heating and, therefore, different geo-neutrino fluxes from the mantle. Seismic tomography reveals features in the deep mantle possibly correlated with radiogenic heating and causing spatial variations in the mantle geo-neutrino flux at the Earth surface. An ocean-based observatory offers the greatest sensitivity to the mantle flux and potential for resolving Earth models and mantle features. Refinements to estimates of the geo-neutrino flux from continental crust reduce uncertainty in measurements of the mantle flux, especially measurements from land-based observatories. These refinements enable the resolution of Earth models using the combined measurements from multiple continental observatories.Comment: 9 pages, 4 figures; Contributed paper TAUP 201