Geostatistical Seismic Inversion with Self-Updating of Local Probability Distributions

Author(s):  
Leonardo Azevedo ◽  
João Narciso ◽  
Rúben Nunes ◽  
Amílcar Soares
Energies ◽  
2022 ◽  
Vol 15 (1) ◽  
pp. 299
Author(s):  
Zhihong Wang ◽  
Tiansheng Chen ◽  
Xun Hu ◽  
Lixin Wang ◽  
Yanshu Yin

In order to solve the problem that elastic parameter constraints are not taken into account in local lithofacies updating in multi-point geostatistical inversion, a new multi-point geostatistical inversion method with local facies updating under seismic elastic constraints is proposed. The main improvement of the method is that the probability of multi-point facies modeling is combined with the facies probability reflected by the optimal elastic parameters retained from the previous inversion to predict and update the current lithofacies model. Constrained by the current lithofacies model, the elastic parameters were obtained via direct sampling based on the statistical relationship between the lithofacies and the elastic parameters. Forward simulation records were generated via convolution and were compared with the actual seismic records to obtain the optimal lithofacies and elastic parameters. The inversion method adopts the internal and external double cycle iteration mechanism, and the internal cycle updates and inverts the local lithofacies. The outer cycle determines whether the correlation between the entire seismic record and the actual seismic record meets the given conditions, and the cycle iterates until the given conditions are met in order to achieve seismic inversion prediction. The theoretical model of the Stanford Center for Reservoir Forecasting and the practical model of the Xinchang gas field in western China were used to test the new method. The results show that the correlation between the synthetic seismic records and the actual seismic records is the best, and the lithofacies matching degree of the inversion is the highest. The results of the conventional multi-point geostatistical inversion are the next best, and the results of the two-point geostatistical inversion are the worst. The results show that the reservoir parameters obtained using the local probability updating of lithofacies method are closer to the actual reservoir parameters. This method is worth popularizing in practical exploration and development.


Geophysics ◽  
2001 ◽  
Vol 66 (4) ◽  
pp. 988-1001 ◽  
Author(s):  
T. Mukerji ◽  
A. Jørstad ◽  
P. Avseth ◽  
G. Mavko ◽  
J. R. Granli

Reliably predicting lithologic and saturation heterogeneities is one of the key problems in reservoir characterization. In this study, we show how statistical rock physics techniques combined with seismic information can be used to classify reservoir lithologies and pore fluids. One of the innovations was to use a seismic impedance attribute (related to the [Formula: see text] ratio) that incorporates far‐offset data, but at the same time can be practically obtained using normal incidence inversion algorithms. The methods were applied to a North Sea turbidite system. We incorporated well log measurements with calibration from core data to estimate the near‐offset and far‐offset reflectivity and impedance attributes. Multivariate probability distributions were estimated from the data to identify the attribute clusters and their separability for different facies and fluid saturations. A training data was set up using Monte Carlo simulations based on the well log—derived probability distributions. Fluid substitution by Gassmann’s equation was used to extend the training data, thus accounting for pore fluid conditions not encountered in the well. Seismic inversion of near‐offset and far‐offset stacks gave us two 3‐D cubes of impedance attributes in the interwell region. The near‐offset stack approximates a zero‐offset section, giving an estimate of the normal incidence acoustic impedance. The far offset stack gives an estimate of a [Formula: see text]‐related elastic impedance attribute that is equivalent to the acoustic impedance for non‐normal incidence. These impedance attributes obtained from seismic inversion were then used with the training probability distribution functions to predict the probability of occurrence of the different lithofacies in the interwell region. Statistical classification techniques, as well as geostatistical indicator simulations were applied on the 3‐D seismic data cube. A Markov‐Bayes technique was used to update the probabilities obtained from the seismic data by taking into account the spatial correlation as estimated from the facies indicator variograms. The final results are spatial 3‐D maps of not only the most likely facies and pore fluids, but also their occurrence probabilities. A key ingredient in this study was the exploitation of physically based seismic‐to‐reservoir property transforms optimally combined with statistical techniques.


2017 ◽  
Vol 5 (4) ◽  
pp. T641-T652 ◽  
Author(s):  
Mark Sams ◽  
Paul Begg ◽  
Timur Manapov

The information within seismic data is band limited and angle limited. Together with the particular physics and geology of carbonate rocks, this imposes limitations on how accurately we can predict the presence of hydrocarbons in carbonates, map the top carbonate, and characterize the porosity distribution through seismic amplitude analysis. Using data for a carbonate reef from the Nam Con Son Basin, Vietnam, the expectations based on rock-physics analysis are that the presence of gas can be predicted only when the porosity at the top of the carbonate is extremely high ([Formula: see text]), but that a fluid contact is unlikely to be observed in the background of significant porosity variations. Mapping the top of the carbonate (except when the top carbonate porosities are low) or a fluid contact requires accurate estimates of changes in [Formula: see text]. The seismic data do not independently support such an accurate estimation of sharp changes in [Formula: see text]. The standard approach of introducing low-frequency models and applying rock-physics constraints during a simultaneous inversion does not resolve the problems: The results are heavily biased by the well control and the initial interpretation of the top carbonate and fluid contact. A facies-based inversion in which the elastic properties are restricted to values consistent with the facies predicted to be present removes the well bias, but it does not completely obviate the need for a reasonably accurate initial interpretation in terms of prior facies probability distributions. Prestack inversion improves the quality of the facies predictions compared with a poststack inversion.


2013 ◽  
Vol 12 (08) ◽  
pp. 1341014 ◽  
Author(s):  
ROBERT J. PETRELLA

Physics-based computational approaches to predicting the structure of macromolecules such as proteins are gaining increased use, but there are remaining challenges. In the current work, it is demonstrated that in energy-based prediction methods, the degree of optimization of the sampled structures can influence the prediction results. In particular, discrepancies in the degree of local sampling can bias the predictions in favor of the oversampled structures by shifting the local probability distributions of the minimum sampled energies. In simple systems, it is shown that the magnitude of the errors can be calculated from the energy surface, and for certain model systems, derived analytically. Further, it is shown that for energy wells whose forms differ only by a randomly assigned energy shift, the optimal accuracy of prediction is achieved when the sampling around each structure is equal. Energy correction terms can be used in cases of unequal sampling to reproduce the total probabilities that would occur under equal sampling, but optimal corrections only partially restore the prediction accuracy lost to unequal sampling. For multiwell systems, the determination of the correction terms is a multibody problem; it is shown that the involved cross-correlation multiple integrals can be reduced to simpler integrals. The possible implications of the current analysis for macromolecular structure prediction are discussed.


Author(s):  
Takuto Omiya ◽  
◽  
Kazuhiro Hotta

In this paper, we perform image labeling based on the probabilistic integration of local and global features. Several conventional methods label pixels or regions using features extracted from local regions and local contextual relationships between neighboring regions. However, labeling results tend to depend on local viewpoints. To overcome this problem, we propose an image labeling method that utilizes both local and global features. We compute the posterior probability distributions of the local and global features independently, and they are integrated by the product. To compute the probability of the global region (entire image), Bag-of-Words is used. In contrast, local cooccurrence between color and texture features is used to compute local probability. In the experiments, we use the MSRC21 dataset. The result demonstrates that the use of global viewpoint significantly improves labeling accuracy.


1997 ◽  
Vol 3 (6) ◽  
pp. 553-568
Author(s):  
L. Schueremans ◽  
D. Van Gemert ◽  
J. Van Dyck

Abstract A probabilistic method to evaluate the reliability of structural masonry elements is presented. Local probability of failure, different failure modes and corresponding limit state functions, probability distributions of basic variables and model uncertainty are discussed. A graphical probability mapping is presented as an easilyaccessible, visual evaluation instrument in the restoration decision process. The proposed methodology is illustrated on tested shear panels, reported in literature.


1997 ◽  
Vol 161 ◽  
pp. 197-201 ◽  
Author(s):  
Duncan Steel

AbstractWhilst lithopanspermia depends upon massive impacts occurring at a speed above some limit, the intact delivery of organic chemicals or other volatiles to a planet requires the impact speed to be below some other limit such that a significant fraction of that material escapes destruction. Thus the two opposite ends of the impact speed distributions are the regions of interest in the bioastronomical context, whereas much modelling work on impacts delivers, or makes use of, only the mean speed. Here the probability distributions of impact speeds upon Mars are calculated for (i) the orbital distribution of known asteroids; and (ii) the expected distribution of near-parabolic cometary orbits. It is found that cometary impacts are far more likely to eject rocks from Mars (over 99 percent of the cometary impacts are at speeds above 20 km/sec, but at most 5 percent of the asteroidal impacts); paradoxically, the objects impacting at speeds low enough to make organic/volatile survival possible (the asteroids) are those which are depleted in such species.


2020 ◽  
Vol 3 (1) ◽  
pp. 10501-1-10501-9
Author(s):  
Christopher W. Tyler

Abstract For the visual world in which we operate, the core issue is to conceptualize how its three-dimensional structure is encoded through the neural computation of multiple depth cues and their integration to a unitary depth structure. One approach to this issue is the full Bayesian model of scene understanding, but this is shown to require selection from the implausibly large number of possible scenes. An alternative approach is to propagate the implied depth structure solution for the scene through the “belief propagation” algorithm on general probability distributions. However, a more efficient model of local slant propagation is developed as an alternative.The overall depth percept must be derived from the combination of all available depth cues, but a simple linear summation rule across, say, a dozen different depth cues, would massively overestimate the perceived depth in the scene in cases where each cue alone provides a close-to-veridical depth estimate. On the other hand, a Bayesian averaging or “modified weak fusion” model for depth cue combination does not provide for the observed enhancement of perceived depth from weak depth cues. Thus, the current models do not account for the empirical properties of perceived depth from multiple depth cues.The present analysis shows that these problems can be addressed by an asymptotic, or hyperbolic Minkowski, approach to cue combination. With appropriate parameters, this first-order rule gives strong summation for a few depth cues, but the effect of an increasing number of cues beyond that remains too weak to account for the available degree of perceived depth magnitude. Finally, an accelerated asymptotic rule is proposed to match the empirical strength of perceived depth as measured, with appropriate behavior for any number of depth cues.


Sign in / Sign up

Export Citation Format

Share Document