Determination of Methane Diffusion Coefficients in Live Oils for Tight Reservoirs at High Pressures

Molecular diffusion plays a critical role in gas injection in tight reservoirs such as liquid-rich shale. Despite recent efforts on measuring diffusion coefficients at high pressures, there is a general lack of the diffusion coefficients in live oil systems at reservoir conditions relevant to the development of these tight reservoirs. The reported diffusion coefficients often differ in orders of magnitude, and there is no consensus on the reliability of the common correlations for liquid phase diffusion coefficients, such as the extended Sigmund correlation. We employed the constant volume diffusion method to measure the high-pressure diffusion coefficients in a newly designed high-pressure tube. The experimental method was first validated using methane + hexadecane and methane + decane, and then used to measure the methane diffusion coefficients in two live oils at reservoir conditions. The obtained data were processed by compositional simulation to determine the diffusion coefficients. The diffusion coefficients measured for methane + hexadecane and methane + decane are in agreement with the existing literature data. For methane + live oil systems, however, the diffusion coefficients estimated by the extended Sigmund correlation are much lower than the measured results. An over ten times adjustment is needed to best fit the pressure decay curves. A further check reveals that for live oil systems, the reduced densities are often in the extrapolated region of the original Sigmund model. The curve in this region of the extended Sigmund correlation has a weak experimental basis, which may be the reason for its large deviation. The estimates from other correlations like Wilke-Chang and Hayduk-Minhas also give very different results. We compared the diffusion coefficients in high-pressure oils reported in the literature, showing a large variation in the reported values. All these indicate the necessity for further study on accurate determination of high-pressure diffusion coefficients in live oils of relevance to shale and other tight reservoirs.

Revisiting the Classical Experiment in Earthquake Control at the Rangely Oil Field, Colorado, 1970, Using a Coupled Flow and Geomechanical Model

ABSTRACT The injection experiment conducted at the Rangely oil field, Colorado, was a pioneering study that showed qualitatively the correlation between reservoir pressure increases and earthquake occurrence. Here, we revisit this field experiment using a mechanistic approach to investigate why and how the earthquakes occurred. Using data collected from decades of field operations, we build a geological model for the Rangely oil field, perform reservoir simulation to history match pore-pressure variations during the experiment, and perform geomechanical simulations to obtain stresses at the main fault, where the earthquakes were sourced. As a viable model, we hypothesize that pressure diffusion occurred through a system of highly permeable fractures, adjacent to the main fault in the field, connecting the injection wells to the area outside of the injection interval where intense seismic activity occurred. We also find that the main fault in the field is characterized by a friction coefficient μ ≈ 0.7—a value that is in good agreement with the classical laboratory estimates conducted by Byerlee for a variety of rock types. Finally, our modeling results suggest that earthquakes outside of the injection interval were released tectonic stresses and thus should be classified as triggered, whereas earthquakes inside the injection interval were driven mostly by anthropogenic pore-pressure changes and thus should be classified as induced.

Seismicity at the Castor gas reservoir driven by pore pressure diffusion and asperities loading

The 2013 seismic sequence at the Castor injection platform offshore Spain, including three earthquakes of magnitude 4.1, occurred during the initial filling of a planned Underground Gas Storage facility. The Castor sequence is one of the most important cases of induced seismicity in Europe and a rare example of seismicity induced by gas injection into a depleted oil field. Here we use advanced seismological techniques applied to an enhanced waveform dataset, to resolve the geometry of the faults, develop a greatly enlarged seismicity catalog and record details of the rupture kinematics. The sequence occurred by progressive fault failure and unlocking, with seismicity initially migrating away from the injection points, triggered by pore pressure diffusion, and then back again, breaking larger asperities loaded to higher stress and producing the largest earthquakes. Seismicity occurred almost exclusively on a secondary fault, located below the reservoir, dipping opposite from the reservoir bounding fault.

A unified modelling and simulation for coupled anomalous transport in porous media and its finite element implementation

This paper presents an unified mathematical and computational framework for mechanics-coupled “anomalous” transport phenomena in porous media. The anomalous diffusion is mainly due to variable fluid flow rates caused by spatially and temporally varying permeability. This type of behaviour is described by a fractional pore pressure diffusion equation. The diffusion transient phenomena is significantly affected by the order of the fractional operators. In order to solve 3D consolidation problems of large scale structures, the fractional pore pressure diffusion equation is implemented in a finite element framework adopting the discretised formulation of fractional derivatives given by Grunwald–Letnikov (GL). Here the fractional pore pressure diffusion equation is implemented in the commercial software Abaqus through an open-source UMATHT subroutine. The similarity between pore pressure, heat and hydrogen transport is also discussed in order to show that it is possible to use the coupled thermal-stress analysis to solve fractional consolidation problems.

Sound absorption of polydisperse heterogeneous porous composites

Sound absorption of polydisperse heterogeneous porous composites is investigated in this paper. The wave equation in polydisperse heterogeneous porous composites is upscaled by using the two-scale method of homogenisation, which allows the material to be modeled as an equivalent fluid with atypical effective parameters. This upscaled model is numerically validated and demonstrates that the dissipation of sound in polydisperse heterogeneous porous composites is due to visco-thermal dissipation in the composite constituents and multiple pressure diffusion in the polydisperse heterogeneous inclusions. Analytical and semi-analytical models are developed for the acoustical effective parameters of polydisperse heterogeneous porous composites with canonical geometry (e.g. porous matrix with cylindrical and spherical inclusions) and with complex geometries. Furthermore, by comparing the sound absorption coefficient of a hard-backed composite layer with that of layers made from the composite constituents alone, it is demonstrated that embedding polydisperse heterogeneous inclusions in a porous matrix can provide a practical way for significantly increasing low frequency sound absorption. The results of this work are expected to serve as a model for the rational design of novel acoustic materials with enhanced sound absorption properties.

Investigations into the opening of fractures during hydraulic testing using a hybrid-dimensional flow formulation

We applied a hybrid-dimensional flow model to pressure transients recorded during pumping experiments conducted at the Reiche Zeche underground research laboratory to study the opening behavior of fractures due to fluid injection. Two distinct types of pressure responses to flow-rate steps were identified that represent radial-symmetric and plane-axisymmetric flow regimes from a conventional pressure-diffusion perspective. We numerically modeled both using a radial-symmetric flow formulation for a fracture that comprises a non-linear constitutive relation for the contact mechanics governing reversible fracture surface interaction. The two types of pressure response can be modeled equally well. A sensitivity study revealed a positive correlation between fracture length and normal fracture stiffness that yield a match between field observations and numerical results. Decomposition of the acting normal stresses into stresses associated with the deformation state of the global fracture geometry and with the local contacts indicates that geometrically induced stresses contribute the more the lower the total effective normal stress and the shorter the fracture. Separating the contributions of the local contact mechanics and the overall fracture geometry to fracture normal stiffness indicates that the geometrical stiffness constitutes a lower bound for total stiffness; its relevance increases with decreasing fracture length. Our study demonstrates that non-linear hydro-mechanical coupling can lead to vastly different hydraulic responses and thus provides an alternative to conventional pressure-diffusion analysis that requires changes in flow regime to cover the full range of observations.

Spatiotemporal Properties of Seismicity and Variations of Shear-Wave Splitting Parameters in the Western Gulf of Corinth (Greece)

The Western Gulf of Corinth (WGoC) exhibits significant seismicity patterns, combining intense microseismic background activity with both seismic swarms and short-lived aftershock sequences. Herein, we present a catalogue of ~9000 events, derived by manual analysis and double-difference relocation, for the seismicity of the WGoC during 2013–2014. The high spatial resolution of the hypocentral distribution permitted the delineation of the activated structures and their relation to major mapped faults on the surface. The spatiotemporal analysis of seismicity revealed a 32-km-long earthquake migration pattern, related to pore-pressure diffusion, triggering moderate mainshock-aftershock sequences, as fluids propagated eastwards in the course of ~15 months. The anisotropic properties of the upper crust were examined through automatic shear-wave splitting (SWS) analysis, with over 2000 SWS measurements at local stations. An average fast shear-wave polarization direction of N98.8° E ± 2.8° was determined, consistent with the direction of the maximum horizontal regional stress. Temporal variations of normalized time-delays between fast and slow shear-waves imply alterations in the level of stress or microcrack fluid saturation during the long-lasting pore-pressure diffusion episode, particularly before major events. The present study provides novel insights regarding seismicity patterns, active fault structures, anisotropic properties of the upper crust and triggering mechanisms of seismicity in the WGoC.

Earthquakes and very deep groundwater perturbation mutually induced

We report unique observations from drilling and hydraulic stimulation at a depth of approximately 4.3 km in two Enhanced Geothermal System (EGS) wells at the Pohang EGS site, South Korea. We surveyed drilling logs and hydraulic stimulation data, simulated pore pressure diffusion around the fault delineated by seismic and drilling log analyses, conducted acoustic image logging through the EGS wells, observed significant water level drops (740 m) in one of the two EGS wells, and obtained hydrochemical and isotopic variation data in conjunction with the microbial community characteristics of the two EGS wells. We discuss the hydraulic and hydrochemical responses of formation pore water to a few key seismic events near the hypocenter. We focused on how the geochemistry of water that flowed back from the geothermal wells changed in association with key seismic events. These were (1) a swarm of small earthquakes that occurred when a significant circulation mud loss occurred during well drilling, (2) the MW 3.2 earthquake during hydraulic stimulation, and (3) the MW 5.5 main shock two months after the end of hydraulic stimulation. This study highlights the value of real-time monitoring and water chemistry analysis, in addition to seismic monitoring during EGS operation.