Fracture Network Model Research Articles

Exploring the impact of fracture interaction on connectivity and flow channelling using grown fracture networks

Quantitative assessment of the flow properties and mechanical stability of naturally fractured rock is frequently practised across the mining, petroleum, geothermal, geological disposal, construction and environmental remediation industries. These fluid and mechanical behaviours are strongly influenced by the connectivity of the fracture system and the size of the intact rock blocks. However, these are amongst the more difficult fracture system properties to characterize and honour in numerical simulations. Nonetheless, they are still the product of interactions between fractures that can be conceptualized as a series of deformation events following geomechanical principles. Generating numerical models of fracture networks by simulating this deformation with a coupled and evolving rock mass and stress field is a significant undertaking. Instead, large-scale fracture network models can be ‘grown’ dynamically according to rules that mimic the underlying mechanical processes and deformation history. This paper explores a computationally efficient rules-based method to generate fracture networks, demonstrates how different types of fracture patterns can be simulated, and illustrates how inclusion of fracture interactions can affect flow and mechanical properties. Relative to methods without fracture interaction and in contrast to some other rules-based approaches, the method described here regularizes and increases fracture connectivity and decreases flow channelling.

Review of discrete fracture network characterization for geothermal energy extraction

Geothermal reservoirs are highly anisotropic and heterogeneous, and thus require a variety of structural geology, geomechanical, remote sensing, geophysical and hydraulic techniques to inform Discrete Fracture Network flow models. Following the Paris Agreement on reduction of carbon emissions, such reservoirs have received more attention and new techniques that support Discrete Fracture Network models were developed. A comprehensive review is therefore needed to merge innovative and traditional technical approaches into a coherent framework to enhance the extraction of geothermal energy from the deep subsurface. Traditionally, statistics extracted from structural scanlines and unmanned aerial vehicle surveys on analogues represent optimum ways to constrain the length of joints, bedding planes, and faults, thereby generating a model of the network of fractures. Combining borehole images with seismic attributes has also proven to be an excellent approach that supports the stochastic generation of Discrete Fracture Network models by detecting the orientation, density, and dominant trends of the fractures in the reservoirs. However, to move forward to flow modelling, computation of transmissivities from pumping tests, and the determination of hydraulically active fractures allow the computation of the hydraulic aperture in permeable sedimentary rocks. The latter parameter is fundamental to simulating flow in a network of discrete fractures. The mechanical aperture can also be estimated based on the characterization of geomechanical parameters (Poisson’s ratio, and Young’s modulus) in Hot Dry Rocks of igneous-metamorphic origin. Compared with previous review studies, this paper will be the first to describe all the geological and hydro-geophysical techniques that inform Discrete Fracture Network development in geothermal frameworks. We therefore envisage that this paper represents a useful and holistic guide for future projects on preparing DFN models.

Differences in the permeability assessment of the fractured reservoir rocks using the conventional and the rough discrete fracture network modeling

Permeability assessment of naturally fractured rocks and fractured rocks after fracturing is critical to the development of oil and gas resources. In this paper, based on the discrete fracture network (DFN) modeling method, the conventional discrete fracture network (C-DFN) and the rough discrete fracture network (R-DFN) models are established. Through the seepage numerical simulation of the fractured rocks under different DFN, the differences in permeability of the fractured rocks under different parameters and their parameter sensitivity are analyzed and discussed. The results show that unconnected and independent fractures in the fracture network may weaken the seepage capacity of the fractured rocks. The fractured rock permeability increases with increase in connectivity and porosity and decreases with increase in maximum branch length and fracture dip. The use of C-DFN to equate the fracture network in the fractured rocks may underestimate the connectivity of the fracture network. For the more realistic R-DFN, the promotion of gas flow by connectivity is dominant when connectivity is high, and the hindrance of gas flow by fracture roughness is dominant when connectivity is low or when it is a single fracture. The permeability of the fractured rocks with R-DFN is more sensitive to the parameters than that of the fractured rocks with C-DFN. The higher the connectivity and porosity of the fractured rocks, the more obvious the difference between the permeability of the fractured rocks evaluated by C-DFN and R-DFN.

Performance evaluation of enhanced geothermal systems with intermittent thermal extraction for sustainable energy production

Efficient extraction of thermal energy from deep geothermal reservoirs using Enhanced Geothermal Systems (EGS) is fundamental for sustainable clean energy generation. This study presents the Intermittent Thermal Extraction (ITE) method, addressing thermal breakthrough challenges encountered during continuous thermal extraction (CTE) within EGS reservoirs. Employing a thermal-hydraulic coupled approach and discrete fracture network modeling, the ITE process for EGS reservoirs is assessed, omitting considerations of stress-induced deformation. The research evaluates the ITE and CTE methods within a typical high-temperature geothermal reservoir in the Himalayan geothermal belt. Analysis includes variations in reservoir temperatures and pressures, alongside an examination of rock matrix permeability and injection fluid temperature impacts on EGS thermal extraction performance. Results demonstrate the ITE method's efficacy in prolonging the EGS reservoir lifespan by preventing thermal breakthroughs and maintaining pressure equilibrium. Under studied conditions, ITE extends the EGS reservoir lifespan by a minimum of 17.7 years, amplifying clean power generation by 13.1% compared to CTE. Additionally, ITE implementation substantially reduces greenhouse gas emissions by an estimated 1.46–5.03 million tons. The rock matrix permeability and injection fluid temperature significantly affect the EGS thermal extraction performance when the fracture aperture is constant. These findings highlight the ITE method's potential to optimize thermal extraction from EGS reservoirs, particularly in remote regions. This study contributes essential insights into advancing geothermal energy extraction techniques, emphasizing ITE as a promising approach for sustainable energy production.

Effect of fracture characteristics on groundwater flow adjacent to high-level radioactive waste repository

In the study, a 2-D numerical model-delineating the high-level radioactive waste repository surrounded by hosting granite-rock was designed to evaluate pressure build-up, elevated temperature, and fracture flow caused by heat generation at the repository. During the early stage (0.1 years), build-up pressure at the repository induced radial groundwater flow. On the other hand, buoyancy-driven vertical flow was dominant during the late stage (2,000 years). The pressure build-up and elevated temperature varied significantly in bentonite, excavation-damaged zone and granite, resulting in different groundwater velocity. Fracture flow in both single and intersecting fractures were also influenced by the pressure and temperature, but their impact varied depending on the geometrical properties (e.g. orientation, elevation, lateral position) of the fracture. In complex discrete fracture network models, build-up pressure in the early stage mainly generated preferential fracture flow through few pathways, while in the late stage, convection-circulating flow induced by the buoyancy was distinct. In the early stage, higher density and connectivity of fracture network did not lead to an increase in average fracture flow, whereas the circulating flow in the late stage became more frequent and larger. Especially, long fractures served as both preferential flow pathways in the early stage and conduits for circulating flow in the late stage. The circulating flow was highly dependent on the spatial distribution of the long fractures, causing a significant uncertainty of average fracture flow in the late stage. The low angle of intersection between fractures resulted in low connectivity that decelerated the average fracture flow in both stages.

Development of an improved three-dimensional rough discrete fracture network model: Method and application

Structure plane is one of the important factors affecting the stability and failure mode of rock mass engineering. Rock mass structure characterization is the basic work of rock mechanics research and the important content of numerical simulation. A new 3-dimensional rough discrete fracture network (RDFN3D) model and its modeling method based on the Weierstrass-Mandelbrot (W-M) function were presented in this paper. The RDFN3D model, which improves and unifies the modelling methods for the complex structural planes, has been realized. The influence of fractal dimension, amplitude, and surface precision on the modeling parameters of RDFN3D was discussed. The reasonable W-M parameters suitable for the roughness coefficient of JRC were proposed, and the relationship between the mathematical model and the joint characterization was established. The RDFN3D together with the smooth 3-dimensional discrete fracture network (DFN3D) models were successfully exported to the drawing exchange format, which will provide a wide application in numerous numerical simulation codes including both the continuous and discontinuous methods. The numerical models were discussed using the COMSOL Multiphysics code and the 3-dimensional particle flow code, respectively. The reliability of the RDFN3D model was preliminarily discussed and analyzed. The roughness and spatial connectivity of the fracture networks have a dominant effect on the fluid flow patterns. The research results can provide a new geological model and analysis model for numerical simulation and engineering analysis of jointed rock mass.

An integrated 3D digital model of stratigraphy, petrophysics and karstified fracture network for the Cristal Cave, NE-Brazil

Digital Outcrop Models (DOMs) are virtual representations of geological features. Although DOMs are widely used tools in geosciences, their integration with other datasets remains relatively underexplored. We combined a DOM, derived from a photogrammetric survey of a carbonate sequence, with lithostratigraphic, petrophysical (porosity, permeability and uniaxial compressive strength), fracture distribution, and karst dissolution information to compose a single integrated three-dimensional digital model. The study site is one of the entrances of the Cristal Cave (São Francisco Craton, Northeastern Brazil), which has been used as a structural and diagenetic outcrop analog for the Brazilian pre-salt carbonate reservoirs. Data from fracture distributions, measured on the exposed surfaces of the cave, were used to build a Discrete Fracture Network, based on the solution of the stereology inverse problem. Fracture apertures were then modified to generate different scenarios of karstification, thus composing Discrete Fracture and Karst Networks. This integrative approach brought relevant insights into the cave development due to dissolution along fracture clusters. Our methodology offers better geological data handling to build static models to be used in a fluid flow modeling environment, contributing to bridge the gap between geophysics/geology and engineering approaches.

Spatial Characterization of Shallow Structures in the Revell Batholith Integrating Seismic Imaging Techniques

The Revell Site, located in Northwestern Ontario within the Canadian Shield, is being assessed as a potential Deep Geological Repository (DGR) for Canada’s used nuclear fuel. Effective DGR establishment requires comprehensive subsurface assessment, particularly in evaluating structural stability, hydrogeological attributes, geological composition, and geochemical properties. Key among these considerations is understanding the three-dimensional characteristics of structural features to ensure the site’s suitability for long-term containment of radioactive materials. This case study focuses on imaging and characterizing structures within the predominantly intact biotite granodiorite-tonalite host rock at the Revell Site. Borehole data reveals these structures as mostly sub-horizontal mafic intrusions, with thicknesses of up to 3.5 m. They appear either as separate, discrete entities or in stacked configurations. Despite their limited thickness, most of these intrusions exhibit discernible attributes in surface seismic images due to their distinct physical properties. Some, however, exhibit velocity variations that decrease their overall reflectivity properties. An integrative approach employing surface seismic, Vertical Seismic Profiles (VSP), and borehole data facilitates the spatial identification of over 30 of these mafic structures. This comprehensive characterization not only lays the foundational framework for future discrete fracture network models but also provides important support for simulations related to fluid flow dynamics, groundwater behavior, contaminant dispersion, and heat transport mechanisms within the Revell Site. This study underscores the critical role of in-depth subsurface characterization in ensuring the secure, long-term management of radioactive materials in DGRs.

Fracture Network Analysis of Karstified Subis Limestone Build-Up in Niah, Sarawak

Understanding complex carbonate fracture networks and karstification at various geological scales is challenging, especially with limited multi-scale datasets. This paper aims to reduce uncertainty in the fracture architecture of Central Luconia karstified reservoirs by narrowing observational gaps between seismic and well data by using the discrete fracture models of exposed limestone outcrops as analogues for the subsurface carbonate reservoir. An outcrop-based fracture network characterisation of a near-surface paleo-karst at Subis Limestone combined with lineament analysis was conducted to extract fracture parameters. The karst structure was first delineated using a digital elevation map and outcrop examination. Then, topology analysis was performed, following the creation of two-dimensional discrete fracture models. Two main fracture sets oriented northeast–southwest and northwest–southeast and 79 potential dolines were identified. Fracture intersections, northeast–southwest major orientations, and drainage systems highly influenced the karst features. The Subis Limestone fracture model revealed that the highest number of fractures and total length of fractures were concentrated in the northern part of the Subis Limestone build-up (X: 250–350, Y: 150–250) and became denser towards the northwest direction of the outcrop (X: 600–800). The fractures in the Subis paleo-karsts appear isolated, with I-nodes ranging from 0.74 to 0.94. Hence, it is crucial to incorporate matrix porosity into multiple scales of fracture network modelling to improve upscaling and the modelling of fracture–vug networks, as well as to minimise the underestimation of discrete fracture networks in fractured and karstified limestone.

A Numerical Method for Evaluating the Collapse of High-Steep Scarp Slopes Based on the Bonded Block Model–Discrete Fracture Network Model

Geotechnical engineering works in deep-incised valleys or open-pit mining areas often encounter high-steep scarp slopes with a slope angle greater than 75°. This type of slope directly threatens the safety of construction personnel, so assessing their stability is essential to ensure construction safety. The natural geometry of high-steep scarp slopes possesses complexity in terms of geometric morphology, structural features of rock mass, and occurrence mechanisms of collapse. There is little research and less emphasis on the evaluation of the collapse risk of high-steep scarp slopes. In particular, the fracture of intact rock or rock bridges is generally ignored in the analysis of collapse processes. A bonded block model (BBM)–discrete fracture network (DFN) coupling characterization model for the high-steep scarp slope is proposed based on a high-steep scarp slope containing dominant joint sets on the left bank of the dam site of the Huangzangsi Water Conservancy Project (Qinghai Province, China). By using the model, the complex geometric forms of the surface of the high-steep scarp slope are quantified, and the fracture process of falling rock masses as well as the controlling effect of dominant joints on the collapse of the scarp slope are revealed. A strength reduction method based on the BBM–DFN model is constructed, and the safety factor of the collapse-prone scarp slope is evaluated. The research results show that (1) the BBM–DFN model can be used to describe the local collapse process; (2) the occurrence of dominant joints plays an important part in controlling the collapse process; (3) there are differences in the safety factor of the scarp slope with different coupling methods; the collapse and failure modes also differ. For safety considerations, the safety factor of the scarp slope on the left bank of the dam site area is determined to be 1.85. The research findings can be used to guide the safety assessment of high-steep scarp slopes and the formulation of both collapse risk prevention and control measures to ensure construction safety in high-steep scarp slope areas.

3D Reconstruction of Pore and Fracture Structures in Crushed Soft Coal with Low Permeability and Its Permeability Analysis: Implications for the Design of CO2 Injection and CH4 Production.

In order to improve the CO2 injection and CH4 production efficiencies during the CO2-ECBM process, it is necessary to clarify the relationship among the complexity of pore and fracture structures, the typicality of the fluid migration path, and the heterogeneity of reservoir permeability. In this study, crushed soft coal with low permeability from Huainan and Huaibei coalfields of China was taken as the research object. First, the three-dimensional (3D) visualization reconstruction of pore and fracture structures was realized. Second, the equivalent pore and fracture network model was constructed. Finally, the permeability evolution and its anisotropy of the coal reservoir were dynamically demonstrated. In this study, the implication of surface porosity on the heterogeneity of pore and fracture structures was first discussed, followed by the implication of coordination number on the anisotropy of fluid flow, and finally, the influence of the anisotropy of fluid flow on the CO2-ECBM process was discussed. The results show that the equivalent pore and fracture network models of the reservoir structure can be constructed based on the digital rock physics technology. The analysis results of porosity, interconnected porosity, typical path of fluid migration, absolute permeability, and surface porosity of each sample have good consistency in characterizing the complexity of pore and fracture structures and the heterogeneity of permeability. The average coordination numbers of RL and LZ samples are 5.99 and 5.78, respectively, and the number of pores and throats is well-balanced, which indicates that LZ and RL collieries are suitable for the development of CO2-ECBM industrial tests. When the interconnected pores and fractures are mainly developed vertically and horizontally, the construction of drilling technology of the CO2-ECBM process should be mainly designed for vertical wells and horizontal wells, respectively. This study has important theoretical and practical significance for the industrial testing and commercialization of CO2-ECBM technology in crushed soft coal with low permeability.

Relationship between box-counting fractal dimension and properties of fracture networks

Due to the capacity to quantify the complexity of a fracture network, fractal dimension (D) is widely applied in analyzing fracture network-related issues, such as connectivity and permeability. While the relationship between D and individual properties of a fracture network has been extensively studied, D is influenced by a combination of multiple attributes of the fracture network. Therefore, this work utilizes multivariate analysis to establish an equation for predicting D, taking into account various properties of the fracture network, namely fracture length, number, and orientation. Monte Carlo simulation is employed to generate a substantial number of fracture network models. Subsequently, relationships between the fractal dimension (D) and various properties are derived for three types of fracture networks with (1) invariant fracture length and random orientation, (2) exponential fracture length and random orientation, and (3) exponential fracture length and von-Mises orientation. The initial analysis focuses on the simplest relationship, wherein the fundamental formula of fractional expression is determined. Then the second and third relationships are obtained through replacing the fixed parameter in the first relationship with the distribution parameters of fracture properties. Correlation analyses between the predicted D and the actual values reveal a remarkably high correlation (>0.99). To validate the established relationships, a fracture network obtained from geological outcrops is utilized. The results demonstrate the validity of the derived relationships. The utilization of these equations enhances the efficiency, practicality, and convenience of estimating fractal dimensions from fracture properties. As a result, the analysis of fracture network-related issues becomes more feasible and accessible.

Investigation of contact thermal resistance at the probe-sample interface in scanning thermal microscopy based on the fractal network model through numerical analysis

Scanning Thermal Microscopy (SThM) enables precise thermal mapping at the nanoscale. Nevertheless, it is challenging to quantitatively determine material thermal properties since the heat conduction at the scanning probe-sample interface is difficult to be measured directly. The heat conduction in the micro-domain must be clarified to enable accurate quantitative testing of thermophysical properties via SThM. This study aims to investigate the contact situation at the probe-sample interface using fractal theory and proposes a numerical model for calculating the contact thermal resistance. To verify the model’s accuracy, various surface geometries are specially designed. Furthermore, the impact of the sample surface fractal characteristics on the test accuracy is discussed through localized thermal signal feedback analysis.

Study on the influence of statistical geometrical characteristics on the permeability tensor of fractured rock masses by the equivalent pipe network method

The permeability of the host rock of the nuclear waste underground repository is one of the most keys for estimating the potential risk for nuclide waste to escape. Therefore, a three-dimensional discrete fracture network (DFN) model using the Monte Carlo method is generated based on the statistical geometrical characteristics of fracture network in Beishan I area which is selected as the site for constructing the nuclear waste repository in China. And the permeability of Beishan I fractured rock masses is calculated through the equivalent pipe network (EPN) method and representative elementary volume (REV) of it is determined. A rotating algorithm is developed to calculate different direction permeability. Furthermore, the influence of fracture statistical geometrical characteristics on REV size is explored by creating a series of DFN models with varying log-normal distribution of side lengths. The results show that the REV size of the Beishan Ⅰ area is 1.66 × 105 m3 and the corresponding principal permeability direction is 111.29°/75.36°. The REV size and permeability both increase with the increment of the mean and variance of the log-normal distribution increases. And with the increase of mean and variance, the growth rate gradually accelerated.

Heat Production Capacity Simulation and Parameter Sensitivity Analysis in the Process of Thermal Reservoir Development

The development of a geothermal system involves changes in the temperature field (T), seepage field (H), stress field (M), and chemical field (C) and the influence among them and injecting the heat extraction working fluid into the injection well that flows (migrating) through natural fractures and exchanges heat with the geothermal high-temperature rock. At the same time, the injection of low-temperature working fluid will induce thermal stress, resulting in changes in the reservoir temperature field and stress field. To study the influence factors and influence degree of heat production performance and mining life under multi-field coupling in the process of thermal reservoir development, based on THMC multi-field coupling numerical simulation software, this paper deeply studies the control differential equations and boundary coupling conditions of rock mass (fracture) deformation, seepage, heat exchange, the chemical reaction, and other processes based on the numerical solution method of the discrete fracture network model, simulating heat production capacity during the deep geothermal resource extraction process. The reservoir geological model analysis and generalization, parameter setting, boundary conditions, initial condition settings, mesh generation, and other steps were carried out in turn. Two different heat extraction working fluids, water, and CO2 were selected for numerical simulation in the mining process. The changes in the thermal reservoir temperature, net heat extraction rate, and SiO2 concentration during the thirty years of systematic mining were compared. The results show that CO2 has a better heat extraction effect. Finally, the reservoir thermal conductivity, heat capacity, well spacing, injection temperature, fracture spacing, fracture permeability, fracture number, fracture length, and other parameters were set, respectively. The parameter variation range was set, and the parameter sensitivity analysis was carried out. The numerical simulation results show that the engineering production conditions (injection temperature, well spacing) have little effect on the thermal efficiency and mining life, and the properties of fractures (fracture permeability, fracture number, fracture length) have a great influence.

Mechanical behaviors and damage characteristics of excavated jointed rocks within 3D printed structures

The presence of structural planes in rock masses reduces their mechanical properties and greatly affects their failure mode, resulting in physical and mechanical anisotropy. This study investigated the mechanical behaviors and damage characteristics of excavated jointed rocks in three-dimensional (3D) printed structures. A digital fracture network model based on the fractal principle was also developed to simulate a complex and rough network of fractures. Using 3D printing technology, rock specimens within a single rough structure surface, different dip angles, and rough fracture network models were established, considering different excavation sizes. Uniaxial compression was performed experimentally to investigate the mechanical response of rock masses. The failure processes of the jointed specimens were analyzed based on the digital image correlation (DIC) technique. The uniaxial compressive strength and elastic modulus were compared considering different excavation sizes. The results show that the geometric characteristics of the joints significantly influence the strength and failure mode of the rock mass. The strength of the layered rock mass shows a significant linear correlation with increasing excavation sizes. This research could provide a reference for studying the mechanical behavior of jointed rock masses using 3D printing technology.

Discontinuity development patterns and the challenges for 3D discrete fracture network modeling on complicated exposed rock surfaces

Natural slopes usually display complicated exposed rock surfaces that are characterized by complex and substantial terrain undulation and ubiquitous undesirable phenomena such as vegetation cover and rockfalls. This study presents a systematic outcrop research of fracture pattern variations in a complicated rock slope, and the qualitative and quantitative study of the complex phenomena impact on three-dimensional (3D) discrete fracture network (DFN) modeling. As the studies of the outcrop fracture pattern have been so far focused on local variations, thus, we put forward a statistical analysis of global variations. The entire outcrop is partitioned into several subzones, and the subzone-scale variability of fracture geometric properties is analyzed (including the orientation, the density, and the trace length). The results reveal significant variations in fracture characteristics (such as the concentrative degree, the average orientation, the density, and the trace length) among different subzones. Moreover, the density of fracture sets, which is approximately parallel to the slope surface, exhibits a notably higher value compared to other fracture sets across all subzones. To improve the accuracy of the DFN modeling, the effects of three common phenomena resulting from vegetation and rockfalls are qualitatively analyzed and the corresponding quantitative data processing solutions are proposed. Subsequently, the 3D fracture geometric parameters are determined for different areas of the high-steep rock slope in terms of the subzone dimensions. The results show significant variations in the same set of 3D fracture parameters across different regions with density differing by up to tenfold and mean trace length exhibiting differences of 3–4 times. The study results present precise geological structural information, improve modeling accuracy, and provide practical solutions for addressing complex outcrop issues.

Permeability evaluation based on Nuclear Magnetic Resonance analysis for gas hydrate reservoir

Permeability (or hydraulic conductivity) is one of the most important parameters for analyzing hydraulic behavior of underground. For example, in the geotechnical engineering, permeability is necessary to calculate water inflow and draw down of water level during excavation, to predict the deformation time due to the consolidation, and to evaluate the barrier of the ground for radioactive waste disposal. In the petroleum industry, permeability is considered an essential parameter which is directly linked to the productivity of oil and natural gas. Therefore, predicting permeability using wireline logging or logging while drilling (LWD) is required during the assessment of the reservoir. Nuclear Magnetic Resonance (NMR) logging is performed at this point. Low-field NMR measurements are sensitive to the hydrogen proton, 1H, in liquid water but not to the 1 H bond in solid form, making it possible to obtain the bulk water volume. In addition, the hydrogen proton NMR- T2 (transverse or spin-spin) relaxation time distribution is used to infer pore-size information and to estimate permeability [1, 2]. The major methods used for estimating permeability, k, are the Schlumberger Doll Research (SDR) model and Timur-Coates (TC) model [3, 4]:
 Where fNMR is the NMR-estimated (effective) porosity. SDR model uses the averaged value of the T2 distribution plotted on log-scale of time and matches it to the experimental reference permeability with the fitting parameter a [mD/(ms)2]. The TC model seeks to differentiate between the fraction of free (mobile) fluid volume (FFV) and bound fluid volume (BFV) and estimates permeability with FFV/BFV ratio and the material constant c [m2/s2]. The permeability in the research and development of methane hydrate, which is expected to be a next-generation resource, has been estimated using above methods. Methane hydrate is an ice-like solid crystalline structure of water and methane molecules. Because hydrate is solid, NMR is not sensitive, the initial (effective) permeability can be estimated by NMR analyzer. On the other hand, A strong disconnect (up to 2 orders of magnitude) between NMR-based estimates of permeability with those determined from a laboratory study of pressure cores from the same reservoir has further cast uncertainty in the application of borehole NMR data in hydrate evaluation [5, 6]. Recently, authors provided comparisons between NMR-derived values and direct fluid flow estimates of permeability in naturally occurring hydrate-bearing sediment [7, 8]. On the basis of the NMR measurement for hydrate-bearing pressure core sediments, which were recovered from the Gulf of Mexico and Alaska North Slope, the following new hydraulic radius model was proposed [7].
 This equation estimates the specific surface based on the NMR signal and has been shown to successfully estimate the initial permeability of a hydrate reservoir by expressing a hydraulic radius model based on Kozeny's equation [9]. However, the proposed equation is not universal and has been found to overestimate for clayey and silty hydrate free sediment. This study summarizes the methods and parameters used to estimate permeability using NMR signals from hydrate-bearing sediment obtained from logging and laboratory tests for the Gulf of Mexico and Alaska North Slope. In addition, we reports parameters for predicting model for muddy sediments which was recovered from overburden of hydrate-reservoir in the Hyuganada sea in 2022 Japan [10, 11]. When using the hydraulic radius model to predict permeability of hydrate-bearing sediment, reasonable estimates can be obtained by setting the shape factor for grains Cs = 1/2 (for a circular grain) and the shape factor for the flow path cn = 4 (for the cylinder network), respectively. Although, for hydrate free muddy sediment, Cs = 1/3 for parallel plate grain shape and cn = 2 for the fracture network model matched its permeability. Furthermore, the conversion parameter from T2 relaxation time to the pore size for hydrate free muddy sediments r2 = 0.009 (mm/s) was confirmed based on the result of mercury injection capillary pressure analysis, which is lower than r2 = 0.0325 (mm/s) of hydrate-bearing sediments.

Influence of Liquid Splitting Behavior at Intersections on Infiltration Dynamics in an Unsaturated Fracture Network

AbstractUnderstanding the macro‐scale flow characteristics in the fractured vadose zone is of great importance for subsurface hydrological and environmental applications. Here we develop an idealized fracture network model composed of a series of linked intersections, aiming to reveal the roles of local fluid flow, storage and splitting behaviors at intersections in controlling macroscopic unsaturated flow in a fracture network. By setting different local flow rules, unsaturated flow dynamics and structure in networks are systematically investigated. Numerical simulations suggest that avalanche infiltration mode, that is, sudden release of a large amount of water from the network, emerges spontaneously, and is modulated by the local splitting behavior. In terms of water infiltration structure, we find that uneven liquid splitting at intersections inevitably leads to preferential flow in networks. The preferential paths strongly rely on the competition between gravitational and capillary forces, relating to network structure and water saturation. When flow is dominated by gravity, preferential paths extend along large‐aperture routes; conversely, when flow is dominated by capillarity, the small‐aperture paths are preferred. In terms of infiltration dynamics, we show that the power spectral density of the water saturation time series in the network follows a power law with an exponent of −2 for all simulations with different structural parameters and local flow rules, suggesting a universal self‐organized criticality behavior for unsaturated flow in fractured rocks. The improved understanding from this study may shed new light on the complex flow dynamics for unsaturated flow in fractured media.

Existence analysis of hydraulic conductivity representative elementary volume in fractured rocks based on three-dimensional discrete fracture network method

The complexity of fractures in rocks results in heterogeneity and anisotropy of their hydraulic properties. The existence of hydraulic conductivity representative elementary volume (KREV) is a fundamental question in comprehending the hydraulic behavior of fractured rock masses. This article introduces an analytical procedure for determining the existence of KREV in fractured rock masses, which considers the distribution of underground apertures using a discrete fracture network (DFN) model. The procedure follows four basic steps: inverse modeling to generate a 3D DFN; calibration of the effective aperture by combining the results of packer tests and flow simulation using finite difference methods; calculation of the equivalent conductivities in 10 domains and 31 different flow directions; analysis of the fit error between the conductivity ellipsoid to determine the existence of KREV. The proposed approach is applied to the fractured rocks surrounding the 3013 exploration adit near the underground powerhouse of the Three Gorges Project as a case study. The results show that most fractures in the study area have apertures ranging from 0.06 to 0.25 mm and follow a lognormal distribution. The KREV size is 14–16 times the fracture spacing. Investigation of this model has significant implications for understanding flow through fracture networks in related fields.