In this work, the application of tensor methodologies for computer-assisted history matching of channelized reservoirs is explored. A tensor-based approach is used for the parameterization of petrophysical parameters to reduce the dimensionality of the parameter estimation problem. Building on the work of Afra and Gildin (2013); Afra et.al. (2014); Afra and Gildin (2016), permeability fields of multiple model realizations are collected in a tensor form which is subsequently decomposed to derive a low-dimensional representation of the dominant spatial structures in the models. This representation then is used to estimate an identifiable reduced set of parameters using an ensemble Kalman filter (EnKF) strategy. This approach is attractive for the parameter estimation of permeabilities because it increases the ability to represent channelized structures in the updates resulting in an improved predictive capacity of the history-matched models. In particular, channel continuity is better preserved than with a Principal Component Analysis (PCA) parameterization. ... Read more

Conflicting objectives are frequently encountered in most real-world problems. When dealing with conflicting objectives, decision makers prefer to obtain a range of possible optimal solutions from which to choose. In theory, methods exists that can produce a range of possible solutions, some of which are “Pareto Optimal”. The application of these methods to solve bi-objective production optimization problems is increasing. A recent paper introduced a method to find points on the boundary of the objective function space by solving a constrained optimization problem using adjoint gradients. In this work, we investigate the applicability of using ensemble optimization (EnOpt) (which relies on approximate ensemble gradients instead of exact adjoint-based gradients) to generate points along a “Pareto” front with acceptable computational effort.. Moreover, we investigate the applicability of this approximate gradient technique to solve constrained optimization problems using the augmented Lagrangian method. Finally, we compare the performance of this bi-objective optimization method to a traditional weighted sum method for bi-objective water flooding optimization of two different synthetic reservoir models. The two objectives used in this work are, undiscounted (0%) net present value (NPV), representing long-term targets and highly discounted (25%) NPV, representing short-term operational targets. The controls are inflow control valve (ICV) settings over time for one model and water injection rate controls for the other. The effect of different starting points and the computational efficiency of the constrained optimization method are also investigated. ... Read more

We consider a technique to estimate an approximate gradient using an ensemble of randomly chosen control vectors, known as Ensemble Optimization (EnOpt) in the oil and gas reservoir simulation community. In particular, we address how to obtain accurate approximate gradients when the underlying numerical mod- els contain uncertain parameters because of geological uncertainties. In that case, ‘robust optimization’ is performed by optimizing the expected value of the objective function over an ensemble of geological mod- els. In earlier publications, based on the pioneering work of Chen et al. (2009), it has been suggested that a straightforward one-to-one combination of random control vectors and random geological models is capa- ble of generating sufficiently accurate approximate gradients. However, this form of EnOpt does not always yield satisfactory results. In a recent article, Fonseca et al. (2015) formulate a modified EnOpt algorithm, referred to here as a Stochastic Simplex Approximate Gradient (StoSAG; in earlier publications referred to as ‘modified robust EnOpt’) and show, via computational experiments, that StoSAG generally yields significantly better gradient approximations than the standard EnOpt algorithm. Here, we provide theoreti- cal arguments to show why StoSAG is superior to EnOpt ... Read more

We consider robust ensemble-based (EnOpt) multiobjective production optimization of on/off inflow-control devices (ICDs) for a sector model inspired by a real-field case. The use of on/off valves as optimization variables leads to a discrete control problem. We propose a reparameterization of such discrete controls in terms of switching times (i.e., we optimize the time at which a particular valve is either open or closed). This transforms the discrete control problem into a continuous control problem that can be efficiently handled with the EnOpt method. In addition, this leads to a significant reduction in the number of controls that is expected to be beneficial for gradient quality when using approximate gradients. We consider an ensemble of sector models where the uncertainty is described by different permeability, porosity, net/gross ratios, and initial water-saturation fields. The controls are the ICD settings over time in the three horizontal injection wells, with approximately 15 ICDs per well. Different optimized strategies resulting from different initial strategies were compared. We achieved a mean 4.2% increase in expected net present value (NPV) at a 10% discount rate compared with a traditional pressure-maintenance strategy. Next, we performed a sequential biobjective optimization and achieved an increase of 9.2% in the secondary objective (25% discounted NPV to emphasize shortterm production gains) for a minimal decrease of 1% in the primary objective (0% discounted NPV to emphasize long-term recovery gains), as averaged over the 100 geological realizations. The work flow was repeated for alternative numbers of ICDs, showing that having fewer control options lowers the expected value for this particular case. The results demonstrate that ensemble-based optimization work flows are able to produce improved robust recovery strategies for realistic field-sector models against acceptable computational cost. ... Read more

With an increase in the number of applications of ensemble optimization (EnOpt) for production optimization, the theoretical understanding of the gradient quality has received little attention. An important factor that influences the quality of the gradient estimate is the number of samples. In this study we use principles from statistical hypothesis testing to quantify the number of samples needed to estimate an ensemble gradient that is comparable in quality to an accurate adjoint gradient. We develop a methodology to estimate the necessary ensemble size to obtain an approximate gradient that is within a predefined angle compared to the adjoint gradient, with a predefined statistical confidence. The method is first applied to the Rosenbrock function (a standard optimization test problem), for a single realization, and subsequently for a case with uncertainty, represented by multiple realizations (robust optimization). The maximum allowed error applied in both experiments is a 10° angle between the directions of the EnOpt gradient and the exact gradient. For the single-realization case we need, depending on the perturbation size, 900, 5 and 3 samples to estimate a "good" gradient with 95% confidence at 50 points in the optimization space for 50 different random sequences. For the robust case, the conventional EnOpt approach is to couple one model realization with one control sample, which leads to a computationally efficient technique to estimate a mean gradient. However, our results show that in order to be 95% confident the original one-to-one model realization to control sample ratio formulation is not sufficient. To achieve the required confidence requires a ratio of 1:1100, i.e. each model realization is paired with 1100 control samples using the original formulation. However, using a modified formulation we need a ratio of 1:10 to stay within the maximum allowed error for 95% of the points in space, though a 1:1 ratio is sufficient for 85% of the points. We also tested our methodology on a reservoir case for deterministic and robust cases, where we observe similar trends in the results. Our results provide insight into the necessary number of samples required for EnOpt, in particular for robust optimization, to achieve a gradient comparable to an adjoint gradient. ... Read more