Numerous fractured hydrocarbon reservoirs exist in the metamorphic basement of the Pannonian Basin in Hungary. Many decades of experience in production have proven that these reservoirs are highly compartmentalised, resulting in a complex mosaic of permeable and impermeable domains situated next to each other. Consequently, in most fields, only a small amount of the total hydrocarbon reserve can be extracted. This paper aims to locate the potential migration pathways inside the most productive basement reservoir of the Pannonian Basin, using a multiscale approach. To achieve this, evaluation well-log data, DFN modelling and a composition analysis of fluid trapped in a vein-filling zeolite phase are combined. Data on a single well are presented as an example. The results of the three approaches indicate the presence of two highly fractured intervals separated by a barely fractured amphibolite. The two zones are probably part of the communicating fracture system inside the single metamorphic mass. The gas analysis further specifies the migrated fluids and indicates hydrocarbons of a composition similar to that of the recently produced oil. Consequently, we conclude that the two zones do not only form an ancient migration pathway but are also members of a more recent hydrocarbon system.

Shale gas reservoirs (SGR) are important replacements for conventional energy resources and have been widely exploited by hydraulic fracturing technologies. On the one hand, due to the inherent ultra-low permeability and porosity, there is stress sensitivity in the reservoirs generally. On the other hand, hydraulic fractures and the stimulated reservoir volume (SRV) generated by the massive hydraulic fracturing operation have contrast properties with the original reservoirs. These two phenomena bring huge challenges in SGR transient pressure analysis. Although some works in the literatures have been done on the transient pressure analysis of multi-fractured horizontal wells in SGR, unfortunately, none of them has taken the stress sensitivity and spatially varying permeability of SRV zone into consideration simultaneously. To fill this gap, this paper first idealizes the SGR to be four linear composite regions. What’s more, SRV zone is further divided into sub-sections on the basis of non-uniform distribution of proppant within SRV zone which easily yields spatially varying permeability away from the main hydraulic fracture. The stress sensitivity is characterized by the varying permeability depended on the pore pressure. By means of perturbation transformation and Laplace transformation, an analytical multi-linear flow model (MLFM) is obtained and validated by the comparison with the previous model. On the basis of our model, the flow regimes are identified and the sensitivity analysis of critical parameters are conducted to further understand the transient pressure behaviors. The research results provided by this work are of significance for well test interpretation and production performance analysis of SGR.

The electrokinetics methods have a great potential to characterize hydrogeological processes in geological media, especially in complex hydrosystems such as fractured formations. In this work, we conceptualize fractured media as a bunch of parallel capillary fractures following the fractal size distribution. This conceptualization permits to obtain analytical models for both the electrical conductivity and the electrokinetic coupling in water saturated fractured media. We explore two different approaches to express the electrokinetic coupling. First, we express the streaming potential coupling coefficient as a function of the zeta potential and then we obtain the effective charge density in terms of macroscopic hydraulic and electrokinetic parameters of porous media. We show that when the surface electrical conductivity is negligible, the proposed models reduces to the previously proposed one based on a bundle of cylindrical capillaries. This model opens up a wide range of applications to monitor the water flow in fractured media.

The study of fractures in the subsurface is very important in unconventional reservoirs since they are the main conduits for hydrocarbon flow. For this reason, a variety of equivalent medium theories have been proposed for the estimation of fracture and fluid properties within reservoir rocks. We experimentally investigated the feasibility of applying the Galvin model in fractured tight stones. For this proposal, three artificial fractured tight sandstone samples with the same background porosity (11.7% ± 1.2%) but different fracture densities of 0.00, 0.0312, and 0.0624 were manufactured. The fracture thickness was 0.06 mm and the fracture diameter was 3 mm in all the fractured samples. Ultrasonic P- and S-wave velocities were measured at 0.5 MHz in a laboratory setting in dry and water-saturated conditions in directions at 0°, 45°, and 90° to the fracture normal. The results were compared with theoretical predictions of the Galvin model. The comparison showed that model predictions significantly underestimated P- and S- wave velocities as well as P-wave anisotropy in water-saturated conditions, but overestimated P-wave anisotropy in dry conditions. By analyzing the differences between the measured results and theoretical predictions, we modified the Galvin model by adding the squirt flow mechanism to it and used the Thomsen model to obtain the elastic moduli in high- and low-frequency limits. The modified model predictions showed good fits with the measured results.

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).

Numerous reservoirs that play a significant role in worldwide petroleum production and reserves contain fractures. Typically, the fractures must form a connected network for a reservoir to be classified as naturally fractured. Characterizing the reservoir with a focus on its fracture network is crucial for modeling and predicting production performance. To simplify the solution, dual continuum modeling techniques are commonly employed. However, to use continuum-scale approaches, properties such as average aperture, permeability, and matrix fracture interaction parameters must be assigned, making it necessary to improve fracture depiction and modeling methods. This study investigates a fractured reservoir with a low matrix permeability and a well-connected fracture network. The focus is on the impact of the hyrachical fracture network on the production performance of gas-based enhanced oil recovery methods. The Discrete Fracture Model (DFN) was utilized to create comprehensive two-dimensional models for three processes: Gas Injection (GI), Water Alternating Gas (WAG), and Foam Assisted Water Alternating Gas (FAWAG). Moreover, dimensionless numbers were employed to establish connections between properties across the entire fracture hierarchy, spanning from minor to major fractures and encompassing the fracture intensity. The results indicate that the FAWAG is more sensitive to fracture types and networks than the WAG and GI processes. Hence, the sensitivity of the individual EOR method to the fracture network requires a respective depth of description of the fracture network. However, other factors, such as reservoir fluids properties and fracture properties, might influence the recovery when the minor fracture networks are excluded. This study has determined that, among the enhanced oil recovery (EOR) techniques examined, the significance of the hyrachical depth of fracture networks diminishes as the ratio of major (primary fracture) aperture to the aperture of medium and minor fractures increases. Additionally, the impact of the assisted-gravity drainage method was greater with increased reservoir height; however, as the intensity ratio increased, the relative importance of medium and minor fracture networks decreased.

Modeling fluid flow in three-dimensional (3D) Discrete Fracture Networks (DFNs) is of relevance in many engineering applications, such as hydraulic fracturing, oil/gas production, geothermal energy extraction, nuclear waste disposal and CO₂ sequestration. A new Boundary Element Method (BEM) technique with discontinuous quadratic elements and a parallel Domain Decomposition Method (DDM) is presented herein for the simulation of the steady-state fluid flow in 3D DFN systems with wellbores, consisting of planar fractures having arbitrary properties and wellbore trajectories. Numerical examples characterized by DFNs of increasing complexity are investigated to evaluate the accuracy and efficiency of the presented technique. The results show that accurate solutions can be obtained with less nodes than with mesh-based methods (e.g. Finite Element Method). In addition, the DDM algorithm used provides a quite fast convergence. The simulation results of the fluid flow around intersections among traces (linear intersections between fractures), intersections between traces and a fracture boundaries, and wellbore intersections is accurate. Source code is available at : https://github.com/BinWang0213/PyDFN3D.