scholarly journals Geological Model of a Storage Complex for a CO2 Storage Operation in a Naturally-Fractured Carbonate Formation

Geosciences ◽  
2018 ◽  
Vol 8 (9) ◽  
pp. 354 ◽  
Author(s):  
Yann Le Gallo ◽  
José de Dios

Investigation into geological storage of CO2 is underway at Hontomín (Spain). The storage reservoir is a deep saline aquifer formed by naturally fractured carbonates with low matrix permeability. Understanding the processes that are involved in CO2 migration within these formations is key to ensure safe operation and reliable plume prediction. A geological model encompassing the whole storage complex was established based upon newly-drilled and legacy wells. The matrix characteristics were mainly obtained from the newly drilled wells with a complete suite of log acquisitions, laboratory works and hydraulic tests. The model major improvement is the integration of the natural fractures. Following a methodology that was developed for naturally fractured hydrocarbon reservoirs, the advanced characterization workflow identified the main sets of fractures and their main characteristics, such as apertures, orientations, and dips. Two main sets of fracture are identified based upon their mean orientation: North-South and East-West with different fracture density for each the facies. The flow capacity of the fracture sets are calibrated on interpreted injection tests by matching their permeability and aperture at the Discrete Fracture Network scale and are subsequently upscaled to the geological model scale. A key new feature of the model is estimated permeability anisotropy induced by the fracture sets.

Author(s):  
Yann Le Gallo ◽  
José Carlos de Dios

Investigation into geological storage of CO2 is underway at the Technology Development Plant (TDP) at Hontomín (Burgos, Spain), the only current onshore injection site in the European Union. The storage reservoir is a deep saline aquifer located within Low Jurassic Formations (Lias and Dogger), formed by fractured carbonates with low matrix permeability. Understanding the processes involved in CO2 migration within this kind of low-primary permeability carbonates influenced by fractures and faults is key to ensure safe operation and reliable plume prediction. During the hydraulic characterization tests, 2300 tons of liquid CO2 and 14000 m3 of synthetic brine were co-injected on site in various sequences to characterize the pressure response of the seal-storage pair [de Dios et al, 2017] The injection tests were analyzed with a compositional dual media model which accounts for both temperature effects (as the CO2 is liquid at the bottom of the wellbore) and multiphase flow hysteresis (to effectively simulate the alternating brine and CO2 injection tests that were performed). The pressure and temperature responses of the storage formation were history-matched mainly through the petrophysical characteristics of the fracture network [Le Gallo et al, 2017]. The dynamic characterization of the fracture network dominates the CO2 migration while the matrix does not appear to significantly contribute to the storage capacity. This initial modeling approach was improved using the characterization workflow developed within the European FP7 CO2ReMove project for sandstone fractured reservoirs [Ringrose et al., 2011; Deflandre et al., 2011]. Fractured reservoirs are challenging to handle because of their high level of heterogeneity that conditions the reservoir behaviour during the injection. In particular, natural fractures have a significant impact on well performance [Ray et al, 2012]. Furthermore, the understanding of the processes involved in CO2 migration within relatively low-permeability storage influenced by fractures and faults is essential for enabling safe storage operation [Iding and Ringrose, 2010]. As part of the European H2020 ENOS project, the site geological model is updated by integration of the recently acquired data such as the image log interpretations from injection and observation wells. The geological model is generated through the analysis and integration of data including borehole images and well test data. Following a methodology developed for naturally fractured hydrocarbon reservoirs [Ray et al., 2012], the image log analysis identified two sets of diffuse fractures. A Discrete Fracture Network [Bourbiaux et al., 2005] was built around both wells which encompass the caprock, storage and underburden formations. The fracture characteristics of the two sets of diffuse fractures, such as orientations, densities and conductivities, are calibrated upon the interpretation of the injection tests [Le Gallo et al, 2017]. For each facies, the DFN characteristics were then upscaled and propagated to the full-field reservoir simulation model as 3D fracture properties (fracture porosity, fracture permeability and equivalent block size).


2019 ◽  
Vol 142 (3) ◽  
Author(s):  
Xiangnan Liu ◽  
Daoyong Yang ◽  
Andrew Chen

Abstract In this paper, pragmatic and robust techniques have been developed to simultaneously interpret absolute permeability and relative permeability together with capillary pressure in a naturally fractured carbonate formation from wireline formation testing (WFT) measurements. By using two sets of pressure and flow rate field data collected by a dual-packer tool, two high-resolution cylindrical near-wellbore numerical models are developed for each dataset on the basis of single- and dual-porosity concepts. Then, simulations and history matchings are performed for both the measured pressure drawdown and buildup profiles, while absolute permeability is determined and relative permeability is interpreted with and without considering capillary pressure. Compared to the experimentally measured relative permeability curves for the same formation collected from the literature, relative permeability interpreted with consideration of capillary pressure has a better match than those without considering capillary pressure. Also, relative permeability obtained from dual-porosity models has similar characteristics to those from single-porosity models especially in the region away from the endpoints, though the computational expenses with dual-porosity models are much larger. Absolute permeabilities in the vertical and the horizontal directions of the upper layer are determined to be 201.0 mD and 86.4 mD, respectively, while those of the lower layer are found to be 342.9 mD and 1.8 mD, respectively. Such a large vertical permeability of the lower layer reflects the contribution of the extensively distributed natural fractures in the vertical direction.


2020 ◽  
Vol 8 (4) ◽  
pp. SP109-SP133 ◽  
Author(s):  
Heloise Bloxsom Lynn ◽  
Bill Goodway

A 3D P-P high-fold full-azimuth full-offset reflection survey was acquired and processed to characterize a naturally fractured carbonate reservoir. The reservoir is a thick carbonate, which will flow commercial oil with a sufficient fracture network. Extensive calibration data include (1) a horizontal borehole’s resistivity image log, (2) the first 24 months cumulative oil produced, by stage, as known from chemical frac tracer data, (3) pre- and postfrac job instantaneous shut-in pressures, (4) microseismic, and (5) wireline log data. We used the cumulative oil production to document the spatially varying amount of aligned vertical porosity (aligned compliance or fracture porosity) connected to the stage borehole location. The stages of high oil production exhibited, for the fracture-perpendicular azimuth, the more positive amplitude variation with angle (AVA) gradients, and dimmer near-angle (6°–15° angles of incidence) amplitudes, compared to the fracture-parallel azimuth. The azimuthal variation of the AVA gradient fit the cos 2θ curve well, indicating the presence of one set of vertical aligned fractures dominating the azimuthal amplitude signature. In a similar fashion, the azimuthal variation of the mathematical intercept, physically the near-angle amplitudes, also fit the cos 2θ curve well. We have constructed crossplots of the azimuthal near-angle amplitude versus the AVA gradient on a bin-by-bin basis: we observed a straight line at bins with elevated oil production (elevated fracture density). A straight line crossplot of the (AVA gradient, mathematical intercept) is the signature of change of the (sensed) porosity, as long as the lithology and pore fluid are held constant. In accord with industry knowledge, we found that porosity affects the P impedance and thus the near-angle amplitudes: the aligned porosity yields azimuthal P impedance (measured at the 6°–15° angles of incidence). Legacy high-fold 3D P-P surveys rich in the 6°–20° angles of incidence should be considered for reprocessing and reinterpretation using these techniques.


2019 ◽  
Vol 3 (2) ◽  
pp. 23 ◽  
Author(s):  
Posadas-Mondragón ◽  
Camacho-Velázquez

In the oil industry, many reservoirs produce from partially penetrated wells, either to postpone the arrival of undesirable fluids or to avoid problems during drilling operations. The majority of these reservoirs are heterogeneous and anisotropic, such as naturally fractured reservoirs. The analysis of pressure-transient tests is a very useful method to dynamically characterize both the heterogeneity and anisotropy existing in the reservoir. In this paper, a new analytical solution for a partially penetrated well based on a fractal approach to capture the distribution and connectivity of the fracture network is presented. This solution represents the complexity of the flow lines better than the traditional Euclidean flow models for single-porosity fractured reservoirs, i.e., for a tight matrix. The proposed solution takes into consideration the variations in fracture density throughout the reservoir, which have a direct influence on the porosity, permeability, and the size distribution of the matrix blocks as a result of the fracturing process. This solution generalizes previous solutions to model the pressure-transient behavior of partially penetrated wells as proposed in the technical literature for the classical Euclidean formulation, which considers a uniform distribution of fractures that are fully connected. Several synthetic cases obtained with the proposed solution are shown to illustrate the influence of different variables, including fractal parameters.


Sign in / Sign up

Export Citation Format

Share Document