8,212 research outputs found
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
Studies in Signal Processing Techniques for Speech Enhancement: A comparative study
Speech enhancement is very essential to suppress the background noise and to increase speech intelligibility and reduce fatigue in hearing. There exist many simple speech enhancement algorithms like spectral subtraction to complex algorithms like Bayesian Magnitude estimators based on Minimum Mean Square Error (MMSE) and its variants. A continuous research is going and new algorithms are emerging to enhance speech signal recorded in the background of environment such as industries, vehicles and aircraft cockpit. In aviation industries speech enhancement plays a vital role to bring crucial information from pilot’s conversation in case of an incident or accident by suppressing engine and other cockpit instrument noises. In this work proposed is a new approach to speech enhancement making use harmonic wavelet transform and Bayesian estimators. The performance indicators, SNR and listening confirms to the fact that newly modified algorithms using harmonic wavelet transform indeed show better results than currently existing methods. Further, the Harmonic Wavelet Transform is computationally efficient and simple to implement due to its inbuilt decimation-interpolation operations compared to those of filter-bank approach to realize sub-bands
Intermittent process analysis with scattering moments
Scattering moments provide nonparametric models of random processes with
stationary increments. They are expected values of random variables computed
with a nonexpansive operator, obtained by iteratively applying wavelet
transforms and modulus nonlinearities, which preserves the variance. First- and
second-order scattering moments are shown to characterize intermittency and
self-similarity properties of multiscale processes. Scattering moments of
Poisson processes, fractional Brownian motions, L\'{e}vy processes and
multifractal random walks are shown to have characteristic decay. The
Generalized Method of Simulated Moments is applied to scattering moments to
estimate data generating models. Numerical applications are shown on financial
time-series and on energy dissipation of turbulent flows.Comment: Published in at http://dx.doi.org/10.1214/14-AOS1276 the Annals of
Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical
Statistics (http://www.imstat.org
Geodesics on the manifold of multivariate generalized Gaussian distributions with an application to multicomponent texture discrimination
We consider the Rao geodesic distance (GD) based on the Fisher information as a similarity measure on the manifold of zero-mean multivariate generalized Gaussian distributions (MGGD). The MGGD is shown to be an adequate model for the heavy-tailed wavelet statistics in multicomponent images, such as color or multispectral images. We discuss the estimation of MGGD parameters using various methods. We apply the GD between MGGDs to color texture discrimination in several classification experiments, taking into account the correlation structure between the spectral bands in the wavelet domain. We compare the performance, both in terms of texture discrimination capability and computational load, of the GD and the Kullback-Leibler divergence (KLD). Likewise, both uni- and multivariate generalized Gaussian models are evaluated, characterized by a fixed or a variable shape parameter. The modeling of the interband correlation significantly improves classification efficiency, while the GD is shown to consistently outperform the KLD as a similarity measure
Parameter Estimation from Time-Series Data with Correlated Errors: A Wavelet-Based Method and its Application to Transit Light Curves
We consider the problem of fitting a parametric model to time-series data
that are afflicted by correlated noise. The noise is represented by a sum of
two stationary Gaussian processes: one that is uncorrelated in time, and
another that has a power spectral density varying as . We present
an accurate and fast [O(N)] algorithm for parameter estimation based on
computing the likelihood in a wavelet basis. The method is illustrated and
tested using simulated time-series photometry of exoplanetary transits, with
particular attention to estimating the midtransit time. We compare our method
to two other methods that have been used in the literature, the time-averaging
method and the residual-permutation method. For noise processes that obey our
assumptions, the algorithm presented here gives more accurate results for
midtransit times and truer estimates of their uncertainties.Comment: Accepted in ApJ. Illustrative code may be found at
http://www.mit.edu/~carterja/code/ . 17 page
Bayesian linear inverse problems in regularity scales
We obtain rates of contraction of posterior distributions in inverse problems
defined by scales of smoothness classes. We derive abstract results for general
priors, with contraction rates determined by Galerkin approximation. The rate
depends on the amount of prior concentration near the true function and the
prior mass of functions with inferior Galerkin approximation. We apply the
general result to non-conjugate series priors, showing that these priors give
near optimal and adaptive recovery in some generality, Gaussian priors, and
mixtures of Gaussian priors, where the latter are also shown to be near optimal
and adaptive. The proofs are based on general testing and approximation
arguments, without explicit calculations on the posterior distribution. We are
thus not restricted to priors based on the singular value decomposition of the
operator. We illustrate the results with examples of inverse problems resulting
from differential equations.Comment: 34 page
Spurious Long Memory in Commodity Futures: Implications for Agribusiness Option Pricing
Long memory, and more precisely fractionally integration, has been put forward as an explanation for the persistence of shocks in a number of economic time series data as well as to reconcile misleading findings of unit roots in data that should be stationary. Recent evidence suggests that long memory characterizes not commodity futures prices but rather price volatility (generally defined as norms of price logreturns). One implication of long memory in volatility is the mispricing of options written on commodity futures, the consequence of which is that fractional Brownian motion should replace geometric Brownian motion as the building block for option pricing solutions. This paper asks whether findings of long memory in volatility might be spurious and caused either by fragile and inaccurate estimation methods and standard errors, by correlated short memory dynamics, or by alternative data generating processes proven to generate the illusion of long memory. We find that for nine out of eleven agricultural commodities for which futures contracts are traded, long memory is spurious but is not caused by the effect of short memory. Alternative explanations are addressed and implications for option pricing are highlighted.Q13, Q14, Marketing, C52, C53, G12, G13,
- …