1
|
Sun L, Li M, Abdelaziz A, Tang X, Liu Q, Grasselli G. An efficient 3D cell-based discrete fracture-matrix flow model for digitally captured fracture networks. INTERNATIONAL JOURNAL OF COAL SCIENCE & TECHNOLOGY 2023; 10:70. [PMID: 37928133 PMCID: PMC10620291 DOI: 10.1007/s40789-023-00625-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/21/2023] [Revised: 05/19/2023] [Accepted: 08/10/2023] [Indexed: 11/07/2023]
Abstract
Complex hydraulic fracture networks are critical for enhancing permeability in unconventional reservoirs and mining industries. However, accurately simulating the fluid flow in realistic fracture networks (compared to the statistical fracture networks) is still challenging due to the fracture complexity and computational burden. This work proposes a simple yet efficient numerical framework for the flow simulation in fractured porous media obtained by 3D high-resolution images, aiming at both computational accuracy and efficiency. The fractured rock with complex fracture geometries is numerically constructed with a cell-based discrete fracture-matrix model (DFM) having implicit fracture apertures. The flow in the complex fractured porous media (including matrix flow, fracture flow, as well as exchange flow) is simulated with a pipe-based cell-centered finite volume method. The performance of this model is validated against analytical/numerical solutions. Then a lab-scale true triaxial hydraulically fractured shale sample is reconstructed, and the fluid flow in this realistic fracture network is simulated. Results suggest that the proposed method achieves a good balance between computational efficiency and accuracy. The complex fracture networks control the fluid flow process, and the opened natural fractures behave as primary fluid pathways. Heterogeneous and anisotropic features of fluid flow are well captured with the present model.
Collapse
Affiliation(s)
- Lei Sun
- Department of Civil and Mineral Engineering, University of Toronto, Toronto, ON M5S 1A4 Canada
| | - Mei Li
- Department of Civil and Mineral Engineering, University of Toronto, Toronto, ON M5S 1A4 Canada
| | - Aly Abdelaziz
- Department of Civil and Mineral Engineering, University of Toronto, Toronto, ON M5S 1A4 Canada
| | - Xuhai Tang
- School of Civil Engineering, Wuhan University, Wuhan, 430072 China
| | - Quansheng Liu
- School of Civil Engineering, Wuhan University, Wuhan, 430072 China
| | - Giovanni Grasselli
- Department of Civil and Mineral Engineering, University of Toronto, Toronto, ON M5S 1A4 Canada
| |
Collapse
|
2
|
The Influence of Lamina Density and Occurrence on the Permeability of Lamellar Shale after Hydration. CRYSTALS 2021. [DOI: 10.3390/cryst11121524] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
The characteristics of laminae in lamellar shale oil reservoirs have important influences on reservoir parameters, especially permeability. In order to explore the influence of lamina density and occurrence on the permeability of lamellar shale after hydration, we studied the lamellar shale in the Chang 7 member of the Yanchang Formation of Triassic in Ordos Basin. By comparing the permeability of bedding shale and lamellar shale with different densities of laminae, it was found that the permeability anisotropy of lamellar shale was stronger. In the direction parallel to the lamina, the permeability increased approximately linearly with an increase in lamina density. The effect of hydration on rock micropore structure and permeability was studied by soaking shale in different fluids. Most of the microfracture in the lamellar shale was parallel to the lamina direction, and hydration led to a widening of the microfracture, which led to the most obvious increase in permeability parallel to the lamina. Collectively, the research results proved that lamina density, occurrence, and hydration have a significant influence on the permeability anisotropy of lamellar shale.
Collapse
|
3
|
The Intensity of Heat Exchange between Rock and Flowing Gas in Terms of Gas-Geodynamic Phenomena. ENTROPY 2021; 23:e23050556. [PMID: 33947172 PMCID: PMC8146585 DOI: 10.3390/e23050556] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/09/2021] [Revised: 04/21/2021] [Accepted: 04/27/2021] [Indexed: 12/03/2022]
Abstract
Gas-induced geodynamic phenomena can occur during underground mining operations if the porous structure of the rock is filled with gas at high pressure. In such cases, the original compact rock structure disintegrates into grains of small dimensions, which are then transported along the mine working space. Such geodynamic events, particularly outbursts of gas and rock, pose a danger both to the life of miners and to the functioning of the mine infrastructure. These incidents are rare in copper ore mining, but they have recently begun to occur, and have not yet been fully investigated. To ensure the safety of mining operations, it is necessary to determine parameters of the rock–gas system for which the energy of the gas will be smaller than the work required to disintegrate and transport the rock. Such a comparison is referred to as an energy balance and serves as a starting point for all engineering analyses. During mining operations, the equilibrium of the rock–gas system is disturbed, and the rapid destruction of the rock is initiated together with sudden decompression of the gas contained in its porous structure. The disintegrated rock is then transported along the mine working space in a stream of released gas. Estimation of the energy of the gas requires investigation of the type of thermodynamic transformation involved in the process. In this case, adiabatic transformation would mean that the gas, cooled in the course of decompression, remains at a temperature significantly lower than that of the surrounding rocks throughout the process. However, if we assume that the transformation is isothermal, then the cooled gas will heat up to the original temperature of the rock in a very short time (<1 s). Because the quantity of energy in the case of isothermal transformation is almost three times as high as in the adiabatic case, obtaining the correct energy balance for gas-induced geodynamic phenomena requires detailed analysis of this question. For this purpose, a unique experimental study was carried out to determine the time required for heat exchange in conditions of very rapid flows of gas around rock grains of different sizes. Numerical simulations reproducing the experiments were also designed. The results of the experiment and the simulation were in good agreement, indicating a very fast rate of heat exchange. Taking account of the parameters of the experiment, the thermodynamic transformation may be considered to be close to isothermal.
Collapse
|
4
|
Houben M, van Eeden JCM, Barnhoorn A, Hangx SJT. Fracture-Induced Permeability in Whitby Mudstone. ENVIRONMENTAL SCIENCE & TECHNOLOGY 2020; 54:9564-9572. [PMID: 32628456 PMCID: PMC7409142 DOI: 10.1021/acs.est.0c00557] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/28/2020] [Revised: 07/02/2020] [Accepted: 07/06/2020] [Indexed: 06/11/2023]
Abstract
Shale host rock and containment potential are largely determined by the connected pore network in the rock, and the connection between the pore network and the naturally present or mechanically induced fracture network together determines the total bulk permeability. Pore connectivity in shales is poorly understood because most of the porosity is present in sub-micrometer-sized pores that are connected through nanometer-sized pore throats. We have used a number of different techniques to investigate the microstructure and permeability of Early Jurassic shales from the UK (Whitby Mudstone), under intact and fractured conditions. Whitby Mudstone is a clay matrix-rich rock (50-70%), with different mineralogical layers on the sub-millimeter scale and very low natural permeability (10-19 to 10-22 m2), representative of many gas shales and caprocks present in Europe. Artificial fracturing of this shale increases its permeability by 2-5 orders of magnitude at low confining pressure (5 MPa). At high confining pressures (30 MPa), permeability changes were more sensitive to the measuring direction with respect to the bedding orientation. Given the distinct lack of well-defined damage zones, most of the permeability increase is controlled by fracture permeability, which is sensitive to the coupled hydro-chemo-mechanical response of the fractures to fluids.
Collapse
Affiliation(s)
- Maartje
E. Houben
- Department
of Earth Sciences, Utrecht University, 3585 CB Utrecht, The Netherlands
| | | | - Auke Barnhoorn
- Department
of Geoscience and Engineering, Delft University
of Technology, 2628 CN Delft, The Netherlands
| | - Suzanne J. T. Hangx
- Department
of Earth Sciences, Utrecht University, 3585 CB Utrecht, The Netherlands
| |
Collapse
|
5
|
The Effect of Hydration on Pores of Shale Oil Reservoirs in the Third Submember of the Triassic Chang 7 Member in Southern Ordos Basin. ENERGIES 2019. [DOI: 10.3390/en12203932] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Shale oil is an unconventional kind of oil and gas resource with great potential. China has huge reserves of shale oil, and shale oil resources are abundant in the third submember of the Triassic Chang 7 member in the southern Ordos Basin. At present, this area is in the initial stage of shale oil exploration and development. The reservoir pore is one of the key factors affecting oil accumulation, drilling safety, and oil production. It is also an important reservoir parameter that must be defined in the exploration stage. In general, the clay content in the shale section is high, and is prone to hydration. In order to study the effect of fluid on the pore type, structure, and distribution of shale oil reservoirs, experiments using X-ray diffraction, a porosity–permeability test, mercury porosimetry, rock casting thin section, and scanning electron microscopy were carried out. The experimental results show that the content of clay and quartz is very high in the studied formation. The pore porosity and permeability of the rock is highly heterogeneous because of the obvious stratigraphic bedding and interbeds. Microstructural observation of rocks shows that the main pore types are intergranular pores, intragranular pores, intercrystalline pores, and organic pores. Crack types are dissolution cracks, contraction cracks of organic matter, and abnormal pressure structural cracks. After hydration, the porosity of rock will increase in varying degrees, and pore size, pore content in different sizes, and pore structure will also change. The results show that the pores of tuff mainly changes at the initial stage of hydration, and the pore change of tuff is the most obvious within 6 hours of soaking in clear water. The influence of hydration on the pore of shale is greater than that of tuff, but the main change stage is later than tuff, and the pore change of shale is the most obvious within 12 to 24 hours of soaking in clear water. The soaking experiment of water-based drilling fluid (WBM-SL) shows that it can plug a certain size of holes and cracks and form a protective layer on the rock surface, thus effectively reducing hydration. In actual construction, multisized solid particles should be allocated in drilling fluid according to the formation pore’s characteristics, and the stability of the protective layer should be guaranteed. This can reduce the accident of well leakage and collapse and is conducive to the efficient and safe development of shale oil.
Collapse
|
6
|
Gomila R, Arancibia G, Mery D, Nehler M, Bracke R, Morata D. Palaeopermeability anisotropy and geometrical properties of sealed-microfractures from micro-CT analyses: An open-source implementation. Micron 2018; 117:29-39. [PMID: 30458300 DOI: 10.1016/j.micron.2018.11.001] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2018] [Revised: 11/05/2018] [Accepted: 11/05/2018] [Indexed: 11/29/2022]
Abstract
Fault zone permeability and the real 3D-spatial distribution of the fault-related fracture networks are critical in the assessment of fault zones behavior for fluids. The study of the real 3D-spatial distribution of the microfracture network, using X-ray micro-computed tomography, is a crucial factor to unravel the real structural permeability conditions of a fault-zone. Despite the availability of several commercial software for rock properties estimation from X-ray micro-computed tomography scanning, their high cost and lack of programmability encourage the use of open-source data treatment. This work presents the implementation of a methodology flow for the quantification of both structural and geometrical parameters (fractures density, fractures aperture, fractures porosity, and fractures surface area), and the modeling of palaeopermeability of fault-related fractured samples, with focus in the proper spatial orientation of both the sample and the results. This is performed with an easy to follow step-by-step implementation, by a combination of open-source software, newly implemented codes, and numerical methods. This approach keeps track of the sample's spatial orientation from the physical to the virtual world, thus assessing any fault-related palaeopermeability anisotropy.
Collapse
Affiliation(s)
- Rodrigo Gomila
- Departamento de Ingeniería Estructural y Geotécnica, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, 7820436, Santiago, Chile; Andean Geothermal Center of Excellence (CEGA, FONDAP-CONICYT), Universidad de Chile, Santiago, Chile.
| | - Gloria Arancibia
- Departamento de Ingeniería Estructural y Geotécnica, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, 7820436, Santiago, Chile; Andean Geothermal Center of Excellence (CEGA, FONDAP-CONICYT), Universidad de Chile, Santiago, Chile; Centro de Investigación en Nanotecnología y Materiales Avanzados (CIEN-UC), Pontificia Universidad Católica de Chile, Santiago, Chile
| | - Domingo Mery
- Andean Geothermal Center of Excellence (CEGA, FONDAP-CONICYT), Universidad de Chile, Santiago, Chile; Departamento de Ciencias de la Computación, Pontificia Universidad Católica de Chile, Santiago, Chile
| | | | - Rolf Bracke
- International Geothermal Centre (GZB), Bochum, Germany
| | - Diego Morata
- Andean Geothermal Center of Excellence (CEGA, FONDAP-CONICYT), Universidad de Chile, Santiago, Chile; Departamento de Geología, Universidad de Chile, Santiago, Chile
| |
Collapse
|
7
|
Experimental Study of the Microstructural Evolution of Glauberite and Its Weakening Mechanism under the Effect of Thermal-Hydrological-Chemical Coupling. Processes (Basel) 2018. [DOI: 10.3390/pr6080099] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
The microstructures of rock gradually evolve with changes in the external environment. This study focused on the microstructure evolution of glauberite and its weakening mechanism under different leaching conditions. The porosity were used as a characteristic index to study the effect of brine temperature and concentration on crack initiation and propagation in glauberite. The research subjects were specimens of ϕ3 × 10 mm cylindrical glauberite core, obtained from a bedded salt deposit buried more than 1000 m underground in the Yunying salt formation, China. The results showed that when the specimens were immersed in solution at low temperature, due to hydration impurities, cracks appeared spontaneously at the centre of the disc and the solution then penetrated the specimens via these cracks and dissolved the minerals around the crack lines. However, with an increase of temperature, the dissolution rate increased greatly, and crack nucleation and dissolved regions appeared simultaneously. When the specimens were immersed in a sodium chloride solution at the same concentration, the porosity s presented gradual upward trends with a rise in temperature, whereas, when the specimens were immersed in the sodium chloride solution at the same temperature, the porosity tended to decrease with the increase of sodium chloride concentration. In the process of leaching, the hydration of illite, montmorillonite, and the residual skeleton of glauberite led to the expansion of the specimen volume, thereby producing the cracks. The diameter expansion rate and the expansion velocity of the specimen increased with temperature increase, whereas, due to the common-ion effect, the porosity of the specimen decreases with the increase of sodium chloride solution concentration.
Collapse
|
8
|
Three-dimensional quantitative fracture analysis of tight gas sandstones using industrial computed tomography. Sci Rep 2017; 7:1825. [PMID: 28500297 PMCID: PMC5431860 DOI: 10.1038/s41598-017-01996-7] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2016] [Accepted: 04/05/2017] [Indexed: 12/09/2022] Open
Abstract
Tight gas sandstone samples are imaged at high resolution industrial X-ray computed tomography (ICT) systems to provide a three-dimensional quantitative characterization of the fracture geometries. Fracture networks are quantitatively analyzed using a combination of 2-D slice analysis and 3-D visualization and counting. The core samples are firstly scanned to produce grayscale slices, and the corresponding fracture area, length, aperture and fracture porosity as well as fracture density were measured. Then the 2-D slices were stacked to create a complete 3-D image using volume-rendering software. The open fractures (vug) are colored cyan whereas the calcite-filled fractures (high density objects) are colored magenta. The surface area and volume of both open fractures and high density fractures are calculated by 3-D counting. Then the fracture porosity and fracture aperture are estimated by 3-D counting. The fracture porosity and aperture from ICT analysis performed at atmospheric pressure are higher than those calculated from image logs at reservoir conditions. At last, the fracture connectivity is determined through comparison of fracture parameters with permeability. Distribution of fracture density and fracture aperture determines the permeability and producibility of tight gas sandstones. ICT has the advantage of performing three dimensional fracture imaging in a non-destructive way.
Collapse
|