Investigation of relative flow characteristics of brine-saturated reservoir formation: A numerical study of the Hawkesbury formation

Abstract Precise knowledge of the relative flow behaviour of CO 2 and brine during CO 2 sequestration in deep saline aquifers and its impact on the sequestration process is required to ensure the safety and efficiency of sequestration projects. This numerical study therefore aims to identify the interaction-induced relative flow behaviour of reservoir rock and its influence on the hydro-mechanical and geochemical phenomena in deep saline aquifers. COMSOL Multiphysics numerical simulator was used to develop a laboratory-scale relative flow model to simulate the CO 2 movement and brine drainage processes in brine-saturated Hawkesbury sandstone samples and finally, was extended it into a field-scale numerical model which can simulate the hydro-mechanical, mineralogical and geochemical behaviours of deep saline aquifers under CO 2 sequestration conditions. A 2-D axisymmetric pore-elastic model was developed using the pore-elastic module available in COMSOL and Buckely-Leverett flow theory was applied to the model using pre-defined partial differential equations. The proposed laboratory-scale model was first validated using experimental permeability data conducted under triaxial drained conditions and the model was then extended to predict relative flow characteristics, such as brine and CO 2 saturation and CO 2 pressure distribution along the sample length under different injection pressures, including both sub- and super-critical conditions and finally, real CO 2 sequestration processes. According to the results, the developed laboratory-scale model simulates the experimental results reasonably well with less than 10% relative error. The numerical results also reveal that there is a considerable effect of CO 2 phase change on the final distribution of the CO 2 and brine saturation profiles. In addition, the CO 2 pressure distribution along the sample length shows a non-linear relationship between the CO 2 pressure and sample length. According to the results of the field-scale model, the long-term interaction of CO 2 causes Hawkesbury formation's pore structure to significantly change.

[1]  A. Lasaga.  Chemical kinetics of water‐rock interactions , 1984 .

[2]  K. Pruess,et al.  TOUGH2 User's Guide Version 2 , 1999 .

[3]  Mohamed Azaroual,et al.  Well injectivity during CO2 storage operations in deep saline aquifers – Part 2: Numerical simulations of drying, salt deposit mechanisms and role of capillary forces , 2014 .

[4]  David Airey,et al.  Numerical simulation of gas flow through porous sandstone and its experimental validation , 2011 .

[5]  B. Chun,et al.  Interfacial tension in high-pressure carbon dioxide mixtures , 1995 .

[6]  Larry W. Lake,et al.  Implications of Coupling Fractional Flow and Geochemistry for CO2 Injection in Aquifers , 2007 .

[7]  Emil O. Frind,et al.  Two‐phase flow in heterogeneous porous media: 2. Model application , 1991 .

[8]  Christopher A. Rochelle,et al.  The impact of chemical reactions on CO2 storage in geological formations: a brief review , 2004, Geological Society, London, Special Publications.

[9]  Ragnhild Korbol,et al.  Sleipner vest CO2 disposal - injection of removed CO2 into the Utsira formation , 1995 .

[10]  David Airey,et al.  Investigation of temperature effect on permeability of naturally fractured black coal for carbon dioxide movement: An experimental and numerical study , 2012 .

[11]  P. Ranjith,et al.  Experimental investigation of geochemical and mineralogical effects of CO2 sequestration on flow characteristics of reservoir rock in deep saline aquifers , 2016, Scientific Reports.

[12]  Nobuo Hirano,et al.  Determination of aperture structure and fluid flow in a rock fracture by high‐resolution numerical modeling on the basis of a flow‐through experiment under confining pressure , 2008 .

[13]  L. V. D. Meer.  Investigations regarding the storage of carbon dioxide in aquifers in the Netherlands , 1992 .

[14]  P. Ranjith,et al.  Determination of effective stress parameters for effective CO2 permeability in deep saline aquifers: An experimental study , 2015 .

[15]  Mohsen Hadian,et al.  Prediction of free flowing porosity and permeability based on conventional well logging data using artificial neural networks optimized by Imperialist competitive algorithm – A case study in the South Pars gas field , 2015 .

[16]  F. Dongyan,et al.  A composite model of hydraulic fractured horizontal well with stimulated reservoir volume in tight oil & gas reservoir , 2015 .

[17]  Long-biao Cheng,et al.  Non-stationary and nonlinear visco-elastic shear creep model for shale , 2011 .

[18]  P. Ranjith,et al.  CO2-induced mechanical behaviour of Hawkesbury sandstone in the Gosford basin: An experimental study , 2015 .

[19]  C. A. Tang,et al.  Coupled analysis of flow, stress and damage (FSD) in rock failure , 2002 .

[20]  James W. Mercer,et al.  A review of immiscible fluids in the subsurface: properties, models, characterization and remediation , 1990 .

[21]  McMillan Burton,et al.  CO2 injectivity into brine aquifers: Why relative permeability matters as much as absolute permeability , 2009 .

[22]  Pathegama Gamage Ranjith,et al.  Influence of CO2-brine co-injection on CO2 storage capacity enhancement in deep saline aquifers: An experimental study on Hawkesbury sandstone formation , 2016 .

[23]  Ruben Juanes,et al.  Impact of relative permeability hysteresis on geological CO2 storage , 2006 .

[24]  Y. Mualem.  A New Model for Predicting the Hydraulic Conductivity , 1976 .

[25]  L. Ayala,et al.  Experimental investigation of shale gas production impairment due to fracturing fluid migration during shut-in time , 2015 .

[26]  M. Azaroual,et al.  Solubility of silica polymorphs in electrolyte solutions, II. Activity of aqueous silica and solid silica polymorphs in deep solutions from the sedimentary Paris Basin , 1997 .

[27]  P. Aagaard,et al.  Thermodynamic and kinetic constraints on reaction rates among minerals and aqueous solutions. II. Rate constants, effective surface area, and the hydrolysis of feldspar , 1984 .

[28]  S. Bauer,et al.  Optimal use of a dome-shaped anticline structure for CO2 storage: a case study in the North German sedimentary basin , 2012, Environmental Earth Sciences.

[29]  Pathegama Gamage Ranjith,et al.  Salinity-dependent strength and stress–strain characteristics of reservoir rocks in deep saline aquifers: An experimental study , 2014 .

[30]  Z. Karpyn,et al.  Investigation of gas flow hindrance due to fracturing fluid leakoff in low permeability sandstones , 2014 .

[31]  Zhangxin Chen.  Formulations and Numerical Methods of the Black Oil Model in Porous Media , 2000, SIAM J. Numer. Anal..

[32]  N. Müller.  Supercritical CO2-Brine Relative Permeability Experiments in Reservoir Rocks—Literature Review and Recommendations , 2011 .

[33]  L. André,et al.  Numerical modeling of fluid–rock chemical interactions at the supercritical CO2–liquid interface during CO2 injection into a carbonate reservoir, the Dogger aquifer (Paris Basin, France) , 2007 .

[34]  C. Doughty,et al.  Modeling Supercritical Carbon Dioxide Injection in Heterogeneous Porous Media , 2004, Vadose Zone Journal.

[35]  Etienne Brosse,et al.  Modelling Fluid-Rock Interaction Induced by the Percolation of CO2-Enriched Solutions in Core Samples: the Role of Reactive Surface Area , 2005 .

[36]  Stephen P. White,et al.  Reservoir storage and containment of greenhouse gases , 1995 .

[37]  G. Dagan.  Statistical Theory of Groundwater Flow and Transport: Pore to Laboratory, Laboratory to Formation, and Formation to Regional Scale , 1986 .

[38]  R. H. Brooks,et al.  Hydraulic properties of porous media , 1963 .

[39]  Mark D. White,et al.  Comparison of relative permeability-saturation-capillary pressure models for simulation of reservoir CO2 injection , 2016 .

[40]  B. Jones,et al.  The Hawkesbury Sandstone South of Sydney, Australia: Triassic Analogue for the Deposit of a Large, Braided River , 1987 .

[41]  S. Bachu,et al.  Hydrogeological and numerical analysis of CO2 disposal in deep aquifers in the Alberta sedimentary basin , 1996 .

[42]  Martin J. Blunt,et al.  Predictive pore‐scale modeling of two‐phase flow in mixed wet media , 2004 .