B. Todd Hoffman and Jef Caers Stanford University, California, USA
|
|
- Derick Warren
- 5 years ago
- Views:
Transcription
1 Sequential Simulation under local non-linear constraints: Application to history matching B. Todd Hoffman and Jef Caers Stanford University, California, USA Introduction Sequential simulation has emerged as a robust and fast method for generating stochastic realizations. The recent development of a sequential simulation method by drawing structures from training images allows generating a wide variety of geological styles in a controlled fashion. Various methods for conditioning realizations to secondary information from geophysical or remote sensing techniques have been developed and applied. However, the current state-of-the-art in sequential simulation lacs the ability to include data that has a strongly non-linear and multiple-point relationship with the unnown. A classical example of such data is a subsurface flow response measured in a well, over time, and may consist of pressure, fractional flow or tracer data. Incorporating such data is of crucial importance in building models with fields such as petroleum engineering and hydrology. Most co-riging or bloc-riging methods cannot incorporate such data since they require a linear or pseudo-linear, often single-point relationship, between the target variable and the data. In this paper we present a method for constraining sequential simulation method to such data. The method is an iterative approach that searches efficiently for realizations that match any type of non-linear data. The approach consists of iteratively perturbing the local conditional probability distribution used to simulate outcomes at each grid node. Theory shows that the proposed perturbation method does not destroy the intended reproduction of geological continuity as provided by a variogram model, or by a training image, nor any conditioning by primary or secondary data. The perturbations are regional, in the sense that various regions of the model can be perturbed by different amounts. Such regional approach provides great flexibility in matching a large variety of local non-linear constraints. The number of regions and their shape can be chosen arbitrarily. The magnitude of model perturbation in each region is optimized using a newly developed parallel 1D optimization method. We show the efficiency and effectiveness of this method using a realistic 3D synthetic example. As application example, we discuss the problem of history matching pressure and flow data in oil and gas fields. Sequential simulation: a recall The aim of sequential simulation, as it was originally construed, it to reproduce the properties of a given multi-variate distribution under some conditioning data. Consider a set of N random variables Z(u i ) for any set of N locations u i, not necessarily on a regular grid. Given a multi-variate law defines on the set of N variables, a decomposition is formed F( u, u, K, u,z,z, K,z ) = F( u,z (n N+ 1)) F( u,z (n N+ 2)) 1 2 N 1 2 N N N N 1 N 1 F( u,z (n + 1)) F( u,z (n)) F(u j,z j (n-j+1) is the conditional cdf (ccdf) of Z(u j ) given a set of data values (n) and the previous (j-1) realizations z ( ) (u i ). This decomposition allows to generate a sample from the multivariate law by sequentially visiting, then sampling at each node along a random path. While the framewor provided by the decomposition loos theoretically appealing, it is limited in its practical extent. The only law for which a practical decomposition is available is the multi-variate Gaussian distribution. Each ccdf is in that case estimated using a linear riging technique, yet the resulting realizations lac geological realism as they reflect the maximum entropy property (disconnected extremes) of the multi-gaussian law, constrained only by the variogram. A more ad-hoc approach to the problem is taen, by considering the fact that any product of a series of conditional distribution results in a specific multi-variate law, represented by its realizations. The multivariate law is never explicitly stated, motivated by the fact that in all practical applications the primary interest lies in the realization and the properties/statistics it exhibits. Drawing at each location along the random path from a series of conditional distributions will result in a given realization. The aim of doing so is to reproduce realizations that reproduce a certain statistic or statistics. Some important examples are: 1
2 Sequential indicator simulation (SISIM, Deutsch and Journel, 1998) relies on indicator riging to determine ccdfs and aims at reproducing the variogram of a series of ranges of values, rather than the overall variogram. Single normal equation simulation (SNESIM, Strebelle, 2002) relies on obtaining the ccdfs from a training image by scanning the image of configurations of the data (n-j+1). In this way, the geological heterogeneity of the training image is reproduced. Any of the above sequential simulation method can be conditioned to so-called hard data or exact sample observations or soft data or indirect observations. The latter require to perform some form of coriging during sequential simulation. The method also allows conditioning to volume-average data as long as the averaging is linear or a power average. Next, we develop a method for conditioning any type of sequential simulation method to non-linear data using a method termed probability perturbation (Caers, 2003) in a regionalized manner. Regional Probability Perturbation Conditioning geostatistical realizations to highly non-linear data requires an iterative technique. An initial realization is created and subsequently changed (perturbed) until it honors the new non-linear data. However, if the realization properties are perturbed directly, other data such as geologic continuity may not be honored. As discussed earlier, a convenient property of simulation is that the realizations are generated from probability models, which are now shortly denoted as P(A B), with A the unnown value to be simulated, and B any hard, soft or linear volume average data and previously simulated values. As D we denote any set of non-linear data. For example A could be the hydraulic conductivity at any location, B some well data and ground penetrating radar data and D some tracer test data. Given a random seed and a series of probability models, P(A B), sequential simulation will generate a unique realization. Consider using one such realization as the initial realization. One way to perturb this realization is the probability model, P(A B) and eep the same random seed. The magnitude of the perturbation of P(A B) will necessarily depend on the difference between the data D and the forward evaluation of the physical law the reproduces that data. For demonstration purposes, we will consider only the case of a binary spatial variable described by an indicator random function model 1 if a given event occurs at u I( u ) = 0 else A={i(u) = 1} could mean that sand occurs at location u, while i(u) = 0 indicates non-sand occurrence. The initial realization containing all locations u will be termed i (0) (u). In the case where the initial realization already honors the non-linear data, there should be no perturbation. Conversely, in the case where the realization does not honor the non-linear data, the realization should be perturbed considerably. To achieve a perturbation that is between no change and considerable change, another probability model, P(A D), that depends on the non-linear data, D, is introduced. A perturbation of P(A B is achieved by combining the conditional probabilities P(A B) and P(A D) into a joint probability model, P(A B,D). The combination of two probability models into a new one is achieved with Journel s conditional independence method (2002). This new probability model is used to populate the next realization. P(A D) is defined as follows: P(A D) = (1-r D ) i (0) (u) + r D P(A) [0,1] where P(A) is the marginal distribution. A set of free parameters, {r D1,,r DK } control how much the model is perturbed in a particular region,. Note that K is the total number of regions and indicates a particular region. Region definition is left to the user, and while they may have any arbitrary shape, they may not overlap. The regions of the realization will be denoted {R 1, R 2,, R K }, and the entire realization is R = (R 1 R 2 R K ). P(A D) is defined for the entire realization, R, but its local value depends on the region definition: when u is located in region R, the perturbation parameter taes on a value of r D. Therefore P(A D) can have different values for different regions of the reservoir. To better understand the relationship between r D and P(A D) and how a perturbation of some initial realization is created consider the two limiting cases all r D =1 and (2) all r D =0. When they all equal zero, P(A D) = i (0) (u) and via Journel s method, P(A B,D) = i (0) (u); hence, the initial realization, i (0) (u), is 2
3 retained in its entirety. When they all are equal to one, P(A D) = P(A); therefore, P(A B,D) = P(A B), and a new realization, i (u), is generated that is equally probable to i (0) (u). Therefore the set of parameters, {r D1,,r DK }, define a perturbation of an initial realization to another equiprobable realization. When parameters are not zero or one but values in between (as they usually are), the resulting realizations will be some combination of the two realizations. Figure 1A shows an initial realization of two facies (or categories), and Figure 1B is the map of P(A D). It is generated using Equation with r D1 =0.80 and r D2 =5. By combining this P(A D) with the same P(A B) used to generate the initial realization, a new realizations, i rd ( u), can be generated (Fig. 1C). Since rd1 is close to one, this region is significantly different from the initial. Conversely, since r D2 is close to zero, this region is very similar to the initial realization. Notice that the facies are in the same location and only the shape has slightly changed. Although the two regions perturbations are quite dissimilar, discontinuities along the border of Regions 1 and 2 are not observed. The reason for not creating artifact discontinuities can be explained by the nature of the sequential simulation algorithm and by the application of the perturbation method. In sequential simulation, each location is simulated based on conditioning data and any previously simulated nodes. The method searches for any such previously simulated nodes in an elliptical search neighborhood. This search neighborhood may (and should) cross the region-boundaries. When simulating a nodal location in one region, the nodal values in any other regions are used to determine P(A B,D), hence creating continuity across the boundaries. Secondly, geological continuity is assured in the perturbation method through the probability P(A B), which is not calculated per region but for all regions together. A similar approach was taen in Hu et al. (2001) in conjunction with the gradual deformation methods that have some degrees of similarity with the probability perturbation method. Yet it should be pointed out that the probability perturbation method does not aim at changes the model gradually as this often leads to slow convergence and the method can be applied to any type of variable. Optimize r D. There may exist a set of r D values, such that i rd ( u), will honor the non-linear data better than the initial realization. Finding the optimum realization, i rdopt ( u), is a problem parameterized by the free parameters, {rd1,,r DK }. Finding the optimum realization is equivalent to finding the optimum set of r D values. {rd1,..., rdk} opt = min{o(rd1,..., rdk ) = Σ KO(rD ) } (2) r D where O(r D1,, r DK ) is the overall objective function, and O(r D ) is the regional objective. O(r D ) is calculated for each region and is defined as some measure of difference between the actual data, D and S the data currently included in the realization, D (r ). D O S = D (r ) D :1, 2Κ K (3) D In general the values of O might depend on all {r D1 r DK }. However, the values in a region are mostly dependent on the properties in that region; thus, we will assume, at least explicitly, that O only depends on r D. Based on this assumption, the best r D for each region can be determined using a one-dimensional optimization routine, but since there are multiple regions, the one-dimensional problems must be solved in parallel (at the same time but without explicit reference to each other). Note that the model is not actually broen into K independent problems where each region would be solved separately. Rather, the realizations are always generated over the entire reservoir, and regions are only used for objective function calculations and perturbation parameter, r D, updating. Although the r D values are updated without explicit reference to each other, dependence is implicitly maintained by creating the models over the entire model. Using the parallel 1D search for a best realization between two equiprobable realizations, one does not expect to converge to a global minimum. Thus, a two-loop optimization routine is necessary where the previous optimum realization, i rdopt ( u), is used as the initial realization in the next step, replacing i (0) (u) in Eq.. This constitutes the outer loop of iteration that will be stopped when the error is below a sufficient tolerance. Also during each outer iteration, the random seed is changed. By doing so, a new equiprobable realization is generated when all r D =1. This allows the method to again search between two equiprobable realizations in the next inner iteration. 3
4 Example In order to obtain accurate predictions of flow in oil and gas reservoirs, a model of the subsurface is created. This model needs to be constrained to a wide variety of data such as hard data from well logs, soft data from seismic, geological conceptual models and many types of dynamic (pressure/flow) data observed from past production or well performance tests. This dynamic data from subsurface porous media is often highly non-linear. An example of how to include and honor all of the various types of data is illustrated. The example uses a synthetic 3D fluvial channel reservoir model as the reference. The reference model was created by Mao (1999) and is based on a North Sea reservoir. It is made up of three distinct facies: channel sand, crevasse splays, and mudstone. The locations of the facies were determined by using an object-based simulation technique. The petrophysical properties were generated independently for each facies using Sequential Gaussian Simulation. The model has aerial dimensions of 10,000 ft by 13,000 ft and is 450 ft thic, and it is divided into 100 gridblocs in the x-direction, 130 in the y-direction, and 30 in the z-direction. The z-direction has three major horizons each with 10 layers. The three horizons have channels of different size, different directions, and different sinuosity (Figure 2). In this example nine vertical wells (six production and three injection) were added to the model for flow simulation. A real-life situation of incomplete information is established in order to generate a reservoir simulation model. The SNESIM algorithm is used to create the geostatistical realizations for the model (Strebelle, 2002). The geological parameters required for the realizations differ from the reference truth as follows: The facies proportions are calculated from the well-logs directly. The proportion of channel sand encountered by the wells is 27% and the proportion of crevasse encountered is 6%. Note that this is different than the reference, which has channel and crevasse proportions equal to 41% and 5%, respectively. The difference in the amount of facies also effects the pore volume. While the reference pore volume is 1.76x10 9 barrels, the generated models are from 1.32 x10 9 to 1.44 x10 9 barrels. Although the reference has different channel shape parameters for different layers, the channel shape parameters for the geostatistical realizations are the same for all layers. Two layers of a five-layer training image are shown in Figure 3. The training image was generated with an object-based technique, and each layer is 200 by 200 gridblocs. The proportions for each facies in the training image come directly from the well-log information discussed in the previous bullet. The porosity and permeability are assumed constant for each facies and obtained by averaging porosity and permeability per facies from wells. The channel porosity is 29%, and its permeability equals 456 md. For the crevasse splay, the porosity is 23%, and the permeability is 253 md, and for the mudstone, the porosity equals 7%, and the permeability is 3 md. The matter of how to define regions within the reservoir remains to be discussed. For this example, streamlines are well suited for the job because they show the direct paths by which fluid will enter a production well. Every production well will have a set of streamlines entering that well. All cells hit by this set of streamlines define the drainage zone for that well, and the various drainage zones will define the region geometry. As the model is perturbed, the facies geometry will change; consequently, the drainage area of the production wells also will change. The region definitions for one stage are displayed in Figure 4 Fractional flow information at the six production wells is the non-linear data that is being incorporated into the model. Although the generated models have different channel shapes than the reference, a match with the reference data is consistently achieved. Figure 5 displays a realization that has included the fractional flow data. The reference data along with the initial and matched data are displayed for the six wells (Figure 6). Although, exact matches are not achieved, significant improvement is observed in all wells, and some wells such as Well 3 and Well 6 show dramatic improvement. The matches obtained would be adequate for most true field applications. Conclusions 1. The proposed regional probability perturbation method can efficiently include highly non-linear data while continuing to honor other pieces of data such as geologic continuity. 2. The method is quite general in the sense that any type of sequential simulation method can be used as well as any type of non-linear data as long as the forward model (e.g. a finite element or difference) is nown. 4
5 3. The method has been shown to wor for the practical example of including production data in a petroleum reservoir model. A good history match was attained on a channel reservoir that has 9 wells, 10 years of production and 50,000 gridblocs. References Caers, J.: Geostatistical history matching under training-image based geological model constraints, paper SPE presented 77 th Annual Fall Meeting of SPE in San Antonio, TX, Sept. 29 Oct. 2, see also SPE 74716, to be printed in SPE Journal, Deutsch, C.V and Journel, A.G. (1998), GSLIB: the geostatistical software library. Oxford University Press. Hu, L-Y., Blanc, G., Noetinger, B., Gradual deformation and iterative calibration of sequential simulations. Mathematical Geology, 33, Journel, A. G., Combining Knowledge from Diverse Sources: an Alternative to Traditional Data Independence Hypothesis, Mathematical Geology, v. 34, No. 1, January Mao, S., Generation of a reference petrophysical/seismic data set: the Stanford V reservoir, Stanford Center for Reservoir Forcasting (SCRF) Report No. 12, May Strebelle, S.: Conditioning Simulation of Complex Structures Multiple-Point Statistics Mathematical Geology, v. 34, No. 1, January Initial Realization (A) P(A D) (B) Perturbed Realization (C) 0.5 P(A) r D1 = 0.80 r D2 = 5 Figure 1: Perturbation of initial realization, (A) to get new realization, (C). Perturbation completed by combining initial probability model, P(A B), with new probability model P(A D), (B). Porosity Layer 6 Layer Layer Figure 2: Reference model for 3D example. Porosity is shown and channels can be distinguished. 5
6 I1 I1 200 Gridblocs 200 Gridblocs Figure 3: Two layers from a five-layer training image. Region Geometry Region 1 Layer 1 Layer 15 Region 2 Region 3 Region 4 Region 5 I1 P2 P1 P3 I2 P4 P6 I3 P5 I1 P2 P1 P3 I2 P4 P6 I3 P5 Region 6 Figure 4: Dynamic region geometry for one stage defined using streamlines. Layer 3 Layer 8 65 Layer Figure 5: History matched realization for 3D example. Water Cut Time (days) 1 ref 1 init 1 match 3 ref 3 init 3 match Water Cut Time (days) 2 ref 2 init 2 match 5 ref 5 init 5 match Water Cut Time (days) 4 ref 4 init 4 match 6 ref 6 init 6 match Figure 6: Water cut reference, initial match and final match of six production wells for one realization. 6
A009 HISTORY MATCHING WITH THE PROBABILITY PERTURBATION METHOD APPLICATION TO A NORTH SEA RESERVOIR
1 A009 HISTORY MATCHING WITH THE PROBABILITY PERTURBATION METHOD APPLICATION TO A NORTH SEA RESERVOIR B. Todd HOFFMAN and Jef CAERS Stanford University, Petroleum Engineering, Stanford CA 94305-2220 USA
More informationHistory matching under training-image based geological model constraints
History matching under training-image based geological model constraints JEF CAERS Stanford University, Department of Petroleum Engineering Stanford, CA 94305-2220 January 2, 2002 Corresponding author
More informationPrograms for MDE Modeling and Conditional Distribution Calculation
Programs for MDE Modeling and Conditional Distribution Calculation Sahyun Hong and Clayton V. Deutsch Improved numerical reservoir models are constructed when all available diverse data sources are accounted
More informationA Parallel, Multiscale Approach to Reservoir Modeling. Omer Inanc Tureyen and Jef Caers Department of Petroleum Engineering Stanford University
A Parallel, Multiscale Approach to Reservoir Modeling Omer Inanc Tureyen and Jef Caers Department of Petroleum Engineering Stanford University 1 Abstract With the advance of CPU power, numerical reservoir
More informationA Geostatistical and Flow Simulation Study on a Real Training Image
A Geostatistical and Flow Simulation Study on a Real Training Image Weishan Ren (wren@ualberta.ca) Department of Civil & Environmental Engineering, University of Alberta Abstract A 12 cm by 18 cm slab
More informationMultiple Point Statistics with Multiple Training Images
Multiple Point Statistics with Multiple Training Images Daniel A. Silva and Clayton V. Deutsch Abstract Characterization of complex geological features and patterns has been one of the main tasks of geostatistics.
More informationIntegration of Geostatistical Modeling with History Matching: Global and Regional Perturbation
Integration of Geostatistical Modeling with History Matching: Global and Regional Perturbation Oliveira, Gonçalo Soares Soares, Amílcar Oliveira (CERENA/IST) Schiozer, Denis José (UNISIM/UNICAMP) Introduction
More informationSampling informative/complex a priori probability distributions using Gibbs sampling assisted by sequential simulation
Sampling informative/complex a priori probability distributions using Gibbs sampling assisted by sequential simulation Thomas Mejer Hansen, Klaus Mosegaard, and Knud Skou Cordua 1 1 Center for Energy Resources
More informationExploring Direct Sampling and Iterative Spatial Resampling in History Matching
Exploring Direct Sampling and Iterative Spatial Resampling in History Matching Matz Haugen, Grergoire Mariethoz and Tapan Mukerji Department of Energy Resources Engineering Stanford University Abstract
More informationShort Note: Some Implementation Aspects of Multiple-Point Simulation
Short Note: Some Implementation Aspects of Multiple-Point Simulation Steven Lyster 1, Clayton V. Deutsch 1, and Julián M. Ortiz 2 1 Department of Civil & Environmental Engineering University of Alberta
More informationRotation and affinity invariance in multiple-point geostatistics
Rotation and ainity invariance in multiple-point geostatistics Tuanfeng Zhang December, 2001 Abstract Multiple-point stochastic simulation of facies spatial distribution requires a training image depicting
More informationCONDITIONING FACIES SIMULATIONS WITH CONNECTIVITY DATA
CONDITIONING FACIES SIMULATIONS WITH CONNECTIVITY DATA PHILIPPE RENARD (1) and JEF CAERS (2) (1) Centre for Hydrogeology, University of Neuchâtel, Switzerland (2) Stanford Center for Reservoir Forecasting,
More informationA workflow to account for uncertainty in well-log data in 3D geostatistical reservoir modeling
A workflow to account for uncertainty in well-log data in 3D geostatistical reservoir Jose Akamine and Jef Caers May, 2007 Stanford Center for Reservoir Forecasting Abstract Traditionally well log data
More informationGeostatistical Reservoir Characterization of McMurray Formation by 2-D Modeling
Geostatistical Reservoir Characterization of McMurray Formation by 2-D Modeling Weishan Ren, Oy Leuangthong and Clayton V. Deutsch Department of Civil & Environmental Engineering, University of Alberta
More informationHigh Resolution Geomodeling, Ranking and Flow Simulation at SAGD Pad Scale
High Resolution Geomodeling, Ranking and Flow Simulation at SAGD Pad Scale Chad T. Neufeld, Clayton V. Deutsch, C. Palmgren and T. B. Boyle Increasing computer power and improved reservoir simulation software
More informationHIERARCHICAL SIMULATION OF MULTIPLE-FACIES RESERVOIRS USING MULTIPLE-POINT GEOSTATISTICS
HIERARCHICAL SIMULATION OF MULTIPLE-FACIES RESERVOIRS USING MULTIPLE-POINT GEOSTATISTICS AREPORT SUBMITTED TO THE DEPARTMENT OF PETROLEUM ENGINEERING OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT OF THE
More informationUsing 3D-DEGA. Omer Inanc Tureyen and Jef Caers Department of Petroleum Engineering Stanford University
Using 3D-DEGA Omer Inanc Tureyen and Jef Caers Department of Petroleum Engineering Stanford University 1 1 Introduction With the advance of CPU power, numerical reservoir models have become an essential
More informationMultiple-point geostatistics: a quantitative vehicle for integrating geologic analogs into multiple reservoir models
Multiple-point geostatistics: a quantitative vehicle for integrating geologic analogs into multiple reservoir models JEF CAERS AND TUANFENG ZHANG Stanford University, Stanford Center for Reservoir Forecasting
More informationRM03 Integrating Petro-elastic Seismic Inversion and Static Model Building
RM03 Integrating Petro-elastic Seismic Inversion and Static Model Building P. Gelderblom* (Shell Global Solutions International BV) SUMMARY This presentation focuses on various aspects of how the results
More informationSimulating Geological Structures Based on Training Images and Pattern Classifications
Simulating Geological Structures Based on Training Images and Pattern Classifications P. Switzer, T. Zhang, A. Journel Department of Geological and Environmental Sciences Stanford University CA, 9435,
More informationA Geomodeling workflow used to model a complex carbonate reservoir with limited well control : modeling facies zones like fluid zones.
A Geomodeling workflow used to model a complex carbonate reservoir with limited well control : modeling facies zones like fluid zones. Thomas Jerome (RPS), Ke Lovan (WesternZagros) and Suzanne Gentile
More informationTensor Based Approaches for LVA Field Inference
Tensor Based Approaches for LVA Field Inference Maksuda Lillah and Jeff Boisvert The importance of locally varying anisotropy (LVA) in model construction can be significant; however, it is often ignored
More informationOn internal consistency, conditioning and models of uncertainty
On internal consistency, conditioning and models of uncertainty Jef Caers, Stanford University Abstract Recent research has been tending towards building models of uncertainty of the Earth, not just building
More informationA 3D code for mp simulation of continuous and
A 3D code for mp simulation of continuous and categorical variables: FILTERSIM Jianbing Wu, Alexandre Boucher & André G. Journel May, 2006 Abstract In most petroleum and geoscience studies, the flow is
More informationConditioning a hybrid geostatistical model to wells and seismic data
Conditioning a hybrid geostatistical model to wells and seismic data Antoine Bertoncello, Gregoire Mariethoz, Tao Sun and Jef Caers ABSTRACT Hybrid geostatistical models imitate a sequence of depositional
More informationSPE Copyright 2002, Society of Petroleum Engineers Inc.
SPE 77958 Reservoir Modelling With Neural Networks And Geostatistics: A Case Study From The Lower Tertiary Of The Shengli Oilfield, East China L. Wang, S. Tyson, Geo Visual Systems Australia Pty Ltd, X.
More informationReservoir Modeling Combining Geostatistics with Markov Chain Monte Carlo Inversion
Reservoir Modeling Combining Geostatistics with Markov Chain Monte Carlo Inversion Andrea Zunino, Katrine Lange, Yulia Melnikova, Thomas Mejer Hansen and Klaus Mosegaard 1 Introduction Reservoir modeling
More informationModeling Uncertainty in the Earth Sciences Jef Caers Stanford University
Modeling spatial continuity Modeling Uncertainty in the Earth Sciences Jef Caers Stanford University Motivation uncertain uncertain certain or uncertain uncertain Spatial Input parameters Spatial Stochastic
More informationHierarchical modeling of multi-scale flow barriers in channelized reservoirs
Hierarchical modeling of multi-scale flow barriers in channelized reservoirs Hongmei Li and Jef Caers Stanford Center for Reservoir Forecasting Stanford University Abstract Channelized reservoirs often
More informationAdaptive spatial resampling as a Markov chain Monte Carlo method for uncertainty quantification in seismic reservoir characterization
1 Adaptive spatial resampling as a Markov chain Monte Carlo method for uncertainty quantification in seismic reservoir characterization Cheolkyun Jeong, Tapan Mukerji, and Gregoire Mariethoz Department
More informationA PARALLEL MODELLING APPROACH TO RESERVOIR CHARACTERIZATION
A PARALLEL MODELLING APPROACH TO RESERVOIR CHARACTERIZATION A DISSERTATION SUBMITTED TO THE DEPARTMENT OF PETROLEUM ENGINEERING AND THE COMMITTEE ON GRADUATE STUDIES OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT
More informationFluid flow modelling with seismic cluster analysis
Fluid flow modelling with seismic cluster analysis Fluid flow modelling with seismic cluster analysis Laurence R. Bentley, Xuri Huang 1 and Claude Laflamme 2 ABSTRACT Cluster analysis is used to construct
More informationAppropriate algorithm method for Petrophysical properties to construct 3D modeling for Mishrif formation in Amara oil field Jawad K.
Appropriate algorithm method for Petrophysical properties to construct 3D modeling for Mishrif formation in Amara oil field Jawad K. Radhy AlBahadily Department of geology, college of science, Baghdad
More informationFast FILTERSIM Simulation with Score-based Distance Function
Fast FILTERSIM Simulation with Score-based Distance Function Jianbing Wu (1), André G. Journel (1) and Tuanfeng Zhang (2) (1) Department of Energy Resources Engineering, Stanford, CA (2) Schlumberger Doll
More informationMarkov Bayes Simulation for Structural Uncertainty Estimation
P - 200 Markov Bayes Simulation for Structural Uncertainty Estimation Samik Sil*, Sanjay Srinivasan and Mrinal K Sen. University of Texas at Austin, samiksil@gmail.com Summary Reservoir models are built
More informationHierarchical Trend Models Based on Architectural Elements
Hierarchical Trend Models Based on Architectural Elements Michael J. Pyrcz (mpyrcz@ualberta.ca) and Clayton V. Deutsch (cdeutsch@ualberta.ca) Department of Civil & Environmental Engineering University
More informationOn Secondary Data Integration
On Secondary Data Integration Sahyun Hong and Clayton V. Deutsch A longstanding problem in geostatistics is the integration of multiple secondary data in the construction of high resolution models. In
More informationWe G Updating the Reservoir Model Using Engineeringconsistent
We G102 09 Updating the Reservoir Model Using Engineeringconsistent 4D Seismic Inversion S. Tian* (Heriot-Watt University), C. MacBeth (Heriot-Watt University) & A. Shams (Heriot-Watt University) SUMMARY
More informationGeostatistics on Stratigraphic Grid
Geostatistics on Stratigraphic Grid Antoine Bertoncello 1, Jef Caers 1, Pierre Biver 2 and Guillaume Caumon 3. 1 ERE department / Stanford University, Stanford CA USA; 2 Total CSTJF, Pau France; 3 CRPG-CNRS
More informationHistory Matching of Structurally Complex Reservoirs Using a Distance-based Model Parameterization
History Matching of Structurally Complex Reservoirs Using a Distance-based Model Parameterization Satomi Suzuki, Guillaume Caumon, Jef Caers S. Suzuki, J. Caers Department of Energy Resources Engineering,
More informationSIMPAT: Stochastic Simulation with Patterns
SIMPAT: Stochastic Simulation with Patterns G. Burc Arpat Stanford Center for Reservoir Forecasting Stanford University, Stanford, CA 94305-2220 April 26, 2004 Abstract Flow in a reservoir is mostly controlled
More informationIterative spatial resampling applied to seismic inverse modeling for lithofacies prediction
Iterative spatial resampling applied to seismic inverse modeling for lithofacies prediction Cheolkyun Jeong, Tapan Mukerji, and Gregoire Mariethoz Department of Energy Resources Engineering Stanford University
More informationJoint quantification of uncertainty on spatial and non-spatial reservoir parameters
Joint quantification of uncertainty on spatial and non-spatial reservoir parameters Comparison between the Method and Distance Kernel Method Céline Scheidt and Jef Caers Stanford Center for Reservoir Forecasting,
More informationVariogram Inversion and Uncertainty Using Dynamic Data. Simultaneouos Inversion with Variogram Updating
Variogram Inversion and Uncertainty Using Dynamic Data Z. A. Reza (zreza@ualberta.ca) and C. V. Deutsch (cdeutsch@civil.ualberta.ca) Department of Civil & Environmental Engineering, University of Alberta
More informationTrend Modeling Techniques and Guidelines
Trend Modeling Techniques and Guidelines Jason A. M c Lennan, Oy Leuangthong, and Clayton V. Deutsch Centre for Computational Geostatistics (CCG) Department of Civil and Environmental Engineering University
More informationDynamic data integration and stochastic inversion of a two-dimensional confined aquifer
Hydrology Days 2013 Dynamic data integration and stochastic inversion of a two-dimensional confined aquifer Dongdong Wang 1 Ye Zhang 1 Juraj Irsa 1 Abstract. Much work has been done in developing and applying
More informationDownscaling saturations for modeling 4D seismic data
Downscaling saturations for modeling 4D seismic data Scarlet A. Castro and Jef Caers Stanford Center for Reservoir Forecasting May 2005 Abstract 4D seismic data is used to monitor the movement of fluids
More informationEnsemble Kalman Filter Predictor Bias Correction Method for Non-Gaussian Geological Facies Detection
Proceedings of the 01 IFAC Worshop on Automatic Control in Offshore Oil and Gas Production, Norwegian University of Science and Technology, Trondheim, Norway, May 31 - June 1, 01 ThC. Ensemble Kalman Filter
More informationAlgorithm-driven and Representation-driven Random Function : A New Formalism for Applied Geostatistics
Algorithm-driven and Representation-driven Random Function : A New Formalism for Applied Geostatistics Alexandre Boucher Dept of Geological and Environmental Sciences Stanford University Abstract This
More informationMPS Simulation with a Gibbs Sampler Algorithm
MPS Simulation with a Gibbs Sampler Algorithm Steve Lyster and Clayton V. Deutsch Complex geologic structure cannot be captured and reproduced by variogram-based geostatistical methods. Multiple-point
More informationSurface-based model conditioning using an hybrid optimization: methodology and application
Surface-based model conditioning using an hybrid optimization: methodology and application Antoine Bertoncello, Jef Caers, Hongmei Li and Tao Sun Department of Energy Resources Engineering Stanford University
More informationE10 Vf~TC~DS FOR HISTORY MATCHING UNDE R ~E ~~AL~JNSTRAINT
E10 Vf~TC~DS FOR HISTORY MATCHING UNDE R ~E ~~AL~JNSTRAINT Jef Caers Stanford University, Petroleum Engineering, Stanford CA 94305-2220, USA Abstrac t Two geostatistical methods for history matching are
More informationUncertainty Quantification Using Distances and Kernel Methods Application to a Deepwater Turbidite Reservoir
Uncertainty Quantification Using Distances and Kernel Methods Application to a Deepwater Turbidite Reservoir Céline Scheidt and Jef Caers Stanford Center for Reservoir Forecasting, Stanford University
More informationB002 DeliveryMassager - Propagating Seismic Inversion Information into Reservoir Flow Models
B2 DeliveryMassager - Propagating Seismic Inversion Information into Reservoir Flow Models J. Gunning* (CSIRO Petroleum) & M.E. Glinsky (BHP Billiton Petroleum) SUMMARY We present a new open-source program
More informationCrosswell Tomographic Inversion with Block Kriging
Crosswell Tomographic Inversion with Block Kriging Yongshe Liu Stanford Center for Reservoir Forecasting Petroleum Engineering Department Stanford University May, Abstract Crosswell tomographic data can
More informationApplication of MPS Simulation with Multiple Training Image (MultiTI-MPS) to the Red Dog Deposit
Application of MPS Simulation with Multiple Training Image (MultiTI-MPS) to the Red Dog Deposit Daniel A. Silva and Clayton V. Deutsch A Multiple Point Statistics simulation based on the mixing of two
More informationCalibration of NFR models with interpreted well-test k.h data. Michel Garcia
Calibration of NFR models with interpreted well-test k.h data Michel Garcia Calibration with interpreted well-test k.h data Intermediate step between Reservoir characterization Static model conditioned
More informationEFFICIENT PRODUCTION OPTIMIZATION USING FLOW NETWORK MODELS. A Thesis PONGSATHORN LERLERTPAKDEE
EFFICIENT PRODUCTION OPTIMIZATION USING FLOW NETWORK MODELS A Thesis by PONGSATHORN LERLERTPAKDEE Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements
More informationA Data estimation Based Approach for Quasi continuous Reservoir Monitoring using Sparse Surface Seismic Data Introduction Figure 1
A Data estimation Based Approach for Quasi continuous Reservoir Monitoring using Sparse Surface Seismic Data Adeyemi Arogunmati* and Jerry M. Harris, Stanford University, California, USA Introduction One
More informationMulti-Objective Stochastic Optimization by Co-Direct Sequential Simulation for History Matching of Oil Reservoirs
Multi-Objective Stochastic Optimization by Co-Direct Sequential Simulation for History Matching of Oil Reservoirs João Daniel Trigo Pereira Carneiro under the supervision of Amílcar de Oliveira Soares
More informationSTOCHASTIC OBJECT-BASED SIMULATION OF CHANNELS CONSTRAINED ABSTRACT INTRODUCTION 1 DATA AND GEOLOGICAL FRAMEWORK BY HIGH RESOLUTION SEISMIC DATA
STOCHASTIC OBJECT-BASED SIMULATION OF CHANNELS CONSTRAINED BY HIGH RESOLUTION SEISMIC DATA S. Viseur, A. Shtuka, J-L. Mallet ABSTRACT Recent progresses in exploration techniques (sonar images, high resolution
More informationDeveloping a Smart Proxy for the SACROC Water-Flooding Numerical Reservoir Simulation Model
SPE-185691-MS Developing a Smart Proxy for the SACROC Water-Flooding Numerical Reservoir Simulation Model Faisal Alenezi and Shahab Mohaghegh, West Virginia University Copyright 2017, Society of Petroleum
More informationIndicator Simulation for Categorical Variables
Reservoir Modeling with GSLIB Indicator Simulation for Categorical Variables Sequential Simulation: the Concept Steps in Sequential Simulation SISIM Program Sequential Simulation: the Concept 2 1 3 1.
More informationSelected Implementation Issues with Sequential Gaussian Simulation
Selected Implementation Issues with Sequential Gaussian Simulation Abstract Stefan Zanon (szanon@ualberta.ca) and Oy Leuangthong (oy@ualberta.ca) Department of Civil & Environmental Engineering University
More informationProbabilistic Graphical Models
School of Computer Science Probabilistic Graphical Models Theory of Variational Inference: Inner and Outer Approximation Eric Xing Lecture 14, February 29, 2016 Reading: W & J Book Chapters Eric Xing @
More informationUniversity of Alberta. Multivariate Analysis of Diverse Data for Improved Geostatistical Reservoir Modeling
University of Alberta Multivariate Analysis of Diverse Data for Improved Geostatistical Reservoir Modeling by Sahyun Hong A thesis submitted to the Faculty of Graduate Studies and Research in partial fulfillment
More informationEnsemble-based decision making for reservoir management present and future outlook. TPD R&T ST MSU DYN and FMU team
Ensemble-based decision making for reservoir management present and future outlook TPD R&T ST MSU DYN and FMU team 11-05-2017 The core Ensemble based Closed Loop Reservoir Management (CLOREM) New paradigm
More informationMaximum A Posteriori Selection with Homotopic Constraint
Maximum A Posteriori Selection with Homotopic Constraint Michael J. Pyrcz and Clayton V. Deutsch Department of Civil & Environmental Engineering University of Alberta Abstract The addition of homotopic
More informationDirect Sequential Co-simulation with Joint Probability Distributions
Math Geosci (2010) 42: 269 292 DOI 10.1007/s11004-010-9265-x Direct Sequential Co-simulation with Joint Probability Distributions Ana Horta Amílcar Soares Received: 13 May 2009 / Accepted: 3 January 2010
More informationPTE 519 Lecture Note Finite Difference Approximation (Model)
PTE 519 Lecture Note 3 3.0 Finite Difference Approximation (Model) In this section of the lecture material, the focus is to define the terminology and to summarize the basic facts. The basic idea of any
More informationModule 1 Lecture Notes 2. Optimization Problem and Model Formulation
Optimization Methods: Introduction and Basic concepts 1 Module 1 Lecture Notes 2 Optimization Problem and Model Formulation Introduction In the previous lecture we studied the evolution of optimization
More information11-Geostatistical Methods for Seismic Inversion. Amílcar Soares CERENA-IST
11-Geostatistical Methods for Seismic Inversion Amílcar Soares CERENA-IST asoares@ist.utl.pt 01 - Introduction Seismic and Log Scale Seismic Data Recap: basic concepts Acoustic Impedance Velocity X Density
More informationThe SPE Foundation through member donations and a contribution from Offshore Europe
Primary funding is provided by The SPE Foundation through member donations and a contribution from Offshore Europe The Society is grateful to those companies that allow their professionals to serve as
More informationGrid-less Simulation of a Fluvio-Deltaic Environment
Grid-less Simulation of a Fluvio-Deltaic Environment R. Mohan Srivastava, FSS Canada Consultants, Toronto, Canada MoSrivastava@fssconsultants.ca and Marko Maucec, Halliburton Consulting and Project Management,
More information11. Petrophysical Modeling
11. Petrophysical Modeling 11.1 Deterministic Modeling When the well logs have been scaled up to the resolution of the cells in the 3D grid, the values for each cell along the well trajectory can be interpolated
More informationTPG4160 Reservoir simulation, Building a reservoir model
TPG4160 Reservoir simulation, Building a reservoir model Per Arne Slotte Week 8 2018 Overview Plan for the lectures The main goal for these lectures is to present the main concepts of reservoir models
More informationCIPC Louis Mattar. Fekete Associates Inc. Analytical Solutions in Well Testing
CIPC 2003 Louis Mattar Fekete Associates Inc Analytical Solutions in Well Testing Well Test Equation 2 P 2 P 1 P + = x 2 y 2 α t Solutions Analytical Semi-Analytical Numerical - Finite Difference Numerical
More informationA Soft Computing-Based Method for the Identification of Best Practices, with Application in the Petroleum Industry
CIMSA 2005 IEEE International Conference on Computational Intelligence for Measurement Systems and Applications Giardini Naxos, Italy, 20-22 July 2005 A Soft Computing-Based Method for the Identification
More informationIndicator Simulation Accounting for Multiple-Point Statistics
Indicator Simulation Accounting for Multiple-Point Statistics Julián M. Ortiz 1 and Clayton V. Deutsch 2 Geostatistical simulation aims at reproducing the variability of the real underlying phenomena.
More informationA family of particle swarm optimizers for reservoir characterization and seismic history matching.
P-487 Summary A family of particle swarm optimizers for reservoir characterization and seismic history matching. Tapan Mukerji*, Amit Suman (Stanford University), Juan Luis Fernández-Martínez (Stanford
More informationA Data-Driven Smart Proxy Model for A Comprehensive Reservoir Simulation
A Data-Driven Smart Proxy Model for A Comprehensive Reservoir Simulation Faisal Alenezi Department of Petroleum and Natural Gas Engineering West Virginia University Email: falenezi@mix.wvu.edu Shahab Mohaghegh
More informationGAS PRODUCTION ANALYSIS:
New Mexico Tech THINKING FOR A NEW MILLENNIUM TIGHT-GAS GAS PRODUCTION ANALYSIS: User Guide for a Type-Curve Matching Spreadsheet Program (TIGPA 2000.1) Her-Yuan Chen Assistant Professor Department of
More informationHistory Matching: Towards Geologically Reasonable Models
Downloaded from orbit.dtu.dk on: Oct 13, 018 History Matching: Towards Geologically Reasonable Models Melnikova, Yulia; Cordua, Knud Skou; Mosegaard, Klaus Publication date: 01 Document Version Publisher's
More informationFlexible Lag Definition for Experimental Variogram Calculation
Flexible Lag Definition for Experimental Variogram Calculation Yupeng Li and Miguel Cuba The inference of the experimental variogram in geostatistics commonly relies on the method of moments approach.
More informationReservoir Characterization with Limited Sample Data using Geostatistics
Reservoir Characterization with Limited Sample Data using Geostatistics By: Sayyed Mojtaba Ghoraishy Submitted to the Department of Chemical and Petroleum Engineering and the Faculty of the Graduate School
More informationFacies Modeling Using a Markov Mesh Model Specification
Math Geosci (2011) 43:611 624 DOI 10.1007/s11004-011-9350-9 Facies Modeling Using a Markov Mesh Model Specification Marita Stien Odd Kolbjørnsen Received: 28 April 2010 / Accepted: 13 February 2011 / Published
More informationarxiv: v1 [stat.ml] 15 Feb 2018
Conditioning of three-dimensional generative adversarial networks for pore and reservoir-scale models arxiv:1802.05622v1 [stat.ml] 15 Feb 2018 Lukas J. Mosser lukas.mosser15@imperial.ac.uk Martin J. Blunt
More information2D Geostatistical Modeling and Volume Estimation of an Important Part of Western Onland Oil Field, India.
and Volume Estimation of an Important Part of Western Onland Oil Field, India Summary Satyajit Mondal*, Liendon Ziete, and B.S.Bisht ( GEOPIC, ONGC) M.A.Z.Mallik (E&D, Directorate, ONGC) Email: mondal_satyajit@ongc.co.in
More informationModeling Uncertainty in the Earth Sciences Jef Caers Stanford University
Modeling response uncertainty Modeling Uncertainty in the Earth Sciences Jef Caers Stanford University Modeling Uncertainty in the Earth Sciences High dimensional Low dimensional uncertain uncertain certain
More informationSPE Distinguished Lecturer Program
SPE Distinguished Lecturer Program Primary funding is provided by The SPE Foundation through member donations and a contribution from Offshore Europe The Society is grateful to those companies that allow
More informationGroundwater in Hydrologic Cycle
Groundwater in Hydrologic Cycle Types of Terrestrial Water Surface Water Soil Moisture Ground water Pores Full of Combination of Air and Water Unsaturated Zone / Zone of Aeration / Vadose (Soil Water)
More informationOur Calibrated Model has No Predictive Value: An Example from the Petroleum Industry
Our Calibrated Model as No Predictive Value: An Example from te Petroleum Industry J.N. Carter a, P.J. Ballester a, Z. Tavassoli a and P.R. King a a Department of Eart Sciences and Engineering, Imperial
More informationSCRF. 22 nd Annual Meeting. April 30-May
SCRF 22 nd Annual Meeting April 30-May 1 2009 1 Research Overview CD annual report with papers Presentations 2 Modeling Uncertainty Distance based modeling of uncertainty - Model selection - Inverse modeling
More informationTechnique to Integrate Production and Static Data in a Self-Consistent Way Jorge L. Landa, SPE, Chevron Petroleum Technology Co.
SPE 71597 Technique to Integrate Production and Static Data in a Self-Consistent Way Jorge L. Landa, SPE, Chevron Petroleum Technology Co. Copyright 2001, Society of Petroleum Engineers Inc. This paper
More informationRubis (NUM) Tutorial #1
Rubis (NUM) Tutorial #1 1. Introduction This example is an introduction to the basic features of Rubis. The exercise is by no means intended to reproduce a realistic scenario. It is assumed that the user
More informationThe Components of Geostatistical Simulation
This file was created by scanning the printed publication. Errors identified by the software have been corrected; however, some errors may remain. The Components of Geostatistical Simulation Carol A. ~otwayl
More informationCONDITIONAL SIMULATION OF TRUNCATED RANDOM FIELDS USING GRADIENT METHODS
CONDITIONAL SIMULATION OF TRUNCATED RANDOM FIELDS USING GRADIENT METHODS Introduction Ning Liu and Dean S. Oliver University of Oklahoma, Norman, Oklahoma, USA; ning@ou.edu The problem of estimating the
More informationIMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL RESERVOIR SIMULATIONS USING THE CIP SCHEME WITH THIRD-ORDER ACCURACY
PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30 - February 1, 2012 SGP-TR-194 IMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL
More information3D MULTIDISCIPLINARY INTEGRATED GEOMECHANICAL FRACTURE SIMULATOR & COMPLETION OPTIMIZATION TOOL
PETROPHYSICS RESERVOIR GEOMECHANICS COMPLETIONS DRILLING PRODUCTION SERVICE 3D MULTIDISCIPLINARY INTEGRATED GEOMECHANICAL FRACTURE SIMULATOR & COMPLETION OPTIMIZATION TOOL INTEGRATED GEOMECHANICAL FRACTURE
More informationElectromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah
Electromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah Summary In this paper we present a new approach to the interpretation
More information