Self-diffusion studies of pore fluids in unconsolidated sediments by PFG NMR

The self-diffusion of water and hexadecane in medium and coarse sands from glacial sand deposits in central Germany were investigated by pulsed field gradient nuclear magnetic resonance (PFG NMR). Due to the restriction of the diffusion path at the pore/grain interface, the measured apparent self-diffusion coefficients (D(Δ)) in the pore space depend on the observation time (Δ) in the PFG NMR experiment. Although the bulk self-diffusion coefficients of water and hexadecane differ by about one order of magnitude, the apparent self-diffusion coefficients in the pore space obey the same characteristic time-behaviour, which depends only on geometrical properties of the pore system. Using the “short-time diffusion” model, surface-to-volume (S/V) ratios and inherent self-diffusion coefficients (D₀) of the pore fluids were extracted from these diffusion measurements. The S/V ratios obtained are independent of the pore fluid used and agree with known geometrical properties of the sand grains. Moreover, the D₀ values are consistent with the corresponding bulk self-diffusion coefficients measured separately. In contrast to these results of PFG NMR, simultaneous investigations of longitudinal (T₁) nuclear magnetic relaxation reveal that the relaxation time of the pore fluid is a less suitable parameter for a quantitative estimation of geometrical properties of the pore/grain interface in these unconsolidated sediments since it depends on chemical properties of the fluid/grain interface.     

A three-scale model for ionic solute transport in swelling clays incorporating ion–ion correlation effects

A new three-scale model is proposed to describe the movement of ionic species of different valences in swelling clays characterized by three separate length scales (nano, micro, and macro) and two levels of porosity (nano- and micropores). At the finest (nano) scale the medium is treated as charged clay particles saturated by aqueous electrolyte solution containing monovalent and divalent ions forming the electrical double layer. A new constitutive law is constructed for the disjoining pressure based on the numerical resolution of non-local problem at the nanoscale which, in contrast to the Poisson–Boltzmann theory for point charge ions, is capable of capturing the short-range interactions between the ions due to their finite size. At the intermediate scale (microscale), the two-phase homogenized particle/electrolyte solution system is represented by swollen clay clusters (or aggregates) with the nanoscale disjoining pressure incorporated in a modified form of Terzaghi’s effective principle. At the macroscale, the electro-chemical–mechanical couplings within clay clusters is homogenized with the ion transport in the bulk fluid lying in the micro pores. The resultant macroscopic picture is governed by a three-scale model wherein ion transport takes place in the bulk solution strongly coupled with the mechanics of the clay clusters which play the role of sources/sinks of mass to the bulk fluid associated with ion adsorption/desorption in the electrical double layer at the nanoscale. Within the context of the quasi-steady version of the multiscale model, wherein the electrolyte solution in the nanopores is assumed at instantaneous thermodynamic equilibrium with the bulk fluid in the micropores, we build-up numerically the ion-adsorption isotherms along with the constitutive law of the retardation coefficients of monovalent and divalent ions. In addition, the constitutive law for the macroscopic swelling pressure is reconstructed numerically showing patterns of attractive forces between particles for bivalent ions for particular ranges of bulk concentrations. The three-scale model is applied to numerically simulate ion diffusion in a compacted clay liner underneath a sanitary landfill. Owing to the distinct constitutive behavior of the swelling pressure and partition coefficient for each ionic species, different compaction regimes and diffusion/adsorption patterns, with totally different characteristic time scales, are observed for sodium and calcium migration in the clay liner.     

A new velocity model for clay-sand mixtures1

None of the standard porosity-velocity models (e.g. the time-average equation, Raymer's equations) is satisfactory for interpreting well-logging data over a broad depth range. Clays in the section are the usual source of the difficulty through the bias and scatter that they introduce into the relationship between porosity and P-wave transit time. Because clays are composed of fine sheet-like particles, they normally form pores with much smaller aspect ratios than those associated with sand grains. This difference in pore geometry provides the key to obtaining more consistent resistivity and sonic log interpretations. A velocity model for Clay–sand mixtures has been developed in terms of the Kuster and Toksöz, effective medium and Gassmann theories. In this model, the total pore space is assumed to consist of two parts: (1) pores associated with sand grains and (2) pores associated with clays (including bound water). The essential feature of the model is the assumption that the geometry of pores associated with sand grains is significantly different from that associated with clays. Because of this, porosity in shales affects elastic compliance differently from porosity in sand-Stones. The predictive power of the model is demonstrated by the agreement between its predictions and laboratory measurements and by its ability to predict sonic logs from other logs over large depth intervals where formations vary from unconsolidated to consolidated sandstones and shales.     

Microstructural Characterization of Granular Cast Iron for Permeable Reactive Barriers

Although intensive research on Fe(0) permeable reactive barriers (PRB) for in situ groundwater remediation has been conducted and multiple applications have been installed in the past two decades, some properties of reactive materials in use have not been fully considered and discussed yet. In the present investigation, a typical granular cast iron has been characterized with different techniques. The grain size distribution not only has an influence on the resulting pore geometry and the surface area but material properties significantly differ between fine and coarse grains. Metallographic analyses revealed large differences in both graphite inclusions and microstructures that likely influence the reactivity. Both graphite and cementite proved to be more resistant toward acidic dissolution compared to Fe⁰. The intrinsic material characteristics described here have not been covered in the existing PRB literature.     

Physical characteristics of high-rank coal reservoirs in different coal-body structures and the mechanism of coalbed methane production

The physical characteristics of coal reservoirs are important for evaluating the potential for gas desorption, diffusion, and seepage during coalbed methane (CBM) production, and influence the performance of CBM wells. Based on data from mercury injection experiments, low-temperature liquid nitrogen adsorption, isothermal adsorption, initial velocity tests of methane diffusion, and gas natural desorption data from a CBM field, herein the physical characteristics of reservoirs of high-rank coals with different coal-body structures are described, including porosity, adsorption/desorption, diffusion, and seepage. Geometric models are constructed for these reservoirs. The modes of diffusion are discussed and a comprehensive diffusion-seepage model is constructed. The following conclusions were obtained. First, the pore distribution of tectonically deformed coal is different from that of normal coal. Compared to normal coal, all types of pore, including micropores (<10 nm), transitional pores (10–100 nm), mesopores (100–1000 nm), and macropores (>1000 nm), are more abundant in tectonically deformed coal, especially mesopores and macropores. The increase in pore abundance is greater with increasing tectonic deformation of coal; in addition, the pore connectivity is altered. These are the key factors causing differences in other reservoir physical characteristics, such as adsorption/desorption and diffusion in coals with different coal-body structures. Second, normal and cataclastic coals mainly contain micropores. The lack of macropores and its bad connectivity limit gas desorption and diffusion during the early stage of CBM production. However, the good connectivity of micropores is favorable for gas desorption and diffusion in later gas production stage. Thus, because of the slow decline in the rate of gas desorption, long-term gas production can easily be obtained from these reservoirs. Third, under natural conditions the adsorption/desorption properties of granulated and mylonitized coal are good, and the diffusion ability is also enhanced. However, for in situ reservoir conditions, the high dependence of reservoir permeability on stress results in a weak seepage of gas; thus, desorption and diffusion is limited. Fourth, during gas production, the pore range in which transitional diffusion takes place always increases, but that for Fick diffusion decreases. This is a reason for the reduction in diffusion capacity, in which micropores and transitional pores are the primary factors limiting gas diffusion. Finally, the proposed comprehensive model of CBM production under in situ reservoir conditions elucidates the key factors limiting gas production, which is helpful for selection of reservoir stimulation methods.     

Inclusion of adsorption in the modelling of early diagenesis

A general one-dimensional diagenetic equation is derived which expresses the effect of the following processes upon the depth and time distribution of a dissolved species in a modern subaqueous sediment: adsorptive or ion exchange equilibrium between pore water and grain surfaces; slower chemical or radiogenic reactions in the pore solution and on the grain surfaces; diffusion in solution and on grain surfaces; advection of grains due to depositional burial; and advection of water due both to burial and to flow past the grains. It is shown that the effects of equilibrium adsorption can be expressed in terms of a lowered diffusion coefficient only in the special case where there is negligible flow of water relative to grains, adsorption is described by a simple linear isotherm, and there is nochemical reaction occurring in the sediment (i.e. only radioactive decay or no reaction at all). Consideration of a commonly used special case of the general equation indicates that if equilibrium adsorption is strong, the depositional term in the resulting diagenetic equation may become quantitatively as important as the terms for diffusion and chemical reaction. Thus, when modelling diagenesis, deposition often cannot be ignored. A coupled pair of differential equations are also derived which express the steady state depth distribution for two ions exhibiting equilibrium ion exchange while undergoing chemical reaction, diffusion, and deposition.     

Breaking Up: Comminution Mechanisms in Sheared Simulated Fault Gouge

The microstructural state and evolution of fault gouge has important implications for the mechanical behaviour, and hence the seismic slip potential of faults. We use 3D discrete element (DEM) simulations to investigate the fragmentation processes operating in fault gouge during shear. Our granular fault gouge models consist of aggregate grains, each composed of several thousand spherical particles stuck together with breakable elastic bonds. The aggregate grains are confined between two blocks of solid material and sheared under a given normal stress. During shear, the grains can fragment in a somewhat realistic way leading to an evolution of grain size, grain shape and overall texture. The ‘breaking up’ of the fault gouge is driven by two distinct comminution mechanisms: grain abrasion and grain splitting. The relative importance of the two mechanisms depends on applied normal stress, boundary wall roughness and accumulated shear strain. If normal stress is sufficiently high, grain splitting contributes significantly to comminution, particularly in the initial stages of the simulations. In contrast, grain abrasion is the dominant mechanism operating in simulations carried out at lower normal stress and is also the main fragmentation mechanism during the later stages of all simulations. Rough boundaries promote relatively more grain splitting whereas smooth boundaries favor grain abrasion. Grain splitting (plus accompanying abrasion) appears to be an efficient mechanism for reducing the mean grain size of the gouge debris and leads rapidly to a power law size distribution with an exponent that increases with strain. Grain abrasion (acting alone) is an effective way to generate excess fine grains and leads to a bimodal distribution of grain sizes. We suggest that these two distinct mechanisms would operate at different stages of a fault’s history. The resulting distributions in grain size and grain shape may significantly affect frictional strength and stability. Our results therefore have implications for the earthquake potential of seismically active faults with accumulations of gouge. They may also be relevant to the susceptibility of rockslides since non-cohesive basal shear zones will evolve in a similar way and potentially control the dynamics of the slide.     

A physical model of the high-frequency seismic signal generated by debris flows

We propose a physical model for the high-frequency (>1 Hz) spectral distribution of seismic power generated by debris flows. The modeled debris flow is assumed to have four regions where the impact rate and impulses are controlled by different mechanisms: the flow body, a coarser-grained snout, a snout lip where particles fall from the snout on the bed, and a dilute front composed of saltating particles. We calculate the seismic power produced by this impact model in two end-member scenarios, a thin-flow and thick-flow limit, which assume that the ratio of grain sizes to flow thicknesses are either near unity or much less than unity. The thin-flow limit is more appropriate for boulder-rich flows that are most likely to generate large seismic signals. As a flow passes a seismic station, the rise phase of the seismic amplitude is generated primarily by the snout while the decay phase is generated first by the snout and then the main flow body. The lip and saltating front generate a negligible seismic signal. When ground properties are known, seismic power depends most strongly on both particle diameter and average flow speed cubed, and also depends on length and width of the flow. The effective particle diameter for producing seismic power is substantially higher than the median grain size and close to the 73rd percentile for a realistic grain size distribution. We discuss how the model can be used to estimate effective particle diameter and average flow speed from an integrated measure of seismic power. © 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd. © 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

Soil creep and convex‐upward velocity profiles: theoretical and experimental investigation of disturbance‐driven sediment transport on hillslopes

The movement of unconsolidated materials near the Earth's surface is often driven by disturbances that occur at a range of spatial and temporal scales. The nature of these disturbances ranges from highly variable, such as tree turnover, to periodic and predictable, such as frost heave or creep. To explore the effect of probabilistic disturbances on surface processes, we formulated a granular creep model with analogy to rate process theory (RPT) used for chemical reactions. According to the theory, individual particles must be energized to a height greater than adjacent particles in order for grain dilation and transport to occur. The height of neighbouring particles (which is akin to activation energy in chemical reactions) varies with slope angle such that energy barriers get smaller in the downslope direction as slopes steepen. When slopes approach the friction‐limited angle of repose, the height of energy barriers approaches zero and grains ?ow in the absence of disturbance. An exponential function is used to describe the probability distribution of particle excitation height although alternative distributions are possible. We tested model predictions of granular dynamics in an experimental sandpile. In the sandpile, acoustic energy serves as the disturbance agent such that grains dilate and shear in response. Particle velocities are controlled by the frequency of energy pulses that result in grain displacement. Using tracer particles, we observed a convex‐upward velocity pro?le near the surface of the sandpile, consistent with predictions of our RPT‐based velocity model. In addition, we depth‐integrated the velocity model to predict how ?ux rates vary with inclination of the sandpile and observed non‐linear ?ux–gradient curves consistent with model predictions. By varying the acoustic energy level in the experimental sandpile, we documented changes in the rate of grain movement; similar changes in modelled velocities were achieved by varying the exponent of the particle excitation probability distribution. The general agreement between observed and modelled granular behaviour in our simple laboratory sandpile supports the utility of RPT‐based methods for modelling transport processes (e.g. soil creep, frost heave, and till deformation), thus enabling us to account for the probabilistic nature of disturbances that liberate sediment in natural landscapes. Copyright © 2004 John Wiley & Sons, Ltd.     

Experimental simulation of plagioclase diagenesis atP-T conditions of 3.5 km burial depth

Dissolution of plagioclase under the physical conditions at shallow to intermediate burial depths is a prime candidate for secondary porosity generation in feldspathic siliciclastic sediments. The diagenetic behavior of granular aggregates of plagioclase feldspar and quartz has been investigated by experimentation performed in a Bridgeman-type pressure vessel. The experiments, each of two weeks duration, simulated pressure-temperature conditions approximating 3.5 km burial depth. By using a double-acting pore-fluid reservoir, solutions of various chemistries were cycled through samples composed of oligoclase or labradorite feldspar and quartz (90:10 wt% respectively).Scanning electron microscope analysis of the post-experiment samples reveals dissolution features and precipitated products. Dissolution voids of 10 microns occur typically in areas of maximum stress such as crack-tips and grain contacts. Dissolution on a larger scale is exemplified by topographical smoothing of grain su faces. The dissolved species are subsequently reprecipitated as Ca-enriched overgrowths (possibly zeolites) and clays. These precipitates are found individually on the scale of 10 microns and collectively as surface coatings on both feldspar and quartz grains. Atomic absorption spectroscopic analyses of the pore fluid suggest that the fluid chemistry is consistent with the observed experimental precipitates.These experiments show that clay coatings are unnecessary precursors to grain surface dissolution and that the diagenetic precipitation is not mineral selective. Also, the mass transfer of the dissolved species appears to be localized because grains displaying both dissolution and precipitation features are commonplace. Volume changes due to mineral transformation/alteration may increase secondary porosity if the dissolved species produced from dissolution are only partially involved in reprecipitation and the remaining dissolved material is flushed out by the pore fluids. However, if the mass transfer is primarily local then permeability would significantly decrease as precipitates may choke the pore throats.     

Effect of pore geometry on Gassmann fluid substitution

Although there is no assumption of pore geometry in derivation of Gassmann's equation, the pore geometry is in close relation with hygroscopic water content and pore fluid communication between the micropores and the macropores. The hygroscopic water content in common reservoir rocks is small, and its effect on elastic properties is ignored in the Gassmann theory. However, the volume of hygroscopic water can be significant in shaly rocks or rocks made of fine particles; therefore, its effect on the elastic properties may be important. If the pore fluids in microspores cannot reach pressure equilibrium with the macropore system, assumption of the Gassmann theory is violated. Therefore, due to pore structure complexity, there may be a significant part of the pore fluids that do not satisfy the assumption of the Gassmann theory. We recommend that this part of pore fluids be accounted for within the solid rock frame and effective porosity be used in Gassmann's equation for fluid substitution. Integrated study of ultrasonic laboratory measurement data, petrographic data, mercury injection capillary pressure data, and nuclear magnetic resonance T₂ data confirms rationality of using effective porosity for Gassmann fluid substitution. The effective porosity for Gassmann's equation should be frequency dependent. Knowing the pore geometry, if an empirical correlation between frequency and the threshold pore‐throat radius or nuclear magnetic resonance T₂ could be set up, Gassmann's equation can be applicable to data measured at different frequencies. Without information of the pore geometry, the irreducible water saturation can be used to estimate the effective porosity.     

模拟天然气水合物的岩石物理特性模型实验

针对水合物沉积的悬浮、颗粒接触和胶结三种微观模式,制作一组微弱胶结非固结高孔隙度人造样品和颗粒填充渐变的三维物理模型. 通过物理模型实验研究天然气水合物物性参数的敏感性. 实验结果表明：在弱颗粒间胶结物和低有效应力的固结差的沉积物中,声波对孔隙流体性质较敏感. 随着温度的降低颗粒胶结,改变原有沉积物的岩石物理特性,速度、弹性模量和频率升高,声波衰减和Vp/VS减小,沉积层内的反射波消隐.     

Image-based relative permeability upscaling from the pore scale

High resolution images acquired from X-ray μ-CT are able to map the internal structure of porous media on which multiphase flow properties can be computed. While the resolution of a few micrometers may be sufficient for capturing the pore space of many sandstones, most carbonates exhibit a large amount of microporosity; pores which are below the image resolution and are not resolved at specific resolution. Neglecting the effect of micropores on fluid flow and transport properties of these rocks can cause erroneous results in particular at partial saturations. Current image-based pore scale models typically only consider macropores for simulating fluid flow. In this paper, we quantify the effect of microporosity on the effective permeability of the wetting phase for heterogeneous model structures with varying amount of micro-to-macro porosity. A multi-scale numerical approach is proposed to couple an average effect of micropores with an explicit representation of macropores. The Brinkman equation is solved using a lattice Boltzmann formulation to facilitate the coupling of Darcy and Stokes equations in micropores and macropores, respectively. The results show good agreement between the fine scale solution and the results of the upscaled models in which microporous regions are homogenised. The paper analyses in particular the choice of the momentum sink parameter at low wetting phase saturations. It is shown that this parameter can be found using either a flux-based calculation of permeability of microporous regions or chosen purely on the basis of the effective permeability of these regions.     

The effect of microporosity on transport properties in porous media

Sizeable amounts of connected microporosity with various origins can have a profound effect on important petrophysical properties of a porous medium such as (absolute/relative) permeability and capillary pressure relationships. We construct pore-throat networks that incorporate both intergranular porosity and microporosity. The latter originates from two separate mechanisms: partial dissolution of grains and pore fillings (e.g. clay). We then use the reconstructed network models to estimate the medium flow properties. In this work, we develop unique network construction algorithms and simulate capillary pressure–saturation and relative permeability–saturation curves for cases with inhomogeneous distributions of pores and micropores. Furthermore, we provide a modeling framework for variable amounts of cement and connectivity of the intergranular porosity and quantifying the conditions under which microporosity dominates transport properties. In the extreme case of a disconnected inter-granular network due to cementation a range of saturations within which neither fluid phase is capable of flowing emerges. To our knowledge, this is the first flexible pore scale model, from first principles, to successfully approach this behavior observed in tight reservoirs.     

Geostatistical modeling of clay spatial distribution in siliciclastic rock samples using the plurigaussian simulation method

In order to implement secondary and enhanced oil recovery processes in complex terrigenous formations as is usual in turbidite deposits, a precise knowledge of the spatial distribution of shale grains is a crucial element for the fluid flow prediction. The reason of this is that the interaction of water with shale grains can significantly modify their size and/or shape, which in turn would cause porous space sealing with the subsequent impact in the flow. In this work, a methodology for stochastic simulations of spatial grains distributions obtained from scanning electron microscopy images of siliciclastic rock samples is proposed. The aim of the methodology is to obtain stochastic models would let us investigate the shale grain behavior under various physico-chemical interactions and flux regimes, which in turn, will help us get effective petrophysical properties (porosity and permeability) at core scale. For stochastic spatial grains simulations a plurigaussian method is applied, which is based on the truncation of several standard Gaussian random functions. This approach is very flexible, since it allows to simultaneously manage the proportions of each grain category in a very general manner and to rigorously handle their spatial dependency relationships in the case of two or more grain categories. The obtained results show that the stochastically simulated porous media using the plurigaussian method adequately reproduces the proportions, basic statistics and sizes of the pore structures present in the studied reference images.     

ULTRASONIC SPECTROSCOPY IN ROCKS: AN EXPERIMENTAL STUDY OF HIGHLY POROUS SYNTHETIC SANDSTONES1

Broad-band ultrasonic impulses and the attenuation spectral ratio technique have been used to investigate the characteristics of synthetic rocks for different saturants, by obtaining the attenuation coefficients in the frequency range of 0.1-1.0 MHz. The general trend of the data indicates that the experimental attenuation coefficient increases with an increase in the grain/pore size. Results show that large grains/pores can increase the attenuation coefficient in all mechanisms, and water-saturated synthetic sandstone rocks result in a higher attenuation coefficient than similar oil-saturated rocks.     

A multiple-point statistics algorithm for 3D pore space reconstruction from 2D images

Fluid flow behavior in a porous medium is a function of the geometry and topology of its pore space. The construction of a three dimensional pore space model of a porous medium is therefore an important first step in characterizing the medium and predicting its flow properties. A stochastic technique for reconstruction of the 3D pore structure of unstructured random porous media from a 2D thin section training image is presented. The proposed technique relies on successive 2D multiple point statistics simulations coupled to a multi-scale conditioning data extraction procedure. The Single Normal Equation Simulation Algorithm (SNESIM), originally developed as a tool for reproduction of long-range, curvilinear features of geological structures, serves as the simulation engine. Various validating criteria such as marginal distributions of pore and grain, directional variograms, multiple-point connectivity curves, single phase effective permeability and two phase relative permeability calculations are used to analyze the results. The method is tested on a sample of Berea sandstone for which a 3D micro-CT scanning image is available. The results confirm that the equi-probable 3D realizations obtained preserve the typical patterns of the pore space that exist in thin sections, reproduce the long-range connectivities, capture the characteristics of anisotropy in both horizontal and vertical directions and have single and two phase flow characteristics consistent with those of the measured 3D micro-CT image.     

骨架导电的混合泥质砂岩通用孔隙结合电阻率模型研究

在水介质中顺序添加分散粘土颗粒、油珠、导电骨架颗粒、层状泥质,并对每一种成分进行连续积分,建立了一种适用骨架导电及含有分散粘土和层状泥质的泥质砂岩通用电阻率模型.通过对该模型的影响因素分析,发现泥质分布形式对模型计算的含水饱和度有很大影响;对应两个不同粘土颗粒电阻率或骨架颗粒电阻率的地层电导率之差,几乎与总含水饱和度无关,而对应两个不同层状泥质电阻率的地层电导率之差,随总含水饱和度增大而增大;骨架胶结指数变化对地层电导率与总含水饱和度关系曲线的影响最大,而粘土胶结指数变化对地层电导率与总含水饱和度关系曲线的影响最小;饱和度指数对地层电导率与总含水饱和度关系曲线的影响随总含水饱和度的增大而减小.通过一组骨架导电的人造岩样的试验,表明当地层水电阻率.小于颗粒电阻率时,该模型可以用于不含粘土的骨架导电的岩石.通过两组分散泥质砂岩岩样实验测量数据和一组层状泥质砂岩测井资料及实际测井资料的计算,表明本文给出的电阻率模型既适用于分散泥质砂岩地层解释又适用于层状泥质砂岩地层解释,同时,还适用于含有分散粘土和层状泥质的混合泥质砂岩地层解释.     

The stress dependent elastic properties of thermally induced microfractures in aeolian Rotliegend sandstone

The impact of thermally induced microfractures on the stress‐sensitive elastic wave properties of aeolian Rotliegend sandstone samples is analysed. It is found that to identify the effects of the microfracture contribution accurately, a correction must first be made to account for water loss (representing a mass loss of 4–6%) from the pore throats and clays due to the heating process, despite care being taken to ensure that the thermally fractured samples re‐adsorb room moisture. Both the original and thermally fractured rocks are stress‐sensitive at the ultrasonic wave frequencies of the laboratory. However, a distinct shift in the estimated distribution of internal rock compliance indicates that the population of thermal microfractures differs in nature from that caused solely by core‐plug extraction damage. In particular, the ratio of normal to tangential compliance is observed to be higher for the thermally generated microfractures than for the broken grain‐grain contacts created by extraction unloading. This can be explained by the intragranular thermal‐fracture surfaces being smoother when compared to the intergranular boundaries. Mechanical hysteresis is observed between the up‐ and downgoing test cycles for both the original and, to a greater extent, the thermally fractured rock. This indicates that there is compaction‐induced movement of the fractures in the samples during application of stress in the laboratory.     

Grain size and chemical controls on the ductile properties of mostly frictional faults at low-temperature hydrothermal conditions

A conceptually simple process which establishes a steady grain size distribution is envisioned to control the ductile creep properties of fault zones that mainly slip by frictional processes. Fracture during earthquakes and aseismic frictional creep tend to reduce grain size. However, sufficiently small grains tend to dissolve so that larger grains grow at their expense, a process called Ostwald ripening. A dynamic stedy state is reached where grain size reduction by fracture is balanced by grain growth from Ostwald ripening. The ductile creep mechanism within fault zones in hard rock is probably pressure solution where the rate is limited by diffusion along load-bearing grain-grain contacts. The diffusion paths that limit Ostwald repening are to a considerable extent the same as those for pressure solution. Active Ostwald ripening thus implies conditions suitable for ductile creep. An analytic theory allows estimation of the steady-state mean grain size and the viscosity for creep implied by this dynamic steady state from material properties and from the width, shear traction, and long-term slip velocity of the fault zone. Numerical models were formulated to compute the steady state grain size distribution. The results indicate that ductile creep, as suggested in the companion paper, is a plausible mechanism for transiently increasing fluid pressure within mostly sealed fault zones so that frictional failure occurs at relatively low shear tractions, 10 MPa. The relevant material properties are too poorly known, however, for the steady state theory (or its extension to a fault that slips in infrequent large earthquakes) to have much predictive value without additional laboratory experiments and studies of exhumed faults.