Although hydraulic fracturing in Liquid-Rich Unconventional Reservoirs (LUR) have become a norm, the recovery factor continues to be low. Use of Enhanced Oil Recovery (EOR) techniques in LUR have recently become more popular to improve the recovery. The objective of this study is to numerically investigate the advantages and disadvantages of the application of CO2 huff-n-puff technique in the LUR formations having complex fracture networks. The study explores the fluid flow mechanisms for oil recovery in the naturally fractured reservoir. A calibrated 3D mechanical earth model with geomechanical and petrophysical property from the Eagle Ford was used for the study. Complex hydraulic fracture model was used to simulate the hydraulic fractures, proppant and fluid distribution around the wellbore. Numerical reservoir simulation on a Perpendicular Bi-section (PEBI) grid was used to capture the permeability, porosity and conductivity distribution due to the proppants in the hydraulic fractures. CO2 huff-n-puff technique using numerical reservoir simulation is used to determine the well performance and recovery factor arising from reservoir fluid viscosity reduction and gas expansion. Effect of fluid thermodynamics to recovery systems in the low permeability reservoir medium is fully captured in approach. Equation of state prepared for simulating the CO2 impact on the oil is prepared with correlating the collected down hole oil sample. Numerical reservoir simulation study coupled with the complex fracture simulation model presents the insights of new means to improve RF in LUR through the injection of CO2. Such EOR method would be critical to increase the long term economic benefits. The study demonstrates that the infill well requirements can be mitigated if the EOR method of Huff-n-puff is utilized in cyclic modes over various time periods of production. Up to 9% extra RF was observed when CO2 Huff-n-puff technique was used as compared to production dependent only on hydraulic fracture stimulation. Parametric sensitivity on job sizes and start timing of EOR in a producing well is used to evaluate the RF. However, the hydraulic fracture geometry and the created footprint along with the time of injection has a larger effect in improving the EOR effectiveness. The methodology provides the demonstration of simulating the EOR methods in unconventional reservoirs for economic assessment. The workflow demonstrates modeling CO2 flooding as an EOR technique on the full wellbore level with complex hydraulic fracture geometry. The approach demonstrated here can be applied to other basins in the unconventional formations to improve the recovery factor.

Although hydraulic fracturing in liquid-rich unconventional reservoirs (LUR) has become a norm, the recovery factor continues to be low. Use of enhanced oil recovery (EOR) techniques in LUR has recently become more popular to improve the recovery. The objective of this study is to numerically investigate the advantages and disadvantages of the application of the CO2 huff-n-puff technique in LUR formations having complex fracture networks. The study explores the fluid flow mechanisms for oil recovery in a naturally fractured reservoir. A calibrated 3D mechanical earth model with geomechanical and petrophysical properties from the Eagle Ford was used for the study. A complex hydraulic fracture model was used to simulate the hydraulic fracture, proppant, and fluid distribution around the wellbore. Numerical reservoir simulation on perpendicular bisection (PEBI) grids was used to capture the permeability, porosity, and conductivity distribution due to the proppants in the hydraulic fractures. The CO2 huff-n-puff technique using numerical reservoir simulation was used to determine the well performance and recovery factor arising from reservoir fluid viscosity reduction and gas expansion. The effect of fluid thermodynamics to recovery systems in the low-permeability reservoir medium was fully captured in this approach. An equation of state prepared for simulating the CO2 impact on the oil was prepared with correlating the collected downhole oil sample. The numerical reservoir simulation study coupled with the complex fracture simulation model presents insights into a new means to improve the recovery factor (RF) in LUR through the injection of CO2. Such EOR method would be critical to increase the long-term economic benefits. The study demonstrates that that infill well requirements can be mitigated if the EOR method of huff-n-puff is utilized in cyclic modes over various time periods of production. Up to 9% extra RF was observed when the CO2 huff-n-puff technique was used as compared to production dependent only on hydraulic fracture stimulation. Parametric sensitivity on job sizes and start timing of EOR in a producing well was used to evaluate the RF. However, the hydraulic fracture geometry and the created footprint along with the time of injection have a larger effect in improving the EOR effectiveness. The methodology demonstrates the simulation of EOR methods in unconventional reservoirs for economic assessment. The workflow demonstrates modeling CO2 flooding as an EOR technique on the full wellbore level with complex hydraulic fracture geometry. The approach can be applied to unconventional formations in other basins to improve the recovery factor.

With increased drilling activity associated with development of unconventional reservoirs, many operators are reporting both stimulation and production interference between wells. Interference between existing production wells (parent wells) and newly completed infill wells (child wells) is often associated with production impairment (Marongiu-Porcu et al. 2015; Ajisafe et al. 2017; Defeu et al. 2018 and Manchanda et al. 2018a). The objective of this work is to develop guidelines concerning infill wells completion scheme to minimize parent-child wells interference in a typical pad, with infill drilling in the Duvernay formation. The area of interest selected within the Duvernay formation consists of three parent wells and two child wells. An integrated mechanical earth model (MEM) was constructed for the area using public databases. The created 3D-geological model included petrophysical and geomechanical properties along with a 2D discrete natural-fracture network representing the distribution of natural fractures in the reservoir. Hydraulic fractures in parent wells were modeled using original stress settings from the 3D-MEM. Then, a dynamic model was constructed for the three parent wells and production simulation was run for five years. Pressure distribution at the time when child wells came into production was extracted and 3D depleted stress distribution was computed using a finite element method that included the effect of pore pressure decrease and principal stress magnitude and orientation changes. Then, hydraulic fracture modeling was performed for the two child wells using the new depleted stress distribution, and finally a five-well dynamic model was created. Sensitivity analyses were performed on the hydraulic fracture parameters of the child wells with the objective of maximizing recovery by accessing more virgin reservoir area between the parent wells. Hydraulic fracture modeling followed by dynamic simulation was done in the pad for multiple cases. Fracture geometry, hydraulic/propped surface area, and fracture conductivity in child wells were extracted and analyzed against production performance of the wells. This study shows a holistic approach in modeling the impact of completion modifications on the child wells performance in an infill drilling scenario. A 3D-geomechanical model coupled with reservoir simulation allowed simulating the propagation of hydraulic fractures in the presence of pressure depleted regions. Results confirmed that the main reason for under-performance of child wells in Duvernay is the stress change induced by the reservoir pressure depletion associated with the parent wells production hence, influencing the child wells hydraulic fractures propagation patterns.

Cyclic gas injection in hydraulically-fractured wells has been successfully applied as an enhanced oil recovery (EOR) method in tight unconventional basins such as the Permian and Eagle Ford. However, displacement processes (continuous gas, solvent or water injection) such as those piloted in the relatively more permeable Bakken formation have not been considered in tighter basins. In this work, we present a novel displacement process that uses alternating injecting and producing hydraulic fractures to flood the inter-fracture region around a horizontal well. We demonstrate that the method is feasible as long as suitable fracture geometries can be generated. A reservoir geomodel of a typical Montney gas condensate reservoir was constructed using publicly available data. Rock mechanics parameters were integrated into the model alongside completion and pumping schedule information to predict hydraulic fracture propagation and geometries representative of the typical stimulated volumes in Montney. A compositional numerical reservoir simulator was then used to test the proposed EOR process in a gas condensate reservoir where we forecast liquid recovery under different frac-to-frac continuous gas injection flooding scenarios. Sensitivities to study the effect of fracture spacing and complexity, solvent composition, starting time of gas injection, and matrix permeability were performed. With simultaneous injection and production from alternating hydraulic fractures, it is possible to flood the volumes between them and consequently avoid the drawbacks of huff ’n’ puff processes. By using a more rigorous fracture description, we can reproduce the interactions between fractures and determine how they affect the conformance of the displacement front. Modeling results showed that frac-to-frac displacement process can significantly improve the condensate recovery compared to primary or even huff n puff EOR process. They also showed that the frac-to-frac EOR process is feasible only if the formation mechanical properties and in-situ stresses are such that the resulting hydraulic fractures exhibit aligned planar geometries. If high-intensity natural fracture networks are present, the hydraulic fractures tend to form complex geometries that negatively affect the conformance of the flooding front. The study also showed that there is an optimal spacing between the injecting and producing fractures that would allow for the efficient utilization of the EOR agent; this spacing was shown to have a strong dependence on matrix permeability. Composition of injected solvent and starting time of gas injection doesn’t seem to have considerable impact on incremental recovery due to frac-to-frac displacement.

Optimized geomodeling and history matching of production data is presented by utilizing an integrated rock and fluid workflow. Facies identification is performed by use of image logs and other geological information. In addition, image logs are used to help define structural geodynamic processes that occurred in the reservoir. Methods of reservoir fluid geodynamics are used to assess the extent of fluid compositional equilibrium, especially the asphaltenes, and thereby the extent of connectivity in these facies. Geochemical determinations are shown to be consistent with measurements of compositional thermodynamic equilibrium. The ability to develop the geo-scenario of the reservoir, the coherent evolution of rock and contained fluids in the reservoir over geologic time, improves the robustness of the geomodel. In particular, the sequence of oil charge, compositional equilibrium, fault block throw, and primary biogenic gas charge are established in this middle Pliocene reservoir with implications for production, field extension,and local basin exploration. History matching of production data prove the accuracy of the geomodel; nevertheless, refinements to the geomodel and improved history matching were obtained by expanded deterministic property estimation from wireline log and other data. Theearly connection of fluid data, both thermodynamic and geochemical, with relevant facies andtheir properties determination enables a more facile method to incorporate this data into the geomodel. Logging data from future wells in the field can be imported into the geomodel allowingdeterministic optimization of this model long after production has commenced. While each reservoir is unique with its own idiosyncrasies, the workflow presented here is generally applicable to all reservoirs and always improves reservoir understanding.