Reservoir-parameter identification using minimum relative entropy-based Bayesian inversion of seismic AVA and marine CSEM data

Geophysics ◽  
2006 ◽  
Vol 71 (6) ◽  
pp. O77-O88 ◽  
Author(s):  
Zhangshuan Hou ◽  
Yoram Rubin ◽  
G. Michael Hoversten ◽  
Don Vasco ◽  
Jinsong Chen

A stochastic joint-inversion approach for estimating reservoir-fluid saturations and porosity is proposed. The approach couples seismic amplitude variation with angle (AVA) and marine controlled-source electromagnetic (CSEM) forward models into a Bayesian framework, which allows for integration of complementary information. To obtain minimally subjective prior probabilities required for the Bayesian approach, the principle of minimum relative entropy (MRE) is employed. Instead of single-value estimates provided by deterministic methods, the approach gives a probability distribution for any unknown parameter of interest, such as reservoir-fluid saturations or porosity at various locations. The distribution means, modes, and confidence intervals can be calculated, providing a more complete understanding of the uncertainty in the parameter estimates. The approach is demonstrated using synthetic and field data sets. Results show that joint inversion using seismic and EM data gives better estimates of reservoir parameters than estimates from either geophysical data set used in isolation. Moreover, a more informative prior leads to much narrower predictive intervals of the target parameters, with mean values of the posterior distributions closer to logged values.

2015 ◽  
Vol 2015 ◽  
pp. 1-12
Author(s):  
Mohammed Alguraibawi ◽  
Habshah Midi ◽  
A. H. M. Rahmatullah Imon

Identification of high leverage point is crucial because it is responsible for inaccurate prediction and invalid inferential statement as it has a larger impact on the computed values of various estimates. It is essential to classify the high leverage points into good and bad leverage points because only the bad leverage points have an undue effect on the parameter estimates. It is now evident that when a group of high leverage points is present in a data set, the existing robust diagnostic plot fails to classify them correctly. This problem is due to the masking and swamping effects. In this paper, we propose a new robust diagnostic plot to correctly classify the good and bad leverage points by reducing both masking and swamping effects. The formulation of the proposed plot is based on the Modified Generalized Studentized Residuals. We investigate the performance of our proposed method by employing a Monte Carlo simulation study and some well-known data sets. The results indicate that the proposed method is able to improve the rate of detection of bad leverage points and also to reduce swamping and masking effects.


2021 ◽  
Author(s):  
Mohammad Shehata ◽  
Hideki Mizunaga

<p>Long-period magnetotelluric and gravity data were acquired to investigate the US cordillera's crustal structure. The magnetotelluric data are being acquired across the continental USA on a quasi-regular grid of ∼70 km spacing as an electromagnetic component of the National Science Foundation EarthScope/USArray Program. International Gravimetreique Bureau compiled gravity Data at high spatial resolution. Due to the difference in data coverage density, the geostatistical joint integration was utilized to map the subsurface structures with adequate resolution. First, a three-dimensional inversion of each data set was applied separately.</p><p>The inversion results of both data sets show a similarity of structure for data structuralizing. The individual result of both data sets is resampled at the same locations using the kriging method by considering each inversion model to estimate the coefficient. Then, the Layer Density Correction (LDC) process's enhanced density distribution was applied to MT data's spatial expansion process. Simple Kriging with varying Local Means (SKLM) was applied to the residual analysis and integration. For this purpose, the varying local means of the resistivity were estimated using the corrected gravity data by the Non-Linear Indicator Transform (NLIT), taking into account the spatial correlation. After that, the spatial expansion analysis of MT data obtained sparsely was attempted using the estimated local mean values and SKLM method at the sections where the MT survey was carried out and for the entire area where density distributions exist. This research presents the integration results and the stand-alone inversion results of three-dimensional gravity and magnetotelluric data.</p>


2006 ◽  
Vol 30 (1) ◽  
pp. 83 ◽  
Author(s):  
Ronald Donato ◽  
Jeffrey Richardson

Diagnosis-based risk adjustment is increasingly seen as an important tool for establishing capitation payments and evaluating appropriateness and efficiency of services provided and has become an important area of research for many countries contemplating health system reform. This paper examines the application of a risk-adjustment method, extensively validated in the United States, known as diagnostic cost groups (DCG), to a large Australian hospital inpatient data set. The data set encompassed hospital inpatient diagnoses and inpatient expenditure for the entire metropolitan population residing in the state of New South Wales. The DCG model was able to explain 34% of individual-level variation in concurrent expenditure and 5.2% in subsequent year expenditure, which is comparable to US studies using inpatient-only data. The degree of stability and internal consistency of the parameter estimates for both the concurrent and prospective models indicate the DCG methodology has face validity in its application to NSW health data sets. Modelling and simulations were conducted which demonstrate the policy applications and significance of risk adjustment model(s) in the Australian context. This study demonstrates the feasibility of using large individual-level data sets for diagnosis-based risk adjustment research in Australia. The results suggest that a research agenda should be established to broaden the options for health system reform.


2001 ◽  
Vol 73 (2) ◽  
pp. 229-240 ◽  
Author(s):  
H. N. Kadarmideen ◽  
R. Rekaya ◽  
D. Gianola

AbstractA Bayesian threshold-liability model with Markov chain Monte Carlo techniques was used to infer genetic parameters for clinical mastitis records collected on Holstein-Friesian cows by one of the United Kingdom’s national recording schemes. Four data sets were created to investigate the effect of data sampling methods on genetic parameter estimates for first and multi-lactation cows, separately. The data sets were: (1) cows with complete first lactations only (8671 cows); (2) all cows, with first lactations whether complete or incomplete (10 967 cows); (3) cows with complete multi-lactations (32 948 records); and (4) all cows with multiple lactations whether complete or incomplete (44 268 records). A Gaussian mixed linear model with sire effects was adopted for liability. Explanatory variables included in the model varied for each data set. Analyses were conducted using Gibbs sampling and estimates were on the liability scale. Posterior means of heritability for clinical mastitis were higher for first lactations (0·11 and 0·10 for data sets 1 and 2, respectively) than for multiple lactations (0·09 and 0·07, for data sets 3 and 4, respectively). For multiple lactations, estimates of permanent environmental variance were higher for complete than incomplete lactations. Repeatability was 0·21 and 0·17 for data sets 3 and 4, respectively. This suggests the existence of effects, other than additive genetic effects, on susceptibility to mastitis that are common to all lactations. In first or multi-lactation data sets, heritability was proportionately 0·10 to 0·19 lower for data sets with all records (in which case the models had days in milk as a covariate) than for data with only complete lactation records (models without days in milk as a covariate). This suggests an effect of data sampling on genetic parameter estimates. The regression of liability on days in milk differed from zero, indicating that the probability of mastitis is higher for longer lactations, as expected. Results also indicated that a regression on days in milk should be included in a model for genetic evaluation of sires for mastitis resistance based on records in progress.


2017 ◽  
Vol 3 (1) ◽  
Author(s):  
Dora Matzke ◽  
Alexander Ly ◽  
Ravi Selker ◽  
Wouter D. Weeda ◽  
Benjamin Scheibehenne ◽  
...  

Whenever parameter estimates are uncertain or observations are contaminated by measurement error, the Pearson correlation coefficient can severely underestimate the true strength of an association. Various approaches exist for inferring the correlation in the presence of estimation uncertainty and measurement error, but none are routinely applied in psychological research. Here we focus on a Bayesian hierarchical model proposed by Behseta, Berdyyeva, Olson, and Kass (2009) that allows researchers to infer the underlying correlation between error-contaminated observations. We show that this approach may be also applied to obtain the underlying correlation between uncertain parameter estimates as well as the correlation between uncertain parameter estimates and noisy observations. We illustrate the Bayesian modeling of correlations with two empirical data sets; in each data set, we first infer the posterior distribution of the underlying correlation and then compute Bayes factors to quantify the evidence that the data provide for the presence of an association.


2020 ◽  
Vol 36 (1) ◽  
pp. 89-115 ◽  
Author(s):  
Harvey Goldstein ◽  
Natalie Shlomo

AbstractThe requirement to anonymise data sets that are to be released for secondary analysis should be balanced by the need to allow their analysis to provide efficient and consistent parameter estimates. The proposal in this article is to integrate the process of anonymisation and data analysis. The first stage uses the addition of random noise with known distributional properties to some or all variables in a released (already pseudonymised) data set, in which the values of some identifying and sensitive variables for data subjects of interest are also available to an external ‘attacker’ who wishes to identify those data subjects in order to interrogate their records in the data set. The second stage of the analysis consists of specifying the model of interest so that parameter estimation accounts for the added noise. Where the characteristics of the noise are made available to the analyst by the data provider, we propose a new method that allows a valid analysis. This is formally a measurement error model and we describe a Bayesian MCMC algorithm that recovers consistent estimates of the true model parameters. A new method for handling categorical data is presented. The article shows how an appropriate noise distribution can be determined.


Geophysics ◽  
2021 ◽  
pp. 1-93
Author(s):  
Joseph Capriotti ◽  
Yaoguo Li

Gravity and gravity gradiometry measurements are commonly used to map density variations in the subsurface. Gravity measurements can characterize gravitational anomalies at both long and short wavelengths effectively, but the cost of collecting a sufficiently spatially dense survey to characterize the short wavelengths can be prohibitive. Gravity gradient data can be quickly collected with short wavelength information at a low noise level, but have decreasing sensitivity to longer wavelengths. We describe a method to jointly invert gravity and gravity gradient data that takes advantage of the differing frequency contents and noise levels of the two methods to create an improved image of the subsurface. Previous work simply treated the inversion as a multiple component gravity inversion, however this can cause unintended errors in the recovered models because each data set is not guaranteed to be fit within its noise level. Our joint inversion methodology ensures that both the gravity and gravity gradient data sets are fit to within their individual noise levels by incorporating a relative weighting parameter, and we describe how to find that parameter. This method can also be used to create an improved broadband gravity anomaly map that has a reduced noise level at long wavelengths using a joint equivalent source reconstruction. We first build a synthetic model using a Minecraft world editor, that has different wavelength anomalies, and show the improvement with joint inversion. These results are also confirmed using a real world example at the R. J. Smith test range in Kauring, Australia.


Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. C1-C13 ◽  
Author(s):  
G. Michael Hoversten ◽  
Florence Cassassuce ◽  
Erika Gasperikova ◽  
Gregory A. Newman ◽  
Jinsong Chen ◽  
...  

Accurately estimating reservoir parameters from geophysical data is vitally important in hydrocarbon exploration and production. We have developed a new joint-inversion algorithm to estimate reservoir parameters directly, using both seismic amplitude variation with angle of incidence (AVA) data and marine controlled-source electromagnetic (CSEM) data. Reservoir parameters are linked to geophysical parameters through a rock-properties model. Errors in the parameters of the rock-properties model introduce errors of comparable size in the reservoir-parameter estimates produced by joint inversion. Tests of joint inversion on synthetic 1D models demonstrate improved fluid saturation and porosity estimates for joint AVA-CSEM data inversion (compared with estimates from AVA or CSEM inversion alone). A comparison of inversions of AVA data, CSEM data, and joint AVA-CSEM data over the North Sea Troll field, at a location for which we have well control, shows that the joint inversion produces estimates of gas saturation, oil saturation, and porosity that are closest (as measured by the rms difference, the [Formula: see text] norm of the difference, and net values over the interval) to the logged values. However, CSEM-only inversion provides the closest estimates of water saturation.


Radiocarbon ◽  
2010 ◽  
Vol 52 (3) ◽  
pp. 901-906 ◽  
Author(s):  
Taiichi Sato ◽  
Hirohisa Sakurai ◽  
Kayo Suzuki ◽  
Yui Takahashi

Radiocarbon ages of single-year tree rings were measured for Kaminoyama wood samples using accelerator mass spectrometry (AMS) in 2 Japanese facilities, MALT and JAEA, in order to investigate the periodic variation of 14C concentrations relating to the 11-yr solar cycle near 26,000 yr BP. Eight sequential measurements of 14C ages were carried out for a set of 13 alternate single-year tree rings covering 26 yr. Averages of the 5 data sets in MALT and 3 data sets in JAEA were 22,146 ± 50 and 22,407 ± 58 14C yr BP, respectively, indicating an offset of 260 ± 77 14C yr. Multiple sequential measurements are advantageous for evaluating offsets. The standard deviation of the residuals of 14C ages from the averages in each data set was 118 14C yr, in contrast to that of 234 14C yr for the combined data sets due to an elimination effect in the offsets. The profiles of weighted mean values for the residuals of 14C ages showed similar enhancements with a width of ∼12 yr for measurements in the 2 facilities. This indicates the reproducibility of the multiple sequential measurements. In the profile for the combined 8 data sets, the 14C enhancement was 73 ± 36 14C yr from the average.


2021 ◽  
Vol 13 (5) ◽  
pp. 2165-2209
Author(s):  
Pavol Zahorec ◽  
Juraj Papčo ◽  
Roman Pašteka ◽  
Miroslav Bielik ◽  
Sylvain Bonvalot ◽  
...  

Abstract. The AlpArray Gravity Research Group (AAGRG), as part of the European AlpArray program, focuses on the compilation of a homogeneous surface-based gravity data set across the Alpine area. In 2017 10 European countries in the Alpine realm agreed to contribute with gravity data for a new compilation of the Alpine gravity field in an area spanning from 2 to 23∘ E and from 41 to 51∘ N. This compilation relies on existing national gravity databases and, for the Ligurian and the Adriatic seas, on shipborne data of the Service Hydrographique et Océanographique de la Marine and of the Bureau Gravimétrique International. Furthermore, for the Ivrea zone in the Western Alps, recently acquired data were added to the database. This first pan-Alpine gravity data map is homogeneous regarding input data sets, applied methods and all corrections, as well as reference frames. Here, the AAGRG presents the data set of the recalculated gravity fields on a 4 km × 4 km grid for public release and a 2 km × 2 km grid for special request. The final products also include calculated values for mass and bathymetry corrections of the measured gravity at each grid point, as well as height. This allows users to use later customized densities for their own calculations of mass corrections. Correction densities used are 2670 kg m−3 for landmasses, 1030 kg m−3 for water masses above the ellipsoid and −1640 kg m−3 for those below the ellipsoid and 1000 kg m−3 for lake water masses. The correction radius was set to the Hayford zone O2 (167 km). The new Bouguer anomaly is station completed (CBA) and compiled according to the most modern criteria and reference frames (both positioning and gravity), including atmospheric corrections. Special emphasis was put on the gravity effect of the numerous lakes in the study area, which can have an effect of up to 5 mGal for gravity stations located at shorelines with steep slopes, e.g., for the rather deep reservoirs in the Alps. The results of an error statistic based on cross validations and/or “interpolation residuals” are provided for the entire database. As an example, the interpolation residuals of the Austrian data set range between about −8 and +8 mGal and the cross-validation residuals between −14 and +10 mGal; standard deviations are well below 1 mGal. The accuracy of the newly compiled gravity database is close to ±5 mGal for most areas. A first interpretation of the new map shows that the resolution of the gravity anomalies is suited for applications ranging from intra-crustal- to crustal-scale modeling to interdisciplinary studies on the regional and continental scales, as well as applications as joint inversion with other data sets. The data are published with the DOI https://doi.org/10.5880/fidgeo.2020.045 (Zahorec et al., 2021) via GFZ Data Services.


Sign in / Sign up

Export Citation Format

Share Document