No abstract

Large numbers of flow simulations are typically required for the determination of optimal well settings. These simulations are often computationally demanding, which poses challenges for the optimizations. In this paper we present a new two-step surrogate treatment (ST) that reduces the computational expense associated with well control optimization. The method is applicable for oil production via waterflood, with well rates optimized at a single control period. The two-step ST entails two separate optimizations, which can both be performed very efficiently. In the first optimization, optimal well-rate ratios (i.e., the fraction of total injection or production associated with each well) are determined such that a measure of velocity variability over the field is minimized, leading to more uniform sweep. In the second step, overall injection and production rates are determined. The flow physics in the first step is highly simplified, while the actual physical system is simulated in the second step. Near-globally-optimal results can be determined in both cases, as the first optimization is posed as a QP problem, and the second step entails just a single optimization variable. Under full parallelization, the overall elapsed time for the ST corresponds to the runtime for 1-2 full-order simulations. Results are presented for multiple well configurations, for 2D and 3D channelized models, and comparisons with formal optimization procedures (mesh adaptive direct search or MADS, and an adjoint-gradient method) are conducted. Three different fluid mobility ratios (M = 1, 3 and 5) are considered. Optimization results demonstrate that the two-step ST provides results in reasonable agreement with those from MADS and adjoint-gradient methods, with speedups of 5× or more. We also show that the ST is applicable in the inner-loop in field development optimization, where it will be especially useful since many different well configurations must be evaluated.

Seismic data incorporation in reservoir simulation models history matching (HM) studies has been continuously growing. 4D seismic data, in contrast with well production data, can provide a very good scenario of fluids arrangement along reservoir. In this work we describe how 3D and 4D seismic data gathered in acquisitions performed in Campos Basin was incorporated in Marlim Sul deep water field geological model reconstruction and in assisted HM (AHM).It is taken advantage of both 3D and 4D seismic data in several stages of the study, for instance, in the construction of a new porosity -most influential in impedance -model by using a methodology based on the inversion of synthetic seismic (calculated by petro-elastic model) in porosity through an optimization process that aims to reduce the difference between observed and synthetic impedance, and when defining influential parameters based on fluids displacement registered by seismic signal, by using a technique based on the creation of transmissibility multipliers parameters regions that considers the fluids displacement shown in 4D signal. Another relevant point is the use of information from reservoir and 3D seismic data when weighting the 4D data in the objective function.Combining the above mentioned techniques with the knowledge of the field -supported by the 3D seismic data -which allowed, for instance, identification of faults -where fault transmissibility multipliers were used as parameters in the HM process -a fairly good agreement on the observed well and seismic production data was achieved. HM studies using AHM tools have been shown a much more time-efficient technique when compared to manual HM. The incorporation of 4D seismic data can considerably improve the HM quality by improving the reservoir description, once it increases the ability of describing fluids arrangement and pressure distribution. The techniques successfully applied in the Marlim Sul field HM support these conclusions.