A Graph-Space Optimal Transport Approach Based on Kaniadakis κ-Gaussian Distribution for Inverse Problems Related to Wave Propagation.

Data-centric inverse problems are a process of inferring physical attributes from indirect measurements. Full-waveform inversion (FWI) is a non-linear inverse problem that attempts to obtain a quantitative physical model by comparing the wave equation solution with observed data, optimizing an objective function. However, the FWI is strenuously dependent on a robust objective function, especially for dealing with cycle-skipping issues and non-Gaussian noises in the dataset. In this work, we present an objective function based on the Kaniadakis κ-Gaussian distribution and the optimal transport (OT) theory to mitigate non-Gaussian noise effects and phase ambiguity concerns that cause cycle skipping. We construct the κ-objective function using the probabilistic maximum likelihood procedure and include it within a well-posed version of the original OT formulation, known as the Kantorovich-Rubinstein metric. We represent the data in the graph space to satisfy the probability axioms required by the Kantorovich-Rubinstein framework. We call our proposal the κ-Graph-Space Optimal Transport FWI (κ-GSOT-FWI). The results suggest that the κ-GSOT-FWI is an effective procedure to circumvent the effects of non-Gaussian noise and cycle-skipping problems. They also show that the Kaniadakis κ-statistics significantly improve the FWI objective function convergence, resulting in higher-resolution models than classical techniques, especially when κ=0.6.

Full-waveform inversion based on generalized Rényi entropy using patched Green's function techniques.

The estimation of physical parameters from data analyses is a crucial process for the description and modeling of many complex systems. Based on Rényi α-Gaussian distribution and patched Green's function (PGF) techniques, we propose a robust framework for data inversion using a wave-equation based methodology named full-waveform inversion (FWI). From the assumption that the residual seismic data (the difference between the modeled and observed data) obeys the Rényi α-Gaussian probability distribution, we introduce an outlier-resistant criterion to deal with erratic measures in the FWI context, in which the classical FWI based on l2-norm is a particular case. The new misfit function arises from the probabilistic maximum-likelihood method associated with the α-Gaussian distribution. The PGF technique works on the forward modeling process by dividing the computational domain into outside target area and target area, where the wave equation is solved only once on the outside target (before FWI). During the FWI processing, Green's functions related only to the target area are computed instead of the entire computational domain, saving computational efforts. We show the effectiveness of our proposed approach by considering two distinct realistic P-wave velocity models, in which the first one is inspired in the Kwanza Basin in Angola and the second in a region of great economic interest in the Brazilian pre-salt field. We call our proposal by the abbreviation α-PGF-FWI. The results reveal that the α-PGF-FWI is robust against additive Gaussian noise and non-Gaussian noise with outliers in the limit α → 2/3, being α the Rényi entropic index.

Statistical Modeling of the Seismic Moments via Mathai Distribution.

Mathai's pathway model is playing an increasingly prominent role in statistical distributions. As a generalization of a great variety of distributions, the pathway model allows the studying of several non-linear dynamics of complex systems. Here, we construct a model, called the Pareto-Mathai distribution, using the fact that the earthquakes' magnitudes of full catalogues are well-modeled by a Mathai distribution. The Pareto-Mathai distribution is used to study artificially induced microseisms in the mining industry. The fitting of a distribution for entire range of magnitudes allow us to calculate the completeness magnitude (Mc). Mathematical properties of the new distribution are studied. In addition, applying this model to data recorded at a Chilean mine, the magnitude Mc is estimated for several mine sectors and also the entire mine.