A numerical model simulating reactive transport and evolution of fracture permeability

Yasuhara, Hajime; Elsworth, Derek

doi:10.1002/nag.513

Cited by 60 publications

(43 citation statements)

References 33 publications

Supporting

Mentioning

Contrasting

Unclassified

Order By: Relevance

“…The fracture aperture is equal to L x . This is similar to what the Yasuhara and Elsworth [] model would predict for these conditions if no critical stress limitations were invoked. The porosity of the control volume decreases with time.…”

Section: Examplessupporting

confidence: 85%

“…Note that the highest excess stress 73 MPa is less than the critical stress, which is 77 MPa for quartz at this temperature. The Yasuhara and Elsworth [] model would have predicted no decrease of the aperture.…”

Section: Examplesmentioning

confidence: 99%

“…Compaction of geological media, especially quartz sand, has been studied experimentally, theoretically, and by modeling, accounting for stress dissolution and reprecipitation [ Revil , , ; Revil et al ., ; Yasuhara et al ., ]. Closing and also opening of fractures due to such phenomena have been studied [ Raj , ; Polak et al ., ; Yasuhara and Elsworth , ; Yasuhara et al ., , ; Liu et al ., ; Elsworth and Yasuhara , , ; Min et al ., ; Taron and Elsworth , ; Taron et al ., ; Zheng and Elsworth , ].…”

Section: Introductionmentioning

confidence: 99%

See 2 more Smart Citations

Stress‐mediated closing of fractures: Impact of matrix diffusion

Neretnieks

2014

JGR Solid Earth

View full text Add to dashboard Cite

Fracture apertures may decrease by several mechanisms when the fractures are subject to stress. This paper considers only stress-enhanced dissolution of the crystals on the stressed surfaces. First, it is argued that the stress-induced dissolution is active already at the smallest difference between effective stress on the stressed surfaces and the unstressed surface of a crystal when in contact with water. This is in contrast to the concept that there exists a critical stress below which, this does not happen, assumed in some earlier studies. Second, and in addition to, but independent of the first argument, it is shown that there is a very strong sink for the stress-enhanced dissolved crystal mass due to diffusion into the porous rock matrix. There, the solute reprecipitates on the crystal surfaces in the matrix, which are subject to lower stress than the crystals bearing the load in the fracture. Diffusion into the porous matrix of the rock has not previously been considered in this context. A simple model that includes this sink is developed, and it is shown that matrix diffusion can be the by far largest sink for the solute and can considerably increase the rate of closure of the fractures. It is further found that under some possibly not uncommon conditions, the dissolution rate of quartz crystals becomes essentially independent of the strength of the sinks for the dissolved silica.

show abstract

Section: Examplessupporting

confidence: 85%

Section: Examplesmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

See 1 more Smart Citation

Stress‐mediated closing of fractures: Impact of matrix diffusion

Neretnieks

2014

JGR Solid Earth

View full text Add to dashboard Cite

show abstract

“…In this section a benchmark test was established in order to verify a laboratory experiment operated by Yasuhara ([48], [49]), in which an artificial fracture sample of novaculite was injected with water for a flow-through experiment. In Yasuhara's experiment, the novaculite sample was split open and the fracture surfaces were profiled.…”

Section: B Yasuhara's Experimentsmentioning

confidence: 99%

Numerical Simulation of Solute Transport in Fractured Rock with Particle Tracking Method

Sun

Su³

et al. 2017

Proceedings of the 2017 International Conference on Applied Mathematics, Modeling and Simulation (AMMS 2017)

View full text Add to dashboard Cite

Abstract-In this work an open source scientific software platform OpenGeoSys (OGS) was utilized for the simulation of fluid flow and solute transport in fractured porous media. A particle based method was proposed which takes use of velocities calculated by finite element method (FEM) to obtain flow patterns, and in which the distribution of the solute mass is represented by particle plume and can be quantitatively obtained by counting the number of particles. Two benchmarks were developed in order to verify the proposed method: one was to simulate solute transport in two dimensional steady flow through infinite fractured porous media with one single crack; the other was to simulate solute transport in an artificial rock fracture in novaculite under laboratory conditions. The simulation results showed that the proposed particle tracking method agreed well with the classical FEM simulation methods. The particle plume could exhibit the flow pattern, as well as the detailed channeling effects and accelerations or decelerations in some regions. This method can be applied as a tool to observe the detailed structure of evolving contaminant plumes in fractured rock.

show abstract

“…The critical stress was then calculated as 79.7 MPa at T = 20 °C for quartz. The dissolution mass flux of pressure solution on contacting asperities was given by Yasuhara and Elsworth as

\frac{d M_{diss}^{P S}}{d t} = \frac{3 V_{m}^{2} ((), σ_{a} - σ_{c}) k_{+} ρ_{g} A_{c}}{R T},

where t is time, k + is the dissolution rate constant of the solid ( k + = 3.09 × 10 −13 mol m −2 s −1 at 20 °C), A c is the contact area of contacting asperity ( A c = π a 2 ), and R is the gas constant ( R = 8.31 J mol −1 k −1 ).…”

Section: Implications On Compaction Mechanism Of Fast Closure Of Fracmentioning

confidence: 99%

Contact mechanism of a rock fracture subjected to normal loading and its impact on fast closure behavior during initial stage of fluid flow experiment

Zhao

Jiang

et al. 2015

Num Anal Meth Geomechanics

View full text Add to dashboard Cite

SUMMARYFast closure of rock fractures has been commonly observed in the initial stage of fluid flow experiments at environmental temperatures under low or moderate normal stresses. To fully understand the mechanisms that drive this fast closure, the evolution of local stresses acting on contacting asperities on the fracture surfaces prior to fluid flow tests needs to be evaluated. In this study, we modeled numerically the asperity deformation and failure processes during initial normal loading, by adopting both elastic and elasticplastic deformation models for the asperities on a real rock fracture with measured surface topography data, and estimated their impact on initial conditions for fluid flow test performed under laboratory conditions. Compared with the previous models that simulate the normal contact of a fracture as the approach of two rigid surfaces without deformations, our models of deformable asperities yielded smaller contact areas and higher stresses on contacting asperities at a given normal stress or normal displacement. The results show that the calculated local stresses were concentrated on the contacts of a few major asperities, resulting in crushing of asperity tips. With these higher contact stresses, however, the predicted closure rates by pressure solution are still several orders of magnitude lower than that of the experimental measurements at the initial stage of fluid flow test. This indicates that single pressure solution may not likely to be the principal compaction mechanism for this fast closure, and that the damages on contacting asperities that occur during the initial normal loading stage may play an important role.

show abstract

A numerical model simulating reactive transport and evolution of fracture permeability

Cited by 60 publications

References 33 publications

Stress‐mediated closing of fractures: Impact of matrix diffusion

Stress‐mediated closing of fractures: Impact of matrix diffusion

Numerical Simulation of Solute Transport in Fractured Rock with Particle Tracking Method

Contact mechanism of a rock fracture subjected to normal loading and its impact on fast closure behavior during initial stage of fluid flow experiment

Contact Info

Product

Resources

About