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

Size: px
Start display at page:

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

Transcription

1 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 the spatial distribution of petrophysical properties of interest, such as facies, porosity, mineral volumes and fluids. Usually this is a two-step approach where the petrophysical properties of interest are derived from inverted elastic models. This sequential approach does not ensure the propagation of the uncertainty related to the seismic inversion problem into the resulting rock property models. This problem can be tackled by inverting the seismic reflection data directly to petrophysical properties (e.g. porosity, volume of shale and water saturation) ensuring the propagation of the uncertainty and measurement errors into the estimated subsurface models. The purpose of this work is to invert seismic reflection data directly to petrophysical properties, to properly propagate the uncertainty related to the seismic inversion problem and measurement errors into the estimated subsurface elastic models. It is presented a novel methodology that combines rock-physics models and a stochastic inversion with global perturbation method, that can quantify the relationship between geologic processes and the corresponding geophysical signatures. This method has been tested on real well-log data and partially stacked seismic data. The application to a real reservoir converges towards the real seismic data and provided realistic petrophysical models and facies volume with the corresponding elastic models retrieved from the rock-physics modeling process. Keywords: Geostatistical Seismic Inversion, Rock-Physics Models, Facies, Reservoir Modeling INTRODUCTION Reservoir s performance is directly related to the natural heterogeneities of the subsurface geology. Within the exploration and production stages, reservoir modeling plays a crucial role in the assessment of the productive zones. The main goal of reservoir characterization is to identify the spatial distribution of the petrophysical properties of interest, such as facies, porosity, mineral volumes and fluid. Frequently, petrophysical properties are derived from inverted elastic models in a two-step approach. When inverting seismic reflection data directly to petrophysical properties, the uncertainty related to the seismic inversion problem and measurement errors, is propagated into the estimated subsurface models. The integration of rock-physics modeling within the inversion loop allows linking the inverted subsurface rock properties with the corresponding elastic response. The inference of petrophysical properties is based on a perturbation technique that performs a stochastic sequential simulation (DSS; Soares 2001) and cosimulation with joint probability distributions (Horta and Soares 2010) of the model parameter space, ensuring the reproduction of the prior probability distributions, honouring the data values at each location, reproducing the original statistics (mean and variance), as well as reproducing the spatial continuity pattern imposed by the variogram model. Contrary to Sequential Gaussian Simulation (SGS; Deutsch and Journel 1998), the use of DSS allows the distribution of the property to be directly simulated without any transforms, as estimated from the experimental data (i.e. well-log data). The selection of the first property to simulate should consider the quality of the available well-log data. As a bestpractice, this property should be the one associated with a larger uncertainty and smoother. Each set of simulated and co-simulated petrophysical properties generate a facies volume with a Bayesian 1 Petroleum Engineering MSc Candidate, Center for Modeling Petroleum Reservoirs, CERENA/DECivil, Instituto Superior Técnico, Universidade Técnica de Lisboa, Lisbon, Portugal. catarina.amaro@tecnico.ulisboa.pt

2 classification. A key input of the inversion for petrophysical properties and facies, are the prior facies proportions. These proportions are normally estimated from the available well-log data before the inversion procedure and propagate the corresponding uncertainty in zones away from the wells. In most cases (e.g. Grana and Della Rossa 2010) there are three main litho-fluid classes: shale, brine sand and oil sand. To constrain even more the reservoir, other sub-categories can be used, such as the stiffness of the mineral material. In simpler cases, it is easy to identify clusters to constrain the prior probabilities, but when dealing with a complex reservoir more properties can helpful. In cases where the training data, constructed from the well-log data is statistically representative of the reservoir conditions, Bayesian classification is a successful application to classify facies. On the other hand, when few wells are available to represent the lithologies and fluid types, a useful method is to increase the training data with Monte Carlo Simulation. The resulting facies models represent the link between the rock properties and the real subsurface geology and are conditioned not only to the available seismic reflection data but also to the existing well-log data. The generated set of models composed by water saturation, porosity, volume of shale and facies volume are used as input in an ensemble of faciesdependent rock-physics models (RPMs) to predict the seismic velocities (upper and lower bounds of seismic velocities) of a rock and/or facies. RPMs are used to link data from different domains, i.e to improve coherency between the subsurface rock properties and elastic properties. They can be represented by a simple regression based on well-log data or a complex physical model, with a number of elastic parameters to be estimated e.g., elastic moduli of matrix and fluid components, critical porosity, aspect ratio, and/or coordination number (Mavko et al. 2009). There are two fundamental tasks for this method to be consistent: a well-log calibration of the rock-physics models and quality control. Then, the main procedure is to ajust a theoretical model to a trend in the data. The solid phase, is the mineral part of a composite made of the mineral frame and the pore fluid (Dvorkin et al. 2014). Granular media models describe the rock as a collection of separate gains that contact between them with a certain stiffness. This model is usually applicable to sandstones (Simm and Bacon 2014) and based on Hertz-Mindlin contact theory (Dvorkin et al. 1994). A very effective approach starts by the definition of the elastic properties the end-members. At zero porosity, the rock must have the properties of the mineral, and at the high porosity, the elastic contact theory determines the elastic properties. The interpolation of this two end-members is based on upper and lower Hashin-Shtrikman bounds. The upper bound is usually associated with the contact cement (stiff-sand model) and the uncemented or low cemented rock (soft-sand model) is represented by the lower bound (Mavko et al. 2009). The effect of pore fluid is accounted by using Gassmanns equation (1951) which is the most common model, within this setting, to predict fluid substitution effects at low seismic frequencies (Mavko et al. 2009). The soft- and stiff-sand models for V P and V S can be computed for a wet rock, by first calculating the models for a room-dry grain pack to posteriorly use them in Gassmanns fluid substitution equation. Densities for the matrix and fluid can be computed using Woods (1955) formula. The facies-dependent rock-physics models allow the calculation of velocities and densities for each facies individually, to posteriorly generate synthetic seismograms following Shueys (1985) 3-terms approximation. The quality of the inversion results are assessed based on the match between the calculated synthetic seismic data and the real seismic reflection data with a correlation coefficient. Within global inversion approaches, the highest correlation coefficient portions of the generated models are selected and used along with the corresponding local correlation coefficients, to produce new sets of rock property models in a co-simulation. At each iteration, it is performed an update of the set of rock properties based on the trace-by-trace match between the real and synthetic seismic data. 2

3 METHODOLOGY Figure 1 - Schematic representation of the general framework of the proposed methodology. The proposed iterative geostatistical seismic inversion methodology (Figure 1) integrates rockphysics modelling contains the following main steps. First, DSS simulates water saturation models using the available well-log as experimental data and a spatial continuity pattern as revealed by a variogram model. Then, co-dss with joint probability distributions is used to co-simulate porosity models given the previously simulated water saturation models and the available well-log data. The simulation of volume of shale models also recurs to co-dss with joint probability distributions, given the water saturation previously inverted model and the available well-log data. Before applying a classification algorithm, a training data is generated from the well-log data. The facies of interest are identified in a petrophysical domain such as water saturation versus porosity and porosity versus volume of shale, resulting in five facies: stiffshale, soft-shale, stiff-brine sand, soft-brine sand and soft-oil sand. Then a Bayesian classification algorithm uses the previously simulated models and the training data to create a facies volume. A facies dependent rock-physics modelling uses the set of three simulated models (water saturation, porosity, volume of shale) and the facies volume to compute V P, V S and density. The resulting P- and S- velocity models along with density are used to compute angledependent reflection coefficients which are then convolved with angle-dependent wavelets to generate synthetic seismic data, following Shueys' linear approximation After generating the partially stacked synthetic seismic reflection data each synthetic 3 seismic trace is individually compared against the corresponding real seismic trace in terms of correlation coefficient. At each iteration, and from the ensemble of rock properties generated during the first step of the proposed algorithm, the portions of these models that ensure the highest correlation coefficient between real and synthetic seismic for all angles are simultaneously selected with the correlation coefficients. The selection procedure is based on cross-over genetic algorithm where the best genes (portions of the petrophysical models from different realizations that ensure the highest correlation coefficient) of each iteration are then used as seed for the generation of a new family of models during the next iterations. Iterate the entire procedure until a given global correlation coefficient between the angle-dependent synthetic and real seismic data is above a certain threshold. REAL DATA EXAMPLE The available dataset comprises four partial angle stacks, with mean reflection angles of 9º, 15º, 21º and 27º and their corresponding dependent wavelets and a set of four well-logs composed by porosity, volume of shale, water saturation and P- and S-Impedance. The inversion grid has 159x419x128 cells in i, j and k directions respectively. The grid cell size was defined to reproduce the original inline and cross-line spacing and the original seismic sampling interval, 2 ms. The joint distributions between water saturation versus porosity and porosity versus volume of shale are used as conditioning data in the direct sequential co-simulation with joint probability distributions (Horta and Soares 2010). In this way, the inverted petro-elastic models mimic the real reservoir conditions by reproducing the relationships between the primary and secondary variables (Figure 3). The spatial continuity pattern of each property individually is imposed by a variogram model, estimated for all properties in study for the vertical direction. These variograms are part of the geostatistical inversion procedure and are used as

4 conditioning data for the stochastic sequential simulation of the petro-elastic properties of interest. As a geostatistical inversion procedure, each model reproduces the well-log data at its locations, the variogram model imposed during the stochastic sequential simulation and the marginal and joint probability distributions as inferred from the available well-log data. The rock property models (water saturation, porosity and volume of shale; Figure 2) agree with the main structures as interpreted from the original seismic data and from previous inversion studies over this reservoir. The interpretation of the mean model of the ensemble of petrophysical models simulated in the last iteration (Figure 2) is one way of interpreting the results of geostatistical seismic inversion. However, it is important to highlight that the individual realizations have higher variability presenting small- and large-scale details. In the inverted models of water saturation, porosity and volume of shale, the well-log data values are reproduced at well locations, as well as the histograms of each property. The proposed methodology also allows the assessment of the uncertainty of each property individually by for example, the variance model at a given iteration. Notice that assessing the uncertainty related with each property individually is of great importance for better reservoir modeling. Usually, these areas are associated with low signal-to-noise ratio, i.e., areas with higher uncertainty, seismic and the inverted models do not match with the seismic. Areas with high variability are related to more uncertainty regarding the model parameters, and usually correspond to areas far from the wells, where fewer data is available. Before the inversion of petrophysical properties (continuous quantities), facies (discrete quantities) were classified in four facies. Each facies must be defined such as it captures the physics of the reservoir s geology. The facies of interest are identified in a petrophysical domain such as water saturation versus porosity and porosity versus volume of shale. First, three main groups are classified based on the well-log data rock properties (Figure 5a and Figure 5c): shale, brine sand and oil sand. This classification is based on assumptions related to the known geology of the reservoir, i.e. a) d) b) c) Figure 3 - Joint distributions from the well-logs (on the left) and from the best-fit models of iteration 6 (on the right) of water saturation versus porosity (a;c) and porosity versus volume of shale (b;d). The joint distributions reproduce the ones estimated from the well-log data. a) Water Saturation b) Porosity c) Volume of Shale Figure 2 - Vertical section mean models and standard deviation computed from the last 32 models retrieved from the last iteration. From top to bottom: water saturation, porosity and volume of shale. 4

5 facies shale has a shale content threshold of 0.3, while the fluid factor was added by water saturation with a threshold of maximum 0.8 for the oil sands. Then, to isolate the sand reservoir, another classification is performed considering rock-physics. This classification is based on the stiffness of the rock and it splits each shale, brine sand and oil sand into two categories, stiff and soft depending on the porosity and shale content. This classification assumes that shale facies belong to the granular media models, such as sands. According to the available well-log data and the geological setting, the reservoir only has one type of oil sand is present, soft oil sand. A Bayesian classification algorithm is then used to classify all four facies for the entire inversion grid considering the previously simulated petrophysical properties. The resulting facies model reproduce the distributions of water saturation, porosity and volume of shale used as training data (Figure 5b and Figure 5d). The resulting facies volume (Figure 4a) is consistent with the petrophysical models (Figure 2), as well as the structures interpreted in the real seismic. The ensemble of facies models generated during the last iteration may be used to derived facies probability for each facies (Figure 4b), showing the abundance of stiff brine sands where the in the same location where the stiff shale has lower probability of occurrence. The reservoir soft oil sand has high probability of occurrence also in areas where the progradation is located but with less continuity when compared to the predominant facies. It is clear the intercalation between sand and shale facies and the predominance of stiff facies with small-scale details Table 1 - Rock-physics model parameters: density ρ, bulk modulus K, and shear modulus μ of the matrix and fluid components for a Φ c of 0.4, n of 4 and Peff of 20 Mpa. ρ (g/cm 3 ) K (GPa) μ (GPa) which is a consequence of constraining the stochastic simulation and co-simulation to the inverted models to posteriorly generate facies volumes. a) b) Figure 4 -Vertical section extracted from the best-fit facies model (top) and the probability of occurence of soft oil sand (bottom). Where facies 1 corresponds to soft shale, facies 2 to stiff shale, facies 3 to soft brine sand, facies 4 to stiff brine sand and facies 5 to soft oil sand. a) b) c) d) Sand Shale Oil n.a. Water n.a. Figure 5 - Training data for facies modeling using the training data from the wells (on the left); joint distributions from the best-fit inverted models (on the right) of water saturation versus porosity and porosity versus volume of shale. Each joint distribution is color coded by facies. 5

6 b) e) c) f) Figure 6 - Vertical sections of the mean elastic models retrieved from the 32 simulations of the last iteratio (on the left) along with the reproduction of the histogram of the inverted model (filled blue) and well-log data (red). From top to bottom, P- (a) and S- wave velocity (b) and density (c). a) d) The elastic models (Figure 6a, Figure 6b, Figure 6c) calculated from the petrophysical models using the facies dependent rock physics models calibrated with the well-log (Table 1) reproduce the well-log data and the histograms are close to the ones computed from the original well-log data(figure 6d, Figure 6e, Figure 6f). The global correlation coefficient (62%, Figure 8) between the synthetic seismic elastic models is higher than the global correlation coefficient achieved by the mean model (Figure 7). The areas associated with The synthetic seismic reflection data computed from the arithmetic mean elastic models (Figure 9) matches successfully the real partially stacked seismic data. a) Nearstack b) Mid-nearstack c) Midstack d) Mid-farstack Figure 7 Correlation coefficients between the real and synthetic seismic data of all partial angle stacks. From left to right: nearstack, mid-nearstack, midstack and mid-farstack. Figure 8 - Global correlation coefficient evolution per iteration. 6

7 It reproduces the location of the main primary reflections and the amplitude variations versus offset. low degree of uncertainty (low standard deviation) are related to areas where the synthetic seismic converged properly towards the real one, associated with high correlation coefficients. On the other hand, areas with low correlation coefficients are possible related with values not considered in the conditioning dataset, as well as seismic data with low signal-tonoise ratio, very common in real datasets. Within this a) b) Figure 9 Real reflection data (a) and synthetic seismic section of the arithmetic mean of 32 simulations of the last iteration, of the partial angle midstack. inverse methodology, the noisy areas are assigned a higher uncertainty throughout the entire procedure. DISCUSSION After 6 iterations, an the inversion procedure reached a global correlation coefficient between the synthetic seismic and real seismic of 62%. The use of partially stacked seismic data allows the inversion of seismic directly for rock properties to better distinguish lithofluid facies, instead of the traditional acoustic models. All models obtained from the entire iterative geostatistical inversion procedure with rock physics integration, reproduce: the values of the conditioning data at its locations (Figure 2); the joint and marginal distributions of water saturation, porosity and volume of shale and the spatial continuity model of each property imposed within the sequential simulation by a variogram. This method is successful to discriminate lithology parameters and detect hydrocarbon probability of occurrence directly from the petrophysical properties. From the petrophysical properties itself, it is possible to detect some geological features such as the shale trend along the presented horizontal slice, with low porosity and high shale content crossing the horizontal time slice ( Figure 10c) which is also present in Figure 10d, classified mainly as stiff shale a) Water Saturation b) Porosity c) Volume of Shale d) Facies Volume Figure 10 Horizontal time-slices (from left to right) of water saturation (a), porosity (b), volume of shale (c) and facies volume (d) (k=120). 7

8 with soft shale in some portions of that area, showing the clear relationship of rock-physics per facies. The variance models (Figure 2) computed between the models generated during the last iteration show lower variability i.e., lower spatial uncertainty in the volume of shale inverted models, following water saturation and porosity. Usually high values of variance are related to areas where the seismic is noisy and the inverted petrophysical models cannot produce synthetic seismic data that fits the real seismic, or the lack of real petrophysical properties. Most of the spatial uncertainty is related with the location of geological boundaries of facies bodies, because of its similarities with the inverted petrophysical models (Figure 2) and the facies volume (Figure 4). Also, it is important to highlight that the interpretation of reservoir models along with their corresponding uncertainty allows better decision making and risk management. CONCLUSIONS The presented methodology was successfully designed, implemented and applied in a real dataset resulting in a good match between the real and the inverted synthetic seismic data. The novel iterative geostatistical seismic inversion methodology that simultaneously integrates seismic reflection data, well-log data and rock physics models can retrieve directly from the seismic data reliable petrophysical models such as water saturation porosity, volume of shale and facies volumes. The results of the presented inversion method are consistent with other seismic inversions applied to this reservoir, but also add value due to the identification of new geological features and proper fluid/lithology characterization. This approach provides a direct connection between the seismic response and the geological (petrophysical) properties, by the application of Rock Physics Models allowing the propagation of the uncertainty related to the seismic inversion in onestep approach. It is an efficient method to guide and improve qualitative interpretation, as well as avoid ambiguities in seismic interpretation related to fluids/lithology, sand/shale and porosity/saturation. This novel method can be applied to all reservoir where the physical link between the elastic and petrophysical properties can be described by a suitable rock physics model, as well as adapt other rock properties and litho-fluid facies. ACKNOWLEDGEMENTS The authors would like to thanks Schlumberger for the donation of the academic licenses of Petrel and CERENA/IST for supporting this work. CA want to thanks Professor Dario Grana from the University of Wyoming for the valuable input to this work and the Department of Geology and Geophysics for the three months stay. APPENDIX A ROCK-PHYSICS MODELS The soft- and stiff-sand models are based on Hertz- Mindlin grain-contact theory and provide estimation of the bulk and shear moduli of a dry rock, assuming a random pack of identical spherical grains under an effective pressure P, with a certain critical porosity (Φ c ) and a coordination number (Mavko et al. 2009). and 3 2 P K HM = n2 (1 Φ c ) 2 μ mat 18π 2 (1 ν) 2 2 P μ HM = 5 4ν 3 (1 Φ c ) 2 μ mat 10 5ν 3n2 2π 2 (1 ν) 2 A - 1 A - 2 where μ mat is the shear modulus of the solid phase and ν is the grain Poissons ratio. The matrix moduli is calculated using Voigt-Reuss- Hill averages for a matrix with sand and clay materials: 8

9 K mat = 1 2 (V 1 ck c + (1 V c )K s + V c + K ) s K c (1 V c ) and A - 3 μ mat = 1 2 (V 1 cμ c + (1 V c )μ s + V c + μ ) A - 4 s μ c (1 V c ) where V c is the volume of clay, K c, μ c, are the bulk and shear moduli, respectively of the clay and and K s, μ s is the bulk and shear moduli of the sand. The bulk (K HM ) and shear (μ HM ) moduli of a room-dry rock is estimated recurring, for example to Hertz-Mindlin grain-contact theory, under the assumption that the sand frame with a random pack of identical spherical grains is under an effective pressure P, with a certain critical porosity ( Φ c ) and a coordination number (Mavko et al. 2009). and 3 2 P K HM = n2 (1 Φ c ) 2 μ mat 18π 2 (1 ν) 2 μ HM = 5 4ν ν 3n 3 (1 Φ c ) 2 μ2 mat P 2π 2 (1 ν) 2 A - 5 A - 6 where μ mat is the shear modulus of the solid phase and ν is the grain Poissons ratio. For the effective porosity values between zero and the critical porosity, this model that best fits the data (modified lower and/or upper Hashin-Shtrikman bound), interpolates the the two end-members, i.e. the solid-phase elastic moduli (K mat and μ mat ) and the elastic moduli of the dry rock (K HM and μ HM ). At any porosity Φ < Φ c, the main point of the soft connector (modified lower Hashin-Shtrikman bound) is given by the following equations: Φ Φ c K soft = [ K HM μ HM Φ Φ c K mat + 4 ] 4 3 μ 3 μ HM HM Φ Φ c μ soft = [ μ HM ξ μ HM 1 6 ξ HM μ HM, ξ HM = 9K HM + 8μ HM K HM + 2μ HM Φ Φ c μ mat + 4 ] 3 ξ μ HM A - 7 A - 8 For for any porosity Φ > Φ c the modified upper Hashin-Shtrikman bound, or the stiff connector is given by: Φ Φ c K stiff = [ K HM μ mat Φ Φ c K stiff = [ K HM μ mat ξ = 9K mat + 8μ mat K mat + 2μ mat Φ Φ c K mat + 4 ] 4 3 μ 3 μ mat mat A Φ Φ c K mat + 4 ] 4 3 μ 3 μ mat mat A - 10 While density (ρ) is simply the arithmetic average of the various solid and fluid components of the rock (weighted according to their volume fractions) velocity is sensitive (Simm and Bacon 2014). Gassmanns equation (1951) is used to model fluid substitution effects at low seismic frequencies (Mavko et al. 2009). P- and S-wave velocities are estimated using matrix and fluid properties: and K sat = K dry + (1 K 2 dry ) K mat Φ + 1 Φ K fl K mat K dry 2 K mat A - 11 μ sat = μ dry A - 12 From the saturated-rock elastic moduli, velocities can be obtained by; 9

10 V P = K sat μ sat, ρ A - 13 and V S = μ sat ρ. A - 14 REFERENCES Dvorkin, J., Nur, A. and Yin H., Effective properties of cemented granular material: Mechanics of Materials, 18, Dvorkin, J., Gutierrez, M.A., Grana, D., Seismic Reflections of Rock Properties. Cambridge University Press. Gassmann, F., Elasticity of porous media: Uber die elastizitat poroser medien, Vierteljahrsschrift der Naturforschenden Gesselschaft, 96, Grana, D. & Della Rossa, E., Probabilistic petrophysical-properties estimation integrating statistical rock physics with seismic inversion. Geophysics, 75(3), p.o21. Horta, A., & Soares, A., Direct sequential Cosimulation with joint probability distributions. Mathematical Geosciences, 42(3), Mavko, G., Mukerji, T., & Dvorkin, J., The Rock Physics Handbook - Tools for Seismic Analysis of Porous Media (2nd edition). Simm, R., Bacon, M., Seismic Amplitude: An Interpreter s Handbook. Cambridge University Press. Soares, A., Direct Sequential Simulation and Cosimulation. Mathematical Geology, 33(8), Shuey, R. T., 1985, A Simplification of the Zoeppritz equations; Geophysics, Vol. 50, No. 4,

ROCK PHYSICS MODELING FOR LITHOLOGY PREDICTION USING HERTZ- MINDLIN THEORY

ROCK PHYSICS MODELING FOR LITHOLOGY PREDICTION USING HERTZ- MINDLIN THEORY ROCK PHYSICS MODELING FOR LITHOLOGY PREDICTION USING HERTZ- MINDLIN THEORY Ida Ayu PURNAMASARI*, Hilfan KHAIRY, Abdelazis Lotfy ABDELDAYEM Geoscience and Petroleum Engineering Department Universiti Teknologi

More information

Uncertainties in rock pore compressibility and effects on time lapse seismic modeling An application to Norne field

Uncertainties in rock pore compressibility and effects on time lapse seismic modeling An application to Norne field Uncertainties in rock pore compressibility and effects on time lapse seismic modeling An application to Norne field Amit Suman and Tapan Mukerji Department of Energy Resources Engineering Stanford University

More information

Calibration of the petro-elastic model (PEM) for 4D seismic studies in multi-mineral rocks Amini, Hamed; Alvarez, Erick Raciel

Calibration of the petro-elastic model (PEM) for 4D seismic studies in multi-mineral rocks Amini, Hamed; Alvarez, Erick Raciel Heriot-Watt University Heriot-Watt University Research Gateway Calibration of the petro-elastic model (PEM) for 4D seismic studies in multi-mineral rocks Amini, Hamed; Alvarez, Erick Raciel DOI: 10.3997/2214-4609.20132136

More information

Th LHR2 08 Towards an Effective Petroelastic Model for Simulator to Seismic Studies

Th LHR2 08 Towards an Effective Petroelastic Model for Simulator to Seismic Studies Th LHR2 08 Towards an Effective Petroelastic Model for Simulator to Seismic Studies A. Briceno* (Heriot-Watt University), C. MacBeth (Heriot-Watt University) & M.D. Mangriotis (Heriot-Watt University)

More information

Reliability of Seismic Data for Hydrocarbon Reservoir Characterization

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.

More information

Integration of seismic interpretation in to geostatistical acoustic impedance Results presentation and discussion

Integration of seismic interpretation in to geostatistical acoustic impedance Results presentation and discussion ANNUAL MEETING MASTER OF PETROLEUM ENGINEERING Integration of seismic interpretation in to geostatistical acoustic impedance Results presentation and discussion SERGIO CRUZ BARDERA 03/May/2016 Instituto

More information

IDENTIFYING PATCHY SATURATION FROM WELL LOGS Short Note. / K s. + K f., G Dry. = G / ρ, (2)

IDENTIFYING PATCHY SATURATION FROM WELL LOGS Short Note. / K s. + K f., G Dry. = G / ρ, (2) IDENTIFYING PATCHY SATURATION FROM WELL LOGS Short Note JACK DVORKIN, DAN MOOS, JAMES PACKWOOD, AND AMOS NUR DEPARTMENT OF GEOPHYSICS, STANFORD UNIVERSITY January 5, 2001 INTRODUCTION Gassmann's (1951)

More information

ROCK PHYSICS DIAGNOSTICS OF NORTH SEA SANDS: LINK BETWEEN MICROSTRUCTURE AND SEISMIC PROPERTIES ABSTRACT

ROCK PHYSICS DIAGNOSTICS OF NORTH SEA SANDS: LINK BETWEEN MICROSTRUCTURE AND SEISMIC PROPERTIES ABSTRACT ROCK PHYSICS DIAGNOSTICS OF NORTH SEA SANDS: LINK BETWEEN MICROSTRUCTURE AND SEISMIC PROPERTIES PER AVSETH, JACK DVORKIN, AND GARY MAVKO Department of Geophysics, Stanford University, CA 94305-2215, USA

More information

GRAIN SORTING, POROSITY, AND ELASTICITY. Jack Dvorkin and Mario A. Gutierrez Geophysics Department, Stanford University ABSTRACT

GRAIN SORTING, POROSITY, AND ELASTICITY. Jack Dvorkin and Mario A. Gutierrez Geophysics Department, Stanford University ABSTRACT GRAIN SORTING, POROSITY, AND ELASTICITY Jack Dvorkin and Mario A. Gutierrez Geophysics Department, Stanford University July 24, 2001 ABSTRACT Grain size distribution (sorting) is determined by deposition.

More information

Statistical Rock Physics

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

More information

We apply a rock physics analysis to well log data from the North-East Gulf of Mexico

We apply a rock physics analysis to well log data from the North-East Gulf of Mexico Rock Physics for Fluid and Porosity Mapping in NE GoM JACK DVORKIN, Stanford University and Rock Solid Images TIM FASNACHT, Anadarko Petroleum Corporation RICHARD UDEN, MAGGIE SMITH, NAUM DERZHI, AND JOEL

More information

23855 Rock Physics Constraints on Seismic Inversion

23855 Rock Physics Constraints on Seismic Inversion 23855 Rock Physics Constraints on Seismic Inversion M. Sams* (Ikon Science Ltd) & D. Saussus (Ikon Science) SUMMARY Seismic data are bandlimited, offset limited and noisy. Consequently interpretation of

More information

Reservoir connectivity uncertainty from stochastic seismic inversion Rémi Moyen* and Philippe M. Doyen (CGGVeritas)

Reservoir connectivity uncertainty from stochastic seismic inversion Rémi Moyen* and Philippe M. Doyen (CGGVeritas) Rémi Moyen* and Philippe M. Doyen (CGGVeritas) Summary Static reservoir connectivity analysis is sometimes based on 3D facies or geobody models defined by combining well data and inverted seismic impedances.

More information

Earth models for early exploration stages

Earth models for early exploration stages ANNUAL MEETING MASTER OF PETROLEUM ENGINEERING Earth models for early exploration stages Ângela Pereira PhD student angela.pereira@tecnico.ulisboa.pt 3/May/2016 Instituto Superior Técnico 1 Outline Motivation

More information

AFI (AVO Fluid Inversion)

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

More information

Seismic reservoir characterization in offshore Nile Delta.

Seismic reservoir characterization in offshore Nile Delta. Seismic reservoir characterization in offshore Nile Delta. Part II: Probabilistic petrophysical-seismic inversion M. Aleardi 1, F. Ciabarri 2, B. Garcea 2, A. Mazzotti 1 1 Earth Sciences Department, University

More information

QUANTITATIVE INTERPRETATION

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

More information

A New AVO Attribute for Hydrocarbon Prediction and Application to the Marmousi II Dataset*

A New AVO Attribute for Hydrocarbon Prediction and Application to the Marmousi II Dataset* A New AVO Attribute for Hydrocarbon Prediction and Application to the Marmousi II Dataset* Changcheng Liu 1 and Prasad Ghosh 2 Search and Discovery Article #41764 (2016) Posted January 25, 2016 *Adapted

More information

Effects of fluid changes on seismic reflections: Predicting amplitudes at gas reservoir directly from amplitudes at wet reservoir

Effects of fluid changes on seismic reflections: Predicting amplitudes at gas reservoir directly from amplitudes at wet reservoir GEOPHYSICS, VOL. 77, NO. 4 (JULY-AUGUST 2012); P. D129 D140, 15 FIGS., 2 TABLES. 10.1190/GEO2011-0331.1 Effects of fluid changes on seismic reflections: Predicting amplitudes at gas reservoir directly

More information

Stochastic vs Deterministic Pre-stack Inversion Methods. Brian Russell

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

More information

Geostatistical History Matching coupled with Adaptive Stochastic Sampling: A zonation-based approach using Direct Sequential Simulation

Geostatistical History Matching coupled with Adaptive Stochastic Sampling: A zonation-based approach using Direct Sequential Simulation Geostatistical History Matching coupled with Adaptive Stochastic Sampling: A zonation-based approach using Direct Sequential Simulation Eduardo Barrela* Instituto Superior Técnico, Av. Rovisco Pais 1,

More information

Net-to-gross from Seismic P and S Impedances: Estimation and Uncertainty Analysis using Bayesian Statistics

Net-to-gross from Seismic P and S Impedances: Estimation and Uncertainty Analysis using Bayesian Statistics Net-to-gross from Seismic P and S Impedances: Estimation and Uncertainty Analysis using Bayesian Statistics Summary Madhumita Sengupta*, Ran Bachrach, Niranjan Banik, esterngeco. Net-to-gross (N/G ) is

More information

Heriot-Watt University

Heriot-Watt University Heriot-Watt University Heriot-Watt University Research Gateway 4D seismic feasibility study for enhanced oil recovery (EOR) with CO2 injection in a mature North Sea field Amini, Hamed; Alvarez, Erick Raciel;

More information

Integrating rock physics modeling, prestack inversion and Bayesian classification. Brian Russell

Integrating rock physics modeling, prestack inversion and Bayesian classification. Brian Russell Integrating rock physics modeling, prestack inversion and Bayesian classification Brian Russell Introduction Today, most geoscientists have an array of tools available to perform seismic reservoir characterization.

More information

Rock Physics Perturbational Modeling: Carbonate case study, an intracratonic basin Northwest/Saharan Africa

Rock Physics Perturbational Modeling: Carbonate case study, an intracratonic basin Northwest/Saharan Africa Rock Physics Perturbational Modeling: Carbonate case study, an intracratonic basin Northwest/Saharan Africa Franklin Ruiz, Carlos Cobos, Marcelo Benabentos, Beatriz Chacon, and Roberto Varade, Luis Gairifo,

More information

Toward an Integrated and Realistic Interpretation of Continuous 4D Seismic Data for a CO 2 EOR and Sequestration Project

Toward an Integrated and Realistic Interpretation of Continuous 4D Seismic Data for a CO 2 EOR and Sequestration Project SPE-183789-MS Toward an Integrated and Realistic Interpretation of Continuous 4D Seismic Data for a CO 2 EOR and Sequestration Project Philippe Nivlet, Robert Smith, Michael A. Jervis, and Andrey Bakulin,

More information

Quantitative interpretation using inverse rock-physics modeling on AVO data

Quantitative interpretation using inverse rock-physics modeling on AVO data Quantitative interpretation using inverse rock-physics modeling on AVO data Erling Hugo Jensen 1, Tor Arne Johansen 2, 3, 4, Per Avseth 5, 6, and Kenneth Bredesen 2, 7 Downloaded 08/16/16 to 129.177.32.62.

More information

Quantitative Interpretation

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

More information

We Density/Porosity Versus Velocity of Overconsolidated Sands Derived from Experimental Compaction SUMMARY

We Density/Porosity Versus Velocity of Overconsolidated Sands Derived from Experimental Compaction SUMMARY We 6 Density/Porosity Versus Velocity of Overconsolidated Sands Derived from Experimental Compaction S. Narongsirikul* (University of Oslo), N.H. Mondol (University of Oslo and Norwegian Geotechnical Inst)

More information

Geological Classification of Seismic-Inversion Data in the Doba Basin of Chad*

Geological Classification of Seismic-Inversion Data in the Doba Basin of Chad* Geological Classification of Seismic-Inversion Data in the Doba Basin of Chad* Carl Reine 1, Chris Szelewski 2, and Chaminda Sandanayake 3 Search and Discovery Article #41899 (2016)** Posted September

More information

Rock physics of a gas hydrate reservoir. Gas hydrates are solids composed of a hydrogen-bonded ROUND TABLE

Rock physics of a gas hydrate reservoir. Gas hydrates are solids composed of a hydrogen-bonded ROUND TABLE ROUND TABLE Rock physics of a gas hydrate reservoir JACK DVORKIN and AMOS NUR, Stanford University, California, U.S. RICHARD UDEN and TURHAN TANER, Rock Solid Images, Houston, Texas, U.S. Gas hydrates

More information

Seismic inversion for reservoir properties combining statistical rock physics and geostatistics: A review

Seismic inversion for reservoir properties combining statistical rock physics and geostatistics: A review GEOPHYSICS, VOL. 75, NO. 5 SEPTEMBER-OCTOBER 2010 ; P. 75A165 75A176, 8 FIGS. 10.1190/1.3478209 Seismic inversion for reservoir properties combining statistical rock physics and geostatistics: A review

More information

GNGTS 2015 Sessione 3.1

GNGTS 2015 Sessione 3.1 seismic reservoir characterization In offshore nile delta. Part I: comparing different methods to derive a reliable rock-physics model M. Aleardi 1, F. Ciabarri 2, F. Peruzzo 2, A. Mazzotti 1 1 Earth Sciences

More information

Sensitivity Analysis of Pre stack Seismic Inversion on Facies Classification using Statistical Rock Physics

Sensitivity Analysis of Pre stack Seismic Inversion on Facies Classification using Statistical Rock Physics Sensitivity Analysis of Pre stack Seismic Inversion on Facies Classification using Statistical Rock Physics Peipei Li 1 and Tapan Mukerji 1,2 1 Department of Energy Resources Engineering 2 Department of

More information

The elastic properties such as velocity, density, impedance,

The elastic properties such as velocity, density, impedance, SPECIAL SECTION: Rr ock Physics physics Lithology and fluid differentiation using rock physics template XIN-GANG CHI AND DE-HUA HAN, University of Houston The elastic properties such as velocity, density,

More information

Downloaded 11/02/16 to Redistribution subject to SEG license or copyright; see Terms of Use at Summary.

Downloaded 11/02/16 to Redistribution subject to SEG license or copyright; see Terms of Use at   Summary. in thin sand reservoirs William Marin* and Paola Vera de Newton, Rock Solid Images, and Mario Di Luca, Pacific Exploración y Producción. Summary Rock Physics Templates (RPTs) are useful tools for well

More information

Integrating reservoir flow simulation with time-lapse seismic inversion in a heavy oil case study

Integrating reservoir flow simulation with time-lapse seismic inversion in a heavy oil case study Integrating reservoir flow simulation with time-lapse seismic inversion in a heavy oil case study Naimeh Riazi*, Larry Lines*, and Brian Russell** Department of Geoscience, University of Calgary **Hampson-Russell

More information

Competing Effect of Pore Fluid and Texture -- Case Study

Competing Effect of Pore Fluid and Texture -- Case Study Competing Effect of Pore Fluid and Texture -- Case Study Depth (m) Sw Sxo. m Poisson's Ratio.. GOC.1 5 7 8 9 P-Impedance OWC 15 GR.. RHOB.5 1 Saturation...5. 1. 1.5 Vs (km/s).. Poisson's Ratio 5 7 P-Impedance

More information

We Simultaneous Joint Inversion of Electromagnetic and Seismic Full-waveform Data - A Sensitivity Analysis to Biot Parameter

We Simultaneous Joint Inversion of Electromagnetic and Seismic Full-waveform Data - A Sensitivity Analysis to Biot Parameter We-09-04 Simultaneous Joint Inversion of Electromagnetic and Seismic Full-waveform Data - A Sensitivity Analysis to Biot Parameter J. Giraud* (WesternGeco Geosolutions), M. De Stefano (WesternGeco Geosolutions)

More information

We LHR3 04 Realistic Uncertainty Quantification in Geostatistical Seismic Reservoir Characterization

We LHR3 04 Realistic Uncertainty Quantification in Geostatistical Seismic Reservoir Characterization We LHR3 04 Realistic Uncertainty Quantification in Geostatistical Seismic Reservoir Characterization A. Moradi Tehrani* (CGG), A. Stallone (Roma Tre University), R. Bornard (CGG) & S. Boudon (CGG) SUMMARY

More information

THE ROCK PHYSICS HANDBOOK

THE ROCK PHYSICS HANDBOOK THE ROCK PHYSICS HANDBOOK TOOLS FOR SEISMIC ANALYSIS IN POROUS MEDIA Gary Mavko Tapan Mukerji Jack Dvorkin Stanford University Stanford University Stanford University CAMBRIDGE UNIVERSITY PRESS CONTENTS

More information

Integrating rock physics and full elastic modeling for reservoir characterization Mosab Nasser and John B. Sinton*, Maersk Oil Houston Inc.

Integrating rock physics and full elastic modeling for reservoir characterization Mosab Nasser and John B. Sinton*, Maersk Oil Houston Inc. Integrating rock physics and full elastic modeling for reservoir characterization Mosab Nasser and John B. Sinton*, Maersk Oil Houston Inc. Summary Rock physics establishes the link between reservoir properties,

More information

Thomas Bayes versus the wedge model: An example inference using a geostatistical prior function

Thomas Bayes versus the wedge model: An example inference using a geostatistical prior function Thomas Bayes versus the wedge model: An example inference using a geostatistical prior function Jason M. McCrank, Gary F. Margrave, and Don C. Lawton ABSTRACT The Bayesian inference is used to estimate

More information

Quantitative Seismic Interpretation An Earth Modeling Perspective

Quantitative Seismic Interpretation An Earth Modeling Perspective Quantitative Seismic Interpretation An Earth Modeling Perspective Damien Thenin*, RPS, Calgary, AB, Canada TheninD@rpsgroup.com Ron Larson, RPS, Calgary, AB, Canada LarsonR@rpsgroup.com Summary Earth models

More information

The Stanford VI reservoir

The Stanford VI reservoir The Stanford VI reservoir Scarlet A. Castro, Jef Caers and Tapan Mukerji Stanford Center for Reservoir Forecasting Stanford Rock Physics and Borehole Geophysics Project May 2005 Abstract An exhaustive

More information

Shaly Sand Rock Physics Analysis and Seismic Inversion Implication

Shaly Sand Rock Physics Analysis and Seismic Inversion Implication Shaly Sand Rock Physics Analysis and Seismic Inversion Implication Adi Widyantoro (IkonScience), Matthew Saul (IkonScience/UWA) Rock physics analysis of reservoir elastic properties often assumes homogeneity

More information

Downloaded 10/02/18 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 10/02/18 to Redistribution subject to SEG license or copyright; see Terms of Use at Multi-scenario, multi-realization seismic inversion for probabilistic seismic reservoir characterization Kester Waters* and Michael Kemper, Ikon Science Ltd. Summary We propose a two tiered inversion strategy

More information

SPE These in turn can be used to estimate mechanical properties.

SPE These in turn can be used to estimate mechanical properties. SPE 96112 Pressure Effects on Porosity-Log Responses Using Rock Physics Modeling: Implications on Geophysical and Engineering Models as Reservoir Pressure Decreases Michael Holmes, SPE, Digital Formation,

More information

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 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

More information

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

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

More information

A021 Petrophysical Seismic Inversion for Porosity and 4D Calibration on the Troll Field

A021 Petrophysical Seismic Inversion for Porosity and 4D Calibration on the Troll Field A021 Petrophysical Seismic Inversion for Porosity and 4D Calibration on the Troll Field T. Coleou* (CGG), A.J. van Wijngaarden (Hydro), A. Norenes Haaland (Hydro), P. Moliere (Hydro), R. Ona (Hydro) &

More information

A look into Gassmann s Equation

A look into Gassmann s Equation A look into Gassmann s Equation Nawras Al-Khateb, CHORUS Heavy Oil Consortium, Department of Geoscience, University of Calgary nawras.alkhateb@ucalgary.ca Summary By describing the influence of the pore

More information

Pressure and Compaction in the Rock Physics Space. Jack Dvorkin

Pressure and Compaction in the Rock Physics Space. Jack Dvorkin Pressure and Compaction in the Rock Physics Space Jack Dvorkin June 2002 0 200 Compaction of Shales Freshly deposited shales and clays may have enormous porosity of ~ 80%. The speed of sound is close to

More information

Lithology prediction and fluid discrimination in Block A6 offshore Myanmar

Lithology prediction and fluid discrimination in Block A6 offshore Myanmar 10 th Biennial International Conference & Exposition P 141 Lithology prediction and fluid discrimination in Block A6 offshore Myanmar Hanumantha Rao. Y *, Loic Michel, Hampson-Russell, Kyaw Myint, Ko Ko,

More information

Seismic reservoir and source-rock analysis using inverse rock-physics modeling: A Norwegian Sea demonstration

Seismic reservoir and source-rock analysis using inverse rock-physics modeling: A Norwegian Sea demonstration 66 Seismic reservoir and source-rock analysis using inverse rock-physics modeling: A Norwegian Sea demonstration Kenneth Bredesen 1, Erling Hugo Jensen 1, 2, Tor Arne Johansen 1, 2, and Per Avseth 3, 4

More information

Measurement of elastic properties of kerogen Fuyong Yan, De-hua Han*, Rock Physics Lab, University of Houston

Measurement of elastic properties of kerogen Fuyong Yan, De-hua Han*, Rock Physics Lab, University of Houston Measurement of elastic properties of kerogen Fuyong Yan, De-hua Han*, Rock Physics Lab, University of Houston Summary To have good understanding of elastic properties of organic shale, it is fundamental

More information

Reservoir Characterization using AVO and Seismic Inversion Techniques

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

More information

Estimating the hydrocarbon volume from elastic and resistivity data: A concept

Estimating the hydrocarbon volume from elastic and resistivity data: A concept INTERPRETER S CORNER Coordinated by Rebecca B. Latimer Estimating the hydrocarbon volume from elastic and resistivity data: A concept CARMEN T. GOMEZ, JACK DVORKIN, and GARY MAVKO, Stanford University,

More information

SEG Houston 2009 International Exposition and Annual Meeting. that the project results can correctly interpreted.

SEG Houston 2009 International Exposition and Annual Meeting. that the project results can correctly interpreted. Calibration of Pre-Stack Simultaneous Impedance Inversion using Rock Physics Scott Singleton and Rob Keirstead, Rock Solid Images Log Conditioning and Rock Physics Modeling Summary Geophysical Well Log

More information

P- and S-Wave Velocity Measurements and Pressure Sensitivity Analysis of AVA Response

P- and S-Wave Velocity Measurements and Pressure Sensitivity Analysis of AVA Response P- and S-Wave Velocity Measurements and Pressure Sensitivity Analysis of AVA Response Tiewei He* University of Alberta, Edmonton, Alberta, Canada tieweihe@phys.ualberta.ca and Douglas Schmitt University

More information

A010 MULTISCALE RESERVOIR CHARACTERIZATION USING

A010 MULTISCALE RESERVOIR CHARACTERIZATION USING 1 A010 MULTISCALE RESERVOIR CHARACTERIZATION USING RODUCTION AND TIME LASE SEISMIC DATA Mokhles MEZGHANI, Alexandre FORNEL, Valérie LANGLAIS, Nathalie LUCET IF, 1 & 4 av de Bois réau, 92852 RUEIL-MALMAISON

More information

SRC software. Rock physics modelling tools for analyzing and predicting geophysical reservoir properties

SRC software. Rock physics modelling tools for analyzing and predicting geophysical reservoir properties SRC software Rock physics modelling tools for analyzing and predicting geophysical reservoir properties Outline About SRC software. Introduction to rock modelling. Rock modelling program structure. Examples

More information

Methane hydrate rock physics models for the Blake Outer Ridge

Methane hydrate rock physics models for the Blake Outer Ridge Stanford Exploration Project, Report 80, May 15, 2001, pages 1 307 Methane hydrate rock physics models for the Blake Outer Ridge Christine Ecker 1 ABSTRACT Seismic analyses of methane hydrate data from

More information

Sequential Simulations of Mixed Discrete-Continuous Properties: Sequential Gaussian Mixture Simulation

Sequential Simulations of Mixed Discrete-Continuous Properties: Sequential Gaussian Mixture Simulation Sequential Simulations of Mixed Discrete-Continuous Properties: Sequential Gaussian Mixture Simulation Dario Grana, Tapan Mukerji, Laura Dovera, and Ernesto Della Rossa Abstract We present here a method

More information

Downloaded 10/25/16 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 10/25/16 to Redistribution subject to SEG license or copyright; see Terms of Use at Facies modeling in unconventional reservoirs using seismic derived facies probabilities Reinaldo J. Michelena*, Omar G. Angola, and Kevin S. Godbey, ireservoir.com, Inc. Summary We present in this paper

More information

Evaluation of Rock Properties from Logs Affected by Deep Invasion A Case Study

Evaluation of Rock Properties from Logs Affected by Deep Invasion A Case Study Evaluation of Rock Properties from Logs Affected by Deep Invasion A Case Study Jahan Zeb a, Reece Murrell b a CGG Australia, 1 Ord Street, West Perth, WA 6005 Contact email: Jahan.Zeb@cgg.com b Esso Australia,

More information

SEISMIC INVERSION OVERVIEW

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,

More information

3D petrophysical modeling - 1-3D Petrophysical Modeling Usning Complex Seismic Attributes and Limited Well Log Data

3D petrophysical modeling - 1-3D Petrophysical Modeling Usning Complex Seismic Attributes and Limited Well Log Data - 1-3D Petrophysical Modeling Usning Complex Seismic Attributes and Limited Well Log Data Mehdi Eftekhari Far, and De-Hua Han, Rock Physics Lab, University of Houston Summary A method for 3D modeling and

More information

Fluid-property discrimination with AVO: A Biot-Gassmann perspective

Fluid-property discrimination with AVO: A Biot-Gassmann perspective Fluid-property discrimination with AVO: A Biot-Gassmann perspective Brian H. Russell, Ken Hedlin 1, Fred J. Hilterman, and Laurence R. Lines ABSTRACT This paper draws together basic rock physics, AVO,

More information

Geostatistical History Matching coupled with Adaptive Stochastic Sampling

Geostatistical History Matching coupled with Adaptive Stochastic Sampling Geostatistical History Matching coupled with Adaptive Stochastic Sampling A Geologically consistent approach using Stochastic Sequential Simulation Eduardo Barrela Nº 79909 Project Thesis Presentation

More information

Rock Physics and Quantitative Wavelet Estimation. for Seismic Interpretation: Tertiary North Sea. R.W.Simm 1, S.Xu 2 and R.E.

Rock Physics and Quantitative Wavelet Estimation. for Seismic Interpretation: Tertiary North Sea. R.W.Simm 1, S.Xu 2 and R.E. Rock Physics and Quantitative Wavelet Estimation for Seismic Interpretation: Tertiary North Sea R.W.Simm 1, S.Xu 2 and R.E.White 2 1. Enterprise Oil plc, Grand Buildings, Trafalgar Square, London WC2N

More information

Time-lapse seismic monitoring and inversion in a heavy oilfield. By: Naimeh Riazi PhD Student, Geophysics

Time-lapse seismic monitoring and inversion in a heavy oilfield. By: Naimeh Riazi PhD Student, Geophysics Time-lapse seismic monitoring and inversion in a heavy oilfield By: Naimeh Riazi PhD Student, Geophysics May 2011 Contents Introduction on time-lapse seismic data Case study Rock-physics Time-Lapse Calibration

More information

Rock Physics Interpretation of microstructure Chapter Jingqiu Huang M.S. Candidate University of Houston

Rock Physics Interpretation of microstructure Chapter Jingqiu Huang M.S. Candidate University of Houston Rock Physics Interpretation of microstructure Chapter2.1 2.2 2.3 Jingqiu Huang M.S. Candidate University of Houston Introduction Theory and models Example in North Sea Introduction Theoretical models Inclusion

More information

Reservoir Characteristics of a Quaternary Channel: Incorporating Rock Physics in Seismic and DC Resistivity Surveys

Reservoir Characteristics of a Quaternary Channel: Incorporating Rock Physics in Seismic and DC Resistivity Surveys Reservoir Characteristics of a Quaternary Channel: Incorporating Rock Physics in Seismic and DC Resistivity Surveys Jawwad Ahmad* University of Alberta, Edmonton, Alberta, Canada jahmad@phys.ualberta.ca

More information

Fluid Substitution Model to Generate Synthetic Seismic Attributes: FluidSub.exe

Fluid Substitution Model to Generate Synthetic Seismic Attributes: FluidSub.exe Fluid Substitution Model to Generate Synthetic Seisic Attributes: FluidSub.exe Sahyun Hong and Clayton V. Deutsch Geostatistical data integration is a ature research field and any of algoriths have been

More information

Downloaded 11/20/12 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 11/20/12 to Redistribution subject to SEG license or copyright; see Terms of Use at AVO crossplot analysis in unconsolidated sediments containing gas hydrate and free gas: Green Canyon 955, Gulf of Mexico Zijian Zhang* 1, Daniel R. McConnell 1, De-hua Han 2 1 Fugro GeoConsulting, Inc.,

More information

Uncertainty analysis for the integration of seismic and CSEM data Myoung Jae Kwon & Roel Snieder, Center for Wave Phenomena, Colorado School of Mines

Uncertainty analysis for the integration of seismic and CSEM data Myoung Jae Kwon & Roel Snieder, Center for Wave Phenomena, Colorado School of Mines Myoung Jae Kwon & Roel Snieder, Center for Wave Phenomena, Colorado School of Mines Summary Geophysical inverse problems consist of three stages: the forward problem, optimization, and appraisal. We study

More information

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

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.

More information

Practical Gassmann fluid substitution in sand/shale sequences

Practical Gassmann fluid substitution in sand/shale sequences first break volume 25, December 2007 tutorial Practical Gassmann fluid substitution in sand/shale sequences Rob Simm * Introduction When performing fluid substitution on log data Gassmann s (1951) model

More information

Fifteenth International Congress of the Brazilian Geophysical Society. Copyright 2017, SBGf - Sociedade Brasileira de Geofísica

Fifteenth International Congress of the Brazilian Geophysical Society. Copyright 2017, SBGf - Sociedade Brasileira de Geofísica Geostatistical Reservoir Characterization in Barracuda Field, Campos Basin: A case study Frank Pereira (CGG)*, Ted Holden (CGG), Mohammed Ibrahim (CGG) and Eduardo Porto (CGG). Copyright 2017, SBGf - Sociedade

More information

Researcher 2015;7(9)

Researcher 2015;7(9) 4D Seismic Feasibility Study using well Logs in Sienna gas Field, West Delta Deep Marine concession, Egypt Helal, A., Shebl, A. 1, ElNaggar, S. 2 and Ezzat, A. 3 1 Faculty of Science, Ain Shams University,

More information

P314 Anisotropic Elastic Modelling for Organic Shales

P314 Anisotropic Elastic Modelling for Organic Shales P314 Anisotropic Elastic Modelling for Organic Shales X. Wu* (British Geological Survey), M. Chapman (British Geological Survey), X.Y. Li (British Geological Survey) & H. Dai (British Geological Survey)

More information

Rock Physics & Formation Evaluation. Special Topic

Rock Physics & Formation Evaluation. Special Topic Volume 30 Issue 5 May 2012 Special Topic Technical Articles Dual representation of multiscale fracture network modelling for UAE carbonate field AVO and spectral decomposition for derisking Palaeogene

More information

Key Laboratory of Geo-detection (China University of Geosciences, Beijing), Ministry of Education, Beijing , China

Key Laboratory of Geo-detection (China University of Geosciences, Beijing), Ministry of Education, Beijing , China Pet.Sci.(20118:43-48 DOI 10.1007/s12182-011-0113-5 43 Experimental study of the relationship between fluid density and saturation and sonic wave velocity of rock samples from the WXS Depression, South

More information

Simultaneous Inversion of Clastic Zubair Reservoir: Case Study from Sabiriyah Field, North Kuwait

Simultaneous Inversion of Clastic Zubair Reservoir: Case Study from Sabiriyah Field, North Kuwait Simultaneous Inversion of Clastic Zubair Reservoir: Case Study from Sabiriyah Field, North Kuwait Osman Khaled, Yousef Al-Zuabi, Hameed Shereef Summary The zone under study is Zubair formation of Cretaceous

More information

Geostatistics for Seismic Data Integration in Earth Models

Geostatistics for Seismic Data Integration in Earth Models 2003 Distinguished Instructor Short Course Distinguished Instructor Series, No. 6 sponsored by the Society of Exploration Geophysicists European Association of Geoscientists & Engineers SUB Gottingen 7

More information

Porosity. Downloaded 09/22/16 to Redistribution subject to SEG license or copyright; see Terms of Use at

Porosity. Downloaded 09/22/16 to Redistribution subject to SEG license or copyright; see Terms of Use at 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

More information

RESEARCH PROPOSAL. Effects of scales and extracting methods on quantifying quality factor Q. Yi Shen

RESEARCH PROPOSAL. Effects of scales and extracting methods on quantifying quality factor Q. Yi Shen RESEARCH PROPOSAL Effects of scales and extracting methods on quantifying quality factor Q Yi Shen 2:30 P.M., Wednesday, November 28th, 2012 Shen 2 Ph.D. Proposal ABSTRACT The attenuation values obtained

More information

Integration of seismic and fluid-flow data: a two-way road linked by rock physics

Integration of seismic and fluid-flow data: a two-way road linked by rock physics Integration of seismic and fluid-flow data: a two-way road linked by rock physics Abstract Yunyue (Elita) Li, Yi Shen, and Peter K. Kang Geologic model building of the subsurface is a complicated and lengthy

More information

Porosity prediction using cokriging with multiple secondary datasets

Porosity prediction using cokriging with multiple secondary datasets Cokriging with Multiple Attributes Porosity prediction using cokriging with multiple secondary datasets Hong Xu, Jian Sun, Brian Russell, Kris Innanen ABSTRACT The prediction of porosity is essential for

More information

IJMGE Int. J. Min. & Geo-Eng. Vol.49, No.1, June 2015, pp

IJMGE Int. J. Min. & Geo-Eng. Vol.49, No.1, June 2015, pp IJMGE Int. J. Min. & Geo-Eng. Vol.49, No.1, June 2015, pp.131-142 Joint Bayesian Stochastic Inversion of Well Logs and Seismic Data for Volumetric Uncertainty Analysis Moslem Moradi 1, Omid Asghari 1,

More information

Rock physics and AVO applications in gas hydrate exploration

Rock physics and AVO applications in gas hydrate exploration Rock physics and AVO applications in gas hydrate exploration ABSTRACT Yong Xu*, Satinder Chopra Core Lab Reservoir Technologies Division, 301,400-3rd Ave SW, Calgary, AB, T2P 4H2 yxu@corelab.ca Summary

More information

3D geostatistical porosity modelling: A case study at the Saint-Flavien CO 2 storage project

3D geostatistical porosity modelling: A case study at the Saint-Flavien CO 2 storage project 3D geostatistical porosity modelling: A case study at the Saint-Flavien CO 2 storage project Maxime Claprood Institut national de la recherche scientifique, Québec, Canada Earth Modelling 2013 October

More information

THE USE OF SEISMIC ATTRIBUTES AND SPECTRAL DECOMPOSITION TO SUPPORT THE DRILLING PLAN OF THE URACOA-BOMBAL FIELDS

THE USE OF SEISMIC ATTRIBUTES AND SPECTRAL DECOMPOSITION TO SUPPORT THE DRILLING PLAN OF THE URACOA-BOMBAL FIELDS THE USE OF SEISMIC ATTRIBUTES AND SPECTRAL DECOMPOSITION TO SUPPORT THE DRILLING PLAN OF THE URACOA-BOMBAL FIELDS Cuesta, Julián* 1, Pérez, Richard 1 ; Hernández, Freddy 1 ; Carrasquel, Williams 1 ; Cabrera,

More information

C002 Petrophysical Seismic Inversion over an Offshore Carbonate Field

C002 Petrophysical Seismic Inversion over an Offshore Carbonate Field C002 Petrophysical Seismic Inversion over an Offshore Carbonate Field T. Coleou* (CGGVeritas), F. Allo (CGGVeritas), O. Colnard (CGGVeritas), I. Machecler (CGGVeritas), L. Dillon (Petrobras), G. Schwedersky

More information

Probabilistic Inversion Technique for Seismic Data

Probabilistic Inversion Technique for Seismic Data Delft University of Technology Faculty of Electrical Engineering, Mathematics and Computer Science Delft Institute of Applied Mathematics Probabilistic Inversion Technique for Seismic Data A thesis submitted

More information

Bayesian Lithology-Fluid Prediction and Simulation based. on a Markov Chain Prior Model

Bayesian Lithology-Fluid Prediction and Simulation based. on a Markov Chain Prior Model Bayesian Lithology-Fluid Prediction and Simulation based on a Markov Chain Prior Model Anne Louise Larsen Formerly Norwegian University of Science and Technology, N-7491 Trondheim, Norway; presently Schlumberger

More information

Oil and Natural Gas Corporation Ltd., VRC(Panvel), WOB, ONGC, Mumbai. 1

Oil and Natural Gas Corporation Ltd., VRC(Panvel), WOB, ONGC, Mumbai. 1 P-259 Summary Data for identification of Porosity Behaviour in Oligocene Lime Stone of D18 Area Of Western Offshore, India V.K. Baid*, P.H. Rao, P.S. Basak, Ravi Kant, V. Vairavan 1, K.M. Sundaram 1, ONGC

More information

Reservoir properties inversion from AVO attributes

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

More information

Hydrogeophysics - Seismics

Hydrogeophysics - Seismics Hydrogeophysics - Seismics Matthias Zillmer EOST-ULP p. 1 Table of contents SH polarized shear waves: Seismic source Case study: porosity of an aquifer Seismic velocities for porous media: The Frenkel-Biot-Gassmann

More information