144 research outputs found
Linearized inverse scattering based on seismic Reverse Time Migration
In this paper we study the linearized inverse problem associated with imaging
of reflection seismic data. We introduce an inverse scattering transform
derived from reverse-time migration (RTM). In the process, the explicit
evaluation of the so-called normal operator is avoided, while other
differential and pseudodifferential operator factors are introduced. We prove
that, under certain conditions, the transform yields a partial inverse, and
support this with numerical simulations. In addition, we explain the recently
discussed 'low-frequency artifacts' in RTM, which are naturally removed by the
new method
Source-Indexed Migration Velocity Analysis with Global Passive Data
The reverse-time migration of global seismic data generated by free-surface multiples is regularly used to constrain the crustal structure, but its accuracy is to a large extent determined by the accuracy of the 3-D background velocity model used for wave propagation. To this improve the velocity model and hence the accuracy of the migrated image, we wish to apply the technique of migration velocity analysis (MVA) to global passive data. Applications of MVA in the active setting typically focus on o ffset- or angle-gather annihilation, a process that takes advantage of data redundancy to form an extended image, and then applies an annihilation operator to determine the success of image formation. Due to the nature of regional-scale passive seismic arrays, it is unlikely that the data in most of these studies will be su cient to form an extended image volume for use in annihilation-based MVA. In order to make use of the sparse and irregular array design of these arrays, we turn towards a shot-pro le moveout scheme for migration velocity analysis introduced by Xie and Yang (2008). In the place of extended image annihilation, we determine the success of the migration velocity model by using a weighted image correlation power norm. We compare pairs of images formed by migrating each teleseismic source by image cross-correlation in the depth direction. We look for a suitable background model by penalizing the amount of correlation power away from zero depth shift. The total weighted correlation power between source-pro le images is then used as the error function and optimized via conjugate gradient. We present the method and a proof-of-concept with 2-D synthetic data
Convergence Rates for Learning Linear Operators from Noisy Data
This paper studies the learning of linear operators between
infinite-dimensional Hilbert spaces. The training data comprises pairs of
random input vectors in a Hilbert space and their noisy images under an unknown
self-adjoint linear operator. Assuming that the operator is diagonalizable in a
known basis, this work solves the equivalent inverse problem of estimating the
operator's eigenvalues given the data. Adopting a Bayesian approach, the
theoretical analysis establishes posterior contraction rates in the infinite
data limit with Gaussian priors that are not directly linked to the forward map
of the inverse problem. The main results also include learning-theoretic
generalization error guarantees for a wide range of distribution shifts. These
convergence rates quantify the effects of data smoothness and true eigenvalue
decay or growth, for compact or unbounded operators, respectively, on sample
complexity. Numerical evidence supports the theory in diagonal and non-diagonal
settings.Comment: To appear in SIAM/ASA Journal on Uncertainty Quantification (JUQ); 34
pages, 5 figures, 2 table
Parameter identification in a semilinear hyperbolic system
We consider the identification of a nonlinear friction law in a
one-dimensional damped wave equation from additional boundary measurements.
Well-posedness of the governing semilinear hyperbolic system is established via
semigroup theory and contraction arguments. We then investigte the inverse
problem of recovering the unknown nonlinear damping law from additional
boundary measurements of the pressure drop along the pipe. This coefficient
inverse problem is shown to be ill-posed and a variational regularization
method is considered for its stable solution. We prove existence of minimizers
for the Tikhonov functional and discuss the convergence of the regularized
solutions under an approximate source condition. The meaning of this condition
and some arguments for its validity are discussed in detail and numerical
results are presented for illustration of the theoretical findings
A global Carleman estimate in a transmission wave equation and application to a one-measurement inverse problem
We consider a transmission wave equation in two embedded domains in ,
where the speed is in the inner domain and in the outer
domain. We prove a global Carleman inequality for this problem under the
hypothesis that the inner domain is strictly convex and . As a
consequence of this inequality, uniqueness and Lip- schitz stability are
obtained for the inverse problem of retrieving a stationary potential for the
wave equation with Dirichlet data and discontinuous principal coefficient from
a single time-dependent Neumann boundary measurement
- …