12,798 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
Full field inversion in photoacoustic tomography with variable sound speed
Recently, a novel measurement setup has been introduced to photoacoustic
tomography, that collects data in the form of projections of the full 3D
acoustic pressure distribution at a certain time instant. Existing imaging
algorithms for this kind of data assume a constant speed of sound. This
assumption is not always met in practice and thus leads to erroneous
reconstructions. In this paper, we present a two-step reconstruction method for
full field detection photoacoustic tomography that takes variable speed of
sound into account. In the first step, by applying the inverse Radon transform,
the pressure distribution at the measurement time is reconstructed point-wise
from the projection data. In the second step, one solves a final time wave
inversion problem where the initial pressure distribution is recovered from the
known pressure distribution at the measurement time. For the latter problem, we
derive an iterative solution approach, compute the required adjoint operator,
and show its uniqueness and stability
Geophysical tomography in engineering geology: an overview
An overview of the tomographic interpretation method in engineering
geophysics is presented, considering the two approaches of the deterministic
tomography inversion, developed for rock elasticity analysis, and the
probability tomography imaging developed in the domain of potential fields
methods. The theoretical basis of both approaches is shortly outlined before
showing a laboratory and a field application.Comment: 4 pages, 6 figures. Invited lecture at the 7th International Congress
of the Brazilian Geophysical Society, Salvador, Bahia, Brazil, 28-31 October
200
Accelerated High-Resolution Photoacoustic Tomography via Compressed Sensing
Current 3D photoacoustic tomography (PAT) systems offer either high image
quality or high frame rates but are not able to deliver high spatial and
temporal resolution simultaneously, which limits their ability to image dynamic
processes in living tissue. A particular example is the planar Fabry-Perot (FP)
scanner, which yields high-resolution images but takes several minutes to
sequentially map the photoacoustic field on the sensor plane, point-by-point.
However, as the spatio-temporal complexity of many absorbing tissue structures
is rather low, the data recorded in such a conventional, regularly sampled
fashion is often highly redundant. We demonstrate that combining variational
image reconstruction methods using spatial sparsity constraints with the
development of novel PAT acquisition systems capable of sub-sampling the
acoustic wave field can dramatically increase the acquisition speed while
maintaining a good spatial resolution: First, we describe and model two general
spatial sub-sampling schemes. Then, we discuss how to implement them using the
FP scanner and demonstrate the potential of these novel compressed sensing PAT
devices through simulated data from a realistic numerical phantom and through
measured data from a dynamic experimental phantom as well as from in-vivo
experiments. Our results show that images with good spatial resolution and
contrast can be obtained from highly sub-sampled PAT data if variational image
reconstruction methods that describe the tissues structures with suitable
sparsity-constraints are used. In particular, we examine the use of total
variation regularization enhanced by Bregman iterations. These novel
reconstruction strategies offer new opportunities to dramatically increase the
acquisition speed of PAT scanners that employ point-by-point sequential
scanning as well as reducing the channel count of parallelized schemes that use
detector arrays.Comment: submitted to "Physics in Medicine and Biology
Nearfield Acoustic Holography using sparsity and compressive sampling principles
Regularization of the inverse problem is a complex issue when using
Near-field Acoustic Holography (NAH) techniques to identify the vibrating
sources. This paper shows that, for convex homogeneous plates with arbitrary
boundary conditions, new regularization schemes can be developed, based on the
sparsity of the normal velocity of the plate in a well-designed basis, i.e. the
possibility to approximate it as a weighted sum of few elementary basis
functions. In particular, these new techniques can handle discontinuities of
the velocity field at the boundaries, which can be problematic with standard
techniques. This comes at the cost of a higher computational complexity to
solve the associated optimization problem, though it remains easily tractable
with out-of-the-box software. Furthermore, this sparsity framework allows us to
take advantage of the concept of Compressive Sampling: under some conditions on
the sampling process (here, the design of a random array, which can be
numerically and experimentally validated), it is possible to reconstruct the
sparse signals with significantly less measurements (i.e., microphones) than
classically required. After introducing the different concepts, this paper
presents numerical and experimental results of NAH with two plate geometries,
and compares the advantages and limitations of these sparsity-based techniques
over standard Tikhonov regularization.Comment: Journal of the Acoustical Society of America (2012
- …