Traditionally, fractured reservoir simulations use Dual-Porosity, Dual-Permeability (DPDK) models that can idealize fractures and misrepresent connectivity. The Embedded Discrete Fracture Modeling (EDFM) approach improves flow predictions by integrating a realistic fracture network grid within a structured matrix grid. However, small fracture cells with high conductivity that pose a challenge for simulators can arise and ad hoc strategies to remove them can alter connectivity or fail for field-scale cases. We present a new gridding algorithm that controls the geometry and topology of the fracture network while enforcing a lower bound on the fracture cell sizes. It honors connectivity and systematically removes cells below a chosen fidelity factor. Furthermore, we implemented a flexible grid coarsening framework based on aggregation and flow-based transmissibility upscaling to convert EDFMs to various coarse representations for simulation speedup. Here, we consider pseudo-DPDK (pDPDK) models to evaluate potential DPDK inaccuracies and the impact of strictly honoring EDFM connectivity via Connected Component within Matrix (CCM) models. We combine these components into a practical workflow that can efficiently generate upscaled EDFMs from stochastic realizations of thousands of geologically realistic natural fractures for ensemble applications. We first consider a simple waterflood example to illustrate our fracture upscaling to obtain coarse (pDPDK and CCM) models. The coarse simulation results show biases consistent with the underlying assumptions (e.g., pDPDK can over-connect fractures). The preservation of fracture connectivity via the CCM aggregation strategy provides better accuracy relative to the fine EDFM forecast while maintaining computational speedup. We then demonstrate the robustness of the proposed EDFM workflow for practical studies through application to an improved oil recovery (IOR) study for a fractured carbonate reservoir. Our automatable workflow enables quick screening of many possibilities since the generation of full-field grids (comprising almost a million cells) and their preprocessing for simulation completes in a few minutes per model. The EDFM simulations, which account for complicated multiphase physics, can be generally performed within hours while coarse simulations are about a few times faster. The comparison of ensemble fine and coarse simulation results shows that on average, a DPDK representation can lead to high upscaling errors in well oil and water production as well as breakthrough time while the use of a more advanced strategy like CCM provides greater accuracy. Finally, we illustrate the use of the Ensemble Smoother with Multiple Data Assimilation (ESMDA) approach to account for field measured data and provide an ensemble of history-matched models with calibrated properties.

We have developed an efficient approach of petroleum reservoir model calibration that integrates 4D seismic surveys together with well-production data. The approach is particularly well-suited for the calibration of high-resolution reservoir properties (permeability) because the field-scale seismic data are areally dense, whereas the production data are effectively averaged over interwell spacing. The joint calibration procedure is performed using streamline-based sensitivities derived from finite-difference flow simulation. The inverted seismic data (i.e., changes in elastic impedance or fluid saturations) are distributed as a 3D high-resolution grid cell property. The sensitivities of the seismic and production surveillance data to perturbations in absolute permeability at individual grid cells are efficiently computed via semianalytical streamline techniques. We generalize previous formulations of streamline-based seismic inversion to incorporate realistic field situations such as changing boundary conditions due to infill drilling, pattern conversion, etc. A commercial finite-difference flow simulator is used for reservoir simulation and to generate the time-dependent velocity fields through which streamlines are traced and the sensitivity coefficients are computed. The commercial simulator allows us to incorporate detailed physical processes including compressibility and nonconvective forces, e.g., capillary pressure effects, while the streamline trajectories provide a rapid evaluation of the sensitivities. The efficacy of our proposed approach was tested with synthetic and field applications. The synthetic example was the Society of Petroleum Engineers benchmark Brugge field case. The field example involves waterflooding of a North Sea reservoir with multiple seismic surveys. In both cases, the advantages of incorporating the time-lapse variations were clearly demonstrated through improved estimation of the permeability heterogeneity, fluid saturation evolution, and swept and drained volumes. The value of the seismic data integration was in particular proven through the identification of the continuity in reservoir sands and barriers, and by the preservation of geologic realism in the calibrated model.

Summary We present the development and field application of a workflow for multiscale reservoir-model calibration that seamlessly integrates production data into the reservoir description from the facies to the grid-cell scale. To start with, the permeability field is parameterized using a novel grid-connectivity-based transformation basis that can be applied with any model geometry, including unstructured and corner-point grids. The parameterization basis functions emerge from spectral decomposition of the grid-connectivity Laplacian and are related to the structural harmonics of the grid. To reconcile data with model resolution during history matching, we first use the coarsest-scale basis functions to identify the large-scale variability. Additional smaller-scale basis elements are then adaptively incorporated to successively refine the model to a level supported by data resolution. During refinement, the inclusion of more detailed basis functions into the parameterization is determined by generic modal frequency when the prior model is unavailable or by using prior information when available. In the final step of the workflow, a streamline-based inversion is performed to locally adjust the reservoir model at grid-cell resolution along preferential-flow paths defined during the coarser-scale parameterization. We demonstrate the suitability and effectiveness of the developed workflow through application to an offshore turbidite reservoir with frequent well intervention, including shut-ins and recompletions. The static model has over 300,000 cells, a complex channelized interpretation with faults, four injector/producer pairs with deviated wells, and over eight years of production history, including water cut and pressure data. The grid-connectivity-based parameterization effectively updates the prior regional permeability at scales and in locations warranted by the data, while preserving the geologic continuity and avoiding ad hoc redefinition of regions given the sparse well pattern. The multiscale calibrated-permeability field indicates flow communication previously unrecognized in static geologic interpretation or manual history matching.

The effects of heterogeneity in Carbon Capture and Storage (CCS) in saline aquifers have been investigated extensively and are known to have important bearings on the storage capacity of the aquifer. In CCS projects, the time-lapse seismic survey has been proposed as a valuable tool for monitoring of CO 2 movement. However, the potential of the time-lapse seismic data for heterogeneity characterization and geologic model updating has not been fully explored. One of the biggest challenges in the quantitative use of time-lapse seismic data during CCS is the complex movement of the CO 2 influenced by compositional effects, geochemical reactions, phase changes and gravity segregation.In this paper, we first introduce compositional streamlines to understand and visualize the flow and transport of CO 2 in the presence of mineral precipitation/dissolution, residual trapping and buoyancy effects. To start with, individual component fluxes are generated by a finite difference fully implicit compositional simulator incorporating all the relevant physics of CO 2 sequestration. The fluxes are then utilized in novel streamline tracing algorithms to generate phase and component streamlines depicting the movement and the trapping of CO 2 in the aquifer. Next, we utilize the compositional streamlines to determine the sensitivity of the time-lapse seismic attributes specifically, interpreted saturation differences, to changes in reservoir properties such as permeability and porosity. The sensitivities are then used in an inverse modeling algorithm to calibrate the geologic model to time-lapse seismic data. The outcome is an improved description of permeability heterogeneity that is consistent with the 4-D seismic response and improved predictions of the CO 2 storage capacity.We have investigated the benefits of time-lapse seismic data integration in improving the performance assessment of CO 2 sequestration using examples involving CO 2 injection under realistic conditions. The first example examines the value of the 4-D seismic data integration in the estimation of storage capacity. The second example systematically studies the impact of viscous to gravity ratio on the performance of time-lapse seismic monitoring and heterogeneity characterization during CCS.