scholarly journals A Comparison of Three Sediment Acoustic Models Using Bayesian Inversion and Model Selection Techniques

2019 ◽  
Vol 11 (5) ◽  
pp. 562 ◽  
Author(s):  
Bo Zou ◽  
Jingsheng Zhai ◽  
Zhanfeng Qi ◽  
Zhaoxing Li

Many geoacoustic models are used to establish the relationship between the physical and acoustic properties of sediments. In this work, Bayesian inversion and model selection techniques are applied to compare combinations of three geoacoustic models and corresponding scattering models—the fluid model with the effective density fluid model (EDFM), the grain-shearing elastic model with the viscosity grain-shearing (VGS(λ)) model, and the poroelastic model with the corrected and reparametrized extended Biot–Stoll (CREB) model. First, the resolution and correlation of parameters for the three models are compared based on estimates of the posterior probability distributions (PPDs), which are obtained by Bayesian inversion using the backscattering strength data. Then, model comparison and selection techniques are utilized to assess the matching degree of model predictions and measurements qualitatively and to ascertain the Bayes factors in favor of each quantitatively. These studies indicate that the fluid and poroelastic models outperform the grain-shearing elastic model, in terms of both parameter resolution and the ability to produce predictions in agreement with measurements for sandy sediments. The poroelastic model is considered to be the best, as the inversion based on it can provide more highly resolved information of sandy sediments. Finally, the attempt to implement geoacoustic inversion with different models provides a relatively feasible remote sensing scheme for various types of sediments under unknown conditions, which needs further validation.

Author(s):  
Daniel Blatter ◽  
Anandaroop Ray ◽  
Kerry Key

Summary Bayesian inversion of electromagnetic data produces crucial uncertainty information on inferred subsurface resistivity. Due to their high computational cost, however, Bayesian inverse methods have largely been restricted to computationally expedient 1D resistivity models. In this study, we successfully demonstrate, for the first time, a fully 2D, trans-dimensional Bayesian inversion of magnetotelluric data. We render this problem tractable from a computational standpoint by using a stochastic interpolation algorithm known as a Gaussian process to achieve a parsimonious parametrization of the model vis-a-vis the dense parameter grids used in numerical forward modeling codes. The Gaussian process links a trans-dimensional, parallel tempered Markov chain Monte Carlo sampler, which explores the parsimonious model space, to MARE2DEM, an adaptive finite element forward solver. MARE2DEM computes the model response using a dense parameter mesh with resistivity assigned via the Gaussian process model. We demonstrate the new trans-dimensional Gaussian process sampler by inverting both synthetic and field magnetotelluric data for 2D models of electrical resistivity, with the field data example converging within 10 days on 148 cores, a non-negligible but tractable computational cost. For a field data inversion, our algorithm achieves a parameter reduction of over 32x compared to the fixed parameter grid used for the MARE2DEM regularized inversion. Resistivity probability distributions computed from the ensemble of models produced by the inversion yield credible intervals and interquartile plots that quantitatively show the non-linear 2D uncertainty in model structure. This uncertainty could then be propagated to other physical properties that impact resistivity including bulk composition, porosity and pore-fluid content.


2019 ◽  
pp. 004912411988245 ◽  
Author(s):  
J. Mulder ◽  
A. E. Raftery

The Schwarz or Bayesian information criterion (BIC) is one of the most widely used tools for model comparison in social science research. The BIC, however, is not suitable for evaluating models with order constraints on the parameters of interest. This article explores two extensions of the BIC for evaluating order-constrained models, one where a truncated unit information prior is used under the order-constrained model and the other where a truncated local unit information prior is used. The first prior is centered on the maximum likelihood estimate, and the latter prior is centered on a null value. Several analyses show that the order-constrained BIC based on the local unit information prior works better as an Occam’s razor for evaluating order-constrained models and results in lower error probabilities. The methodology based on the local unit information prior is implemented in the R package “BICpack” which allows researchers to easily apply the method for order-constrained model selection. The usefulness of the methodology is illustrated using data from the European Values Study.


2017 ◽  
Vol 141 (1) ◽  
pp. EL32-EL37 ◽  
Author(s):  
Guangying Zheng ◽  
Yiwang Huang ◽  
Jian Hua ◽  
Xiuyu Xu ◽  
Fei Wang

2021 ◽  
Author(s):  
Philipp Kaniuth ◽  
Martin N. Hebart

AbstractRepresentational Similarity Analysis (RSA) has emerged as a popular method for relating representational spaces from human brain activity, behavioral data, and computational models. RSA is based on the comparison of representational dissimilarity matrices (RDM), which characterize the pairwise dissimilarities of all conditions across all features (e.g. fMRI voxels or units of a model). However, classical RSA treats each feature as equally important. This ‘equal weights’ assumption contrasts with the flexibility of multivariate decoding, which reweights individual features for predicting a target variable. As a consequence, classical RSA may lead researchers to underestimate the correspondence between a model and a brain region and, for model comparison, it may lead to selecting the inferior model. While previous work has suggested that reweighting can improve model selection in RSA, it has remained unclear to what extent these results generalize across datasets and data modalities. To fill this gap, the aim of this work is twofold: First, utilizing a range of publicly available datasets and three popular deep neural networks (DNNs), we seek to broadly test feature-reweighted RSA (FR-RSA) applied to computational models and reveal the extent to which reweighting model features improves RDM correspondence and affects model selection. Second, we propose voxel-reweighted RSA, a novel use case of FR-RSA that reweights fMRI voxels, mirroring the rationale of multivariate decoding of optimally combining voxel activity patterns. We find that reweighting individual model units (1) markedly improves the fit between model RDMs and target RDMs derived from several fMRI and behavioral datasets and (2) affects model selection, highlighting the importance of considering FR-RSA. For voxel-reweighted RSA, improvements in RDM correspondence were even more pronounced, demonstrating the utility of this novel approach. We additionally demonstrate that classical noise ceilings can be exceeded when FR-RSA is applied and propose an updated approach for their computation. Taken together, our results broadly validate the use of FR-RSA for improving the fit between computational models, brain and behavioral data, possibly allowing us to better adjudicate between competing computational models. Further, our results suggest that FR-RSA applied to brain measurement channels could become an important new method to assess the match between representational spaces.


2015 ◽  
Vol 137 (4) ◽  
pp. 2391-2391 ◽  
Author(s):  
Anthony L. Bonomo ◽  
Nicholas P. Chotiros ◽  
Marcia J. Isakson

Sign in / Sign up

Export Citation Format

Share Document