Facilitating analysis of Monte Carlo dense matrix inversion algorithm scaling behaviour through simulation

2013 ◽  
Vol 4 (6) ◽  
pp. 473-479 ◽  
Author(s):  
Janko Straßburg ◽  
Vassil N. Alexandrov
Geophysics ◽  
2015 ◽  
Vol 80 (1) ◽  
pp. R31-R41 ◽  
Author(s):  
Andrea Zunino ◽  
Klaus Mosegaard ◽  
Katrine Lange ◽  
Yulia Melnikova ◽  
Thomas Mejer Hansen

Determination of a petroleum reservoir structure and rock bulk properties relies extensively on inference from reflection seismology. However, classic deterministic methods to invert seismic data for reservoir properties suffer from some limitations, among which are the difficulty of handling complex, possibly nonlinear forward models, and the lack of robust uncertainty estimations. To overcome these limitations, we studied a methodology to invert seismic reflection data in the framework of the probabilistic approach to inverse problems, using a Markov chain Monte Carlo (McMC) algorithm with the goal to directly infer the rock facies and porosity of a target reservoir zone. We thus combined a rock-physics model with seismic data in a single inversion algorithm. For large data sets, the McMC method may become computationally impractical, so we relied on multiple-point-based a priori information to quantify geologically plausible models. We tested this methodology on a synthetic reservoir model. The solution of the inverse problem was then represented by a collection of facies and porosity reservoir models, which were samples of the posterior distribution. The final product included probability maps of the reservoir properties in obtained by performing statistical analysis on the collection of solutions.


1950 ◽  
Vol 4 (31) ◽  
pp. 127-127 ◽  
Author(s):  
George E. Forsythe ◽  
Richard A. Leibler

Author(s):  
Andre Bannwart Perina ◽  
Paulo Matias ◽  
Eduardo Marques ◽  
Vanderlei Bonato ◽  
Joao Miguel Gago Pontes de Brito Lima

1995 ◽  
Vol 06 (01) ◽  
pp. 25-45
Author(s):  
STEFANO ANTONELLI ◽  
MARCO BELLACCI ◽  
ANDREA DONINI ◽  
RENATA SARNO

We present the first tests and results from a study of QCD with two flavours of dynamical Wilson fermions using the Hybrid Monte Carlo Algorithm (HMCA) on APE100 machines. The simulations have been performed on 64 lattice for the pure gauge HMCA and on 84, 123×32 lattices for full QCD configurations. We discuss the inversion algorithm for the fermionic operator, the methods used to overcome the problems arising using a 32 bit machine and the implementation of a new random number generator for APE100 machines. We propose different scenarios for the simulation of physical observables, with respect to the memory capacity and speed of different APE100 configurations.


Sign in / Sign up

Export Citation Format

Share Document