# Nonlinear Bayesian joint inversion of seismic reflection

Save this PDF as:

Size: px
Start display at page:

Download "Nonlinear Bayesian joint inversion of seismic reflection"

## Transcription

1 Geophys. J. Int. (2) 142, Nonlinear Bayesian joint inversion of seismic reflection coefficients Tor Erik Rabben 1,, Håkon Tjelmeland 2, and Bjørn Ursin 1 1 Department of Petroleum Engineering and Applied Geophysics, Norwegian University of Science and Technology, 7491 Trondheim, Norway 2 Department of Mathematical Sciences, Norwegian University of Science and Technology, 7491 Trondheim, Norway Accepted 1999 November 11. Received 1999 October 6; in original form 27 February 9 SUMMARY Inversion of the seismic reflection coefficients are formulated in a Bayesian framework. Measured reflection coefficients and model parameters are assigned statistical distributions based on information known prior to the inversion, and together with the forward model uncertainties are propagated into the final result. This enables a quantification of the reliability of the inversion. A quadratic approximation to the Zoeppritz equations is used as the forward model. Compared with the linear approximation the bias is reduced and the uncertainty estimates are more reliable. The differences when using the quadratic approximations and the exact expressions are minor. The solution algorithm is sampling based and because of the nonlinear forward model the Metropolis-Hastings algorithm is used. To achieve convergence it is important to keep strict control of the accept probability in the algorithm. Joint inversion using information from both reflected PP-waves and converted PS-waves yield smaller bias and reduced uncertainties compared to using only reflected PP-waves. Key words: Seismic inversion statistical methods

2 2 T.E. Rabben, H. Tjelmeland, and B. Ursin 1 INTRODUCTION The seismic reflection coefficients contain information about elastic parameters in the subsurface. In an amplitude versus angle (AVA) inversion, the main objective is to estimate elastic parameters from the reflection coefficients. Elastic parameters can in turn be used for sediment classification and extraction of fluid properties. Given a model it is easy to find an answer to the inverse problem, the real challenge is to find a reasonable answer (Hampson, 1991). Mathematically, the main problem is nonuniqueness. Analytical expressions for the reflection coefficients are given by the Zoeppritz equations. These equations are complicated and highly nonlinear. In the situation of reflections between two isotropic media the equations involves 5 parameters. Ursin & Tjåland (1996) showed that in practice only up to three parameters can be estimated from pre-critical PP reflection coefficients. To help overcome this, linear approximations to the Zoeppritz equations have been derived (e.g. Aki & Richards (198) and Smith & Gidlow (1987)), both for their simplicity and the favourable reduction from 5 to 3 parameters. One other possible solution is to include PS reflections in the inversion. Because pressure waves and shear waves sense different rock and pore-fluid properties, joint PP and PS data can provide superior lithologic discrimination (Margrave et al., 21). Using a simplified model and incorporating all available data might not be enough to overcome the nonuniqueness problem and produce a reasonable answer. In most cases regularisation is also necessary. Regularising the inverse problem means finding a physical meaningful stable solution (Tenorio, 21). A different approach to handle nonuniqueness is Bayesian formulation and uncertainty estimation. In the Bayesian statistics the model parameters and AVA data will no longer be treated as deterministic constants but will have statistical distributions assigned to them. The end result will not only consist of an estimate of the model parameters but also an estimate of the uncertainty in the parameters. It is very important to note that a Bayesian approach does not remove the nonuniqueness, it only helps identifying it.

3 Nonlinear Bayesian inversion 3 Previous work with inversion of reflection coefficients involves different models and solution algorithms. Smith & Gidlow (1987) used a linear approximation to the PP reflection coefficient and a least square approach to solve the inversion. Stewart (199) and Veire & Landrø (26) extended this to joint PP and PS inversion in a least square setting. Margrave et al. (21) gave a nice introduction to joint inversion and compared the results with only PP inversion. For a tutorial of least square inversion and how to regularise it see Lines & Treitel (1984) and for a more comprehensive treatment of regularisation Tenorio (21). Buland & Omre (23a; 23b) formulated the PP inversion using linear approximations, in the Bayesian framework. In the latter they not only estimate the elastic parameters but also the wavelet and noise level in the AVA data, together with uncertainties. For details on Bayesian modeling from a geophysics standpoint see Sen & Stoffa (1996) and from a statistical standpoint, Robert & Casella (1999) and Liu (21). Our approach to the inversion of reflection coefficients is to use an isotropic, quadratic approximation to the Zoeppritz equations instead of the linear ones used earlier. Stovas & Ursin (21; 23) gave implicit, second order expressions for reflection and transmission coefficients in both isotropic and transversely isotropic media. They showed that quadratic approximations are superior to the linear ones for intermediate pre-critical reflection angles, but the number of parameters are still three as for the linear approximations. However, for a comparison, we also perform the inversion using a linear approximation and the exact Zoeppritz equations. We formulate the inversion in a Bayesian framework following Buland & Omre (23b) and test both PP and joint PP and PS inversion. The use of statistical distributions enable us to impose spatial correlation and correlations between model parameters and reflection angles in a very natural way. The solution algorithm is based on sampling and because of the nonlinear model we use the Metropolis-Hastings algorithm (Robert & Casella, 1999; Liu, 21) following closely the work of Tjelmeland & Eidsvik (25). We define our computational domain to be a two dimensional surface, e.g. the top reservoir. The reason is to avoid the wavelet estimation and the convolution in the modeling and instead focus on the

4 4 T.E. Rabben, H. Tjelmeland, and B. Ursin nonlinearity. Our parameters to invert for are m, σe, 2 and σm. 2 The first one is the three elastic parameters in the Zoeppritz approximation. The two last parameters are data driven scalars and their main objective is to stabilise the inversion algorithm, but can to some extent quantify the noise level in the reflection coefficients. In the following sections we first give explicit expressions for quadratic approximations of the reflection coefficients. We then define the Bayesian model and describe the statistical inversion algorithm. In the numerical examples we define four synthetic test cases, show the inversion results, and conclude with a discussion of the results. 2 MODEL The parametrization we are using for the reflection coefficients is in P-wave and S-wave impedance and density, as suggested by Dȩbski & Tarantola (1995). Stovas & Ursin (23) derived implicit second order expressions for reflections between two transversely isotropic media. Explicit expressions for PP and PS-reflections simplified for two isotropic media, read 1 I α r PP = 2 cos 2 θ p I 4 sin 2 I β θ s 1 ( ) α Ī β 2 tan2 θ p 1 4γ 2 cos 2 ρ θ p ( ) 2 + tan θ p tan θ s [4γ 2 (1 (1 + γ 2 ) sin 2 Iβ θ p ) ( 4γ 2 (1 ( γ2 ) sin 2 Iβ θ p ) + Ī β ( γ 2 (1 (2 + γ 2 ) sin 2 θ p ) 1 4 Ī β ) ρ ) ] 2 ) ( ρ r PS = tan θ p tan θ s {[(1 cos θ s (cos θ s + γ cos θ p )) [ (1 cos θ s (cos θ s γ cos θ p )) [ ( 1 I α 2 cos 2 θ p I + α + ( 2 I β Ī β 1 2 cos 2 θ s 8 sin 2 θ s ( 4 sin 2 θ s 1 2 (tan2 θ p + tan 2 θ s ) ( 2 I β Ī β ρ ) 1 2 ) Iβ Ī β ) ]} ρ ρ ) 1 2 ] ρ ] ρ (1) (2)

5 Nonlinear Bayesian inversion 5 where γ = β/ᾱ is the background v S /v P -ratio, θ p is the angle of the incoming P-wave (and also the reflected P-wave because of isotropic medium) and θ s is the angle of the reflected S-wave. I α is P-wave impedance, I β is S-wave impedance and ρ is density. I α is difference between the lower and upper medium and I α denotes the average, equal definitions for I β and ρ. Appendix A shows how to derive the of reflection coefficients between two isotropic, elastic media for other parametrizations. In (1) and (2) the normalization is with respect to vertical energy flux. The relation to the more common amplitude normalized expressions is only a scalar factor. For r PP this constant is equal to one and the linear terms, after a changes of parameters, are equal to the amplitude normalized ones found in Aki & Richards (198). The variables to invert for are denoted m, σ 2 e, and σ 2 m. The first one is the obvious and consists of parameters found in (1) and (2), m = { Iα I α dimensional spatial grid with dimensions n y n x,, I β I β, ρ }. It is defined on a two m = {m ij R Dm ; i = 1..n y, j = 1..n x }. (3) D m is the number of unknown in each node and equal to three for the isotropic approximations (1) and (2). The two last parameters, σ 2 e and σ 2 m, are scalars quantifying variance levels. In our Bayesian formulation two covariance matrices have to be specified, and the choice made will influence the solution. We have therefore assumed these two matrices to be known only up to a multiplicative variance factor and let the estimation of the two factors be a part of the inversion procedure. In theory, the full covariance matrices could also be estimated this way but we have choosen to do this only for the variance levels, denoted σ 2 e for the data errors and σ 2 m for the parameter errors. The input to the inversion, the measured data, is denoted d. It is defined over the same spatial grid as for m, d = {d ij R D d ; i = 1..n y, j = 1..n x }. (4) D d is the sum measured PP and PS angles, d = {r PP (θ), r PS (θ)} in each node. It is important

6 6 T.E. Rabben, H. Tjelmeland, and B. Ursin that we have more measurements than parameters, i.e. D d D m, otherwise the inversion will be underdetermined. The forward model is the link from m to d. In the most exact case with full Zoeppritz it is written d = f z (m) + e z, (5) where f z is the Zoeppritz equations. The term e z is the observation errors related to the always present noise in measurements. In theory it will also consist of modelling errors but these are neglected under the assumption that the approximations in Zoeppritz are valid. Our focus will be on the quadratic approximations and the corresponding forward model is d = f(m) + e, (6) where f is (1) and (2) in the case of both PP and PS reflections. The error term can be approximated with e e z as long as the quadratic approximations are good. We can also formulate the linear forward model d = F m + e, (7) where F is a matrix containing only the linear part of (1) and (2). In general, e is not approximating e z as good as e does. A common way of inverting for m is to minimise e in (6) in an iterative nonlinear least squares algorithm, but this will not produce uncertainty measures. Therefore, we will reformulate the problem in a Bayesian framework. In such a formulation the parameters and measured data are no longer treated as being fixed, but have statistical distributions assigned to them. These distributions can be divided into three groups: prior, likelihood, and posterior. Prior distributions reflect the information known before the data is measured. This information can be theoretically founded, e.g. velocities can not be negative, or case dependent, e.g. information from the geological setting. The purpose is to restrict the search space to

7 Nonlinear Bayesian inversion 7 only reasonable values. However, it has also a down side, prior distributions can greatly influence the final solution. If the prior distribution is chosen to have very small variance, the solution will be close to it, regardless of the information in the measured data. The parameter σm 2 is introduced in the Bayesian model to help overcome this problem by adjusting the variance in the prior. In our formulation the prior distribution is denoted π(m, σe, 2 σm), 2 where π represent any distribution. Likelihood distributions are the second group. This is the analog to the forward model. In our Bayesian setting it is expressed π(d m, σe) 2 and describes the distribution of d given the parameters mand σ 2 e. The parameter σ 2 e is introduced to the likelihood for the same reason as σm 2 was introduced to the prior. The last distribution is the posterior and is the distribution we would like assess in the inversion procedure. It is the reverse of the likelihood, i.e. given d, what is the distribution of m, σe, 2 and σm, 2 and is expressed π(m, σe, 2 σm d). 2 Relating the posterior to the prior and likelihood is done by applying the Bayes rule to the posterior distribution, and the final answer is π(m, σe, 2 σm d) 2 π(d m, σe) 2 π(m, σe, 2 σm). 2 (8) The factor π(d) becomes an unknown constant, that we skip, and hence the proportionality in (8). It can be simplified by assuming some of the parameters to be independent. The details of the Bayesian formulation and choices of distributions are found in Appendix B. 3 INVERSION ALGORITHM Because of the nonlinearity in the likelihood, no analytical solution to the posterior is known. To approximate it we have to produce samples based on the prior and likelihood information, see equation (8), to build the posterior distribution. It is sampled with a Metropolis-Hastings algorithm (Robert & Casella, 1999; Liu, 21) which sequentially updates m, σe, 2 and σm, 2 keeping the two others fixed when updating. A Metropolis-Hastings algorithm consists of two steps:

8 8 T.E. Rabben, H. Tjelmeland, and B. Ursin (i) Propose a new sample (ii) Accept the sample with a probability p In the first step essentially any distribution can be used, the acceptance probability in the second step will correct for the errors introduced in the proposal distribution. If one proposes new values from the full conditionals (the distribution of parameter to propose given all other parameters) all proposals are accepted. There is, however, several reasons why this proposal distribution it not always used. It can be nonlinear and impossible to generate sample from, or very complicated or high dimensional and thereby very expensive to sample from. Back to our problem, we will sample the posterior by sequentially update m, σe, 2 and σm, 2 keeping the two others fixed. That is, instead of sampling π(m, σe, 2 σm d) 2 we will sample π(m d, σe, 2 σm), 2 π(σe d, 2 m, σm), 2 and π(σm d, 2 m, σe) 2 iteratively. After an update of all three we will have new sample from the posterior (8). In the process of sampling the three distributions we again use the Metropolis-Hastings algorithm. When sampling m we can not use the full conditional since the likelihood contains the nonlinear forward model. We therefore have to use a different distribution and the natural choice is a linearised posterior distribution involving the linear forward model (7) in the likelihood. This also results in an acceptance probability less than 1. For the two last posterior distributions the situation is different. Here we are able to produce samples from π(σe d, 2 m, σm) 2 and π(σm d, 2 m, σe). 2 The sampling of σe 2 and σm 2 are therefore examples of Gibbs updates. Details about the sampling procedure are found in Appendix C and the linear proposal distribution for m is derived in Appendix D. 4 NUMERICAL MODEL We will use a synthetic model to test the inversion algorithm. From a chosen true m we can use the exact Zoeppritz equations to generate synthetic measurements d. Since the isotropic Zoeppritz equations have 5 parameters we have to fix two variables in addition to the three in m. We have chosen the P-wave velocity in the upper medium and the background v P /v S

9 Nonlinear Bayesian inversion 9 ratio to be constant. A total of four cases will be generated, that is only PP data and both PP and PS data, each of them with and without noise added. For all these four cases we will compare inversions with three different models: linear, quadratic, and exact Zoeppritz. The truth we have chosen is, as discussed earlier, parametrized in contrasts in P-wave and S-wave impedance and density and they are all ranging from.2 to.5 as shown in Fig. 1. These contrasts are very strong and the purpose is to test the nonlinear inversion algorithm. The size of the computational grid is n x = n y = 1 with a spacing of 25m in each direction. In the Metropolis-Hastings algorithm the linear model is used to propose an update of the posterior, and it is accepted with a certain probability. To control this probability, and hence the acceptance rate, is very important. If it is equal to 1 all the proposals will be accepted and the result is linear inversion. On the contrary, if it is close to it will produce very few updates which leads to poor convergence and run time problems. The main reason for introducing the two scalars σe 2 and σm 2 and their distributions is to help overcome this problem. They are both data driven and will help control the acceptance rate and stabilise the inversion algorithm. To some extent σe 2 will also quantify the noise level, but it is always with respect to the forward model used. More specifically, it quantifies how well the model can reproduce d from the posterior m, but this does not necessary quantify the bias. The expected value of the prior of m in all cases is half of the true m in Fig. 1. The reason for this is to see if the inversion relies too heavily on prior information, if this is the case it will show as a large bias. Our first set of measurements is PP reflections without any noise. Fig. 2 shows the reflection coefficients and the corresponding bias in both the linear and quadratic approximations. We use four angles from to 55 and the reason for skipping the bias plots for θ = is that with P-wave impedance as one of the parameters both approximations are exact at normal incidence. For the three other angles we see that the quadratic approximations are better than the linear, but when critical angle is approached, in the two lower corners, even the quadratic will have large bias.

10 1 T.E. Rabben, H. Tjelmeland, and B. Ursin The second example we will test is PP inversion when noise is included. In the prior of m we add spatially correlated normal distributed noise to break the smoothness. The noise added to the measurements d is normally distributed and here it is both spatially correlated and correlated between the angles. Fig. 3 shows d when the noise is included. The variance is constant but is most visible for the normal incidence reflection where the range of values is smallest. After inverting PP reflections we will turn to joint PP and PS inversion. We will use the four coefficients in Fig. 2 and include three PS reflections with incoming P-waves between 2 and 55. The idea is to see what additional information the PS reflections will bring. Fig. 4 displays the exact Zoeppritz reflection coefficients together with the bias in the approximations, similar to the PP case. Here it is even more clear how superior the quadratic approximations are. The last synthetic example is joint inversion of reflection coefficients including noise. The prior used is the same as in the previous case with noise. Also for the measurements the situations is almost equal. The variance is still constant, the noise is spatially correlated and we have correlations between the angles, but the correlations between P-wave and S-wave reflections are set to be zero. Fig. 5 shows the measurements generated for joint inversion of PP and PS reflections including noise. 5 PP INVERSION The result after the inversion is the posterior distribution of m, σe, 2 and σm. 2 The two last ones are scalars and the distributions are simple to visualise, but this is not the case for m since it is multivariate. We therefore calculate the mean of the distribution and subtracts the truth to produce the bias. Fig. 6 shows the absolute value of this bias in the case without noise. Second, we have plotted the standard deviation of the distribution in Fig. 7. The two last scalar distributions are displayed i Fig. 8. Case number two that we tested was PP data including noise. The absolute value of the

11 Nonlinear Bayesian inversion 11 bias of the mean is displayed in Fig. 9, the standard deviation is in Fig. 1, and the scalar distributions are in Fig JOINT PP AND PS INVERSION The two last synthetic examples are joint inversion with and without noise included. As for the two PP cases we display the absolute value of the bias and the standard deviation for m and the full distributions for σe 2 and σm. 2 Figs 12 and 13 show the absolute value of the bias and the standard deviation in the case without noise and Fig. 14 the corresponding scalar distributions. Figs 15, 16, and 17 show the same plots for the last case, joint inversion of PP and PS data with noise added. 7 DISCUSSION In this section we will compare the results from all the four synthetic test cases. We will start by looking at m in the two cases without noise followed by the two with noise and conclude by looking at σe 2 and σm. 2 The absolute value of the bias in PP inversion and joint inversion, both without noise, found in Figs 6 and 12, have interesting common features and differences. Because of the parametrization in P-wave impedance we have in practice no bias for this parameter, independent of forward model and synthetic case example. For the two other parameters the bias is several orders of magnitude higher. In the PP case the performance of the three forward models are fairly equal for these two parameters, with the nonlinear ones slightly better as expected. In the joint inversion the situation is different. The linear inversion is actually worse, while the nonlinear ones have improved substantially. In the S-wave impedance we see some effects of too low acceptance rate in the areas to the right. Looking at Fig. 4 we see that this coincides with the areas where the linear model has large bias. Next, we keep these previous results in mind and compare the corresponding standard deviations in Figs 7 and 13. First of all we note that the additional PS information reduces

12 12 T.E. Rabben, H. Tjelmeland, and B. Ursin the standard deviation for all three parameters. Second, the standard deviation for the P- wave impedance are much smaller compared with the two other parameters, as for the bias. For the linear models the deviations are almost constant. This is because we have chosen space invariant noise and prior distributions and the linear likelihood. The reason for not being exactly constant is the sampling based algorithm we are using. In Fig. 7 the linear inversion shows a lower standard deviations for S-wave impedance than for the density. Comparing these two with the corresponding bias we see that this is a misleading result. The uncertainty in the nonlinear inversions yields a far more realistic result. This is an effect of using a linear forward model which is incorrect. The fit to the data is very good, but the parameter values are incorrect (they have large bias). Turning to the two last synthetic cases, the inversion of data with noise added, we first look at the bias in Figs 9 and 15. One important change is the scale of the P-wave impedance, the bias is now of the same magnitude as for the S-wave impedance and the density, but the parameter that is most difficult to resolve still seems to be the S-wave impedance. This effect is less pronounced in the joint inversion case. When comparing the three models we see that the differences are minor among them. The amount of noise added is such that the model used is less important. Looking at the standard deviation in Figs 1 and 16 we see the same feature as in the two first cases, the introduction of PS information greatly reduces the uncertainty for all three models and all three parameters. In addition we note that in this case with noise added, P-wave impedance has the lowest uncertainty while S-wave impedance has the highest. Last we will briefly comment on the two last scalar parameters, σe 2 and σm. 2 As mentioned earlier, the main reason for including these are to help stabilising the inversion. They are relative numbers and their value will depend on the prior distributions and forward models used. The reason for displaying them is for completeness since they are a part of the Bayesian model.

13 Nonlinear Bayesian inversion 13 8 CONCLUSION In our approach to inversion of reflection coefficients we have reformulated the problem in a Bayesian framework. The use of prior and likelihood distributions enables us to assess the uncertainty in the inversion result, the posterior. At the same time it enables us easily to impose correlations and covariances in the modeling. We have also used new quadratic approximations to the Zoeppritz equations and compared them with both the linearised and exact equations. Because of the nonlinear forward models we have tested the performance of the Metropolis-Hastings algorithm to sample form the nonlinear likelihood. Last, we have also included inversion of joint PP and PS and compared with the more common PP inversion. Our work shows that the Metropolis-Hastings algorithm works for this type of application, but it is crucial to control the acceptance rate to achieve proper convergence. The quadratic approximations outperformed the linear ones when the reflection data have a low noise level, otherwise the inversion will yield more or less equal results. On the other hand, the differences between exact Zoeppritz equations and the quadratic approximations are minor, even in the case without noise. Last, by including PS reflection data in the inversion the bias and the uncertainties are greatly reduced. ACKNOWLEDGMENT We wish to thank BP, Hydro, Schlumberger, Statoil, and The Research Council of Norway for their support through the Uncertainty in Reservoir Evaluation (URE) project. REFERENCES Aki, K. & Richards, P. G., 198. Quantitative Seismology: Theory and Methods, W. H. Freeman and Co. Buland, A. & Omre, H., 23a. Bayesian linearized AVO inversion, Geophysics, 68, Buland, A. & Omre, H., 23b. Joint AVO inversion, wavelet estimation and noise-level estimation using a spatially coupled hierarchical bayesian model, Geophysical Prospecting, 51,

14 14 T.E. Rabben, H. Tjelmeland, and B. Ursin Dȩbski, W. & Tarantola, A., Information on elastic parameters obtained from the amplitudes of reflected waves, Geophysics, 6, Hampson, D., AVO inversion, theory and practice, The Leading Edge. Lines, L. R. & Treitel, S., A review of least-squares inversion and its application to geophysical problems, Geophysical Prospecting, 32, Liu, J. S., 21. Monte Carlo strategies in scientific computing, Springer. Margrave, G. F., Stewart, R. R., & Larsen, J. A., 21. Joint PP and PS seismic inversion, The Leading Edge. Robert, C. P. & Casella, G., Monte Carlo statistical methods, Springer. Sen, M. K. & Stoffa, P. L., Bayesian inference, Gibbs sampler and uncertainty estimation in geophysical inversion, Geophysical Prospecting, 44, Smith, G. C. & Gidlow, P. M., Weighted stacking for rock property estimation and detection of gas, Geophysical Prospecting, 35, Stewart, R. R., 199. Joint P and P-SV seismic inversion, Tech. rep., CREWES. Stovas, A. & Ursin, B., 21. Second-order approximations of the reflection and transmission coefficients between two visco-elastic isotropic media, Journal of Seismic Exploration, 9, Stovas, A. & Ursin, B., 23. Reflection and transmission responses of layered transversely isotropic viscoelastic media, Geophysical Prospecting, 51, Tenorio, L., 21. Statistical regularization of inverse problems, SIAM Review, 43, Tjelmeland, H. & Eidsvik, J., 25. Directional Metropolis-Hastings update for posteriors with nonlinear likelihoods, in Geostatistics Banff 24, edited by O. Leuangthong & C. V. Deutsch. Ursin, B. & Tjåland, E., The information content of the elastic relfection matrix, Geophysical Journal International, 125, Veire, H. H. & Landrø, M., 26. Simultaneous inversion of PP and PS seismic data, Geophysics, 71, R1 R1.

15 Nonlinear Bayesian inversion 15 APPENDIX A: SECOND ORDER REFLECTION COEFFICIENTS The complete set of reflection coefficients between two transversely isotropic media for a downgoing wave, derived by Stovas & Ursin (21), reads R D = r PP r SP r PS r SS g 11 g 12 f g f(g 11 g 22 ). g f(g 2 11 g 22 ) g 22 + g 12 f (A1) The parametrization we have chosen to express f, g 11, g 22, and g 12 in is shear and plane wave modulus together with density. This choice results in the most compact expressions when simplified to isotropic media: 1 g 11 4 cos 2 θ p 2 sin 2 1 θ s (1 4 tan2 θ p ) g 22 2 sin 2 θ s 1 4 cos = 2 θ s 1(1 4 tan2 θ s ) 1 g 12 k(cos θ s (cos θ s + γ cos θ p ) 1) 2 k 1 f k(cos θ s (cos θ s γ cos θ p ) 1) k 2 M M µ µ ρ (A2) with the plane wave modulus, M = ρα 2 ; the shear modulus, µ = ρβ 2 ; and density, ρ. The difference is M = M 2 M 1 = ρ 2 α 2 2 ρ 1 α 2 1 and the average is M = (ρ 2 α ρ 1 α 2 1)/2 with equal definitions for µ and ρ. The background v S /v P ratio is γ = β/ᾱ and the constant k is tan θp tan θ s. See Fig. A1 for the definition of the angles and parameters. r PP = For a downgoing P-wave reflected to P and S-waves the explicit expressions are 1 M 4 cos 2 θ p M 2 µ sin2 θ s µ (1 tan2 θ p ) ρ ( ) 2 µ + tan θ p tan θ s [γ 2 (1 (1 + γ 2 ) sin 2 θ p ) 1 4 ( ρ ) 2 + sin 2 θ s ( µ µ ρ µ (A3) ) ]

16 16 T.E. Rabben, H. Tjelmeland, and B. Ursin and r PS = {[ tan θ p tan θ s (1 cos θ s (cos θ s + γ cos θ p )) µ µ 1 2 [ ] + 1 (1 cos θ s (cos θ s γ cos θ p )) µ 2 µ 1 ρ 2 [ ( ) 1 M 4 cos 2 θ p M + 1 µ 4 sin 2 θ 4 cos 2 s θ s µ ] ρ ( ) ]} ρ 2 tan 2 θ p tan 2 θ s. To express in other parametrizations, e.g. (1) and (2), the following relations can be used M M = (1 4 3 γ2 ) K K + 4 µ 3γ2 µ M M M M = 2(1 γ2 ) 2 M M M M µ = (1 2γ2 ) λ λ + 2γ2 µ µ σ + µ γ 2 µ = 2 α ᾱ + ρ = 2 I α I α µ = 2 β β µ µ = 2 I β Ī β ρ + ρ ρ, (A4) (A5) where K is bulk modulus, M = K µ; λ is Lamé constant, K = λ + 2µ; σ is Poisson s ratio, σ = 1 2 (M 2µ)/(M µ); and I is impedance, I α = ρα and I β = ρβ. APPENDIX B: BAYESIAN FORMULATION The joint posterior posterior distribution (8) can be simplified to π(m, σ 2 e, σ 2 m d) π(d m, σ 2 e) π(m σ 2 m) π(σ 2 e) π(σ 2 m), (B1) by assuming the prior of m and likelihood of d to be independent of σ 2 e and σ 2 m respectively and σ 2 e and σ 2 m to be a priori independent. We are assuming the error to be multivariate Gaussian distributed with zero mean and

17 Nonlinear Bayesian inversion 17 covariance σ 2 eσ e, π(e σ 2 e) = N (e;, σ 2 eσ e ). (B2) The prior of m is also assumed to be a multivariate Gaussian distribution and the fact that the forward model is deterministic together with (B2) makes also the likelihood of d multivariate Gaussian π(m σ 2 m) = N (m; µ m, σ 2 mσ m ) π(d m, σ 2 e) = N (d; f(m), σ 2 eσ e ). (B3) Here f(m) is the nonlinear forward model in (6), µ m is the prior expected value of m, and Σ e and Σ m the covariance matrices excluding the level factors σe 2 and σm. 2 The Σ s allow for specifying a priori spatial correlations, correlations between parameters within a grid point, and correlated noise for neighbouring reflection angles, while the level factors, σm 2 and σe, 2 will be explored as a part of the inversion. For details about the multivariate Gaussian distribution see Appendix E. A graphical representation of the prior of m, the likelihood, and the deterministic relation is shown in the graph in Fig. B1. The prior of σe 2 and σm 2 are assigned inverse gamma distributions π(σ 2 e) = IG(σ 2 e; α e, β e ) π(σ 2 m) = IG(σ 2 m; α m, β m ), (B4) with corresponding parameters α and β. For details about the inverse gamma distribution see Appendix E. APPENDIX C: SAMPLING Analytical calculations of the posterior (B1) is not possible. Instead, the joint posterior distributions for m, σe, 2 and σm 2 will be sampled iteratively in a Metropolis-Hastings algorithm. In a simple case of sampling x from the posterior π(x y), in our case x could be m, σe 2 or σm 2 and y be d, the algorithm consists of the two steps (i) Propose a new sample x u q( y)

18 18 T.E. Rabben, H. Tjelmeland, and B. Ursin (ii) Accept x u with a probability p(x u x) = min where q is the proposal distribution. { } 1, π(xu y) q(x y) π(x y) q(x u y) First, since σ 2 e and σ 2 m are assumed to be a posteriori independent they can be sampled in parallel. Their full conditional are π(σ 2 e d, m) π(d m, σ 2 e) π(σ 2 e) π(σ 2 m m) π(m σ 2 m) π(σ 2 m), (C1) and with the choice of likelihood and prior, (B3) and (B4), the posterior also become inverse gamma distributed only with modified parameters. A proof of this is shown in Appendix E. The explicit posterior expressions read ( π(σe d, 2 m) IG σe 2 α + n ) e 2, β + n e s2 e 2 (C2) with s 2 e = 1 n e (d f(m)) T Σ 1 e (d f(m)) (C3) and ( π(σm m) 2 IG with σ 2 m α + n m 2, β + n m s2 m 2 s 2 m = 1 n m (m µ m ) T Σ 1 m (m µ m ) ) (C4) (C5) The scalars are, from (3) and (4), n e = n x n y D d and n m = n x n y D m. Since the posterior distributions are simply inverse gamma distributed, this will also become our proposal distribution and the acceptance probability will be identical to one. The updates of σe 2 and σm 2 are therefore examples of Gibbs updates. Second, sampling the posterior of m is not as simple as for the σ s. It is not possible to choose the full conditional as the proposal distribution because of the nonlinear likelihood, there is no algorithm to sample from it. The natural choice is then to use the linearised likelihood and the posterior then becomes π lin (m d, σ 2 e, σ 2 m) π lin (d m, σ 2 e) π(m σ 2 m), (C6)

19 Nonlinear Bayesian inversion 19 where π lin (d m, σe) 2 is the linearised likelihood using (7). The use of a simplified proposal distribution also implies that the acceptance rate will be different from zero. In fact, if we try to update the complete m, the acceptance rate will be impractically low - if not zero. We therefore have to update only blocks of m and repeat this several times to produce an update of m in the joint posterior distribution. More specifically, the proposal distribution is conditioned on d and m in two regions A and B together with the scalars σe 2 and σm, 2 m u A π lin ( m B, d A, d B, σ 2 e, σ 2 m). (C7) Fig. C1 shows the different blocks. The reason for including B is to impose continuity of the update with respect to the surrounding grid points. Block C does not contribute to the distribution of m u A. For details on the linear update distribution see Appendix D. The acceptance probability reads { } p(m u m) = min 1, π(mu d, σe, 2 σm) 2 π(m d, σe, 2 σm) πlin (m A m B, d A, d B, σe, 2 σm) 2 2 π lin (m u A m, (C8) B, d A, d B, σe, 2 σm) 2 where the first fraction is calculated using the exact, nonlinear version of (C6). APPENDIX D: LINEAR PROPOSAL DISTRIBUTION To sample m u A from the distribution (C7) we need the linear forward model (7). For block A and B it reads d A = F A m A + e A d B = F B m B + e B. Partitioning the prior of m and the noise model yield m A N µ A, σ2 mσ m,aa m B µ B σmσ 2 m,ba and e A N, σ2 eσ e,aa σeσ 2 e,ab σeσ 2 e,ba σeσ 2 e,bb. e B σmσ 2 m,ab σmσ 2 m,bb (D1) (D2) (D3)

20 2 T.E. Rabben, H. Tjelmeland, and B. Ursin We are now able to calculate the mean and correlations between m A, m B, d A, and d B that we need. For instance E(d A ) = E(F A m A + e A ) = F A µ A Cov(m A, d B ) = Cov(m A, F B m B + ɛ B ) = σ 2 mσ m,ab F T B (D4) Cov(d A, d B ) = Cov(F A m A + e A, F B m B + e B ) = Cov(F A m A, F B m B ) Cov(e A, e B ) = σ 2 mf A Σ m,ab F T B + σ 2 eσ e,ab. The full distribution becomes m A µ A σ mσ 2 m,aa m B µ B σmσ 2 m,ba = N, d A F A µ A σ mf 2 A Σ m,aa d B F B µ B σmf 2 B Σ m,ba σmσ 2 m,ab σmσ 2 m,aa FA T σmσ 2 m,ab FB T σmσ 2 m,bb σmσ 2 m,ba FA T σmσ 2 m,bb FB T, σmf 2 A Σ m,ab σmf 2 A Σ m,aa FA T + σ2 eσ e,aa σmf 2 A Σ m,ab FB T + σ2 eσ e,ab σmf 2 B Σ m,bb σmf 2 B Σ m,ba FA T + σ2 eσ e,ba σmf 2 B Σ m,bb FB T + σ2 eσ e,bb (D5) which we conveniently group and define by m A x R = N µ A µ R, Σ AA Σ RA Σ AR Σ RR. (D6)

21 Nonlinear Bayesian inversion 21 Using the expression for conditional normal distribution, equation (E3), the proposal distribution becomes N (m u A ; µ u, Σ u ) with µ u = µ A + Σ AR Σ 1 RR (x R µ R ) Σ u = Σ AA Σ AR Σ 1 RR Σ RA. (D7) APPENDIX E: STATISTICAL DISTRIBUTIONS A multivariate Gaussian variable x with expectation vector µ and covariance matrix Σ has the probability function { 1 N (x; µ, Σ) = exp 1 } (2π) n/2 Σ 1/2 2 (x µ)t Σ 1 (x µ), (E1) where n is the dimension of x. For two multivariate Gaussian variables x 1 and x 2 with joint distribution x 1 N 1, Σ 11 Σ 12, Σ 21 Σ 22 x 2 µ 2 (E2) the conditional expectation and variance are µ 1 2 = µ 1 + Σ 12 Σ 1 22 (x 2 µ 2 ) Σ 1 2 = Σ 11 Σ 12 Σ 1 22 Σ 21. The inverse gamma probability function is (E3) IG(x; α, β) = βα Γ(α) ( 1 x where x, α >, and β >. ) α+1 e β/x (E4) Given the prior distribution σ 2 IG(α, β) and measurements x N (µ, σ 2 Σ), the

22 22 T.E. Rabben, H. Tjelmeland, and B. Ursin posterior distribution of σ 2 is π(σ 2 x) π(x σ 2 ) π(σ 2 ) where = N (µ, σ 2 Σ) IG(α, β) { 1 exp 1 } (σ 2 ) n/2 Σ 1/2 2 σ 2 (x µ) T Σ 1 (x µ) ( ) α+1 1 exp { βσ } 2 σ 2 ( ) α+1+n/2 1 exp { σ 2 (β + s 2 n2 } ) IG σ 2 ( σ 2 α n + 2, β + n ) s2 2 (E5) s 2 = 1 n (x µ)t Σ 1 (x µ) (E6) and n is the dimension of x. Clearly, the posterior is also inverse gamma but with modified parameters. This paper has been produced using the Blackwell Publishing GJI L A TEX2e class file.

23 Nonlinear Bayesian inversion 23 I α I α.5 I β I β.5 Distance [m] Distance [m].2 Distance [m] Distance [m].2 ρ.5 Distance [m] Distance [m].2 Figure 1. The true contrasts in m. Upper left is P-wave impedance, upper right is S-wave impedance, and lower is density

24 24 T.E. Rabben, H. Tjelmeland, and B. Ursin r PP θ =.2.1 Linear Quadratic θ = θ = θ = Figure 2. To the left is PP reflection coefficients from the Zoeppritz model for 4 different incidence angles. The two right columns show the bias in the linear and quadratic approximations, relative to the Zoeppritz model, for the nonzero angles. For θ = the bias is zero.

25 Nonlinear Bayesian inversion 25 θ = θ = θ = 4 θ = Figure 3. PP reflection coefficients from the Zoeppritz model and including multivariate normal distributed noise.

26 26 T.E. Rabben, H. Tjelmeland, and B. Ursin r PS.4 Linear Quadratic.3 θ = θ = θ = Figure 4. To the left is PS reflection coefficients from the Zoeppritz model for 3 nonzero incidence angles. The two right columns show the bias in the linear and quadratic approximations, relative to the Zoeppritz model.

27 Nonlinear Bayesian inversion 27 r PP, θ = r PS, θ = r PP, θ = 2 r PS, θ = r PP, θ = 4 r PS, θ = r PP, θ = Figure 5. PP and PS reflection coefficients from the Zoeppritz model and including multivariate normal distributed noise.

28 28 T.E. Rabben, H. Tjelmeland, and B. Ursin Linear Quadratic Zoeppritz Iα Iα.5.6 Iβ Iα.3.6 ρ.3 Figure 6. Bias in the posterior distribution of m from PP inversion without noise. Each column is the result of three inversions using three different models; linear, quadratic, and exact Zoeppritz. The rows displays the three different parameters of m; contrasts in P-wave impedance, S-wave impedance, and density.

29 Nonlinear Bayesian inversion 29 Linear Quadratic Zoeppritz Iα Iα Iβ Iα ρ.12 Figure 7. Standard deviation in the posterior distribution of m from PP inversion without noise.

30 3 T.E. Rabben, H. Tjelmeland, and B. Ursin Probability Linear Quadratic Zoeppritz σ 2 m σ 2 e Figure 8. Posterior distribution of σ 2 m and σ 2 e from PP inversion without noise

31 Linear Quadratic Zoeppritz Nonlinear Bayesian inversion Iα Iα Iβ Iα ρ.13 Figure 9. Bias in the posterior distribution of m from PP inversion including multivariate Gaussian distributed noise.

32 32 T.E. Rabben, H. Tjelmeland, and B. Ursin Linear Quadratic Zoeppritz.1 Iα Iα.5.1 Iβ Iα.5.1 ρ.5 Figure 1. Standard deviation in the posterior distribution of m from PP inversion including multivariate Gaussian distributed noise.

33 Nonlinear Bayesian inversion 33 Probability Linear Quadratic Zoeppritz σ 2 m σ 2 e Figure 11. Posterior distribution of σ 2 m and σ 2 e from PP inversion including multivariate Gaussian distributed noise

34 34 T.E. Rabben, H. Tjelmeland, and B. Ursin Linear Quadratic Zoeppritz Iα Iα.5.6 Iβ Iα.3.6 ρ.3 Figure 12. Bias in the posterior distribution of m from joint PP and PS inversion without noise.

35 Nonlinear Bayesian inversion 35 Linear Quadratic Zoeppritz Iα Iα Iβ Iα ρ.12 Figure 13. Standard deviation in the posterior distribution of m from joint PP and PS inversion without noise.

36 36 T.E. Rabben, H. Tjelmeland, and B. Ursin.2.15 Linear Quadratic Zoeppritz.15 Probability σ 2 m σ 2 e Figure 14. Posterior distribution of σ 2 m and σ 2 e from joint PP and PS inversion without noise

37 Linear Quadratic Zoeppritz Nonlinear Bayesian inversion Iα Iα Iβ Iα ρ.13 Figure 15. Bias in the posterior distribution of m from joint PP and PS inversion including multivariate Gaussian distributed noise.

38 38 T.E. Rabben, H. Tjelmeland, and B. Ursin Linear Quadratic Zoeppritz.1 Iα Iα.5.1 Iβ Iα.5.1 ρ.5 Figure 16. Standard deviation in the posterior distribution of m from joint PP and PS inversion including multivariate Gaussian distributed noise.

39 Nonlinear Bayesian inversion Linear Quadratic Zoeppritz.15 Probability σ 2 m σ 2 e Figure 17. Posterior distribution of σ 2 m and σ 2 e from joint PP and PS inversion including multivariate Gaussian distributed noise

40 4 T.E. Rabben, H. Tjelmeland, and B. Ursin θ p θ p θ s α 1, β 1, ρ 1 α 2, β 2, ρ 2 Figure A1. Relation between the downgoing P-wave θ p, reflected P-wave θ p, and reflected S-wave θ s and the definition of velocities and density in upper and lower medium.

41 Nonlinear Bayesian inversion 41 µ m σ 2 mσ m m f(m) σ 2 eσ e d Figure B1. Directed acyclic graph. The upper square is the prior model, the lower is the likelihood. The relation between m and d is the deterministic forward model.

42 42 T.E. Rabben, H. Tjelmeland, and B. Ursin A B C Figure C1. Block A is the block to be updated, B is a boundary zone of limited thickness and C is the rest.

### Bayesian lithology/fluid inversion comparison of two algorithms

Comput Geosci () 14:357 367 DOI.07/s596-009-9155-9 REVIEW PAPER Bayesian lithology/fluid inversion comparison of two algorithms Marit Ulvmoen Hugo Hammer Received: 2 April 09 / Accepted: 17 August 09 /

### NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET. Directional Metropolis Hastings updates for posteriors with nonlinear likelihoods

NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET Directional Metropolis Hastings updates for posteriors with nonlinear likelihoods by Håkon Tjelmeland and Jo Eidsvik PREPRINT STATISTICS NO. 5/2004 NORWEGIAN

### Stochastic vs Deterministic Pre-stack Inversion Methods. Brian Russell

Stochastic vs Deterministic Pre-stack Inversion Methods Brian Russell Introduction Seismic reservoir analysis techniques utilize the fact that seismic amplitudes contain information about the geological

### Linearized AVO and Poroelasticity for HRS9. Brian Russell, Dan Hampson and David Gray 2011

Linearized AO and oroelasticity for HR9 Brian Russell, Dan Hampson and David Gray 0 Introduction In this talk, we combine the linearized Amplitude ariations with Offset (AO) technique with the Biot-Gassmann

### Reservoir Characterization using AVO and Seismic Inversion Techniques

P-205 Reservoir Characterization using AVO and Summary *Abhinav Kumar Dubey, IIT Kharagpur Reservoir characterization is one of the most important components of seismic data interpretation. Conventional

### Correspondence between the low- and high-frequency limits for anisotropic parameters in a layered medium

GEOPHYSICS VOL. 74 NO. MARCH-APRIL 009 ; P. WA5 WA33 7 FIGS. 10.1190/1.3075143 Correspondence between the low- and high-frequency limits for anisotropic parameters in a layered medium Mercia Betania Costa

### Reservoir properties inversion from AVO attributes

Reservoir properties inversion from AVO attributes Xin-gang Chi* and De-hua Han, University of Houston Summary A new rock physics model based inversion method is put forward where the shaly-sand mixture

### AVAZ inversion for fracture orientation and intensity: a physical modeling study

AVAZ inversion for fracture orientation and intensity: a physical modeling study Faranak Mahmoudian*, Gary F. Margrave, and Joe Wong, University of Calgary. CREWES fmahmoud@ucalgary.ca Summary We present

### SEISMIC INVERSION OVERVIEW

DHI CONSORTIUM SEISMIC INVERSION OVERVIEW Rocky Roden September 2011 NOTE: Terminology for inversion varies, depending on the different contractors and service providers, emphasis on certain approaches,

### Inversion of seismic AVA data for porosity and saturation

Inversion of seismic AVA data for porosity and saturation Brikt Apeland Thesis for the degree Master of Science Department of Earth Science University of Bergen 27th June 213 2 Abstract Rock physics serves

### AFI (AVO Fluid Inversion)

AFI (AVO Fluid Inversion) Uncertainty in AVO: How can we measure it? Dan Hampson, Brian Russell Hampson-Russell Software, Calgary Last Updated: April 2005 Authors: Dan Hampson, Brian Russell 1 Overview

### Statistical Rock Physics

Statistical - Introduction Book review 3.1-3.3 Min Sun March. 13, 2009 Outline. What is Statistical. Why we need Statistical. How Statistical works Statistical Rock physics Information theory Statistics

### AVO attribute inversion for finely layered reservoirs

GEOPHYSICS, VOL. 71, NO. 3 MAY-JUNE 006 ; P. C5 C36, 16 FIGS., TABLES. 10.1190/1.197487 Downloaded 09/1/14 to 84.15.159.8. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/

### Approximations to the Zoeppritz equations and their use in AVO analysis

GEOPHYSICS, VOL. 64, NO. 6 (NOVEMBER-DECEMBER 1999; P. 190 197, 4 FIGS., 1 TABLE. Approximations to the Zoeppritz equations their use in AVO analysis Yanghua Wang ABSTRACT To efficiently invert seismic

### A simple algorithm for band-limited impedance inversion

Impedance inversion algorithm A simple algorithm for band-limited impedance inversion Robert J. Ferguson and Gary F. Margrave ABSTRACT This paper describes a seismic inversion method which has been cast

### Bayesian inversion (I):

Bayesian inversion (I): Advanced topics Wojciech Dȩbski Instytut Geofizyki PAN debski@igf.edu.pl Wydział Fizyki UW, 13.10.2004 Wydział Fizyki UW Warszawa, 13.10.2004 (1) Plan of the talk I. INTRODUCTION

### CPSC 540: Machine Learning

CPSC 540: Machine Learning MCMC and Non-Parametric Bayes Mark Schmidt University of British Columbia Winter 2016 Admin I went through project proposals: Some of you got a message on Piazza. No news is

### Gaussian processes. Chuong B. Do (updated by Honglak Lee) November 22, 2008

Gaussian processes Chuong B Do (updated by Honglak Lee) November 22, 2008 Many of the classical machine learning algorithms that we talked about during the first half of this course fit the following pattern:

### Probabilistic Graphical Models

Probabilistic Graphical Models Brown University CSCI 2950-P, Spring 2013 Prof. Erik Sudderth Lecture 12: Gaussian Belief Propagation, State Space Models and Kalman Filters Guest Kalman Filter Lecture by

### An investigation of the free surface effect

An investigation of the free surface effect Nasser S. Hamarbitan and Gary F. Margrave, An investigation of the free surface effect ABSTRACT When P and S seismic waves are incident on a solid-air interface

### Probabilistic Graphical Models

2016 Robert Nowak Probabilistic Graphical Models 1 Introduction We have focused mainly on linear models for signals, in particular the subspace model x = Uθ, where U is a n k matrix and θ R k is a vector

### Introduction to Bayesian methods in inverse problems

Introduction to Bayesian methods in inverse problems Ville Kolehmainen 1 1 Department of Applied Physics, University of Eastern Finland, Kuopio, Finland March 4 2013 Manchester, UK. Contents Introduction

### Stat260: Bayesian Modeling and Inference Lecture Date: February 10th, Jeffreys priors. exp 1 ) p 2

Stat260: Bayesian Modeling and Inference Lecture Date: February 10th, 2010 Jeffreys priors Lecturer: Michael I. Jordan Scribe: Timothy Hunter 1 Priors for the multivariate Gaussian Consider a multivariate

### Integration of Rock Physics Models in a Geostatistical Seismic Inversion for Reservoir Rock Properties

Integration of Rock Physics Models in a Geostatistical Seismic Inversion for Reservoir Rock Properties Amaro C. 1 Abstract: The main goal of reservoir modeling and characterization is the inference of

### Delivery: an open-source Bayesian seismic inversion tool. James Gunning, CSIRO Petroleum Michael Glinsky,, BHP Billiton

Delivery: an open-source Bayesian seismic inversion tool James Gunning, CSIRO Petroleum Michael Glinsky,, BHP Billiton Application areas Development, appraisal and exploration situations Requisite information

### When one of things that you don t know is the number of things that you don t know

When one of things that you don t know is the number of things that you don t know Malcolm Sambridge Research School of Earth Sciences, Australian National University Canberra, Australia. Sambridge. M.,

### Research - Recent work:

Research - Recent work: Henning Omre - www.math.ntnu.no/~omre Professor in Statistics Department of Mathematical Sciences Norwegian University of Science & Technology Trondheim Norway December 2016 Current

### Quantitative Interpretation

Quantitative Interpretation The aim of quantitative interpretation (QI) is, through the use of amplitude analysis, to predict lithology and fluid content away from the well bore. This process should make

### Analysis of a prior model for a blocky inversion of seismic amplitude versus offset data

Analysis of a prior model for a blocky inversion of seismic amplitude versus offset data Short title: Blocky inversion of seismic AVO data Ulrich Theune* StatoilHydro Research Center, 7005 Trondheim, NORWAY,

### Unphysical negative values of the anelastic SH plane wave energybased transmission coefficient

Shahin Moradi and Edward S. Krebes Anelastic energy-based transmission coefficient Unphysical negative values of the anelastic SH plane wave energybased transmission coefficient ABSTRACT Computing reflection

### A Case Study on Simulation of Seismic Reflections for 4C Ocean Bottom Seismometer Data in Anisotropic Media Using Gas Hydrate Model

A Case Study on Simulation of Seismic Reflections for 4C Ocean Bottom Seismometer Data in Anisotropic Media Using Gas Hydrate Model Summary P. Prasada Rao*, N. K. Thakur 1, Sanjeev Rajput 2 National Geophysical

### P-wave impedance, S-wave impedance and density from linear AVO inversion: Application to VSP data from Alberta

Linear AVO inversion P-wave impedance, S-wave impedance and density from linear AVO inversion: Application to VSP data from Alberta Faranak Mahmoudian and Gary F. Margrave ABSTRACT In AVO (amplitude variation

### The prediction of reservoir

Risk Reduction in Gas Reservoir Exploration Using Joint Seismic-EM Inversion GAS EXPLORATION By Yoram Rubin, G. Michael Hoversten, Zhangshuan Hou and Jinsong Chen, University of California, Berkeley A

### Elastic impedance inversion from robust regression method

Elastic impedance inversion from robust regression method Charles Prisca Samba 1, Liu Jiangping 1 1 Institute of Geophysics and Geomatics,China University of Geosciences, Wuhan, 430074 Hubei, PR China

### AVAZ and VVAZ practical analysis to estimate anisotropic properties

AVAZ and VVAZ practical analysis to estimate anisotropic properties Yexin Liu*, SoftMirrors Ltd., Calgary, Alberta, Canada yexinliu@softmirrors.com Summary Seismic anisotropic properties, such as orientation

### Parameter estimation: A new approach to weighting a priori information

c de Gruyter 27 J. Inv. Ill-Posed Problems 5 (27), 2 DOI.55 / JIP.27.xxx Parameter estimation: A new approach to weighting a priori information J.L. Mead Communicated by Abstract. We propose a new approach

### Rock physics and AVO analysis for lithofacies and pore fluid prediction in a North Sea oil field

Rock physics and AVO analysis for lithofacies and pore fluid prediction in a North Sea oil field Downloaded 09/12/14 to 84.215.159.82. Redistribution subject to SEG license or copyright; see Terms of Use

### Bayesian inference: what it means and why we care

Bayesian inference: what it means and why we care Robin J. Ryder Centre de Recherche en Mathématiques de la Décision Université Paris-Dauphine 6 November 2017 Mathematical Coffees Robin Ryder (Dauphine)

### The Bayesian Approach to Multi-equation Econometric Model Estimation

Journal of Statistical and Econometric Methods, vol.3, no.1, 2014, 85-96 ISSN: 2241-0384 (print), 2241-0376 (online) Scienpress Ltd, 2014 The Bayesian Approach to Multi-equation Econometric Model Estimation

### Comparative Study of AVO attributes for Reservoir Facies Discrimination and Porosity Prediction

5th Conference & Exposition on Petroleum Geophysics, Hyderabad-004, India PP 498-50 Comparative Study of AVO attributes for Reservoir Facies Discrimination and Porosity Prediction Y. Hanumantha Rao & A.K.

### An introduction to Sequential Monte Carlo

An introduction to Sequential Monte Carlo Thang Bui Jes Frellsen Department of Engineering University of Cambridge Research and Communication Club 6 February 2014 1 Sequential Monte Carlo (SMC) methods

### COM336: Neural Computing

COM336: Neural Computing http://www.dcs.shef.ac.uk/ sjr/com336/ Lecture 2: Density Estimation Steve Renals Department of Computer Science University of Sheffield Sheffield S1 4DP UK email: s.renals@dcs.shef.ac.uk

### Towards Interactive QI Workflows Laurie Weston Bellman*

Laurie Weston Bellman* Summary Quantitative interpretation (QI) is an analysis approach that is widely applied (Aki and Richards, 1980, Verm and Hilterman, 1995, Avseth et al., 2005, Weston Bellman and

### Bayesian Linear Models

Bayesian Linear Models Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department of Forestry & Department

### 1 Using standard errors when comparing estimated values

MLPR Assignment Part : General comments Below are comments on some recurring issues I came across when marking the second part of the assignment, which I thought it would help to explain in more detail

### Large Scale Modeling by Bayesian Updating Techniques

Large Scale Modeling by Bayesian Updating Techniques Weishan Ren Centre for Computational Geostatistics Department of Civil and Environmental Engineering University of Alberta Large scale models are useful

### Dynamic Data Driven Simulations in Stochastic Environments

Computing 77, 321 333 (26) Digital Object Identifier (DOI) 1.17/s67-6-165-3 Dynamic Data Driven Simulations in Stochastic Environments C. Douglas, Lexington, Y. Efendiev, R. Ewing, College Station, V.

### Bayesian Monte Carlo Filtering for Stochastic Volatility Models

Bayesian Monte Carlo Filtering for Stochastic Volatility Models Roberto Casarin CEREMADE University Paris IX (Dauphine) and Dept. of Economics University Ca Foscari, Venice Abstract Modelling of the financial

### Bayesian Linear Models

Bayesian Linear Models Sudipto Banerjee 1 and Andrew O. Finley 2 1 Department of Forestry & Department of Geography, Michigan State University, Lansing Michigan, U.S.A. 2 Biostatistics, School of Public

### Inference in Bayesian Networks

Andrea Passerini passerini@disi.unitn.it Machine Learning Inference in graphical models Description Assume we have evidence e on the state of a subset of variables E in the model (i.e. Bayesian Network)

### Transdimensional Markov Chain Monte Carlo Methods. Jesse Kolb, Vedran Lekić (Univ. of MD) Supervisor: Kris Innanen

Transdimensional Markov Chain Monte Carlo Methods Jesse Kolb, Vedran Lekić (Univ. of MD) Supervisor: Kris Innanen Motivation for Different Inversion Technique Inversion techniques typically provide a single

### PATTERN RECOGNITION AND MACHINE LEARNING CHAPTER 13: SEQUENTIAL DATA

PATTERN RECOGNITION AND MACHINE LEARNING CHAPTER 13: SEQUENTIAL DATA Contents in latter part Linear Dynamical Systems What is different from HMM? Kalman filter Its strength and limitation Particle Filter

### Bayesian Prediction of Code Output. ASA Albuquerque Chapter Short Course October 2014

Bayesian Prediction of Code Output ASA Albuquerque Chapter Short Course October 2014 Abstract This presentation summarizes Bayesian prediction methodology for the Gaussian process (GP) surrogate representation

### Reliability of Seismic Data for Hydrocarbon Reservoir Characterization

Reliability of Seismic Data for Hydrocarbon Reservoir Characterization Geetartha Dutta (gdutta@stanford.edu) December 10, 2015 Abstract Seismic data helps in better characterization of hydrocarbon reservoirs.

### Parametric Techniques

Parametric Techniques Jason J. Corso SUNY at Buffalo J. Corso (SUNY at Buffalo) Parametric Techniques 1 / 39 Introduction When covering Bayesian Decision Theory, we assumed the full probabilistic structure

### Pre-Stack Seismic Inversion and Amplitude Versus Angle Modeling Reduces the Risk in Hydrocarbon Prospect Evaluation

Advances in Petroleum Exploration and Development Vol. 7, No. 2, 2014, pp. 30-39 DOI:10.3968/5170 ISSN 1925-542X [Print] ISSN 1925-5438 [Online] www.cscanada.net www.cscanada.org Pre-Stack Seismic Inversion

### Stochastic inversion of seismic PP and PS data for reservoir parameter. estimation

Stochastic inversion of seismic PP and PS data for reservoir parameter estimation Shortened Title: Stochastic inversion of PP and PS data Jinsong Chen (jchen@lbl.gov) Lawrence Berkeley National Laboratory,

### Stat 535 C - Statistical Computing & Monte Carlo Methods. Arnaud Doucet.

Stat 535 C - Statistical Computing & Monte Carlo Methods Arnaud Doucet Email: arnaud@cs.ubc.ca 1 1.1 Outline Introduction to Markov chain Monte Carlo The Gibbs Sampler Examples Overview of the Lecture

### Mark Gales October y (x) x 1. x 2 y (x) Inputs. Outputs. x d. y (x) Second Output layer layer. layer.

University of Cambridge Engineering Part IIB & EIST Part II Paper I0: Advanced Pattern Processing Handouts 4 & 5: Multi-Layer Perceptron: Introduction and Training x y (x) Inputs x 2 y (x) 2 Outputs x

### Earthscope Imaging Science & CIG Seismology Workshop

Earthscope Imaging Science & CIG Seismology Introduction to Direct Imaging Methods Alan Levander Department of Earth Science Rice University 1 Two classes of scattered wave imaging systems 1. Incoherent

### Log Ties Seismic to Ground Truth

26 GEOPHYSICALCORNER Log Ties Seismic to Ground Truth The Geophysical Corner is a regular column in the EXPLORER, edited by R. Randy Ray. This month s column is the first of a two-part series titled Seismic

### MCMC algorithms for fitting Bayesian models

MCMC algorithms for fitting Bayesian models p. 1/1 MCMC algorithms for fitting Bayesian models Sudipto Banerjee sudiptob@biostat.umn.edu University of Minnesota MCMC algorithms for fitting Bayesian models

### Fundamental Issues in Bayesian Functional Data Analysis. Dennis D. Cox Rice University

Fundamental Issues in Bayesian Functional Data Analysis Dennis D. Cox Rice University 1 Introduction Question: What are functional data? Answer: Data that are functions of a continuous variable.... say

### V003 How Reliable Is Statistical Wavelet Estimation?

V003 How Reliable Is Statistical Wavelet Estimation? J.A. Edgar* (BG Group) & M. van der Baan (University of Alberta) SUMMARY Well logs are often used for the estimation of seismic wavelets. The phase

### Bayesian Paradigm. Maximum A Posteriori Estimation

Bayesian Paradigm Maximum A Posteriori Estimation Simple acquisition model noise + degradation Constraint minimization or Equivalent formulation Constraint minimization Lagrangian (unconstraint minimization)

### The Origin of Deep Learning. Lili Mou Jan, 2015

The Origin of Deep Learning Lili Mou Jan, 2015 Acknowledgment Most of the materials come from G. E. Hinton s online course. Outline Introduction Preliminary Boltzmann Machines and RBMs Deep Belief Nets

### Variational Principal Components

Variational Principal Components Christopher M. Bishop Microsoft Research 7 J. J. Thomson Avenue, Cambridge, CB3 0FB, U.K. cmbishop@microsoft.com http://research.microsoft.com/ cmbishop In Proceedings

### TIV Contrast Source Inversion of mcsem data

TIV Contrast ource Inversion of mcem data T. Wiik ( ), L. O. Løseth ( ), B. Ursin ( ), K. Hokstad (, ) ( )epartment for Petroleum Technology and Applied Geophysics, Norwegian University of cience and Technology

### Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications

Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications Dongbin Xiu Department of Mathematics, Purdue University Support: AFOSR FA955-8-1-353 (Computational Math) SF CAREER DMS-64535

### Bandlimited impedance inversion: using well logs to fill low frequency information in a non-homogenous model

Bandlimited impedance inversion: using well logs to fill low frequency information in a non-homogenous model Heather J.E. Lloyd and Gary F. Margrave ABSTRACT An acoustic bandlimited impedance inversion

### P-wave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry

GEOPHYSICS, VOL. 62, NO. 3 (MAY-JUNE 1997); P. 713 722, 7 FIGS., 1 TABLE. P-wave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry Andreas Rüger ABSTRACT

### Nonparametric Drift Estimation for Stochastic Differential Equations

Nonparametric Drift Estimation for Stochastic Differential Equations Gareth Roberts 1 Department of Statistics University of Warwick Brazilian Bayesian meeting, March 2010 Joint work with O. Papaspiliopoulos,

### Tools for Parameter Estimation and Propagation of Uncertainty

Tools for Parameter Estimation and Propagation of Uncertainty Brian Borchers Department of Mathematics New Mexico Tech Socorro, NM 87801 borchers@nmt.edu Outline Models, parameters, parameter estimation,

### The Basic Properties of Surface Waves

The Basic Properties of Surface Waves Lapo Boschi lapo@erdw.ethz.ch April 24, 202 Love and Rayleigh Waves Whenever an elastic medium is bounded by a free surface, coherent waves arise that travel along

### Estimation of Operational Risk Capital Charge under Parameter Uncertainty

Estimation of Operational Risk Capital Charge under Parameter Uncertainty Pavel V. Shevchenko Principal Research Scientist, CSIRO Mathematical and Information Sciences, Sydney, Locked Bag 17, North Ryde,

Data Using a Facies Based Bayesian Seismic Inversion, Forties Field, UKCS Kester Waters* (Ikon Science Ltd), Ana Somoza (Ikon Science Ltd), Grant Byerley (Apache Corp), Phil Rose (Apache UK) Summary The

### The Kalman Filter ImPr Talk

The Kalman Filter ImPr Talk Ged Ridgway Centre for Medical Image Computing November, 2006 Outline What is the Kalman Filter? State Space Models Kalman Filter Overview Bayesian Updating of Estimates Kalman

### Multilevel Statistical Models: 3 rd edition, 2003 Contents

Multilevel Statistical Models: 3 rd edition, 2003 Contents Preface Acknowledgements Notation Two and three level models. A general classification notation and diagram Glossary Chapter 1 An introduction

### Computer Practical: Metropolis-Hastings-based MCMC

Computer Practical: Metropolis-Hastings-based MCMC Andrea Arnold and Franz Hamilton North Carolina State University July 30, 2016 A. Arnold / F. Hamilton (NCSU) MH-based MCMC July 30, 2016 1 / 19 Markov

### x. Figure 1: Examples of univariate Gaussian pdfs N (x; µ, σ 2 ).

.8.6 µ =, σ = 1 µ = 1, σ = 1 / µ =, σ =.. 3 1 1 3 x Figure 1: Examples of univariate Gaussian pdfs N (x; µ, σ ). The Gaussian distribution Probably the most-important distribution in all of statistics

### Frequency dependent AVO analysis

Frequency dependent AVO analysis of P-, S-, and C-wave elastic and anelastic reflection data Kris Innanen k.innanen@ucalgary.ca CREWES 23 rd Annual Sponsor s Meeting, Banff AB, Nov 30-Dec 2, 2011 Outline

### An empirical study of hydrocarbon indicators

An empirical study of hydrocarbon indicators Brian Russell 1, Hong Feng, and John Bancroft An empirical study of hydrocarbon indicators 1 Hampson-Russell, A CGGVeritas Company, Calgary, Alberta, brian.russell@cggveritas.com

### Application of the Ensemble Kalman Filter to History Matching

Application of the Ensemble Kalman Filter to History Matching Presented at Texas A&M, November 16,2010 Outline Philosophy EnKF for Data Assimilation Field History Match Using EnKF with Covariance Localization

### CHARACTERIZATION OF SATURATED POROUS ROCKS WITH OBLIQUELY DIPPING FRACTURES. Jiao Xue and Robert H. Tatham

CHARACTERIZATION OF SATURATED POROUS ROCS WITH OBLIQUELY DIPPING FRACTURES Jiao Xue and Robert H. Tatham Department of Geological Sciences The University of Texas at Austin ABSTRACT Elastic properties,

### The Monte Carlo Method: Bayesian Networks

The Method: Bayesian Networks Dieter W. Heermann Methods 2009 Dieter W. Heermann ( Methods)The Method: Bayesian Networks 2009 1 / 18 Outline 1 Bayesian Networks 2 Gene Expression Data 3 Bayesian Networks

### Amplitude analysis of isotropic P-wave reflections

GEOPHYSICS, VOL. 71, NO. 6 NOVEMBER-DECEMBER 2006; P. C93 C103, 7 FIGS., 2 TABLES. 10.1190/1.2335877 Amplitude analysis of isotropic P-wave reflections Mirko van der Baan 1 and Dirk Smit 2 ABSTRACT The

Geostatistical Reservoir Characterization of Deepwater Channel, Offshore Malaysia Trisakti Kurniawan* and Jahan Zeb, Petronas Carigali Sdn Bhd, Jimmy Ting and Lee Chung Shen, CGG Summary A quantitative

### STA 4273H: Statistical Machine Learning

STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! http://www.utstat.utoronto.ca/~rsalakhu/ Sidney Smith Hall, Room 6002 Lecture 7 Approximate

### Comparing Non-informative Priors for Estimation and Prediction in Spatial Models

Environmentrics 00, 1 12 DOI: 10.1002/env.XXXX Comparing Non-informative Priors for Estimation and Prediction in Spatial Models Regina Wu a and Cari G. Kaufman a Summary: Fitting a Bayesian model to spatial

### Using Estimating Equations for Spatially Correlated A

Using Estimating Equations for Spatially Correlated Areal Data December 8, 2009 Introduction GEEs Spatial Estimating Equations Implementation Simulation Conclusion Typical Problem Assess the relationship

### Bayesian Inference. Chapter 4: Regression and Hierarchical Models

Bayesian Inference Chapter 4: Regression and Hierarchical Models Conchi Ausín and Mike Wiper Department of Statistics Universidad Carlos III de Madrid Master in Business Administration and Quantitative

### Numerical Modeling for Different Types of Fractures

umerical Modeling for Different Types of Fractures Xiaoqin Cui* CREWES Department of Geoscience University of Calgary Canada xicui@ucalgary.ca and Laurence R. Lines Edward S. Krebes Department of Geoscience

### Dimension-Independent likelihood-informed (DILI) MCMC

Dimension-Independent likelihood-informed (DILI) MCMC Tiangang Cui, Kody Law 2, Youssef Marzouk Massachusetts Institute of Technology 2 Oak Ridge National Laboratory 2 August 25 TC, KL, YM DILI MCMC USC

### DAG models and Markov Chain Monte Carlo methods a short overview

DAG models and Markov Chain Monte Carlo methods a short overview Søren Højsgaard Institute of Genetics and Biotechnology University of Aarhus August 18, 2008 Printed: August 18, 2008 File: DAGMC-Lecture.tex

### Impact of Rock Physics Depth Trends and Markov Random. Fields on Hierarchical Bayesian Lithology/Fluid Prediction.

Impact of Rock Physics Depth Trends and Markov Random Fields on Hierarchical Bayesian Lithology/Fluid Prediction. Kjartan Rimstad 1 & Henning Omre 1 1 Norwegian University of Science and Technology, Trondheim,

### Bayesian Inference: Principles and Practice 3. Sparse Bayesian Models and the Relevance Vector Machine

Bayesian Inference: Principles and Practice 3. Sparse Bayesian Models and the Relevance Vector Machine Mike Tipping Gaussian prior Marginal prior: single α Independent α Cambridge, UK Lecture 3: Overview

### Bertrand Six, Olivier Colnard, Jean-Philippe Coulon and Yasmine Aziez CGGVeritas Frédéric Cailly, Total

4-D Seismic Inversion: A Case Study Offshore Congo Bertrand Six, Olivier Colnard, Jean-Philippe Coulon and Yasmine Aziez CGGVeritas Frédéric Cailly, Total Summary The first 4D seismic survey in Congo was

### Multivariate spatial modeling

Multivariate spatial modeling Point-referenced spatial data often come as multivariate measurements at each location Chapter 7: Multivariate Spatial Modeling p. 1/21 Multivariate spatial modeling Point-referenced